We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.

We present our analysis of the magnetic field structures from 6000 au to 100 au scales in the Class 0 protostar B335 inferred from our JCMT POL-2 observations and the ALMA archival polarimetric data. To interpret the observational results, we perform a series of (non-)ideal MHD simulations of the collapse of a rotating non-turbulent dense core, whose initial conditions are adopted to be the same as observed in B335, and generate synthetic polarization maps. The comparison of our JCMT and simulation results suggests that the magnetic field on a 6000 au scale in B335 is pinched and well aligned with the bipolar outflow along the east-west direction. Among all our simulations, the ALMA polarimetric results are best explained with weak magnetic field models having an initial mass-to-flux ratio of 9.6. However, we find that with the weak magnetic field, the rotational velocity on a 100 au scale and the disk size in our simulations are larger than the observational estimates by a factor of several. An independent comparison of our simulations and the gas kinematics in B335 observed with the SMA and ALMA favors strong magnetic field models with an initial mass-to-flux ratio smaller than 4.8. We discuss two possibilities resulting in the different magnetic field strengths inferred from the polarimetric and molecular-line observations, (1) overestimated rotational-to-gravitational energy in B335 and (2) additional contributions in the polarized intensity due to scattering on a 100 au scale.

Distinguishing diagnostic observational signatures produced by MHD models is essential in understanding the physics for the formation of protostellar disks in the ALMA era. Developing suitable tools along with time evolution will facilitate better identification of diagnostic features. With a ray-tracing based radiative transfer code Perspective, we explore time evolution of MHD models carried out in Li, Krasnopolsky & Shang (2013) - most of which have $90^\circ$ misalignment between the rotational axis and the magnetic field. Four visible object types can be characterized, origins of which are dependent on the initial conditions. Our results show complex spiraling density, velocity and polarization structures. The systems are under constant change, but many of those distinctive features are present already early on, and they grow more visible in time, but most could not be identified from the data without examining their change in time. The results suggest that spiraling pseudodisk structures could function as an effective observation signature of the formation process, and we witness accretion in the disk with eccentric orbits which appear as spiral-like perturbation from simple circular Keplerian orbits. Magnetically aligned polarization appears purely azimuthal in the disk and magnetic field can lead to precession of the disk.

Previous axisymmetric investigations in two dimensions (2D) have shown that rings and gaps develop naturally in non-ideal magnetohydrodynamic (MHD) disk-wind systems, especially in the presence of ambipolar diffusion (AD). Here we extend these 2D simulations to three dimensions (3D) and find that rings and gaps are still formed in the presence of a moderately strong ambipolar diffusion. The rings and gaps form from the same basic mechanism that was identified in the 2D simulations, namely, the redistribution of the poloidal magnetic flux relative to the disk material as a result of the reconnection of a sharply pinched poloidal magnetic field lines. Thus, the less dense gaps are more strongly magnetized with a large poloidal magnetic field compared to the less magnetized (dense) rings. The rings and gaps start out rather smoothly in 3D simulations that have axisymmetric initial conditions. Non-axisymmetric variations arise spontaneously at later times, but they do not grow to such an extent as to disrupt the rings and gaps. These disk substructures persist to the end of the simulations, lasting up to 3000 orbital periods at the inner edge of the simulated disk. The longevity of the perturbed yet still coherent rings make them attractive sites for trapping large grains that would otherwise be lost to rapid radial migration due to gas drag. As the ambipolar diffusivity decreases, both the disk and the wind become increasingly turbulent, driven by the magnetorotational instability, with tightly-wound spiral arms becoming more prominent in the disk.

Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal MHD effects proposed to enable the formation of Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare gas kinematics traced by the H13CO+ (3-2) and C18O (2-1) lines in the Class 0 protostar B335 obtained with our ALMA observations. A central compact (~1"-2") component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. We constructed kinematical models to fit the observed velocity structures and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. The infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85+/-0.2 km/s at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km/s at a radius of 100 au. The Hall parameter for H13CO+ is estimated to be >>1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between ions and neutrals could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0.1") are needed to examine this possibility and ambipolar diffusion on a smaller scale.

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded.

