Protostellar Disk Formation Research Articles

The ALMA Legacy Survey of Class 0/I Disks in Corona australis, Aquila, chaMaeleon, oPhiuchus north, Ophiuchus, Serpens (CAMPOS). I. Evolution of Protostellar Disk Radii

We surveyed nearly all the embedded protostars in seven nearby clouds (Corona Australis, Aquila, Chamaeleon I and II, Ophiuchus North, Ophiuchus, Serpens) with the Atacama Large Millimeter/submillimeter Array at 1.3 mm observations with a resolution of 0.″1. This survey detected 184 protostellar disks, 90 of which were observed at a resolution of 14–18 au, making it one of the most comprehensive high-resolution disk samples across various protostellar evolutionary stages to date. Our key findings include the detection of new annular substructures in two Class I and two flat-spectrum sources, while 21 embedded protostars exhibit distinct asymmetries or substructures in their disks. We find that protostellar disks have a substantially large variability in their radii across all evolutionary classes. In particular, the fraction of large disks with sizes above 60 au decreases as the protostar evolves from Class 0 to Class I. Compiling the literature data, we discovered an increasing trend of the gas disk radii to dust disk radii ratio (R gas,Kep/R mm) with increasing bolometric temperature (T bol). Our results indicate that the dust and gas disk radii decouple during the early Class I stage. However, in the Class 0 stage, the dust and gas disk sizes are similar, which allows for a direct comparison between models and observational data at the earliest stages of protostellar evolution. We show that the distribution of radii in the 52 Class 0 disks in our sample is in high tension with various disk formation models, indicating that protostellar disk formation remains an unsolved question.

The First Estimation of the Ambipolar Diffusivity Coefficient from Multi-scale Observations of the Class 0/I Protostar, HOPS-370

Protostars are born in magnetized environments. As a consequence, the formation of protostellar disks can be suppressed by the magnetic field, efficiently removing the angular momentum of the infalling material. Nonideal MHD effects are proposed as one way to allow protostellar disks to form. Thus, it is important to understand their contributions to observations of protostellar systems. We derive an analytical equation to estimate the ambipolar diffusivity coefficient at the edge of the protostellar disk in the Class 0/I protostar, HOPS-370, for the first time, under the assumption that the disk radius is set by ambipolar diffusion. Using previous results of the protostellar mass, disk mass, disk radius, density and temperature profiles, and magnetic field strength, we estimate the ambipolar diffusivity coefficient to be 1.7−1.4+1.5×1019cm2s−1 . We quantify the contribution of ambipolar diffusion by estimating its dimensionless Elsässer number to be ∼1.7−1.0+1.0 , indicating its dynamical importance in this region. We compare our results to those of the chemical calculations of the ambipolar diffusivity coefficient using the Non-Ideal Magnetohydrodynamics Coefficients and Ionization Library, which are consistent with our results. In addition, we compare our derived ambipolar diffusivity coefficient to the diffusivity coefficients for ohmic dissipation and the Hall effect, and find ambipolar diffusion is dominant in our density regime. These results demonstrate a new methodology to understand nonideal MHD effects in observations of protostellar disks. More detailed modeling of the magnetic field, envelope, and microphysics, along with a larger sample of protostellar systems, is needed to further understand the contributions of nonideal MHD.

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot, and high-pressure regions that formed during the gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas and dust interaction and resolve potential zones of dust concentration during star and disc formation stages. We explore whether the thermal and dynamical conditions that developed during protostellar disc formation can generate gas flows that efficiently mix and transport the well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a two-dimensional axisymmetric geometry and followed the azimuthal component of the velocity. The dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10, and 100 micron-sized neutral, spherical dust grains. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 $M_ odot $ and explore the effects of different initial rotation rates and cloud core sizes. Our study underlines mechanisms for the early transport of dust from the inner hot disc regions via the occurrence of two transient gas motions, namely meridional flow and outflow. The vortical flow fosters dynamical mixing and retention of dust, while the thermal pressure driven outflow replenishes dust in the outer disc. Notably, these phenomena occur regardless of the initial cloud core mass, size, and rotation rate. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to the outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse.

