Search

A big question in the field of star and planet formation is the time at which substantial dust grain growth occurs. The observed properties of dust emission across different wavelength ranges have been used as an indication that millimeter-sized grains are already present in the envelopes of young protostars. However, this interpretation is in tension with results from coagulation simulations, which are not able to produce such large grains in these conditions. In this work, we show analytically that the production of millimeter-sized grains in protostellar envelopes is impossible under the standard assumptions about the coagulation process. We discuss several possibilities that may serve to explain the observed dust emission in the absence of in-situ grain growth to millimeter sizes., Accepted to ApJ

Massive, young stars are the main source of energy that maintains multiphase structure and turbulence in the interstellar medium (ISM), and without this "feedback" the star formation rate (SFR) would be much higher than is observed. Rapid energy loss in the ISM and efficient energy recovery by stellar feedback lead to co-regulation of SFRs and the ISM state. Realistic approaches to this problem should solve the dynamical evolution of the ISM, including star formation, and the input of feedback energy self-consistently and accurately. Here, we present the TIGRESS-NCR numerical framework, in which UV radiation, supernovae, cooling and heating processes, and gravitational collapse are modeled explicitly. We use an adaptive ray tracing method for UV radiation transfer from star clusters represented by sink particles, accounting for attenuation by dust and gas. We solve photon-driven chemical equations to determine the abundances of H (time-dependent) and C/O-bearing species (steady-state), which then set cooling and heating rates self-consistently. Applying these methods, we present high-resolution magnetohydrodynamics simulations of differentially rotating local galactic disks representing typical conditions of nearby star-forming galaxies. We analyze ISM properties and phase distributions and show good agreement with existing multiwavelength galactic observations. We measure midplane pressure components (turbulent, thermal, and magnetic) and the weight, demonstrating that vertical dynamical equilibrium holds. We quantify the ratios of pressure components to the SFR surface density, which we call the feedback yields. The TIGRESS-NCR framework will allow for a wide range of parameter exploration, including low metallicity system., Comment: ApJ accepted