Radial substructures in circumstellar disks are now routinely observed by ALMA. There is also growing evidence that disk winds drive accretion in such disks. We show through 2D (axisymmetric) simulations that rings and gaps develop naturally in magnetically coupled disk-wind systems on the scale of tens of au, where ambipolar diffusion (AD) is the dominant non-ideal MHD effect. In simulations where the magnetic field and matter are moderately coupled, the disk remains relatively laminar with the radial electric current steepened by AD into a thin layer near the midplane. The toroidal magnetic field sharply reverses polarity in this layer, generating a large magnetic torque that drives fast accretion, which drags the poloidal field into a highly pinched radial configuration. The reconnection of this pinched field creates magnetic loops where the net poloidal magnetic flux (and thus the accretion rate) is reduced, yielding dense rings. Neighbouring regions with stronger poloidal magnetic fields accrete faster, forming gaps. In better magnetically coupled simulations, the so-called `avalanche accretion streams' develop continuously near the disk surface, rendering the disk-wind system more chaotic. Nevertheless, prominent rings and gaps are still produced, at least in part, by reconnection, which again enables the segregation of the poloidal field and the disk material similar to the more diffusive disks. However, the reconnection is now driven by the non-linear growth of MRI channel flows. The formation of rings and gaps in rapidly accreting yet laminar disks has interesting implications for dust settling and trapping, grain growth, and planet formation.

Efficient magnetic braking is a formidable obstacle to the formation of rotationally supported disks (RSDs) around protostars in magnetized dense cores. We have previously shown, through 2D (axisymmetric) non-ideal MHD simulations, that removing very small grains (VSGs: ~10 \AA$~$to few 100 \AA) can greatly enhance ambipolar diffusion and enable the formation of RSDs. Here we extend the simulations of disk formation enabled by VSG removal to 3D. We find that the key to this scenario of disk formation is that the drift velocity of the magnetic field almost cancels out the infall velocity of the neutrals in the $10^2$-$10^3$ AU-scale "pseudo-disk" where the field lines are most severely pinched and most of protostellar envelope mass infall occurs. As a result, the bulk neutral envelope matter can collapse without dragging much magnetic flux into the disk-forming region, which lowers the magnetic braking efficiency. We find that the initial disks enabled by VSG removal tend to be Toomre-unstable, which leads to the formation of prominent spiral structures that function as centrifugal barriers. The piling-up of infall material near the centrifugal barrier often produces dense fragments of tens of Jupiter masses, especially in cores that are not too strongly magnetized. Some fragments accrete onto the central stellar object, producing bursts in mass accretion rate. Others are longer lived, although whether they can survive long-term to produce multiple systems remains to be ascertained. Our results highlight the importance of dust grain evolution in determining the formation and properties of protostellar disks and potentially multiple systems.

We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.

Rings and gaps are being observed in an increasing number of disks around young stellar objects. We illustrate the formation of such radial structures through idealized, 2D (axisymmetric) resistive MHD simulations of coupled disk-wind systems threaded by a relatively weak poloidal magnetic field (plasma-$\beta \sim 10^3$). We find two distinct modes of accretion depending on the resistivity and field strength. A small resistivity or high field strength promotes the development of rapidly infalling `avalanche accretion streams' in a vertically extended disk envelope that dominates the dynamics of the system, especially the mass accretion. The streams are suppressed in simulations with larger resistivities or lower field strengths, where most of the accretion instead occurs through a laminar disk. In these simulations, the disk accretion is driven mainly by a slow wind that is typically accelerated by the pressure gradient from a predominantly toroidal magnetic field. Both wind-dominated and stream-dominated modes of accretion create prominent features in the surface density distribution of the disk, including rings and gaps, with a strong spatial variation of the magnetic flux relative to the mass. Regions with low mass-to-flux ratios accrete quickly, leading to the development of gaps, whereas regions with higher mass-to-flux ratios tend to accrete more slowly, allowing matter to accumulate and form dense rings. In some cases, avalanche accretion streams are observed to produce dense rings directly through continuous feeding. We discuss the implications of ring and gap formation driven by winds and streams on grain growth and planet formation.