Polarization Holes as an Indicator of Magnetic Field−Angular Momentum Alignment. I. Initial Tests

The formation of protostellar disks is still a mystery, largely due to the difficulties in observations that can constrain theories. For example, the 3D alignment between the rotation of the disk and the magnetic fields (B-fields) in the formation environment is critical in some models, but so far it is impossible to observe. Here we study the possibility of probing the alignment between B-field and disk rotation using “polarization holes” (PHs). PHs are widely observed and are caused by unresolved B-field structures. With ideal magnetohydrodynamic simulations, we demonstrate that different initial alignments between B-field and angular momentum (AM) can result in B-field structures that are distinct enough to produce distinguishable PHs. Thus, PHs can potentially serve as probes for alignments between B-field and AM in disk formation.

Multiple power-law tails in the density and column-density distribution in contracting star-forming clumps

ABSTRACT We present a numerical study of the evolution of power-law tails (PLTs) in the (column-)density distributions (N-PDF, ρ-PDF) in contracting star-forming clumps in primordial gas, without and with some initial rotational and/or turbulent support. In all considered runs multiple PLTs emerge shortly after the formation of the first protostar. The first PLT (PLT 1) in the ρ-PDF is a stable feature with slope q1 ≃ −1.3 which corresponds – under the condition of preserved spherical symmetry – to the outer envelope of the protostellar object with density profile ρ ∝ l−2 in the classical Larson–Penston collapse model, where l is the radius. The second PLT (PLT 2) in the ρ-PDF is stable in the pure-infall runs but fluctuates significantly in the runs with initial support against gravity as dozens of protostars form and their mutual tidal forces change the density structure. Its mean slope, 〈q2〉 ≃ −2, corresponds to a density profile of ρ ∝ l−3/2 which describes a core in free fall in the classical Larson–Penston collapse model or an attractor solution at scales with dominating protostellar gravity. PLT 1 and PLT 2 in the N-PDFs are generally consistent with the observational data of Galactic low-mass star-forming regions from Herschel data. In the runs with initial support against gravity a third PLT (PLT 3) in the ρ-PDFs appears simultaneously with or after the emergence of PLT 2. It is very shallow, with mean slope of 〈q3〉 ≃ −1, and is associated with the formation of thin protostellar accretion discs.

Possible episodic infall towards a compact disk in B335

Context. Previous observations of the isolated Class 0 source B335 have presented evidence of ongoing infall in various molecular lines, such as HCO+, HCN, and CO. There have been no confirmed observations of a rotationally supported disk on scales greater than ~12 au. Aims. The presence of an outflow in B335 suggests that a disk is also expected to be present or undergoing formation. To constrain the earliest stages of protostellar evolution and disk formation, we aim to map the region where gas falls inwards and observationally constrain its kinematics. Furthermore, we aim to put strong limits on the size and orientation of any disk-like structure in B335. Methods. We used high angular resolution 13CO data from the Atacama Large Millimeter/submillimeter Array (ALMA) and combined it with shorter-baseline archival data to produce a high-fidelity image of the infall in B335. We also revisited the imaging of high-angular resolution Band 6 continuum data to study the dust distribution in the immediate vicinity of B335. Results. Continuum emission shows an elliptical structure (10 by 7 au) with a position angle 5 degrees east of north, consistent with the expectation for a forming disk in B335. A map of the infall velocity (as estimated from the 13CO emission), shows evidence of asymmetric infall, predominantly from the north and south. Close to the protostar, infall velocities appear to exceed free-fall velocities. Three-dimensional (3D) radiative transfer models, where the infall velocity is allowed to vary within the infall region, may explain the observed kinematics. Conclusions. The data suggest that a disk has started to form in B335 and that gas is falling towards that disk. However, kinematically-resolved line data towards the disk itself is needed to confirm the presence of a rotationally supported disk around this young protostar. The high infall velocities we measured are not easily reconcilable with a magnetic braking scenario, suggesting that there is a pressure gradient that allows the infall velocity to vary in the region.

