We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

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.

Yen et al. (2012) advanced a direct approach for the calculation of self-gravitational force to second order accuracy based on uniform grid discretization. This method improves the accuracy of N-body calculation by using exact integration of kernel functions and employing the Fast Fourier Transform (FFT) to reduce complexity of computation to nearly linear. This direct approach is free of artificial boundary conditions, however, the applicability is limited by the uniform discretization of grids. We report here an advancement in the direct method with the implementation of adaptive mesh refinement (AMR) and maintaining second-order accuracy, which breaks the barrier set by uniform grid discretization. The adoption of graphic process units (GPUs) can significantly speed up the computation and make application of this method possible to astrophysical systems of gaseous disk galaxies and protoplanetary disks.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

The high-velocity molecular jet driven by Class 0 protostar IRAS 04166+2706 exhibits a unique saw-tooth velocity pattern. It consists of a series of well-aligned symmetric knots with similar averaged speeds, whose speeds at peaks of emission decreases roughly linearly away from the origin. Recent ALMA observations of knots R6 and B6 reveal kinematic behavior with expansion velocity increasing linearly from the axis to the edge. This pattern can be formed by a spherically expanding wind with axial density concentration. In this picture, the diverging velocity profile naturally possesses an increasing expansion velocity away from the axis, resulting in a tooth-like feature on the position-velocity diagram through projection. Such geometric picture predicts a correspondence between the slopes of the teeth and the outflow inclination angles, and the same inclination angle of 52$^\circ$ of the IRAS 04166+2706 can generally explain the whole pattern. Aided by numerical simulations in the framework of unified wind model by Shang et al. (2006), the observed velocity pattern can indeed be generated. A proper geometrical distribution of the jet and wind material is essential to the reconstruction the ejection history of the system.

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.