We present ALMA 1.3 mm continuum, 12CO, C18O, and SO data for the Class 0 protostars, Lupus 3 MMS, IRAS 15398-3559, and IRAS 16253-2429 at resolutions of ~100 AU. By measuring a rotational profile in C18O, a 100 AU Keplerian disk around a 0.3 Msun protostar is observed in Lupus 3 MMS. No 100 AU Keplerian disks are observed in IRAS 15398-3559 and IRAS 16253-2429. Nevertheless, embedded compact (<30 AU) continuum components are detected. The C18O emission in IRAS 15398-3559 shows signatures of infall with a constant angular momentum. IRAS 16253-2429 exhibits signatures of infall and rotation, but its rotational profile is unresolved. By fitting the C18O data with our kinematic models, the protostellar masses and the disk radii are inferred to be 0.01 Msun and 20 AU in IRAS 15398-3559, and 0.03 Msun and 6 AU in IRAS 16253-2429. By comparing the specific angular momentum profiles from 10,000 to 100 AU in 8 Class 0 and I protostars, we find that the evolution of envelope rotation can be described with conventional inside-out collapse models. In comparison with a sample of 18 protostars with known disk radii, our results reveal signs of disk growth, with the disk radius increasing as M*^{0.8+/-0.14} or t^{1.09+\-0.37} in the Class 0 stage, where M* is the protostellar mass and t is the age. The disk growth rate slows down in the Class I stage. Besides, we find a hint that the mass accretion rate declines as t^{-0.26+\-0.04} from the Class 0 to I stages.

It has been shown that a realistic level of magnetization of dense molecular cloud cores can suppress the formation of a rotationally supported disk (RSD) through catastrophic magnetic braking in the axisymmetric ideal MHD limit. In this study, we present conditions for the formation of RSDs through non-ideal MHD effects computed self-consistently from an equilibrium chemical network. We find that removing from the standard MRN distribution the large population of very small grains (VSGs) of ~10 $\AA$ to few 100 $\AA$ that dominate the coupling of the bulk neutral matter to the magnetic field increases the ambipolar diffusivity by ~1--2 orders of magnitude at densities below 10$^{10}$ cm$^{-3}$. The enhanced ambipolar diffusion (AD) in the envelope reduces the amount of magnetic flux dragged by the collapse into the circumstellar disk-forming region. Therefore, magnetic braking is weakened and more angular momentum can be retained. With continuous high angular momentum inflow, RSDs of tens of AU are able to form, survive, and even grow in size, depending on other parameters including cosmic-ray ionization rate, magnetic field strength, and rotation speed. Some disks become self-gravitating and evolve into rings in our 2D (axisymmetric) simulations, which have the potential to fragment into (close) multiple systems in 3D. We conclude that disk formation in magnetized cores is highly sensitive to chemistry, especially to grain sizes. A moderate grain coagulation/growth to remove the large population of VSGs, either in the prestellar phase or during free-fall collapse, can greatly promote AD and help formation of tens of AU RSDs.

We explore kinematics and morphologies of molecular outflows driven by young protostars using magnetohydrodynamic simulations in the context of the unified wind model of Shang et al. The model explains the observed high-velocity jet and low-velocity shell features. In this work we investigate how these characteristics are affected by the underlying temperature and magnetic field strength. We study the problem of a warm wind running into a cold ambient toroid by using a tracer field that keeps track of the wind material. While an isothermal equation of state is adopted, the effective temperature is determined locally based on the wind mass fraction. In the unified wind model, the density of the wind is cylindrically stratified and highly concentrated toward the outflow axis. Our simulations show that for a sufficiently magnetized wind, the jet identity can be well maintained even at high temperatures. However, for a high temperature wind with low magnetization, the thermal pressure of the wind gas can drive material away from the axis, making the jet less collimated as it propagates. We also study the role of the poloidal magnetic field of the toroid. It is shown that the wind-ambient interface becomes more resistant to corrugation when the poloidal field is present, and the poloidal field that bunches up within the toroid prevents the swept-up material from being compressed into a thin layer. This suggests that the ambient poloidal field may play a role in producing a smoother and thicker swept-up shell structure in the molecular outflow.

We have conducted ALMA cycle 2 observations in the 1.3 mm continuum and in the C18O (2-1) and SO (5_6-4_5) lines at a resolution of ~0.3" toward the Class 0 protostar B335. The 1.3 mm continuum, C18O, and SO emission all show central compact components with sizes of ~40-180 AU within more extended components. The C18O component shows signs of infalling and rotational motion. By fitting simple kinematic models to the C18O data, the protostellar mass is estimated to be 0.05 Msun. The specific angular momentum, on a 100 AU scale, is ~4.3E-5 km/s*pc. A similar specific angular momentum, ~3E-5 to 5E-5 km/s*pc, is measured on a 10 AU scale from the velocity gradient observed in the central SO component, and there is no clear sign of an infalling motion in the SO emission. By comparing the infalling and rotational motion, our ALMA results suggest that the observed rotational motion has not yet reached Keplerian velocity neither on a 100 AU nor even on a 10 AU scale. Consequently, the radius of the Keplerian disk in B335 (if present) is expected to be 1-3 AU. The expected disk radius in B335 is one to two orders of magnitude smaller than those of observed Keplerian disks around other Class 0 protostars. Based on the observed infalling and rotational motion from 0.1 pc to inner 100 AU scales, there are two possible scenarios to explain the presence of such a small Keplerian disk in B335: magnetic braking and young age. If our finding is the consequence of magnetic braking, ~50% of the angular momentum of the infalling material within a 1000 AU scale might have been removed, and the magnetic field strength on a 1000 AU scale is estimated to be ~200 uG. If it is young age, the infalling radius in B335 is estimated to be ~2700 AU, corresponding to a collapsing time scale of ~5E4 yr.

