-
Star Formation in a Strongly Magnetized Cloud
Authors:
Yu Cheng,
Jonathan C. Tan,
Paola Caselli,
Laura Fissel,
Hector G. Arce,
Francesco Fontani,
Matthew D. Goodson,
Mengyao Liu,
Nicholas Galitzki
Abstract:
We study star formation in the Center Ridge 1 (CR1) clump in the Vela C giant molecular cloud, selected as a high column density region that shows the lowest level of dust continuum polarization angle dispersion, likely indicating that the magnetic field is relatively strong. We observe the source with the ALMA 7m-array at 1.05~mm and 1.3~mm wavelengths, which enable measurements of dust temperatu…
▽ More
We study star formation in the Center Ridge 1 (CR1) clump in the Vela C giant molecular cloud, selected as a high column density region that shows the lowest level of dust continuum polarization angle dispersion, likely indicating that the magnetic field is relatively strong. We observe the source with the ALMA 7m-array at 1.05~mm and 1.3~mm wavelengths, which enable measurements of dust temperature, core mass and astrochemical deuteration. A relatively modest number of eleven dense cores are identified via their dust continuum emission, with masses spanning from 0.17 to 6.7 Msun. Overall CR1 has a relatively low compact dense gas fraction compared with other typical clouds with similar column densities, which may be a result of the strong magnetic field and/or the very early evolutionary stage of this region. The deuteration ratios, Dfrac, of the cores, measured with N2H+(3-2) and N2D+(3-2) lines, span from 0.011 to 0.85, with the latter being one of the highest values yet detected. The level of deuteration appears to decrease with evolution from prestellar to protostellar phase. A linear filament, running approximately parallel with the large scale magnetic field orientation, is seen connecting the two most massive cores, each having CO bipolar outflows aligned orthogonally to the filament. The filament contains the most deuterated core, likely to be prestellar and located midway between the protostars. The observations permit measurement of the full deuteration structure of the filament along its length, which we present. We also discuss the kinematics and dynamics of this structure, as well as of the dense core population.
△ Less
Submitted 25 May, 2021; v1 submitted 4 January, 2021;
originally announced January 2021.
-
Deuterium Chemodynamics of Massive Pre-Stellar Cores
Authors:
Chia-Jung Hsu,
Jonathan C. Tan,
Matthew D. Goodson,
Paola Caselli,
Bastian Körtgen,
Yu Cheng
Abstract:
High levels of deuterium fractionation of $\rm N_2H^+$ (i.e., $\rm D_{frac}^{N_2H^+} \gtrsim 0.1$) are often observed in pre-stellar cores (PSCs) and detection of $\rm N_2D^+$ is a promising method to identify elusive massive PSCs. However, the physical and chemical conditions required to reach such high levels of deuteration are still uncertain, as is the diagnostic utility of $\rm N_2H^+$ and…
▽ More
High levels of deuterium fractionation of $\rm N_2H^+$ (i.e., $\rm D_{frac}^{N_2H^+} \gtrsim 0.1$) are often observed in pre-stellar cores (PSCs) and detection of $\rm N_2D^+$ is a promising method to identify elusive massive PSCs. However, the physical and chemical conditions required to reach such high levels of deuteration are still uncertain, as is the diagnostic utility of $\rm N_2H^+$ and $\rm N_2D^+$ observations of PSCs. We perform 3D magnetohydrodynamics simulations of a massive, turbulent, magnetised PSC, coupled with a sophisticated deuteration astrochemical network. Although the core has some magnetic/turbulent support, it collapses under gravity in about one freefall time, which marks the end of the simulations. Our fiducial model achieves relatively low $\rm D_{frac}^{N_2H^+} \sim 0.002$ during this time. We then investigate effects of initial ortho-para ratio of $\rm H_2$ ($\rm OPR^{H_2}$), temperature, cosmic ray (CR) ionization rate, CO and N-species depletion factors and prior PSC chemical evolution. We find that high CR ionization rates and high depletion factors allow the simulated $\rm D_{frac}^{N_2H^+}$ and absolute abundances to match observational values within one freefall time. For $\rm OPR^{H_2}$, while a lower initial value helps the growth of $\rm D_{frac}^{N_2H^+}$, the spatial structure of deuteration is too widespread compared to observed systems. For an example model with elevated CR ionization rates and significant heavy element depletion, we then study the kinematic and dynamic properties of the core as traced by its $\rm N_2D^+$ emission. The core, undergoing quite rapid collapse, exhibits disturbed kinematics in its average velocity map. Still, because of magnetic support, the core often appears kinematically sub-virial based on its $\rm N_2D^+$ velocity dispersion.
△ Less
Submitted 5 January, 2021; v1 submitted 19 October, 2020;
originally announced October 2020.
-
A systematic comparison of two-equation RANS turbulence models applied to shock-cloud interactions
Authors:
Matthew D. Goodson,
Fabian Heitsch,
Karl Eklund,
Virginia A. Williams
Abstract:
Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixin…
▽ More
Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.
△ Less
Submitted 25 March, 2017;
originally announced March 2017.
-
Structure, Dynamics and Deuterium Fractionation of Massive Pre-Stellar Cores
Authors:
Matthew D. Goodson,
Shuo Kong,
Jonathan C. Tan,
Fabian Heitsch,
Paola Caselli
Abstract:
High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with…
▽ More
High levels of deuterium fraction in N$_2$H$^+$ are observed in some pre-stellar cores. Single-zone chemical models find that the timescale required to reach observed values ($D_{\rm frac}^{{\rm N}_2{\rm H}^+} \equiv {\rm N}_2{\rm D}^+/{\rm N}_2{\rm H}^+ \gtrsim 0.1$) is longer than the free-fall time, possibly ten times longer. Here, we explore the deuteration of turbulent, magnetized cores with 3D magnetohydrodynamics simulations. We use an approximate chemical model to follow the growth in abundances of N$_2$H$^+$ and N$_2$D$^+$. We then examine the dynamics of the core using each tracer for comparison to observations. We find that the velocity dispersion of the core as traced by N$_2$D$^+$ appears slightly sub-virial compared to predictions of the Turbulent Core Model of McKee & Tan, except at late times just before the onset of protostar formation. By varying the initial mass surface density, the magnetic energy, the chemical age, and the ortho-to-para ratio of H$_2$, we also determine the physical and temporal properties required for high deuteration. We find that low initial ortho-to-para ratios ($\lesssim 0.01$) and/or multiple free-fall times ($\gtrsim 3$) of prior chemical evolution are necessary to reach the observed values of deuterium fraction in pre-stellar cores.
△ Less
Submitted 30 October, 2016; v1 submitted 22 September, 2016;
originally announced September 2016.
-
Chemical enrichment of the pre-solar cloud by supernova dust grains
Authors:
Matthew D. Goodson,
Ian Luebbers,
Fabian Heitsch,
Christopher C. Frazer
Abstract:
The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternati…
▽ More
The presence of short-lived radioisotopes (SLRs) in solar system meteorites has been interpreted as evidence that the solar system was exposed to a supernova shortly before or during its formation. Yet results from hydrodynamical models of SLR injection into the proto-solar cloud or disc suggest that gas-phase mixing may not be efficient enough to reproduce the observed abundances. As an alternative, we explore the injection of SLRs via dust grains as a way to overcome the mixing barrier. We numerically model the interaction of a supernova remnant containing SLR-rich dust grains with a nearby molecular cloud. The dust grains are subject to drag forces and both thermal and non-thermal sputtering. We confirm that the expanding gas shell stalls upon impact with the dense cloud and that gas-phase SLR injection occurs slowly due to hydrodynamical instabilities at the cloud surface. In contrast, dust grains of sufficient size (> 1 micron) decouple from the gas and penetrate into the cloud within 0.1 Myr. Once inside the cloud, the dust grains are destroyed by sputtering, releasing SLRs and rapidly enriching the dense (potentially star-forming) regions. Our results suggest that SLR transport on dust grains is a viable mechanism to explain SLR enrichment.
△ Less
Submitted 26 August, 2016; v1 submitted 22 July, 2016;
originally announced July 2016.