Predicting Stellar Mass Accretion: An Optimized Echo State Network Approach in Time Series Modeling

Modeling the dynamics of the formation and evolution of protostellar disks as well as the history of stellar mass accretion typically involve the numerical solution of complex systems of coupled differential equations. The resulting mass accretion history of protostars is known to be highly episodic due to recurrent instabilities and also exhibits short timescale flickering. By leveraging the strong predictive abilities of neural networks, we extract some of the critical temporal dynamics experienced during the mass accretion including periods of instability. Particularly, we utilize a novel form of the Echo-State Neural Network (ESN), which has been shown to efficiently deal with data having inherent nonlinearity. We introduce the use of Optimized-ESN (Opt-ESN) to make model-independent time series forecasting of mass accretion rate in the evolution of protostellar disks. We apply the network to multiple hydrodynamic simulations with different initial conditions and exhibiting a variety of temporal dynamics to demonstrate the predictability of the Opt-ESN model. The model is trained on simulation data of $\sim 1-2$ Myr, and achieves predictions with a low normalized mean square error ($\sim 10^{-5}$ to $10^{-3}$) for forecasts ranging between 100 and 3800 yr. This result shows the promise of the application of machine learning based models to time-domain astronomy.

Testing a New Model of Embedded Protostellar Disks against Observations: The Majority of Orion Class 0/I Disks Are Likely Warm, Massive, and Gravitationally Unstable

We formulate a parameterized model of embedded protostellar disks and test its ability to estimate disk properties by fitting dust-continuum observations. The main physical assumptions of our model are motivated by a recent theoretical study of protostellar disk formation; these assumptions include that the disk should be marginally gravitationally unstable, and that the dominant dust heating mechanism is internal accretion heating instead of external protostellar irradiation. These assumptions allow our model to estimate reliably the disk mass even when the observed emission is optically thick and to determine self-consistent disk (dust) temperatures. Using our model to fit multiwavelength observations of 163 disks in the VANDAM Orion survey, we find that the majority (57%) of this sample can be fit well by our model. Using our model, we produce new estimates of Orion protostellar disk properties. We find that these disks are generally warm and massive, with a typical star-to-disk mass ratio in Class 0/I. We also discuss why our estimates differ from those in previous studies and the implications of our results on disk evolution and fragmentation.

Grain growth during protostellar disc formation

ABSTRACT Recent observations indicate that mm/cm-sized grains may exist in the embedded protostellar discs. How such large grains grow from the micron size (or less) in the earliest phase of star formation remains relatively unexplored. In this study, we take a first step to model the grain growth in the protostellar environment, using 2D (axisymmetric) radiation hydrodynamic and grain growth simulations. We show that the grain growth calculations can be greatly simplified by the ‘terminal velocity approximation’, where the dust drift velocity relative to the gas is proportional to its stopping time, which is proportional to the grain size. We find that the grain–grain collision from size-dependent terminal velocity alone is too slow to convert a significant fraction of the initially micron-sized grains into mm/cm sizes during the deeply embedded Class 0 phase. Substantial grain growth is achieved when the grain–grain collision speed is enhanced by a factor of 4. The dust growth above and below the disc midplane enables the grains to settle faster towards the midplane, which increases the local dust-to-gas ratio, which, in turn, speeds up further growth there. How this needed enhancement can be achieved is unclear, although turbulence is a strong possibility that deserves further exploration.

Centrifugal barrier and super-Keplerian rotation in protostellar disc formation