A two dimensional hydrochemical hybrid code, KM2, is constructed to deal with astrophysical problems that would require coupled hydrodynamical and chemical evolution. The code assumes axisymmetry in cylindrical coordinate system, and consists of two modules: a hydrodynamics module and a chemistry module. The hydrodynamics module solves hydrodynamics using a Godunov-type finite volume scheme and treats included chemical species as passively advected scalars. The chemistry module implicitly solves non-equilibrium chemistry and change of the energy due to thermal processes with transfer of external ultraviolet radiation. Self-shielding effects on photodissociation of CO and H$_2$ are included. In this introductory paper, the adopted numerical method is presented, along with code verifications using the hydrodynamics modules, and a benchmark on the chemistry module with reactions specific to a photon-dominated region (PDR). Finally, as an example of the expected capability, the hydrochemical evolution of a PDR is presented based on the PDR benchmark.

The formation of rotationally supported protostellar disks is suppressed in ideal MHD in non-turbulent cores with aligned magnetic field and rotation axis. A promising way to resolve this so-called "magnetic braking catastrophe" is through turbulence. The reason for the turbulence-enabled disk formation is usually attributed to the turbulence-induced magnetic reconnection, which is thought to reduce the magnetic flux accumulated in the disk-forming region. We advance an alternative interpretation, based on magnetic decoupling-triggered reconnection of severely pinched field lines close to the central protostar and turbulence-induced warping of the pseudodisk of Galli and Shu. Such reconnection weakens the central split magnetic monopole that lies at the heart of the magnetic braking catastrophe under flux freezing. We show, through idealized numerical experiments, that the pseudodisk can be strongly warped, but not completely destroyed, by a subsonic or sonic turbulence. The warping decreases the rates of angular momentum removal from the pseudodisk by both magnetic torque and outflow, making it easier to form a rotationally supported disk. More importantly, the warping of the pseudodisk out of the disk-forming, equatorial plane greatly reduces the amount of magnetic flux threading the circumstellar, disk-forming region, further promoting disk formation. The beneficial effects of pseudodisk warping can also be achieved by a misalignment between the magnetic field and rotation axis. These two mechanisms of disk formation, enabled by turbulence and field-rotation misalignment respectively, are thus unified. We find that the disks formed in turbulent magnetized cores are rather thick and significantly magnetized. Implications of these findings, particularly for the thick young disk inferred in L1527, are briefly discussed.

HH 212 is a nearby (400 pc) Class 0 protostellar system showing several components that can be compared with theoretical models of core collapse. We have mapped it in 350 GHz continuum and HCO+ J=4-3 emission with ALMA at up to ~ 0.4" resolution. A flattened envelope and a compact disk are seen in continuum around the central source, as seen before. The HCO+ kinematics shows that the flattened envelope is infalling with small rotation (i.e., spiraling) into the central source, and thus can be identified as a pseudodisk in the models of magnetized core collapse. Also, the HCO+ kinematics shows that the disk is rotating and can be rotationally supported. In addition, to account for the missing HCO+ emission at low-redshifted velocity, an extended infalling envelope is required, with its material flowing roughly parallel to the jet axis toward the pseudodisk. This is expected if it is magnetized with an hourglass B-field morphology. We have modeled the continuum and HCO+ emission of the flattened envelope and disk simultaneously. We find that a jump in density is required across the interface between the pseudodisk and the disk. A jet is seen in HCO+ extending out to ~ 500 AU away from the central source, with the peaks upstream of those seen before in SiO. The broad velocity range and high HCO+ abundance indicate that the HCO+ emission traces internal shocks in the jet.

