No Arabic abstract
We present 2.5-D global, ideal MHD simulations of magnetically and rotationally driven protostellar jets from Keplerian accretion discs, wherein only the initial magnetic field strength at the inner radius of the disc, $B_{rm i}$, is varied. Using the AMR-MHD code AZEUS, we self-consistently follow the jet evolution into the observational regime ($>10^3,mathrm{AU}$) with a spatial dynamic range of $sim6.5times10^5$. The simulations reveal a three-component outflow: 1) A hot, dense, super-fast and highly magnetised jet core; 2) a cold, rarefied, trans-fast and highly magnetised sheath surrounding the jet core and extending to a tangential discontinuity; and 3) a warm, dense, trans-slow and weakly magnetised shocked ambient medium entrained by the advancing bow shock. The simulations reveal power-law relationships between $B_{rm i}$ and the jet advance speed, $v_{rm jet}$, the average jet rotation speed, $langle v_varphirangle$, as well as fluxes of mass, momentum, and kinetic energy. Quantities that do not depend on $B_{rm i}$ include the plasma-$beta$ of the transported material which, in all cases, seems to asymptote to order unity. Jets are launched by a combination of the magnetic tower and bead-on-a-wire mechanisms, with the former accounting for most of the jet acceleration---even for strong fields---and continuing well beyond the fast magnetosonic point. At no time does the leading bow shock leave the domain and, as such, these simulations generate large-scale jets that reproduce many of the observed properties of protostellar jets including their characteristic speeds and transported fluxes.
We aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts Deltatau_mod, caused by dense clumps spiralling on the protostar, with the differences of dynamical timescales between the knots Deltatau_obs in CARMA 7. We found that the time spacing between the bursts have a bi-modal distribution caused by isolated and clustered luminosity bursts. The former are characterized by long quiescent periods between the bursts with Deltatau_mod = a few * (10^3-10^4) yr, whereas the latter occur in small groups with time spacing between the bursts Deltatau_mod= a few * (10-10^2) yr. For the clustered bursts, the distribution of Deltatau_mod in our models can be fit reasonably well to the distribution of Deltatau_obs in the protostellar jet of CARMA 7, if a certain correction for the (yet unknown) inclination angle with respect to the line of sight is applied. The K-S test on the model and observational data sets suggests the best-fit values for the inclination angles of 55-80 deg., which become narrower (75-80 deg.) if only strong luminosity bursts are considered. The dynamical timescales of the knots in the jet of CARMA 7 are too short for a meaningful comparison with the long time spacings between isolated bursts in our models. The exact sequences of time spacings between the luminosity bursts in our models and knots in the jet of CARMA 7 were found difficult to match. (abridged)
We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.
Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic position-velocity cuts are constructed from self-similar MHD DWs over a broad range of parameters, and three different methods are applied for estimating the specific angular momentum. We find that the launch radius inferred using the well-known relation from Anderson et al. (2006) can markedly differ from the true outermost launch radius $r_{out}$ of the DW. The double-peak separation and flow width methods provide only a strict lower limit to $r_{out}$. This bias is independent of angular resolution and can reach a factor ten. In contrast, the rotation curve method gives a good estimate of $r_{out}$ when the flow is well resolved, and an upper limit otherwise. The magnetic lever arm is always underestimated. Only comparison with synthetic predictions can take into account properly all observational effects. As an application, we present a comparison with ALMA observations of HH212 at resolutions from 250 au to 16 au, which represents the most stringent observational test of MHD DW to date. This comparison confirms our predicted biases for the double-peak separation method, and the large $r_{out}sim40~$au and small magnetic lever arm first suggested by Tabone et al. (2017). We also derive the first accurate analytical expression for the fraction of disk angular momentum extracted by an MHD disk wind of given radial extent, magnetic lever arm, and mass flux. Application to HH212 confirms that MHD DWs are serious candidates for the steady angular momentum extraction process in young disks.
(Abridged) The purpose of this paper is to explore and compare the physical and chemical structure of Class I low-mass protostellar sources on protoplanetary disc scales. We present a study of the dust and gas emission towards a representative sample of 12 Class I protostars from the Ophiuchus molecular cloud with the Atacama Large Millimeter/submillimeter Array (ALMA). The continuum at 0.87 mm and molecular transitions from C17O, C34S, H13CO+, CH3OH, SO2 , and C2H were observed at high angular resolution (0.4, ~60 au diameter) towards each source. Disc and stellar masses are estimated from the continuum flux and position-velocity diagrams, and six of the sources show disc-like structures. Towards the more luminous sources, compact emission and large line widths are seen for transitions of SO2 that probe warm gas (Eu ~200 K). In contrast, C17O emission is detected towards the least evolved and less luminous systems. No emission of CH3OH is detected towards any of the continuum peaks, indicating an absence of warm CH3OH gas towards these sources. A power-law relation is seen between the stellar mass and the bolometric luminosity, corresponding to a mass accretion rate of (2.4 +/- 0.6) x 10^-7 Msun/year for the Class I sources. This mass accretion rate is lower than the expected value if the accretion is constant in time and rather points to a scenario of accretion occurring in bursts. The differentiation between C17O and SO2 suggests that they trace different physical components: C17O traces the densest and colder regions of the disc-envelope system, while SO2 may be associated with regions of higher temperature, such as accretion shocks. The lack of warm CH3OH emission suggests that there is no hot-core-like region around any of the sources and that the CH3OH column density averaged over the disc is low.
X-ray observations of protostellar jets show evidence of strong shocks heating the plasma up to temperatures of a few million degrees. In some cases, the shocked features appear to be stationary. They are interpreted as shock diamonds. We aim at investigating the physics that guides the formation of X-ray emitting stationary shocks in protostellar jets, the role of the magnetic field in determining the location, stability, and detectability in X-rays of these shocks, and the physical properties of the shocked plasma. We performed a set of 2.5-dimensional magnetohydrodynamic numerical simulations modelling supersonic jets ramming into a magnetized medium and explored different configurations of the magnetic field. The model takes into account the most relevant physical effects, namely thermal conduction and radiative losses. We compared the model results with observations, via the emission measure and the X-ray luminosity synthesized from the simulations. Our model explains the formation of X-ray emitting stationary shocks in a natural way. The magnetic field collimates the plasma at the base of the jet and forms there a magnetic nozzle. After an initial transient, the nozzle leads to the formation of a shock diamond at its exit which is stationary over the time covered by the simulations (~ 40 - 60 yr; comparable with time scales of the observations). The shock generates a point-like X-ray source located close to the base of the jet with luminosity comparable with that inferred from X-ray observations of protostellar jets. For the range of parameters explored, the evolution of the post-shock plasma is dominated by the radiative cooling, whereas the thermal conduction slightly affects the structure of the shock.