ABSTRACT With the advent of ALMA, it is now possible to observationally constrain how discs form around deeply embedded protostars. In particular, the recent ALMA C3H2 line observations of the nearby protostar L1527 have been interpreted as evidence for the so-called ‘centrifugal barrier,’ where the protostellar envelope infall is gradually decelerated to a stop by the centrifugal force in a region of super-Keplerian rotation. To test the concept of centrifugal barrier, which was originally based on angular momentum conserving-collapse of a rotating test particle around a fixed point mass, we carry out simple axisymmetric hydrodynamic simulations of protostellar disc formation including a minimum set of ingredients: self-gravity, rotation, and a prescribed viscosity that enables the disc to accrete. We find that a super-Keplerian region can indeed exist when the viscosity is relatively large but, unlike the classic picture of centrifugal barrier, the infalling envelope material is not decelerated solely by the centrifugal force. The region has more specific angular momentum than its surrounding envelope material, which points to an origin in outward angular momentum transport in the disc (subject to the constraint of disc expansion by the infalling envelope), rather than the spin-up of the envelope material envisioned in the classic picture as it falls closer to the centre in order to conserve angular momentum. For smaller viscosities, the super-Keplerian rotation is weaker or non-existing. We conclude that, despite the existence of super-Keplerian rotation in some parameter regime, the classic picture of centrifugal barrier is not supported by our simulations.

Influence of the magnetic field on the formation of protostellar disks

Abstract We numerically model the collapse of magnetic rotating protostellar clouds with mass of 10 M ⊙ {M}_{\odot } . The simulations are carried out with the help of 2D MHD code Enlil. The structure of the cloud at the isothermal stage of the collapse is investigated for the cases of weak, moderate, and strong initial magnetic field. Simulations reveal the universal hierarchical structure of collapsing protostellar clouds, consisting of the flattened envelope with the qausi-magnetostatc disk inside and the first core in its center. The size of the primary disk increases with the initial magnetic energy of the cloud. The magnetic braking efficiently transports the angular momentum from the primary disk into the envelope in the case, when the initial magnetic energy of the cloud is more than 20% of its gravitational energy. The intensity of the outflows launched from the region near the boundary of the first core increases with initial magnetic energy. The “dead” zone with small ionization fraction, x < 1 0 − 11 x\lt 1{0}^{-11} , forms inside the first hydrostatic core and at the base of the outflow. Ohmic dissipation and ambipolar diffusion determine conditions for further formation of the protostellar disk in this region.

Magnetic Spirals in Accretion Flows Originated from Misaligned Magnetic Fields

Abstract Misalignment between rotation and magnetic fields has been suggested to be one type of physical mechanism that can ease the effects of magnetic braking during the collapse of cloud cores leading to the formation of protostellar disks. However, its essential factors are poorly understood. Therefore, we perform a more detailed analysis of the physics involved. We analyze existing simulation data to measure the system torques, mass accretion rates, and Toomre Q parameters. We also examine the presence of shocks in the system. While advective torques are generally the strongest, we find that magnetic and gravitational torques can play substantial roles in how angular momentum is transferred during the disk formation process. Magnetic torques can shape the accretion flows, creating two-armed magnetized inflow spirals aligned with the magnetic field. We find evidence of an accretion shock that is aligned according to the spiral structure of the system. Inclusion of ambipolar diffusion as explored in this work has shown a slight influence in the small-scale structures but not in the main morphology. We discuss potential candidate systems where some of these phenomena could be present.

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

ABSTRACT We use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

The interplay between ambipolar diffusion and Hall effect on magnetic field decoupling and protostellar disc formation

ABSTRACT Non-ideal magnetohydrodynamic (MHD) effects have been shown recently as a robust mechanism of averting the magnetic braking ‘catastrophe’ and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay between ambipolar diffusion (AD) and Hall effect during the protostellar collapse. We find that removing the smallest grain population (≲10 nm) from the standard MRN size distribution is sufficient for enabling disc formation. Further varying the grain sizes can result in either a Hall-dominated or an AD-dominated collapse; both form discs of tens of au in size regardless of the magnetic field polarity. The direction of disc rotation is bimodal in the Hall-dominated collapse but unimodal in the AD-dominated collapse. We also find that AD and Hall effect can operate either with or against each other in both radial and azimuthal directions, yet the combined effect of AD and Hall is to move the magnetic field radially outward relative to the infalling envelope matter. In addition, microphysics and magnetic field polarity can leave profound imprints both on observables (e.g. outflow morphology, disc to stellar mass ratio) and on the magnetic field characteristics of protoplanetary discs. Including Hall effect relaxes the requirements on microphysics for disc formation, so that prestellar cores with cosmic ray ionization rate of ≲2–3 × 10−16 s−1 can still form small discs of ≲10 au radius. We conclude that disc formation should be relatively common for typical prestellar core conditions, and that microphysics in the protostellar envelope is essential to not only disc formation, but also protoplanetary disc evolution.