(Abridged) In this review we focus on the observations and theory of the formation of early disks and outflows, and their connections with the first phases of planet formation. Large rotationally supported circumstellar disks, although common around more evolved young stellar objects, are rarely detected during the earliest, "Class 0" phase; however, a few excellent candidates have been discovered recently around both low and high mass protostars. In this early phase, prominent outflows are ubiquitously observed; they are expected to be associated with at least small magnetized disks. Disk formation - once thought to be a simple consequence of the conservation of angular momentum during hydrodynamic core collapse - is far more subtle in magnetized gas. In this case, the rotation can be strongly magnetically braked. Indeed, both analytic arguments and numerical simulations have shown that disk formation is suppressed in the strict ideal magnetohydrodynamic (MHD) limit for the observed level of core magnetization. We review what is known about this "magnetic braking catastrophe", possible ways to resolve it, and the current status of early disk observations. Outflows are also intimately linked to disk formation; they are a natural product of magnetic fields and rotation and are important signposts of star formation. We review new developments on early outflow generation since PPV. The properties of early disks and outflows are a key component of planet formation in its early stages and we review these major connections.

Stars form in dense cores of molecular clouds that are observed to be significantly magnetized. In the simplest case of a laminar (non-turbulent) core with the magnetic field aligned with the rotation axis, both analytic considerations and numerical simulations have shown that the formation of a large, $10^2\au$-scale, rotationally supported protostellar disk is suppressed by magnetic braking in the ideal MHD limit for a realistic level of core magnetization. This theoretical difficulty in forming protostellar disks is termed "magnetic braking catastrophe". A possible resolution to this problem, proposed by \citeauthor{HennebelleCiardi2009} and \citeauthor{Joos+2012}, is that misalignment between the magnetic field and rotation axis may weaken the magnetic braking enough to enable disk formation. We evaluate this possibility quantitatively through numerical simulations. We confirm the basic result of \citeauthor{Joos+2012} that the misalignment is indeed conducive to disk formation. In relatively weakly magnetized cores with dimensionless mass-to-flux ratio $\gtrsim 5$, it enabled the formation of rotationally supported disks that would otherwise be suppressed if the magnetic field and rotation axis are aligned. For more strongly magnetized cores, disk formation remains suppressed, however, even for the maximum tilt angle of $90\degree$. If dense cores are as strongly magnetized as indicated by OH Zeeman observations (with a mean dimensionless mass-to-flux ratio $\sim 2$), it would be difficult for the misalignment alone to enable disk formation in the majority of them. We conclude that, while beneficial to disk formation, especially for the relatively weak field case, the misalignment does not completely solve the problem of catastrophic magnetic braking in general.

Magnetic flux redistribution lies at the heart of the problem of star formation in dense cores of molecular clouds that are magnetized to a realistic level. If all of the magnetic flux of a typical core were to be dragged into the central star, the stellar field strength would be orders of magnitude higher than the observed values. This well-known "magnetic flux problem" can in principle be resolved through non-ideal MHD effects. Two dimensional (axisymmetric) calculations have shown that ambipolar diffusion, in particular, can transport magnetic flux outward relative to matter, allowing material to enter the central object without dragging the field lines along. We show through simulations that such axisymmetric protostellar accretion flows are unstable in three dimensions to magnetic interchange instability in the azimuthal direction. The instability is driven by the magnetic flux redistributed from the matter that enters the central object. It typically starts to develop during the transition from the prestellar phase of star formation to the protostellar mass accretion phase. In the latter phase, the magnetic flux is transported outward mainly through advection, by strongly magnetized low-density regions that expand against the collapsing inflow. The tussle between the gravity-driven infall and magnetically driven expansion leads to a filamentary inner accretion flow, more disordered than previously pictured. The efficient outward transport of magnetic flux by advection lowers the field strength at small radii, making the magnetic braking less efficient and the formation of rotationally supported disks easier in principle. However, we find no evidence for such disks in any of our rotating collapse simulations. We conclude that the inner protostellar accretion flow is shaped to a large extent by this magnetic interchange instability. How disks form in such an environment is unclear.

