-
Metallicity Dependence of Pressure-Regulated Feedback-Modulated Star Formation in the TIGRESS-NCR Simulation Suite
Authors:
Chang-Goo Kim,
Eve C. Ostriker,
Jeong-Gyu Kim,
Munan Gong,
Greg L. Bryan,
Drummond B. Fielding,
Sultan Hassan,
Matthew Ho,
Sarah M. R. Jeffreson,
Rachel S. Somerville,
Ulrich P. Steinwandel
Abstract:
We present a new simulation suite for the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are: (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating…
▽ More
We present a new simulation suite for the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are: (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including metallicity over $Z'\equiv Z/Z_\odot\sim0.1-3$, gas surface density $Σ_{\rm gas}\sim5-150 M_{\odot}{\rm pc^{-2}}$, and stellar surface density $Σ_{\rm star}\sim 1-50 M_{\odot}{\rm pc^{-2}}$. The range of emergent star formation rate surface density is $Σ_{\rm SFR}\sim 10^{-4}-0.5 M_{\odot}{\rm kpc^{-2}yr^{-1}}$ and ISM total midplane pressure is $P_{\rm tot}/k_B=10^3-10^6{\rm cm^{-3}K}$, with $P_{\rm tot}$ equal to the ISM weight $W$. For given $Σ_{\rm gas}$ and $Σ_{\rm star}$, we find $Σ_{\rm SFR} \propto Z'^{0.3}$. We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. We characterize feedback modulation in terms of the yield $Υ$, defined as the ratio of each stress to $Σ_{\rm SFR}$. The thermal feedback yield varies sensitively with both weight and metallicity as $Υ_{\rm th}\propto W^{-0.46}Z'^{-0.53}$, while the combined turbulent and magnetic feedback yield shows weaker dependence $Υ_{\rm turb+mag}\propto W^{-0.22}Z'^{-0.18}$. The reduction in $Σ_{\rm SFR}$ at low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.
△ Less
Submitted 6 June, 2024; v1 submitted 29 May, 2024;
originally announced May 2024.
-
Coagulation-Fragmentation Equilibrium for Charged Dust: Abundance of Submicron Grains Increases Dramatically in Protoplanetary Disks
Authors:
Vitaly Akimkin,
Alexei V. Ivlev,
Paola Caselli,
Munan Gong,
Kedron Silsbee
Abstract:
Dust coagulation in protoplanetary disks is not straightforward and is subject to several slow-down mechanisms, such as bouncing, fragmentation and radial drift to the star. Furthermore, dust grains in UV-shielded disk regions are negatively charged due to collisions with the surrounding electrons and ions, which leads to their electrostatic repulsion. For typical disk conditions, the relative vel…
▽ More
Dust coagulation in protoplanetary disks is not straightforward and is subject to several slow-down mechanisms, such as bouncing, fragmentation and radial drift to the star. Furthermore, dust grains in UV-shielded disk regions are negatively charged due to collisions with the surrounding electrons and ions, which leads to their electrostatic repulsion. For typical disk conditions, the relative velocities between micron-size grains are small and their collisions are strongly affected by the repulsion. On the other hand, collisions between pebble-size grains can be too energetic, leading to grain fragmentation. The aim of the present paper is to study a combined effect of the electrostatic and fragmentation barriers on dust evolution. We numerically solve the Smoluchowski coagulation-fragmentation equation for grains whose charging occurs under conditions typical for the inner disk regions, where thermal ionization operates. We find that dust fragmentation efficiently resupplies the population of small grains under the electrostatic barrier. As a result, the equilibrium abundance of sub-micron grains is enhanced by several orders of magnitude compared to the case of neutral dust. For some conditions with fragmentation velocities $\sim 1$ m s$^{-1}$, macroscopic grains are completely destroyed.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Implementation of chemistry in the Athena++ code
Authors:
Munan Gong,
Ka-Wai Ho,
James M. Stone,
Eve C. Ostriker,
Paola Caselli,
Tommaso Grassi,
Chang-Goo Kim,
Jeong-Gyu Kim,
Goni Halevi
Abstract:
Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating…
▽ More
Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magneto-hydrodynamic code Athena++. We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing the user to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ODEs with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photo-dissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++, aiming to provide a robust and flexible code for the astrochemical simulation community.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
The Physical Drivers and Observational Tracers of CO-to-H2 Conversion Factor Variations in Nearby Barred Galaxy Centers
Authors:
Yu-Hsuan Teng,
Karin M. Sandstrom,
Jiayi Sun,
Munan Gong,
Alberto D. Bolatto,
I-Da Chiang,
Adam K. Leroy,
Antonio Usero,
Simon C. O. Glover,
Ralf S. Klessen,
Daizhong Liu,
Miguel Querejeta,
Eva Schinnerer,
Frank Bigiel,
Yixian Cao,
Melanie Chevance,
Cosima Eibensteiner,
Kathryn Grasha,
Frank P. Israel,
Eric J. Murphy,
Lukas Neumann,
Hsi-An Pan,
Francesca Pinna,
Mattia C. Sormani,
J. D. T. Smith
, et al. (2 additional authors not shown)
Abstract:
The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NG…
▽ More
The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $α_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $α_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $α_\rm{CO}$ with low-J $^{12}$CO optical depths ($τ_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $τ_\rm{CO}$ correlation explains most of the $α_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $α_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $τ_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $α_\rm{CO}$ variation. We also test current simulation-based $α_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations.
△ Less
Submitted 10 April, 2023;
originally announced April 2023.
-
Dust hot spots at 10 au scales around the Class 0 binary IRAS 16293-2422 A: a departure from the passive irradiation model
Authors:
M. J. Maureira,
M. Gong,
J. E. Pineda,
H. B. Liu,
K. Silsbee,
P. Caselli,
J. Zamponi,
D. Segura-Cox,
A. Schmiedeke
Abstract:
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 <100 au) down to ~10 au scales. The circumbinary material's spectral index…
▽ More
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 <100 au) down to ~10 au scales. The circumbinary material's spectral index, alpha, has a median of 3.1 and a dispersion of ~0.2, providing no firm evidence of mm-sizes grains therein. Continuum substructures with brightness temperature peaks of T_b~60-80 K at 1.3 mm are observed near the disks at both wavelengths. These peaks do not overlap with strong variations of alpha, indicating they trace high-temperature spots instead of regions with significant optical depth variations. The lower limits to the inferred dust temperature in the hot spots are 122, 87 and 49 K. Depending on the assumed dust opacity index, these values can be several times higher. They overlap with high gas temperatures and enhanced complex organic molecular (COM) emission. This newly resolved dust temperature distribution is in better agreement with the expectations from mechanical instead of the most commonly assumed radiative heating. In particular, we find that the temperatures agree with shock heating predictions. This evidence and recent studies highlighting accretion heating in Class 0 disks suggest that mechanical heating (shocks, dissipation powered by accretion, etc.) is important during the early stages and should be considered when modeling and measuring properties of deeply embedded protostars and disks.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Introducing TIGRESS-NCR: I. Co-Regulation of the Multiphase Interstellar Medium and Star Formation Rates
Authors:
Chang-Goo Kim,
Jeong-Gyu Kim,
Munan Gong,
Eve C. Ostriker
Abstract:
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 sh…
▽ More
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.
△ Less
Submitted 27 March, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations
Authors:
Jeong-Gyu Kim,
Munan Gong,
Chang-Goo Kim,
Eve C. Ostriker
Abstract:
We present an efficient heating/cooling method coupled with chemistry and ultraviolet (UV) radiative transfer, which can be applied to numerical simulations of the interstellar medium (ISM). We follow the time-dependent evolution of hydrogen species (H$_2$, H, H$^+$), assume carbon/oxygen species (C, C$^+$, CO, O, and O$^+$) are in formation-destruction balance given the non-steady hydrogen abunda…
▽ More
We present an efficient heating/cooling method coupled with chemistry and ultraviolet (UV) radiative transfer, which can be applied to numerical simulations of the interstellar medium (ISM). We follow the time-dependent evolution of hydrogen species (H$_2$, H, H$^+$), assume carbon/oxygen species (C, C$^+$, CO, O, and O$^+$) are in formation-destruction balance given the non-steady hydrogen abundances, and include essential heating/cooling processes needed to capture thermodynamics of all ISM phases. UV radiation from discrete point sources and the diffuse background is followed through adaptive ray tracing and a six-ray approximation, respectively, allowing for H$_2$ self-shielding; cosmic ray (CR) heating and ionization are also included. To validate our methods and demonstrate their application for a range of density, metallicity, and radiation field, we conduct a series of tests, including the equilibrium curves of thermal pressure vs. density, the chemical and thermal structure in photo-dissociation regions, H I-to-H$_2$ transitions, and the expansion of H II regions and radiative supernova remnants. Careful treatment of photochemistry and CR ionization is essential for many aspects of ISM physics, including identifying the thermal pressure at which cold and warm neutral phases co-exist. We caution that many current heating and cooling treatments used in galaxy formation simulations do not reproduce the correct thermal pressure and ionization fraction in the neutral ISM. Our new model is implemented in the MHD code Athena and incorporated in the TIGRESS simulation framework, for use in studying the star-forming ISM in a wide range of environments.
△ Less
Submitted 14 October, 2022;
originally announced October 2022.
-
Dust grains cannot grow to millimeter sizes in protostellar envelopes
Authors:
Kedron Silsbee,
Vitaly Akimkin,
Alexei V. Ivlev,
Leonardo Testi,
Munan Gong,
Paola Caselli
Abstract:
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, whi…
▽ More
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.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
The role of neutral hydrogen in setting the abundances of molecular species in the Milky Way's diffuse interstellar medium. II. Comparison between observations and theoretical models
Authors:
Daniel R. Rybarczyk,
Munan Gong,
Snezana Stanimirovic,
Brian Babler,
Claire E. Murray,
Jan Martin Winters,
Gan Luo,
T. M. Dame,
Lucille Steffes
Abstract:
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 coa…
▽ More
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_s<40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\lesssim10^{12}~\mathrm{cm^{-2}}$, are consistent with modest density, FUV radiation field, and CRIR models, typical of the diffuse molecular ISM. Structures with spin temperature $T_s>40~\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.
△ Less
Submitted 10 December, 2021;
originally announced December 2021.
-
The role of neutral hydrogen in setting the abundances of molecular species in the Milky Way's diffuse interstellar medium. I. Observational constraints from ALMA and NOEMA
Authors:
Daniel R. Rybarczyk,
Snezana Stanimirovic,
Munan Gong,
Brian Babler,
Claire E. Murray,
Maryvonne Gerin,
Jan Martin Winters,
Gan Luo,
T. M. Dame,
Lucille Steffes
Abstract:
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 f…
▽ More
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
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
Resolving the formation of cold HI filaments in the high velocity cloud complex C
Authors:
Antoine Marchal,
Peter G. Martin,
Munan Gong
Abstract:
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 th…
▽ More
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)
△ Less
Submitted 23 June, 2021;
originally announced June 2021.
-
Impact of magneto-rotational instability on grain growth in protoplanetary disks: II. Increased grain collisional velocities
Authors:
Munan Gong,
Alexei V. Ivlev,
Vitaly Akimkin,
Paola Caselli
Abstract:
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 turbulen…
▽ More
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.
△ Less
Submitted 17 June, 2021;
originally announced June 2021.
-
Thermal damping of Weak Magnetosonic Turbulence in the Interstellar Medium
Authors:
Kedron Silsbee,
Alexei V. Ivlev,
Munan Gong
Abstract:
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 approximati…
▽ More
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.
△ Less
Submitted 11 November, 2021; v1 submitted 27 November, 2020;
originally announced November 2020.
-
The environmental dependence of the X_CO conversion factor
Authors:
Munan Gong,
Eve C. Ostriker,
Chang-Goo Kim,
Jeong-Gyu Kim
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…
▽ More
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.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
Impact of magneto-rotational instability on grain growth in protoplanetary disks: I. Relevant turbulence properties
Authors:
Munan Gong,
Alexei V. Ivlev,
Bo Zhao,
Paola Caselli
Abstract:
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 nume…
▽ More
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 $α$-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.
△ Less
Submitted 12 February, 2020;
originally announced February 2020.
-
Chondrule Formation by the Jovian Sweeping Secular Resonance
Authors:
Munan Gong,
Xiaochen Zheng,
Douglas N. C. Lin,
Kedron Silsbee,
Clement Baruteau,
Shude Mao
Abstract:
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 protoplanet…
▽ More
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 $τ_\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.
△ Less
Submitted 26 August, 2019;
originally announced August 2019.
-
The X_CO conversion factor from galactic multiphase ISM simulations
Authors:
Munan Gong,
Eve C. Ostriker,
Chang-Goo Kim
Abstract:
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 galact…
▽ More
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.
△ Less
Submitted 26 March, 2018;
originally announced March 2018.
-
A Simple and Accurate Network for Hydrogen and Carbon Chemistry in the ISM
Authors:
Munan Gong,
Eve C. Ostriker,
Mark G. Wolfire
Abstract:
Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in nume…
▽ More
Chemistry plays an important role in the interstellar medium (ISM), regulating heating and cooling of the gas, and determining abundances of molecular species that trace gas properties in observations. Although solving the time-dependent equations is necessary for accurate abundances and temperature in the dynamic ISM, a full chemical network is too computationally expensive to incorporate in numerical simulations. In this paper, we propose a new simplified chemical network for hydrogen and carbon chemistry in the atomic and molecular ISM. We compare results from our chemical network in detail with results from a full photo-dissociation region (PDR) code, and also with the Nelson & Langer (1999) (NL99) network previously adopted in the simulation literature. We show that our chemical network gives similar results to the PDR code in the equilibrium abundances of all species over a wide range of densities, temperature, and metallicities, whereas the NL99 network shows significant disagreement. Applying our network in 1D models, we find that the $\mathrm{CO}$-dominated regime delimits the coldest gas and that the corresponding temperature tracks the cosmic ray ionization rate in molecular clouds. We provide a simple fit for the locus of $\mathrm{CO}$ dominated regions as a function of gas density and column. We also compare with observations of diffuse and translucent clouds. We find that the $\mathrm{CO}$, $\mathrm{CHx}$ and $\mathrm{OHx}$ abundances are consistent with equilibrium predictions for densities $n=100-1000~\mathrm{cm^{-3}}$, but the predicted equilibrium $\mathrm{C}$ abundance is higher than observations, signaling the potential importance of non-equilibrium/dynamical effects.
△ Less
Submitted 4 April, 2019; v1 submitted 27 October, 2016;
originally announced October 2016.
-
Prestellar Core Formation, Evolution, and Accretion from Gravitational Fragmentation in Turbulent Converging Flows
Authors:
Munan Gong,
Eve C. Ostriker
Abstract:
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 numbe…
▽ More
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.
△ Less
Submitted 8 April, 2015;
originally announced April 2015.
-
Herschel/PACS view of disks around low-mass stars and brown dwarfs in the TW Hya association
Authors:
Yao Liu,
Gregory J. Herczeg,
Munan Gong,
Katelyn N. Allers,
Joanna M. Brown,
Adam L. Kraus,
Michael C. Liu,
Evgenya L. Shkolnik,
Ewine F. van Dishoeck
Abstract:
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\,μ{\rm{m}}$ and $100\,μ{\rm{m}}$ and three targets at $160\,μ{\rm{m}}$. Our observations, combined with previous photometry from 2MASS, WISE, and SCUBA-2, ena…
▽ More
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\,μ{\rm{m}}$ and $100\,μ{\rm{m}}$ and three targets at $160\,μ{\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.
△ Less
Submitted 7 November, 2014;
originally announced November 2014.
-
The Supernova Triggered Formation and Enrichment of Our Solar System
Authors:
M. Gritschneder,
D. N. C. Lin,
S. D. Murray,
Q. -Z. Yin,
M. -N. Gong
Abstract:
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 sy…
▽ More
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.
△ Less
Submitted 31 October, 2011;
originally announced November 2011.