Do we need non-ideal magnetohydrodynamic to model protostellar discs?

ABSTRACT We investigate and discuss protostellar discs in terms of where the various non-ideal magnetohydrodynamics (MHD) processes are important. We find that the traditional picture of a magnetized disc (where Ohmic resistivity is dominant near the mid-plane, surrounded by a region dominated by the Hall effect, with the remainder of the disc dominated by ambipolar diffusion) is a great oversimplification. In simple parametrized discs, we find that the Hall effect is typically the dominant term throughout the majority of the disc. More importantly, we find that in much of our parametrized discs, at least two non-ideal processes have coefficients within a factor of 10 of one another, indicating that both are important and that naming a dominant term underplays the importance of the other terms. Discs that were self-consistently formed in our previous studies are also dominated by the Hall effect, and the ratio of ambipolar diffusion and Hall coefficients is typically less than 10, suggesting that both terms are equally important and listing a dominant term is misleading. These conclusions become more robust once the magnetic field geometry is taken into account. In agreement with the literature we review, we conclude that non-ideal MHD processes are important for the formation and evolution of protostellar discs. Ignoring any of the non-ideal processes, especially ambipolar diffusion and the Hall effect, yields an incorrect description of disc evolution.

Gravitoviscous protoplanetary discs with a dust component – IV. Disc outer edges, spectral indices, and opacity gaps

ABSTRACT The crucial initial step in planet formation is the agglomeration of micron-sized dust into macroscopic aggregates. This phase is likely to happen very early during the protostellar disc formation, which is characterized by active gas dynamics. We present numerical simulations of protostellar/protoplanetary disc long-term evolution, which includes gas dynamics with self-gravity in the thin-disc limit, and bidisperse dust grain evolution due to coagulation, fragmentation, and drift through the gas. We show that the decrease of the grain size to the disc periphery leads to sharp outer edges in dust millimetre emission, which are explained by a drop in dust opacity coefficient rather than by dust surface density variations. These visible outer edges are at the location where average grain size ≈λ/2$\pi$, where λ is the observational wavelength, so discs typically look more compact at longer wavelengths if dust size decreases outwards. This allows a simple recipe for reconstructing grain sizes in disc outer regions. Discs may look larger at longer wavelengths if grain size does not reach λ/2$\pi$ for some wavelength. Disc visible sizes evolve non-monotonically over the first million years and differ from dust and gas physical sizes by factor of a few. We compare our model with recent observation data on gas and dust disc sizes, far-infrared fluxes, and spectral indices of protoplanetary discs in Lupus. We also show that non-monotonic variations of the grain size in radial direction can cause wavelength-dependent opacity gaps, which are not associated with any physical gaps in the dust density distribution.

Non-ideal magnetohydrodynamics versus turbulence – I. Which is the dominant process in protostellar disc formation?

ABSTRACT Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1 M⊙ gas cloud that is threaded with a magnetic field and imposed with both rotational and turbulent velocities. We investigate magnetic fields that are parallel/antiparallel and perpendicular to the rotation axis, two rotation rates, and four Mach numbers. Disc formation occurs preferentially in the models that include non-ideal MHD where the magnetic field is antiparallel or perpendicular to the rotation axis. This is independent of the initial rotation rate and level of turbulence, suggesting that subsonic turbulence plays a minimal role in influencing the formation of discs. Aside from first core outflows that are influenced by the initial level of turbulence, non-ideal MHD processes are more important than turbulent processes during the formation of discs around low-mass stars.

Protostellar disk formation by a nonrotating, nonaxisymmetric collapsing cloud: model and comparison with observations