Dense, star-forming, cores of molecular clouds are observed to be significantly magnetized. A realistic magnetic field of moderate strength has been shown to suppress, through catastrophic magnetic braking, the formation of a rotationally supported disk during the protostellar accretion phase of low-mass star formation in the ideal MHD limit. We address, through 2D (axisymmetric) simulations, the question of whether realistic levels of nonideal effects, computed with a simplified chemical network including dust grains, can weaken the magnetic braking enough to enable a rotationally supported disk to form. We find that ambipolar diffusion, the dominant nonideal MHD effect over most of the density range relevant to disk formation, does not enable disk formation, at least in 2D. The reason is that ambipolar diffusion allows the magnetic flux that would be dragged into the central stellar object in the ideal MHD limit to pile up instead in a small circumstellar region, where the magnetic field strength (and thus the braking efficiency) is greatly enhanced. We also find that, on the scale of tens of AU or more, a realistic level of Ohmic dissipation does not weaken the magnetic braking enough for a rotationally supported disk to form, either by itself or in combination with ambipolar diffusion. The Hall effect, the least explored of these three nonideal MHD effects, can spin up the material close to the central object to a significant, supersonic rotation speed, even when the core is initially non-rotating, although the spun-up material remains too sub-Keplerian to form a rotationally supported disk. The problem of catastrophic magnetic braking that prevents disk formation in dense cores magnetized to realistic levels remains unresolved. Possible resolutions of this problem are discussed.

Stars form in dense cores of magnetized molecular clouds. If the magnetic flux threading the cores is dragged into the stars, the stellar field would be orders of magnitude stronger than observed. This well-known "magnetic flux problem" demands that most of the core magnetic flux be decoupled from the matter that enters the star. We carry out the first exploration of what happens to the decoupled magnetic flux in 3D, using an MHD version of the ENZO adaptive mesh refinement code. The field-matter decoupling is achieved through a sink particle treatment, which is needed to follow the protostellar accretion phase of star formation. We find that the accumulation of the decoupled flux near the accreting protostar leads to a magnetic pressure buildup. The high pressure is released anisotropically, along the path of least resistance. It drives a low-density expanding region in which the decoupled magnetic flux is expelled. This decoupling-enabled magnetic structure has never been seen before in 3D MHD simulations of star formation. It generates a strong asymmetry in the protostellar accretion flow, potentially giving a kick to the star. In the presence of an initial core rotation, the structure presents an obstacle to the formation of a rotationally supported disk, in addition to magnetic braking, by acting as a rigid magnetic wall that prevents the rotating gas from completing a full orbit around the central object. We conclude that the decoupled magnetic flux from the stellar matter can strongly affect the protostellar collapse dynamics.

Stars form in dense cores of molecular clouds that are observed to be significantly magnetized. A dynamically important magnetic field presents a significant obstacle to the formation of protostellar disks. Recent studies have shown that magnetic braking is strong enough to suppress the formation of rotationally supported disks in the ideal MHD limit. Whether non-ideal MHD effects can enable disk formation remains unsettled. We carry out a first study on how disk formation in magnetic clouds is modified by the Hall effect, the least explored of the three non-ideal MHD effects in star formation (the other two being ambipolar diffusion and Ohmic dissipation). For illustrative purposes, we consider a simplified problem of a non-self-gravitating, magnetized envelope collapsing onto a central protostar of fixed mass. We find that the Hall effect can spin up the inner part of the collapsing flow to Keplerian speed, producing a rotationally supported disk. The disk is generated through a Hall-induced magnetic torque. Disk formation occurs even when the envelope is initially non-rotating, provided that the Hall coefficient is large enough. When the magnetic field orientation is flipped, the direction of disk rotation is reversed as well. The implication is that the Hall effect can in principle produce both regularly rotating and counter-rotating disks around protostars. We conclude that the Hall effect is an important factor to consider in studying the angular momentum evolution of magnetized star formation in general and disk formation in particular.

Disk formation in magnetized cloud cores is hindered by magnetic braking. Previous work has shown that for realistic levels of core magnetization, the magnetic field suppresses the formation of rotationally supported disks during the protostellar mass accretion phase of low-mass star formation both in the ideal MHD limit and in the presence of ambipolar diffusion for typical rates of cosmic ray ionization. Additional effects, such as ohmic dissipation, the Hall effect, and protostellar outflow, are needed to weaken the magnetic braking and enable the formation of persistent, rotationally supported, protostellar disks. In this paper, we first demonstrate that the classic microscopic resistivity is not large enough to enable disk formation by itself. We then experiment with a set of enhanced values for the resistivity in the range $\eta=10^{17}$--$10^{22}$ cm^2/s. We find that a value of order $10^{19}$ cm^2/s is needed to enable the formation of a 100 AU-scale Keplerian disk; the value depends somewhat on the degree of core magnetization. The required resistivity is a few orders of magnitude larger than the classic microscopic values. Whether it can be achieved naturally during protostellar collapse remains to be determined.