Characterizing the physical conditions at disk scales in Class 0 sources is crucial for constraining the protostellar accretion process and the initial conditions for planet formation. We use ALMA 1.3 mm and 3 mm observations to investigate the physical conditions of the dust around the Class 0 binary IRAS 16293-2422 A (sep, 18 pages, 7 figures. Published in ApJL

We have complemented existing observations of HI absorption with new observations of HCO$^+$, C$_2$H, HCN, and HNC absorption from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the direction of 20 background radio continuum sources with $4^\circ \leq |b| \leq 81^\circ$ to constrain the atomic gas conditions that are suitable for the formation of diffuse molecular gas. We find that these molecular species form along sightlines where $A_V \gtrsim 0.25$, consistent with the threshold for the HI-to-H$_2$ transition at solar metallicity. Moreover, we find that molecular gas is associated only with structures that have an HI optical depth $> 0.1$, a spin temperature $< 80$ K, and a turbulent Mach number $\gtrsim 2$. We also identify a broad, faint component to the HCO$^+$ absorption in a majority of sightlines. Compared to the velocities where strong, narrow HCO$^+$ absorption is observed, the HI at these velocities has a lower cold neutral medium (CNM) fraction and negligible CO emission. The relative column densities and linewidths of the different molecular species observed here are similar to those observed in previous experiments over a range of Galactic latitudes, suggesting that gas in the solar neighborhood and gas in the Galactic plane are chemically similar. For a select sample of previously-observed sightlines, we show that the absorption line profiles of HCO$^+$, HCN, HNC, and C$_2$H are stable over periods of $\sim 3$ years and $\sim 25$ years, likely indicating that molecular gas structures in these directions are at least $\gtrsim 100$ AU in size, 31 pages, 12 figures, 8 tables. Accepted to ApJ

The physical properties of galactic halo gas have a profound impact on the life cycle of galaxies. As gas travels through a galactic halo, it undergoes dynamical interactions, influencing its impact on star formation and the chemical evolution of the galactic disk. In the Milky-Way halo, considerable effort has been made to understand the spatial distribution of neutral gas, which are mostly in the form of large complexes. However, the internal variations of their physical properties remains unclear. In this study, we investigate the thermal and dynamical state of the neutral gas in HVCs. High-resolution observations (1.'1) of the 21 cm line emission in the EN field of the DHIGLS HI survey are used to analyze the physical properties of the bright concentration C I B located at an edge of complex C. We use the Gaussian decomposition code ROHSA to model its multiphase content, and perform a power spectrum analysis to analyze its multi-scale structure. Physical properties of some 200 structures extracted using dendrograms are examined. We identify two distinct regions, one of which has a prominent protrusion extending from the edge of complex C that exhibits an ongoing phase transition from warm diffuse gas to cold dense gas and filaments. The scale at which the warm gas becomes unstable and undergoes a thermal condensation is about 15 pc, corresponding to a cooling time about 1.5 Myr. We find that a transition from subsonic to trans-sonic turbulence is associated with the thermal condensation. A large scale perspective of complex C suggests that hydrodynamic instabilities are involved in creating the structured concentration C I B and the phase transition therein. However, the details of the dynamical and thermal processes remain unclear and will require further investigation, through both observations and numerical simulations. (Shortened for arxiv), 32 pages, 31 figures, 7 tables. Accepted for publication in ApJ

Grain growth in protoplanetary disks is the first step towards planet formation. One of the most important pieces in the grain growth model is calculating the collisional velocity between two grains induced by the turbulent motion of the gas. The collisional velocities in previous works are obtained based on the assumption that the turbulence is hydrodynamic with the Kolmogorov power spectrum. However, realistic protoplanetary disks are magnetized, and turbulent motions can be induced by the magneto-rotational instabilities (MRI). To understand the impact of MRI turbulence on grain growth, we carry out 3D MHD simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. Figure 1 illustrates the magnetic field structure in one of the MRI simulations. We investigate two fundamental parameters of the turbulence that can have large impacts on grain collision velocities, the kinetic energy spectrum and the turbulence autocorrelation time. We find that the convergence of the turbulence α-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum with a -4/3 slope, shallower than the -5/3 slope in the Kolmogorov turbulence (Figure 2). The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop at small scales, with a slope close to -1, steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence can have a significant impact on the grain collisional velocities. Figure 3 shows the collisional velocities in Kolmogorov turbulence and turbulence observed in the simulations. For small grains that are well coupled to the gas, the collisional velocities can vary by orders of magnitude across different turbulence models. To understand the dust growth in protoplanetary disks, better computational techniques are needed to achieve the large dynamical range necessary for resolving the turbulence spectrum in the future. Reference: Gong M., Ivlev A. V., Zhao B., Caselli P., 2020, ApJ, 891, 172 https://ui.adsabs.harvard.edu/abs/2020ApJ...891..172G/abstract

CO is the most widely used observational tracer of molecular gas. The observable CO luminosity is translated to H_2 mass via a conversion factor, X_CO, which is a source of uncertainty and bias. Despite variations in X_CO, the empirically-determined solar neighborhood value is often applied across different galactic environments. To improve understanding of X_CO, we employ 3D magnetohydrodynamics simulations of the interstellar medium (ISM) in galactic disks with a large range of gas surface densities, allowing for varying metallicity, far-ultraviolet (FUV) radiation, and cosmic ray ionization rate (CRIR). With the TIGRESS simulation framework we model the three-phase ISM with self-consistent star formation and feedback, and post-process outputs with chemistry and radiation transfer to generate synthetic CO(1--0) and (2--1) maps. Our models reproduce the observed CO excitation temperatures, line-widths, and line ratios in nearby disk galaxies. X_CO decreases with increasing metallicity, with a power-law slope of -0.8 for the (1--0) line and -0.5 for the (2--1) line. X_CO also decreases at higher CRIR, and is insensitive to the FUV radiation. As density increases, X_CO first decreases due to increasing excitation temperature, and then increases when the emission is fully saturated. We provide fits between X_CO and observable quantities such as the line ratio, peak antenna temperature, and line brightness, which probe local gas conditions. These fits, which allow for varying beam size, may be used in observations to calibrate out systematic biases. We also provide estimates of the CO-dark H_2 fraction at different gas surface densities, observational sensitivities, and beam sizes., Comment: Accepted by ApJ

We compare observations of HI from the Very Large Array (VLA) and the Arecibo Observatory and observations of HCO$^+$ from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the diffuse ($A_V\lesssim1$) interstellar medium (ISM) to predictions from a photodissociation region (PDR) chemical model and multi-phase ISM simulations. Using a coarse grid of PDR models, we estimate the density, FUV radiation field, and cosmic ray ionization rate (CRIR) for each structure identified in HCO$^+$ and HI absorption. These structures fall into two categories. Structures with $T_s40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\gtrsim10^{12}~\mathrm{cm^{-2}}$, are consistent with high density, FUV radiation field, and CRIR models, characteristic of environments close to massive star formation. The latter are also found in directions with a significant fraction of thermally unstable HI. In at least one case, we rule out the PDR model parameters, suggesting that alternative mechanisms (e.g., non-equilibrium processes like turbulent dissipation and/or shocks) are required to explain the observed HCO$^+$ in this direction. Similarly, while our observations and simulations of the turbulent, multi-phase ISM agree that HCO$^+$ formation occurs along sightlines with $N(\mathrm{HI})\gtrsim10^{21}~\mathrm{cm^{-2}}$, the simulated data fail to explain HCO$^+$ column densities $\gtrsim\rm{few}\times10^{12}~\mathrm{cm^{-2}}$. Since a majority of our sightlines with HCO$^+$ had such high column densities, this likely indicates that non-equilibrium chemistry is important for these lines of sight., 15 pages, 3 figures, 1 table. Accepted to ApJ

We present a generic mechanism for the thermal damping of compressive waves in the interstellar medium (ISM), occurring due to radiative cooling. We solve for the dispersion relation of magnetosonic waves in a two-fluid (ion-neutral) system in which density- and temperature-dependent heating and cooling mechanisms are present. We use this dispersion relation, in addition to an analytic approximation for the nonlinear turbulent cascade, to model dissipation of weak magnetosonic turbulence. We show that in some ISM conditions, the cutoff wavelength for magnetosonic turbulence becomes tens to hundreds of times larger when the thermal damping is added to the regular ion-neutral damping. We also run numerical simulations which confirm that this effect has a dramatic impact on cascade of compressive wave modes., Comment: Accepted to ApJ

Turbulence is the dominant source of collisional velocities for grains with a wide range of sizes in protoplanetary disks. So far, only Kolmogorov turbulence has been considered for calculating grain collisional velocities, despite the evidence that turbulence in protoplanetary disks may be non-Kolmogorov. In this work, we present calculations of grain collisional velocities for arbitrary turbulence models characterized by power-law spectra and determined by three dimensionless parameters: the slope of the kinetic energy spectrum, the slope of the auto-correlation time, and the Reynolds number. The implications of our results are illustrated by numerical simulations of the grain size evolution for different turbulence models. We find that for the modeled cases of the Iroshnikov-Kraichnan turbulence and the turbulence induced by the magneto-rotational instabilities, collisional velocities of small grains are much larger than those for the standard Kolmogorov turbulence. This leads to faster grain coagulation in the outer regions of protoplanetary disks, resulting in rapid increase of dust opacity in mm-wavelength and possibly promoting planet formation in very young disks., Comment: Accepted by ApJ

Chondrules are silicate spheroids found in meteorites, serving as important fossil records of the early solar system. In order to form chondrules, chondrule precursors must be heated to temperatures much higher than the typical conditions in the current asteroid belt. One proposed mechanism for chondrule heating is the passage through bow shocks of highly eccentric planetesimals in the protoplanetary disk in the early solar system. However, it is difficult for planetesimals to gain and maintain such high eccentricities. In this paper, we present a new scenario in which planetesimals in the asteroid belt region are excited to high eccentricities by the Jovian sweeping secular resonance in a depleting disk, leading to efficient formation of chondrules. We study the orbital evolution of planetesimals in the disk using semi-analytic models and numerical simulations. We investigate the dependence of eccentricity excitation on the planetesimal's size as well as the physical environment, and calculate the probability for chondrule formation. We find that 50 - 2000 km planetesimals can obtain eccentricities larger than 0.6 and cause effective chondrule heating. Most chondrules form in high velocity shocks, in low density gas, and in the inner disk. The fraction of chondrule precursors which become chondrules is about 4 - 9 % between 1.5 - 3 AU. Our model implies that the disk depletion timescale is $\tau_\mathrm{dep}\approx 1~\mathrm{Myr}$, comparable to the age spread of chondrules; and that Jupiter formed before chondrules, no more than 0.7 Myr after the formation of the CAIs., Comment: Accepted for publication by ApJ

Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out numerical simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. We find that the convergence of the turbulence $\alpha$-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum of $k^{-4/3}$. The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop close to $k^{-1}$ at small scales, with a slope steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence with the injection scale at the disk scale height can potentially have a significant impact on the grain collisional velocities., Comment: Accepted by ApJ

CO(J=1-0) line emission is a widely used observational tracer of molecular gas, rendering essential the X_CO factor, which is applied to convert CO luminosity to H_2 mass. We use numerical simulations to study how X_CO depends on numerical resolution, non-steady-state chemistry, physical environment, and observational beam size. Our study employs 3D magnetohydrodynamics (MHD) simulations of galactic disks with solar neighborhood conditions, where star formation and the three-phase interstellar medium (ISM) are self-consistently regulated by gravity and stellar feedback. Synthetic CO maps are obtained by post-processing the MHD simulations with chemistry and radiation transfer. We find that CO is only an approximate tracer of H_2. On parsec scales, W_CO is more fundamentally a measure of mass-weighted volume density, rather than H_2 column density. Nevertheless, $\langle X_\mathrm{CO} \rangle=0.7-1.0\times10^{20}~\mathrm{cm^{-2}K^{-1}km^{-1}s}$ consistent with observations, insensitive to the evolutionary ISM state or radiation field strength if steady-state chemistry is assumed. Due to non-steady-state chemistry, younger molecular clouds have slightly lower X_CO and flatter profiles of X_CO versus extinction than older ones. The CO-dark H_2 fraction is 26-79 %, anti-correlated with the average extinction. As the observational beam size increases from 1 pc to 100 pc, X_CO increases by a factor of ~ 2. Under solar neighborhood conditions, X_CO in molecular clouds is converged at a numerical resolution of 2 pc. However, the total CO abundance and luminosity are not converged even at the numerical resolution of 1 pc. Our simulations successfully reproduce the observed variations of X_CO on parsec scales, as well as the dependence of X_CO on extinction and the CO excitation temperature., Comment: accepted by ApJ

We investigate prestellar core formation and accretion based on three-dimensional hydrodynamic simulations. Our simulations represent local $\sim 1$pc regions within giant molecular clouds where a supersonic turbulent flow converges, triggering star formation in the post-shock layer. We include turbulence and self-gravity, applying sink particle techniques, and explore a range of inflow Mach number ${\cal M}=2-16$. Two sets of cores are identified and compared: $t_1$-cores are identified of a time snapshot in each simulation, representing dense structures in a single cloud map; $t_\mathrm{coll}$-cores are identified at their individual time of collapse, representing the initial mass reservoir for accretion. We find that cores and filaments form and evolve at the same time. At the stage of core collapse, there is a well-defined, converged characteristic mass for isothermal fragmentation that is comparable to the critical Bonner-Ebert mass at the post-shock pressure. The core mass functions (CMFs) of $t_\mathrm{coll}$-cores show a deficit of high-mass cores ($\gtrsim 7M_\odot$) compared to the observed stellar initial mass function (IMF). However, the CMFs of $t_1$-cores are similar to the observed CMFs and include many low-mass cores that are gravitationally stable. The difference between $t_1$-cores and $t_\mathrm{coll}$-cores suggests that the full sample from observed CMFs may not evolve into protostars. Individual sink particles accrete at a roughly constant rate throughout the simulations, gaining one $t_\mathrm{coll}$-core mass per free-fall time even after the initial mass reservoir is accreted. High-mass sinks gain proportionally more mass at late times than low-mass sinks. There are outbursts in accretion rates, resulting from clumpy density structures falling into the sinks.

We conducted Herschel/PACS observations of five very low-mass stars or brown dwarfs located in the TW Hya association with the goal of characterizing the properties of disks in the low stellar mass regime. We detected all five targets at $70\,\mu{\rm{m}}$ and $100\,\mu{\rm{m}}$ and three targets at $160\,\mu{\rm{m}}$. Our observations, combined with previous photometry from 2MASS, WISE, and SCUBA-2, enabled us to construct SEDs with extended wavelength coverage. Using sophisticated radiative transfer models, we analyzed the observed SEDs of the five detected objects with a hybrid fitting strategy that combines the model grids and the simulated annealing algorithm and evaluated the constraints on the disk properties via the Bayesian inference method. The modelling suggests that disks around low-mass stars and brown dwarfs are generally flatter than their higher mass counterparts, but the range of disk mass extends to well below the value found in T Tauri stars, and the disk scale heights are comparable in both groups. The inferred disk properties (i.e., disk mass, flaring, and scale height) in the low stellar mass regime are consistent with previous findings from large samples of brown dwarfs and very low-mass stars. We discuss the dependence of disk properties on their host stellar parameters and find a significant correlation between the Herschel far-IR fluxes and the stellar effective temperatures, probably indicating that the scaling between the stellar and disk masses (i.e., $M_{\rm{disk}} \propto M_{\star}$) observed mainly in low-mass stars may extend down to the brown dwarf regime., Comment: Accepted for publication in A&A

We investigate the enrichment of the pre-solar cloud core with short lived radionuclides (SLRs), especially 26Al. The homogeneity and the surprisingly small spread in the ratio 26Al/27Al observed in the overwhelming majority of calcium-aluminium-rich inclusions (CAIs) in a vast variety of primitive chondritic meteorites places strong constraints on the formation of the the solar system. Freshly synthesized radioactive 26Al has to be included and well mixed within 20kyr. After discussing various scenarios including X-winds, AGB stars and Wolf-Rayet stars, we come to the conclusion that triggering the collapse of a cold cloud core by a nearby supernova is the most promising scenario. We then narrow down the vast parameter space by considering the pre-explosion survivability of such a clump as well as the cross-section necessary for sufficient enrichment. We employ numerical simulations to address the mixing of the radioactively enriched SN gas with the pre-existing gas and the forced collapse within 20kyr. We show that a cold clump of 10Msun at a distance of 5pc can be sufficiently enriched in 26Al and triggered into collapse fast enough - within 18kyr after encountering the supernova shock - for a range of different metallicities and progenitor masses, even if the enriched material is assumed to be distributed homogeneously in the entire supernova bubble. In summary, we envision an environment for the birth place of the Solar System 4.567Gyr ago similar to the situation of the pillars in M16 nowadays, where molecular cloud cores adjacent to an HII region will be hit by a supernova explosion in the future. We show that the triggered collapse and formation of the Solar System as well as the required enrichment with radioactive 26Al are possible in this scenario., 12 pages, 8 figures, accepted for publication in ApJ. Resolution of most figures degraded to fit within arXiv size limits. A full resolution version is available at http://www.usm.uni-muenchen.de/~gritschm/Gritschneder_2011_sun.pdf