HH 211-mms is one of the youngest Class 0 protostellar systems in Perseus at ~ 235 pc away. We have mapped its central region at up to ~ 7 AU (0.03") resolution. A dusty disk is seen deeply embedded in a flattened envelope, with an intensity jump in dust continuum at ~ 350 GHz. It is nearly edge-on and is almost exactly perpendicular to the jet axis. It has a size of ~ 30 au along the major axis. It is geometrically thick, indicating that the (sub)millimeter light emitting grains have yet to settle to the midplane. Its inner part is expected to have transformed into a Keplerian rotating disk with a radius of ~ 10 au. A rotating disk atmosphere and a compact rotating bipolar outflow are detected in SO. The outflow fans out from the inner disk surfaces and is rotating in the same direction as the flattened envelope, and hence could trace a disk wind carrying away angular momentum from the inner disk. From the rotation of the disk atmosphere, the protostellar mass is estimated to be <~ 50 M_Jup. Together with results from the literature, our result favors a model where the disk radius grows linearly with the protostellar mass, as predicted by models of pre-stellar dense core evolution that asymptotes to an $r^{-1}$ radial profile for both the column density and angular velocity.

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.

HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there.

The central problem in forming a star is the angular momentum in the circumstellar disk which prevents material from falling into the central stellar core. An attractive solution to the "angular momentum problem" appears to be the ubiquitous (low-velocity and poorly-collimated) molecular outflows and (high-velocity and highly-collimated) protostellar jets accompanying the earliest phase of star formation that remove angular momentum at a range of disk radii. Previous observations suggested that outflowing material carries away the excess angular momentum via magneto-centrifugally driven winds from the surfaces of circumstellar disks down to ~ 10 AU scales, allowing the material in the outer disk to transport to the inner disk. Here we show that highly collimated protostellar jets remove the residual angular momenta at the ~ 0.05 AU scale, enabling the material in the innermost region of the disk to accrete toward the central protostar. This is supported by the rotation of the jet measured down to ~ 10 AU from the protostar in the HH 212 protostellar system. The measurement implies a jet launching radius of ~ 0.05_{-0.02}^{+0.05} AU on the disk, based on the magneto-centrifugal theory of jet production, which connects the properties of the jet measured at large distances to those at its base through energy and angular momentum conservation.

HH 212 is a nearby (400 pc) Class 0 protostellar system recently found to host a "hamburger"-shaped dusty disk with a radius of ~ 60 AU, deeply embedded in an infalling-rotating flattened envelope. We have spatially resolved this envelope-disk system with the Atacama Large Millimeter/submillimeter Array at up to ~ 16 AU (0.04") resolution. The envelope is detected in HCO+ J=4-3 down to the dusty disk. Complex organic molecules (COMs) and doubly deuterated formaldehyde (D2CO) are detected above and below the dusty disk within ~ 40 AU of the central protostar. The COMs are methanol (CH3OH), deuterated methanol (CH2DOH), methyl mercaptan (CH3SH), and formamide (NH2CHO, a prebiotic precursor). We have modeled the gas kinematics in HCO+ and COMs, and found a centrifugal barrier at a radius of ~ 44 AU, within which a Keplerian rotating disk is formed. This indicates that HCO+ traces the infalling-rotating envelope down to centrifugal barrier and COMs trace the atmosphere of a Keplerian rotating disk within the centrifugal barrier. The COMs are spatially resolved for the first time, both radially and vertically, in the atmosphere of a disk in the earliest, Class 0 phase of star formation. Our spatially resolved observations of COMs favor their formation in the disk rather than a rapidly infalling (warm) inner envelope. The abundances and spatial distributions of the COMs provide strong constraints on models of their formation and transport in low-mass star formation.

In the earliest (so-called "Class 0") phase of sunlike (low-mass) star formation, circumstellar disks are expected to form, feeding the protostars. However, such disks are difficult to resolve spatially because of their small sizes. Moreover, there are theoretical difficulties in producing such disks in the earliest phase, due to the retarding effects of magnetic fields on the rotating, collapsing material (so-called "magnetic braking"). With the Atacama Large Millimeter/submillimeter Array (ALMA), it becomes possible to uncover such disks and study them in detail. HH 212 is a very young protostellar system. With ALMA, we not only detect but also spatially resolve its disk in dust emission at submillimeter wavelength. The disk is nearly edge-on and has a radius of ~ 60 AU. Interestingly, it shows a prominent equatorial dark lane sandwiched between two brighter features, due to relatively low temperature and high optical depth near the disk midplane. For the first time, this dark lane is seen at submillimeter wavelength, producing a "hamburger"-shaped appearance that is reminiscent of the scattered-light image of an edge-on disk in optical and near infrared. Our observations open up an exciting possibility of directly detecting and characterizing small disks around the youngest protostars through high-resolution imaging with ALMA, which provides strong constraints on theories of disk formation.

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.

Forbidden neon emission from jets of low-mass young stars can be used to probe the underlying high-energy processes in these systems. We analyze spectra of the jet of DG Tau obtained with the Very Large Telescope/X-Shooter spectrograph in 2010. [Ne III] $\lambda$3869 is clearly detected in the innermost 3" microjet and the outer knot located at $\sim$6".5. The velocity structure of the inner microjet can be decomposed into the low-velocity component (LVC) at $\sim -70$ km/s and the high-velocity component (HVC) at $\sim -180$ km/s. Based on the observed [Ne III] flux and its spatial extent, we suggest the origins of the [Ne III] emission regions and their relation with known X-ray sources along the jet. The flares from the hard X-ray source close to the star may be the main ionization source of the innermost microjet. The fainter soft X-ray source at 0".2 from the star may provide sufficient heating to help to sustain the ionization fraction against the recombination in the flow. The outer knot may be reionized by shocks faster than 100 km/s such that [Ne III] emission reappears and that the soft X-ray emission at 5".5 is produced. Velocity decomposition of the archival Hubble Space Telescope spectra obtained in 1999 shows that the HVC had been faster, with a velocity centroid of $\sim -260$ km/s. Such a decrease in velocity may potentially be explained by the expansion of the stellar magnetosphere, changing the truncation radius and thus the launching speed of the jet. The energy released by magnetic reconnections during relaxation of the transition can heat the gas up to several tens of megakelvin and provide the explanation for on-source keV X-ray flares that ionize the neon microjet.

$Aims.$ We study the relation between the jet and the outflow in the IRAS 04166+2706 protostar. This Taurus protostar drives a molecular jet that contains multiple emission peaks symmetrically located from the central source. The protostar also drives a wide-angle outflow consisting of two conical shells. $Methods.$ We have used the Atacama Large Millimeter/submillimeter Array (ALMA) interferometer to observe two fields along the IRAS 04166+2706 jet. The fields were centered on a pair of emission peaks that correspond to the same ejection event, and were observed in CO(2-1), SiO(5-4), and SO(65-54). $ Results.$ Both ALMA fields present spatial distributions that are approximately elliptical and have their minor axes aligned with the jet direction. As the velocity increases, the emission in each field moves gradually across the elliptical region. This systematic pattern indicates that the emitting gas in each field lies in a disk-like structure that is perpendicular to the jet axis and is expanding away from the jet. A small degree of curvature in the first-moment maps indicates that the disks are slightly curved in the manner expected for bow shocks moving away from the IRAS source. A simple geometrical model confirms that this scenario fits the main emission features. $Conclusions.$ The emission peaks in the IRAS 04166+2706 jet likely represent internal bow shocks where material is being ejected laterally away from the jet axis. While the linear momentum of the ejected gas is dominated by the component in the jet direction, the sideways component is not negligible, and can potentially affect the distribution of gas in the surrounding outflow and core.

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.

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.

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.

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.

The regular satellites found around Neptune ($\approx 17~M_{\Earth}$) and Uranus ($\approx 14.5~M_{\Earth}$) suggest that past gaseous circumplanetary disks may have co-existed with solids around sub-Neptune-mass protoplanets ($< 17~M_{\Earth}$). These disks have been shown to be cool, optically thin, quiescent, with low surface density and low viscosity. Numerical studies of the formation are difficult and technically challenging. As an introductory attempt, three-dimensional global simulations are performed to explore the formation of circumplanetary disks around sub-Neptune-mass protoplanets embedded within an isothermal protoplanetary disk at the inviscid limit of the fluid in the absence of self-gravity. Under such conditions, a sub-Neptune-mass protoplanet can reasonably have a rotationally supported circumplanetary disk. The size of the circumplanetary disk is found to be roughly one-tenth of the corresponding Hill radius, which is consistent with the orbital radii of irregular satellites found for Uranus. The protoplanetary gas accretes onto the circumplanetary disk vertically from high altitude and returns to the protoplanetary disk again near the midplane. This implies an open system in which the circumplanetary disk constantly exchanges angular momentum and material with its surrounding prenatal protoplanetary gas.

Neon emission lines are good indicators of high-excitation regions close to a young stellar system because of their high ionization potentials and large critical densities. We have discovered [Ne III]{\lambda}3869 emission from the microjets of Sz 102, a low-mass young star in Lupus III. Spectroastrometric analyses of two-dimensional [Ne III] spectra obtained from archival high-dispersion ($R\approx 33,000$) Very Large Telescope/UVES data suggest that the emission consists of two velocity components spatially separated by ~ 0."3, or a projected distance of ~ 60 AU. The stronger redshifted component is centered at ~ +21 km/s with a line width of ~ 140 km/s, and the weaker blueshifted component at ~ -90 km/s with a line width of ~ 190 km/s. The two components trace velocity centroids of the known microjets and show large line widths that extend across the systemic velocity, suggesting their potential origins in wide-angle winds that may eventually collimate into jets. Optical line ratios indicate that the microjets are hot ($T\lesssim1.6\times10^4$ K) and ionized ($n_e\gtrsim5.7\times10^4$ cm$^{-3}$). The blueshifted component has ~ 13% higher temperature and ~ 46% higher electron density than the redshifted counterpart, forming a system of asymmetric pair of jets. The detection of the [Ne III]{\lambda}3869 line with the distinct velocity profile suggests that the emission originates in flows that may have been strongly ionized by deeply embedded hard X-ray sources, most likely generated by magnetic processes. The discovery of [Ne III]{\lambda}3869 emission along with other optical forbidden lines from Sz 102 support the picture of wide-angle winds surrounding magnetic loops in the close vicinity of the young star. Future high sensitivity X-ray imaging and high angular-resolution optical spectroscopy may help confirm the picture proposed.

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.

The molecular outflow from IRAS 04166+2706 was mapped with the Submillimeter Array (SMA) at 350 GHz continuum and CO J = 3$-$2 at an angular resolution of ~1 arcsec. The field of view covers the central arc-minute, which contains the inner four pairs of knots of the molecular jet. On the channel map, conical structures are clearly present in the low velocity range (|V$-$V$_0$|$<$10 km s$^{-1}$), and the highly collimated knots appear in the Extremely High Velocity range (EHV, 50$>$|V$-$V$_0$|$>$30 km s$^{-1}$). The higher angular resolution of ~1 arcsec reveals the first blue-shifted knot (B1) that was missing in previous PdBI observation of Sant\'iago-Garc\'ia et al. (2009) at an offset of ~6 arcsec to the North-East of the central source. This identification completes the symmetric sequence of knots in both the blue- and red-shifted lobes of the outflow. The innermost knots R1 and B1 have the highest velocities within the sequence. Although the general features appear to be similar to previous CO J = 2$-$1 images in Sant\'iago-Garc\'ia et al. (2009), the emission in CO J = 3$-$2 almost always peaks further away from the central source than that of CO J = 2$-$1 in the red-shifted lobe of the channel maps. This gives rise to a gradient in the line-ratio map of CO J = 3$-$2/J = 2$-$1 from head to tail within a knot. A large velocity gradient (LVG) analysis suggests that the differences may reflect a higher gas kinetic temperature at the head. We also explore possible constraints imposed by the non-detection of SiO J = 8$-$7.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

We present optical spectrophotometric monitoring of four active T Tauri stars (DG Tau, RY Tau, XZ Tau, RW Aur A) at high spectral resolution ($R \ga 1 \times 10^4$), to investigate the correlation between time variable mass ejection seen in the jet/wind structure of the driving source and time variable mass accretion probed by optical emission lines. This may allow us to constrain the understanding of the jet/wind launching mechanism, the location of the launching region, and the physical link with magnetospheric mass accretion. In 2010, observations were made at six different epochs to investigate how daily and monthly variability might affect such a study. We perform comparisons between the line profiles we observed and those in the literature over a period of decades and confirm the presence of time variability separate from the daily and monthly variability during our observations. This is so far consistent with the idea that these line profiles have a long term variability (3-20 years) related to episodic mass ejection suggested by the structures in the extended flow components. We also investigate the correlations between equivalent widths and between luminosities for different lines. We find that these correlations are consistent with the present paradigm of steady magnetospheric mass accretion and emission line regions that are close to the star.

The ejection of matter in the close vicinity of a young stellar object is investigated, treating the accretion disk as a gravitationally bound reservoir of matter. By solving the resistive MHD equations in 2D axisymmetry using our version of the Zeus-3D code with newly implemented resistivity, we study the effect of magnetic diffusivity in the magnetospheric accretion-ejection mechanism. Physical resistivity was included in the whole computational domain so that reconnection is enabled by the physical as well as the numerical resistivity. We show, for the first time, that quasi-stationary fast ejecta of matter, which we call {\em micro-ejections}, of small mass and angular momentum fluxes, can be launched from a purely resistive magnetosphere. They are produced by a combination of pressure gradient and magnetic forces, in presence of ongoing magnetic reconnection along the boundary layer between the star and the disk, where a current sheet is formed. Mass flux of micro-ejection increases with increasing magnetic field strength and stellar rotation rate, and is not dependent on the disk to corona density ratio and amount of resistivity.

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.

The bright-rimmed cloud SFO 22 was observed with the 45 m telescope of Nobeyama Radio Observatory in the ^{12}CO (J = 1-0), ^{13}CO (J = 1-0), and C^{18}O (J = 1-0) lines, where well-developed head-tail structure and small line widths were found. Such features were predicted by radiation-driven implosion models, suggesting that SFO 22 may be in a quasi-stationary equilibrium state. We compare the observed properties with those from numerical models of a photo-evaporating cloud, which include effects of magnetic pressure and heating due to strong far-ultraviolet (FUV) radiation from an exciting star. The magnetic pressure may play a more important role in the density structures of bright-rimmed clouds, than the thermal pressure that is enhanced by the FUV radiation. The FUV radiation can heat the cloud surface to near 30 K, however, its effect is not enough to reproduce the observed density structure of SFO 22. An initial magnetic field of 5 \mu G in our numerical models produces the best agreement with the observations, and its direction can affect the structures of bright-rimmed clouds.

We present a set of compelling science cases for the ALMA Band 1 receiver suite. For these cases, we assume in tandem the updated nominal Band 1 frequency range of 35-50 GHz with a likely extension up to 52 GHz; together these frequencies optimize the Band 1 science return. The scope of the science cases ranges from nearby stars to the re-ionization edge of the Universe. Two cases provide additional leverage on the present ALMA Level One Science Goals and are seen as particularly powerful motivations for building the Band 1 Receiver suite: (1) detailing the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires continuum observations out to ~35 GHz (~9mm); and (2) detecting CO 3-2 line emission from galaxies like the Milky Way during the epoch of re-ionization, i.e., 6 < z < 10, also requires Band 1 receiver coverage. The range of Band 1 science is wide, however, and includes studies of very small dust grains in the ISM, pulsar wind nebulae, radio supernovae, X-ray binaries, the Galactic Center (i.e., Sgr A*), dense cloud cores, complex carbon-chain molecules, masers, magnetic fields in the dense ISM, jets and outflows from young stars, distant galaxies, and galaxy clusters (i.e., the Sunyaev-Zel'dovich Effect). A comparison of ALMA and the Jansky VLA (JVLA) at the same frequencies of Band 1 finds similar sensitivity performance at 40-50 GHz, with a slight edge for ALMA at higher frequencies (e.g., within a factor of 2 for continuum observations). With its larger number of instantaneous baselines, however, ALMA Band 1data will have greater fidelity than those from the JVLA at similar frequencies.

In this paper, the RW Aur A microjet is studied from the point of view of X-wind models. The archived HST/STIS spectra of optical forbidden lines [O I], [S II], and [N II] from RW Aur A, taken in Cycle 8 with seven parallel slits along the jet axis, spaced at 0".07 apart, were analyzed. Images, position-velocity diagrams, and line ratios among the species were constructed, and compared with synthetic observations generated by selected solutions of the X-wind. Prominent features arising in a steady state X-wind could be identified within the convolved images, full-widths at half maxima and high-velocity peaks on both of the redshifted and blueshifted jets. The well-known asymmetric velocity profiles of the opposite jets are built into the selected models. We discuss model selections within the existing uncertainties of stellar parameters and inclination angle of the system. In this framework, the mass-loss rates that were inferred to be decreasing along the jet axis in the literature are the results of slowly decreasing excitation conditions and electron density profiles. Despite the apparent asymmetry in terminal velocities, line intensities and mass-loss rates, the average linear momenta from the opposite sides of the jet are actually balanced. These previously hard-to-explain features of the asymmetric RW Aur A jet system now find a different but self-consistent interpretation within the X-wind framework.

The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity can affect properties of the circumplanetary disk considerably when the mass of the protoplanet is $M_p \lesssim 33M_\oplus$, where $M_\oplus$ is the Earth mass. However, effects of viscosity on the circumplanetary disk are negligibly small when the mass of the protoplanet $M_p \gtrsim 33M_\oplus$. We find that when $M_p \lesssim 33M_\oplus$, viscosity can disrupt the spiral structure of the gas around the planet considerably and make the gas smoothly distributed, which makes the torques exerted on the protoplanet weaker. Thus, viscosity can make the migration speed of a protoplanet lower. After including viscosity, size of the circumplanetary disk can be decreased by a factor of $\gtrsim 20%$. Viscosity helps to transport gas into the circumplanetary disk from the differentially rotating circumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 50% after viscosity is taken into account when $M_p \lesssim 33M_\oplus$. Effects of viscosity on the formation of planets and satellites are briefly discussed.

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.

FS Tau B is one of the few T Tauri stars that possess a jet and a counterjet as well as an optically-visible cavity wall. We obtained images and spectra of its jet-cavity system in the near-infrared H and K bands using Subaru/IRCS and detected the jet and the counterjet in the [Fe II] 1.644 \mu m line for the first time. Within the inner 2" the blueshifted jet is brighter, whereas beyond ~ 5" the redshifted counterjet dominates the [Fe II] emission. The innermost blueshifted knot is spectrally resolved to have a large line width of ~ 110 km/s, while the innermost redshifted knot appears spectrally unresolved. The velocity ratio of the jet to the counterjet is ~ 1.34, which suggests that FS Tau B is driving an asymmetric jet, similar to those found in several T Tauri Stars. Combining with optical observations in the literature, we showed that the blueshifted jet has lower density and higher excitation than the redshifted counterjet. We suggest that the asymmetry in brightness and velocity is the manifestation of a bipolar outflow driving at different mass-loss rates, while maintaining balance of linear momentum. A full explanation to the asymmetry in the FS Tau B system awaits detail modeling and further investigation of the kinematic structure of the wind-associated cavity walls.

We perform numerical simulations in the close vicinity of a slowly rotating young stellar object. Using our own resistive MHD Zeus347 code in 2D axisymmetry, magnetospheric interaction experiences four robust stages in evolution. Quasi-stationary axial and conical streams of outflowing matter can last many orbital periods as results of the resistivity-facilitated magnetic reconnections. The shape of the magnetic field depends on resistivity in the magnetosphere.

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.

The protostellar jet driven by L1448C was observed in the SiO J=8-7 and CO J=3-2 lines and 350 GHz dust continuum at ~1" resolution with the Submillimeter Array (SMA). A narrow jet from the northern source L1448C(N) was observed in the SiO and the high-velocity CO. The jet consists of a chain of emission knots with an inter-knot spacing of ~2" (500 AU) and a semi-periodic velocity variation. The innermost pair of knots, which are significant in the SiO map but barely seen in the CO, are located at ~1" (250 AU) from the central source, L1448C(N). Since the dynamical time scale for the innermost pair is only ~10 yr, SiO may have been formed in the protostellar wind through the gas-phase reaction, or been formed on the dust grain and directly released into the gas phase by means of shocks. It is found that the jet is extremely active with a mechanical luminosity of ~7 L_sun, which is comparable to the bolometric luminosity of the central source (7.5 L_sun). The mass accretion rate onto the protostar derived from the mass-loss rate is ~10^{-5} M_sun/yr. Such a high mass accretion rate suggests that the mass and the age of the central star are 0.03-0.09 M_sun and (4-12)x10^3 yr, respectively, implying that the central star is in the very early stage of protostellar evolution. The low-velocity CO emission delineates two V-shaped shells with a common apex at L1448C(N). The kinematics of these shells are reproduced by the model of a wide opening angle wind. The co-existence of the highly-collimated jets and the wide-opening angle shells can be explained by the unified X-wind model" in which highly-collimated jet components correspond to the on-axis density enhancement of the wide-opening angle wind. The CO $J$=3--2 map also revealed the second outflow driven by the southern source L1448C(S) located at ~8.3" (2000 AU) from L1448C(N).

The flux and line shape of the fine-structure transitions of \NeII\ and \NeIII\ at 12.8 and 15.55\,$\mu$m and of the forbidden transitions of \OI\ $\lambda6300$ are calculated for young stellar objects with a range of mass-loss rates and X-ray luminosities using the X-wind model of jets and the associated wide-angle winds. For moderate and high accretion rates, the calculated \NeII\ line luminosity is comparable to or much larger than produced in X-ray irradiated disk models. All of the line luminosities correlate well with the main parameter in the X-wind model, the mass-loss rate, and also with the assumed X-ray luminosity --- and with one another. The line shapes of an approaching jet are broad and have strong blue-shifted peaks near the effective terminal velocity of the jet. They serve as a characteristic and testable aspect of jet production of the neon fine-structure lines and the \OI\ forbidden transitions.

HH 211 is a highly collimated jet originating from a nearby young Class 0 protostar. Here is a follow-up study of the jet with our previous observations at unprecedented resolution up to ~ 0.3" in SiO (J=8-7), CO (J=3-2), and SO (N_J=8_9-7_8). SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks in the jet. Although the emissions of these molecules show roughly the same morphology of the jet, there are detailed differences. In particular, the CO emission traces the jet closer to the source than the SiO and SO emissions. In addition, in the better resolved internal shocks, both the CO and SO emission are seen slightly ahead of the SiO emission. The jet is clearly seen on both sides of the source with more than one cycle of wiggle. The wiggle is reflection-symmetric about the source and can be reasonably fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself could be a very low-mass protobinary with a total mass of ~ 60 M_Jup and a binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase are found to be highly enhanced in the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the source where the dynamical time scale is <10 yrs. The abundance enhancements of these molecules are closely related to the internal shocks. The detected SiO is either the consequence of the release of Si-bearing material from dust grains or of its formation via gas chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the shocks.

Stellar magnetosphere and accretion disk interact, and a result should be outflow launched from the innermost vicinity of a protostellar object. We simulated physical conditions in this region by resistive MHD simulations. Outflows resembling the observed ones do not happen in the closest vicinity, except for quasi-stationary funnel flows onto the star, but could occur at few tens of stellar radii above the star. Numerical simulations we performed using our resistive version of ZEUS-3D code, ZEUS347.

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

Studying the earliest stages in the birth of stars is crucial for understanding how they form. Brown dwarfs with masses between that of stars and planets are not massive enough to maintain stable hydrogen-burning fusion reactions during most of their lifetime. Their origins are subject to much debate in recent literature because their masses are far below the typical mass where core collapse is expected to occur. We present the first confirmed evidence that brown dwarfs undergo a phase of molecular outflow that is typical of young stars. Using the Submillimeter Array, we have obtained a map of a bipolar molecular outflow from a young brown dwarf. We estimate an outflow mass of 1.6 x 10^-4 M_Sun and a mass-loss rate of 1.4 x 10^-9 M_Sun. These values are over two orders of magnitude smaller than the typical ones for T Tauri stars. From our millimiter continuum data and our own analysis of Spitzer infrared photometry, we estimate that the brown dwarf has a disk with a mass of 8 x 10^-3 M_Sun and an outer disk radius of 80 AU. Our results demonstrate that the bipolar molecular outflow operates down to planetary masses, occurring in brown dwarfs as a scaled-down version of the universal process seen in young stars.

The interaction of accretion disks with the magnetospheres of young stars can produce X-winds and funnel flows. With the assumption of axial symmetry and steady state flow, the problem can be formulated in terms of quantities that are conserved along streamlines, such as the Bernoulli integral (BI), plus a partial differential equation (PDE), called the Grad-Shafranov equation (GSE), that governs the distribution of streamlines in the meridional plane. The GSE plus BI yields a PDE of mixed type, elliptic before critical surfaces where the flow speed equals certain characteristic wave speeds are crossed and hyperbolic afterward. The computational difficulties are exacerbated by the locations of the critical surfaces not being known in advance. To overcome these obstacles, we consider a variational principle by which the GSE can be attacked by extremizing an action integral, with all other conserved quantities of the problem explicitly included as part of the overall formulation. To simplify actual applications we adopt the cold limit of a negligibly small ratio of the sound speed to the speed of Keplerian rotation in the disk where the X-wind is launched. We also ignore the obstructing effects of any magnetic fields that might thread a disk approximated to be infinitesimally thin. We then introduce trial functions with adjustable coefficients to minimize the variations that give the GSE. We tabulate the resulting coefficients so that other workers can have analytic forms to reconstruct X-wind solutions for various astronomical, cosmochemical, and meteoritical applications.

We have mapped the protostellar jet HH 211 in 342 GHz continuum, SiO ($J=8-7$), and CO ($J=3-2$) emission at $\sim$ \arcs{1} resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum at the center of the jet, tracing an envelope and a possible optically thick compact disk (with a size $<$ 130 AU) around the protostar. A knotty jet is seen in CO and SiO as in \H2{}, but extending closer to the protostar. It consists of a chain of knots on each side of the protostar, with an interknot spacing of $\sim$ \arcs{2}$-$\arcs{3} or 600$-$900 AU and the innermost pair of knots at only $\sim$ \arcsa{1}{7} or 535 AU from the protostar. These knots likely trace unresolved internal (bow) shocks (i.e., working surfaces) in the jet, with a velocity range up to $\sim$ 25 \vkm{}. The two-sided mass-loss rate of the jet is estimated to be $\sim (0.7-2.8)\times 10^{-6}$ \solarmass{} yr$^{-1}$. The jet is episodic, precessing, and bending. A velocity gradient is seen consistently across two bright SiO knots (BK3 and RK2) perpendicular to the jet axis, with $\sim$ 1.5$\pm$0.8 \vkm{} at $\sim$ 30$\pm$15 AU, suggesting a presence of a jet rotation. The launching radius of the jet, derived from the potential jet rotation, is $\sim$ 0.15$-$0.06 AU in the inner disk.

Context. Molecular clouds near the H II regions tend to harbor more luminous protostars. Aims. Our aim in this paper is to investigate whether or not radiation-driven implosion mechanism enhances luminosity of protostars near regions of high-ionizing fluxes. Methods. We performed numerical simulations to model collapse of cores exposed to UV radiation from O stars. We investigated dependence of mass loss rates on the initial density profiles of cores and variation of UV fluxes. We derived simple analytic estimates of accretion rates and final masses of protostars. Results. Radiation-driven implosion mechanism can increase accretion rates of protostars by 1-2 orders of magnitude. On the other hand, mass loss due to photo-evaporation is not large enough to have a significant impact on the luminosity. The increase of accretion rate makes luminosity 1-2 orders higher than those of protostars that form without external triggering. Conclusions. Radiation-driven implosion can help explain the observed higher luminosity of protostars in molecular clouds near H II regions.

HH 212 is a nearby (460 pc) protostellar jet discovered in H$_2$ powered by a Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. It is highly collimated and symmetric with matched pairs of bow shocks on either side of the source. We have mapped it in 850 $\mu$m continuum, SiO ($J=8-7$), CO ($J=3-2$), SO ($N_J=8_9-7_8$), HCO$^+$ ($J = 4-3$), and H$^{13}$CO$^+$ ($J = 4-3$) emission simultaneously at $\sim$ 1$''$ resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum around the source, mainly arising from an inner envelope (i.e., the inner part of a previously seen flattened envelope) and a possible disk. The inner envelope is also seen with rotation in CO, HCO$^+$, and probably SO. Like H$_2$ emission, CO and SiO emission are seen along the jet axis but extending closer to the source, tracing the bow shocks with a broad range of velocities and the continuous structures in between. SO emission is seen only around the source, forming a jetlike structure extending along the jet axis from the source, likely tracing the jet near the launching region. The jet is episodic and bending. It may also be slightly precessing as the jetlike SO structure shows a slight S-shaped symmetry about the source. A hint of jet rotation is also seen across the jet axis. Internal outflow shells are seen in CO and HCO$^+$, associated with the bow shocks in the inner part of the jet. The bases of the HCO$^+$ shells are seen with a hint of rotation similar to that seen in the inner envelope, probably consisted mainly of the material extended from the inner envelope and even the possible disk.

Two explanations exist for the short-lived radionuclides present in the solar system when the CAIs first formed. They originated either from the ejecta of a supernova or by the in situ irradiation of nebular dust by energetic particles. With a half-life of only 53 days, Beryllium-7 is then the key discriminant, since it can be made only by irradiation. We calculate the yield of Be-7. Within model uncertainties associated mainly with nuclear cross sections, we obtain agreement with the experimental value. Moreover, if Be-7 and Be-10 have the same origin, the irradiation time must be short. The x-wind model provides a natural astrophysical setting that gives the requisite conditions. The decoupling of the Al-26 and Be-10 observed in some rare CAIs receives a quantitative explanation when rare gradual events are considered. Finally, we show that the presence of supernova-produced Fe-60 in the solar accretion disk does not necessarily mean that other short-lived radionuclides have a stellar origin.

We observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.

We have mapped the SiO J=5-4 line at 217GHz from the HH211 molecular outflow with the Submillimeter Array (SMA). The high resolution map (1.6''x0.9'') shows that the SiO J=5-4 emission comes from the central narrow jet along the outflow axis with a width of ~0.8'' (~250 AU) FWHM. The SiO jet consists of a chain of knots separated by 3-4'' (~1000 AU) and most of the SiO knots have counterparts in shocked H_2 emission seen in a new, deep VLT near-infrared image of the outflow. A new, innermost pair of knots are discovered at just +/-2'' from the central star. The line ratio between the SiO J=5-4 data and upper limits from the SiO J=1-0 data of Chandler & Richer (2001) suggests that these knots have a temperature in excess of 300-500 K and a density of (0.5-1) x10^7 cm^{-3}. The radial velocity measured for these knots is ~30 km/s, comparable to the maximum velocity seen in the entire jet. The high temperature, high density, and velocity structure observed in this pair of SiO knots suggest that they are closely related to the primary jet launched close to the protostar.

HH 212 is a highly collimated jet discovered in H2 powered by a young Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. We have mapped around it in 1.33 mm continuum, 12CO ($J=2-1$), 13CO ($J=2-1$), C18O ($J=2-1$), and SO ($J_K = 6_5-5_4$) emission at $\sim$ \arcs{2.5} resolution with the Submillimeter Array. A dust core is seen in the continuum around the source. A flattened envelope is seen in C18O around the source in the equator perpendicular to the jet axis, with its inner part seen in 13CO. The structure and kinematics of the envelope can be roughly reproduced by a simple edge-on disk model with both infall and rotation. In this model, the density of the disk is assumed to have a power-law index of $p=-1.5$ or -2, as found in other low-mass envelopes. The envelope seems dynamically infalling toward the source with slow rotation because the kinematics is found to be roughly consistent with a free fall toward the source plus a rotation of a constant specific angular momentum. A 12CO outflow is seen surrounding the H2 jet, with a narrow waist around the source. Jetlike structures are also seen in 12CO near the source aligned with the H2 jet at high velocities. The morphological relationship between the H2 jet and the 12CO outflow, and the kinematics of the 12CO outflow along the jet axis are both consistent with those seen in a jet-driven bow shock model. SO emission is seen around the source and the H2 knotty shocks in the south, tracing shocked emission around them.

In January 2003, the Chandra Orion Ultradeep Project (COUP) detected about 1400 young stars during a 13.2 day observation of the Orion Nebula Cluster (ONC). This paper studies a well-defined sample of 28 solar-mass COUP sources to characterize the magnetic activity of analogs of the young Sun and thereby to improve understanding of the effects of solar X-rays on the solar nebula during the era of planet formation. We find that active young Suns spend 70% of their time in a characteristic state with relatively constant flux and magnetically confined plasma with temperatures kT_2 = 2.1 * kT_1. During characteristic periods, the 0.5-8 keV X-ray luminosity is about 0.03% of the bolometric luminosity. One or two powerful flares per week with peak luminosities logL_x ~ 30-32 erg/s are typically superposed on this characteristic emission accompanied by heating of the hot plasma component from ~2.4 keV to ~7 keV at the flare peak. The energy distribution of flares superposed on the characteristic emission level follows the relationship dN/dE ~ E^-1.7. The flare rates are consistent with the production of sufficiently energetic protons to spawn a spallogenic origin of some important short-lived radionuclides found in ancient meteorites. The X-rays can ionize gas in the circumstellar disk at a rate of 6 10^-9 ionizations per second at 1 AU from the central star, orders of magnitude above cosmic ray ionization rates. The estimated energetic particle fluences are sufficient to a account for many isotopic anomalies observed in meteoritic inclusions.

We have analyzed a number of intense X-ray flares observed in the Chandra Orion Ultradeep Project (COUP), a 13 days observation of the Orion Nebula Cluster (ONC). Analysis of the flare decay allows to determine the size, peak density and magnetic field of the flaring structure. A total of 32 events (the most powerful 1% of COUP flares), have sufficient statistics for the analysis. A broad range of decay times (from 10 to 400 ks) are present in the sample. Peak flare temperatures are often very high, with half of the flares in the sample showing temperatures in excess of 100 MK. Significant sustained heating is present in the majority of the flares. The magnetic structures which are found, are in a number of cases very long, with semi-lengths up to 10^12 cm, implying the presence of magnetic fields of hundreds of G extending to comparable distance from the stellar photosphere. These very large sizes for the flaring structures ($ >> R_*) are not found in more evolved stars, where, almost invariably, the same type of analysis results in structures with L <= R_*. As the majority of young stars in the ONC are surrounded by disks, we speculate that the large magnetic structures which confine the flaring plasma are actually the same type of structures which channel the plasma in the magnetospheric accretion paradigm, connecting the star's photosphere with the accretion disk.

We present the result of a systematic search for the iron Kalpha fluorescent line at 6.4 keV among 1616 X-ray sources detected by ultra-deep Chandra observations of the Orion Nebula Cluster and the obscured Orion Molecular Cloud 1 population as part of the Chandra Orion Ultra-deep Project (COUP). Seven sources are identified to have an excess emission at 6.4 keV among 127 control sample sources with significant counts in the 6.0-9.0 keV band. These seven sources are young stellar objects (YSOs) characterized by intense flare-like flux variations, thermal spectra, and near-infrared (NIR) counterparts. The observed equivalent widths of the line cannot be attributed to the fluorescence by interstellar or circumstellar matter along the line of sight. The X-ray spectral fits and NIR colors of the 6.4 keV sources show that these sources have X-ray absorption of > 1x10^22 cm^(-2) and NIR excess emission, which is not expected when the fluorescence occurs at the stellar photosphere. We therefore conclude that the iron fluorescent line of YSOs arises from reflection off of circumstellar disks, which are irradiated by the hard X-ray continuum emission of magnetic reconnection flares.

We have constructed the foundations to a series of theoretical diagnostic methods to probe the jet phenomenon in young stars as observed at various optical forbidden lines. We calculate and model in a self-consistent manner the physical and radiative processes which arise within an inner disk-wind driven magnetocentrifugally from the circumstellar accretion disk of a young sun-like star. Comparing with real data taken at high angular resolution, our approach will provide the basis of systematic diagnostics for jets and their related young stellar objects, to attest the emission mechanisms of such phenomena. This work can help bring first-principle theoretical predictions to confront actual multi-wavelength observations, and will bridge the link between many very sophiscated numerical simulations and observational data. Analysis methods discussed here are immediately applicable to new high-resolution data obtained with HST and Adaptic Optics.

In order to compare the x-wind with observations, one needs to be able to calculate its thermal and ionization properties. We formulate the physical basis for the streamline-by-streamline integration of the ionization and heat equations of the steady x-wind. In addition to the well-known processes associated with the interaction of stellar and accretion-funnel hot-spot radiation with the wind, we include X-ray heating and ionization, mechanical heating, and a revised calculation of ambipolar diffusion heating. The mechanical heating arises from fluctuations produced by star-disk interactions of the time dependent x-wind that are carried by the wind to large distances where they are dissipated in shocks, MHD waves, and turbulent cascades. We model the time-averaged heating by the scale-free volumetric heating rate, $\Gamma_{\rm mech} = \alpha \rho v^3 s^{-1}$, where $\rho$ and $v$ are the local mass density and wind speed, respectively, $s$ is the distance from the origin, and $\alpha$ is a phenomenological constant. When we consider a partially-revealed but active young stellar object, we find that choosing $\alpha \sim 10^{-3}$ in our numerical calculations produces temperatures and electron fractions that are high enough for the x-wind jet to radiate in the optical forbidden lines at the level and on the spatial scales that are observed. We also discuss a variety of applications of our thermal-chemical calculations that can lead to further observational checks of x-wind theory.