Outflows can be loaded and accelerated to high speeds along rapidly rotating, open magnetic field lines by centrifugal forces. Whether such magnetocentrifugally driven winds are stable is a longstanding theoretical problem. As a step towards addressing this problem, we perform the first large-scale 3D MHD simulations that extend to a distance $\sim 10^2$ times beyond the launching region, starting from steady 2D (axisymmetric) solutions. In an attempt to drive the wind unstable, we increase the mass loading on one half of the launching surface by a factor of $\sqrt{10}$, and reduce it by the same factor on the other half. The evolution of the perturbed wind is followed numerically. We find no evidence for any rapidly growing instability that could disrupt the wind during the launching and initial phase of propagation, even when the magnetic field of the magnetocentrifugal wind is toroidally dominated all the way to the launching surface. The strongly perturbed wind settles into a new steady state, with a highly asymmetric mass distribution. The distribution of magnetic field strength is, in contrast, much more symmetric. We discuss possible reasons for the apparent stability, including stabilization by an axial poloidal magnetic field, which is required to bend field lines away from the vertical direction and produce a magnetocentrifugal wind in the first place.

Dynamically significant magnetic fields are routinely observed in molecular clouds, with mass-to-flux ratio lambda = (2 pi sqrt{G}) (Sigma/B) ~ 1 (here Sigma is the total column density and B is the field strength). It is widely believed that ``subcritical'' clouds with lambda < 1 cannot collapse, based on virial arguments by Mestel and Spitzer and a linear stability analysis by Nakano and Nakamura. Here we confirm, using high resolution numerical models that begin with a strongly supersonic velocity dispersion, that this criterion is a fully nonlinear stability condition. All the high-resolution models with lambda <= 0.95 form ``Spitzer sheets'' but collapse no further. All models with lambda >= 1.02 collapse to the maximum numerically resolvable density. We also investigate other factors determining the collapse time for supercritical models. We show that there is a strong stochastic element in the collapse time: models that differ only in details of their initial conditions can have collapse times that vary by as much as a factor of 3. The collapse time cannot be determined from just the velocity dispersion; it depends also on its distribution. Finally, we discuss the astrophysical implications of our results.

We present the results of a series of time-dependent numerical simulations of cold, magnetocentrifugally launched winds from accretion disks. Our simulations span four and half decades of mass loading; in the context of a disk with a launching region from $0.1\AU$ to $1.0\AU$ around a $1\solarmass$ star and a field strength of about $20\gauss$ at the inner disk edge, this amounts to mass loss rates of $1\times 10^{-9}$ -- $3\times 10^{-5}\solarmassyear$ from each side of the disk. We find that the degree of collimation of the wind increases with mass loading; however even the ``lightest'' wind simulated is significantly collimated compared with the force-free magnetic configuration of the same magnetic flux distribution. The implication is that for flows from young stellar objects a radial field approximation is inappropriate. Surprisingly, the terminal velocity of the wind and the magnetic lever arm are still well-described by the analytical solutions for a radial field geometry. We also find that the isodensity contours and Alfv\'en surface are very nearly self-similar in mass loading. The wind becomes unsteady above some critical mass loading rate. For a small enough injection speed, we are able to obtain the first examples of a class of heavily-loaded magnetocentrifugal winds with magnetic fields completely dominated by the toroidal component all the way to the launching surface.

We follow numerically the time evolution of axisymmetric outflows driven magneto-centrifugally from the inner portion of accretion disks, from their launching surface to large, observable distances. Special attention is paid to the collimation of part of the outflow into a dense, narrow jet around the rotation axis, after a steady state has been reached. For parameters typical of T Tauri stars, we define a fiducial ``jet'' as outlined by the contour of constant density at 10^4 cm^{-3}. We find that the jet, so defined, appears nearly cylindrical well above the disk, in agreement with previous asymptotic analyses. Closer to the equatorial plane, the density contour can either bulge outwards or pinch inwards, depending on the conditions at the launching surface, particularly the mass flux distribution. We find that even though a dense, jet-like feature is always formed around the axis, there is no guarantee that the high-density axial jet would dominate the more tenuous, wide-angle part of the wind. Specifically, on the 100 AU scale, resolvable by HST and ground-based adaptive optics for nearby T Tauri winds, the fraction of the wind mass flux enclosed by the fiducial jet can vary substantially, again depending on the launching conditions. We show two examples in which the fraction is ~20% and ~45%. These dependences may provide a way to constrain the conditions at the launching surface, which are poorly known at present.