Context. Planet-forming disks are fundamental objects that are thought to be inherited from large scale rotation through the conservation of angular momentum during the collapse of a prestellar dense core. Aims. We investigate the possibility for a protostellar disk to be formed from a motionless dense core that contains nonaxisymmetric density fluctuations. The rotation is thus generated locally by the asymmetry of the collapse. Methods. We study the evolution of the angular momentum in a nonaxisymmetric collapse of a dense core from an analytical point of view. To test the theory, we performed three-dimensional simulations of a collapsing prestellar dense core using adaptative mesh refinement. We started from a nonaxisymmetrical situation, considering a dense core with random density perturbations that follow a turbulence spectrum. We analyzed the emerging disk by comparing the angular momentum it contains with the one expected from our analytic development. We studied the velocity gradients at different scales in the simulation as is done with observations. Results. We show that the angular momentum in the frame of a stellar object, which is not located at the center of mass of the core, is not conserved due to inertial forces. Our simulations of such nonaxisymmetrical collapse quickly produce accretion disks at the small scales in the core. The analysis of the kinematics at different scales in the simulated core reveals projected velocity gradients of amplitudes similar to the ones observed in protostellar cores and for which directions vary, sometimes even reversing when small and large scales are compared. These complex kinematics patterns appear in recent observations and could be a discriminating feature with models where rotation is inherited from large scales. Our results from simulations without initial rotation are more consistent with these recent observations than when solid-body rotation is initially imprinted. Lastly, we show that the disks that formed in this scenario of nonaxisymmetrical gravitational collapse grow to reach sizes larger than those that are observed, and then fragment. We show that including a magnetic field in these simulations reduces the size of the outcoming disks and it prevents them from fragmenting, as is shown by previous studies. Conclusions. We show that in a nonaxisymmetrical collapse, the formation of a disk can be induced by small perturbations of the initial density field in the core, even in the absence of global large-scale rotation of the core. In this scenario, large disks are generic features that are natural consequences of the hydrodynamical fluid interactions and self-gravity. Since recent observations have shown that most disks are significantly smaller and have a size of a few tens of astronomical units, our study suggests that magnetic braking is the most likely explanation. The kinematics of our model are consistent with typically observed values of velocity gradients and specific angular momentum in protostellar cores. These results open a new avenue in which our understanding of the early phases of disk formation can be explored since they suggest that a fraction of the protostellar disks could be the product of nonaxisymmetrical collapse, rather than directly resulting from the conservation of preexisting large scale angular momentum in rotating cores.

Hall effect in protostellar disc formation and evolution

ABSTRACT The Hall effect is recently shown to be efficient in magnetized dense molecular cores and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and ohmic dissipation, with realistic magnetic diffusivities computed from our equilibrium chemical network. We find that Hall effect only becomes efficient when the large population of very small grains (VSGs: ≲100 Å) is removed from the standard Mathis–Rumpl–Nordsieck size distribution. With such an enhanced Hall effect, however, the bimodality of disc formation does not continue into the main accretion phase. The outer part of the initial ∼40 au disc formed in the anti-aligned configuration ($\boldsymbol {\Omega \cdot B}\lt 0$) flattens into a thin rotationally supported Hall current sheet as Hall effect moves the poloidal magnetic field radially inward relative to matter, leaving only the inner ≲10–20 au RSD. In the aligned configuration ($\boldsymbol {\Omega \cdot B}\gt 0$), disc formation is suppressed initially but a counter-rotating disc forms subsequently due to efficient azimuthal Hall drift. The counter-rotating disc first grows to ∼30 au as Hall effect moves the magnetic field radially outward, but only the inner ≲10 au RSD is long lived like in the anti-aligned case. Besides removing VSGs, cosmic ray ionization rate should be below a few 10−16 s−1 for Hall effect to be efficient in disc formation. We conclude that Hall effect produces small ≲10–20 au discs regardless of the polarity of the magnetic field, and that radially outward diffusion of magnetic fields remains crucial for disc formation and growth.

Investigating the complex velocity structures within dense molecular cloud cores with GBT-Argus

ABSTRACT We present the first results of high-spectral resolution (0.023 km s−1) N2H+ observations of dense gas dynamics at core scales (∼0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position–velocity (PV) space for all five targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming discs nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disc formation.