It is widely believed that T Tauri winds are driven magnetocentrifugally from accretion disks close to the central stars. The exact launching conditions are uncertain. We show that a general relation exists between the poloidal and toroidal velocity components of a magnetocentrifugal wind at large distances and the rotation rate of the launching surface, independent of the uncertain launching conditions. We discuss the physical basis of this relation and verify it using a set of numerically-determined large-scale wind solutions. Both velocity components are in principle measurable from spatially resolved spectra, as has been done for the extended low-velocity component (LVC) of the DG Tau wind by Bacciotti et al. For this particular source, we infer that the spatially resolved LVC originates from a region on the disk extending from $\sim 0.3$ to $\sim 4.0\AU$ from the star, which is consistent with, and a refinement over, the previous rough estimate of Bacciotti et al.

We present self-similar solutions that describe the gravitational collapse of rotating, isothermal, magnetic molecular-cloud cores, relevant to the formation of rotationally supported protostellar disks. This work focuses on the evolution after a point mass first forms at the center and generalizes previous results by Contopoulos, Ciolek, & Konigl that did not include rotation. Our model includes ambipolar diffusion and magnetic braking and allows us to examine the full range of expected behaviors and their dependence on the physical parameters. For typical parameter values, the inflow first passes through an ambipolar-diffusion shock (at a radius r_a, where the magnetic flux decouples from the matter), and later through a centrifugal shock at r=r_c, the outer edge of a rotationally supported disk of mass M_d. By the time (~10^5 yr) that the central mass M_c grows to ~1 M_\sun, r_a may be larger than 1000 AU, r_c larger than 100 AU, and M_d/M_c smaller than 10%. Disk properties are consistent with data on T Tauri systems, and our results imply that protostellar disks may well be Keplerian also during earlier phases. We show that the disk is likely to drive centrifugal outflows transporting angular momentum and mass, and we incorporate these effects into the model. We verify that gravitational torques and magnetorotational instability-induced turbulence typically do not play an important role in the angular momentum transport. We also present solutions for the limiting cases of fast rotation (where collapse results in a massive disk with such a large outer radius that it traps the ambipolar diffusion front) and strong braking (where no disk is formed and the collapse resembles that of a nonrotating core at small radii), as well as solutions for the rotational collapse of ideal-MHD and nonmagnetic model cores.

Protostellar jets and winds are probably driven magnetocentrifugally from the surface of accretion disks close to the central stellar objects. The exact launching conditions on the disk, such as the distributions of magnetic flux and mass ejection rate, are poorly unknown. They could be constrained from observations at large distances, provided that a robust model is available to link the observable properties of the jets and winds at the large distances to the conditions at the base of the flow. We discuss the difficulties in constructing such large-scale wind models, and describe a novel technique which enables us to numerically follow the acceleration and propagation of the wind from the disk surface to arbitrarily large distances and the collimation of part of the wind into a dense, narrow ``jet'' around the rotation axis. Special attention is paid to the shape of the jet and its mass flux relative to that of the whole wind. The mass flux ratio is a measure of the jet formation efficiency.

The magneto-centrifugal model for jet formation is studied by time-dependent simulations reaching steady state in a cold gas with negligible fluid pressure, in an axisymmetric geometry, using a modification of the Zeus3D code adapted to parallel computers. The number of boundary conditions imposed at the coronal base takes into account the existence of the fast and Alfvenic critical surfaces, avoiding over-determination of the flow. The size and shape of the computational box is chosen to include these critical surfaces, reducing the influence of the outer boundary conditions. As there is a region, near the origin, where the inclination of field lines to the axis is too small to drive a centrifugal wind, we inject a thin, axial jet, expected to form electromagnetically near black holes. Acceleration and collimation appear for wide generic conditions. A reference run is shown in detail, with a wind leaving the computational volume in the axial direction with a poloidal velocity equal to 4 times the poloidal Alfven speed, collimated inside 11 degrees. Finally, the critical surfaces, fieldlines, thrust, energy, torque and mass discharge of the outgoing wind are shown for simulations with various profiles of mass and magnetic flux at the base of the corona.