Photochemistry and Heating/Cooling of the Multiphase Interstellar Medium with UV Radiative Transfer for Magnetohydrodynamic Simulations

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 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.

57 pages, 22 figures; accepted for publication in ApJS

References in corpus (74)

• Theory of Star Formation
• Athena: A New Code for Astrophysical MHD
• The effect of photo-ionization on the cooling rates of enriched, astrophysical plasmas
• The Effects of Stellar Rotation. II. A Comprehensive Set of Starburst99 Models
• Momentum Injection by Supernovae in the Interstellar Medium
• Time-Dependent Ionization in Radiatively Cooling Gas
• Escape fraction of ionizing photons during reionization: effects due to supernova feedback and runaway OB stars
• Grackle: a Chemistry and Cooling Library for Astrophysics
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• Turbulent Structure of a Stratified Supernova-Driven Interstellar Medium
• Simulating the formation of molecular clouds. II. Rapid formation from turbulent initial conditions
• H3+ in Diffuse Interstellar Clouds: a Tracer for the Cosmic-Ray Ionization Rate
• The Dynamics of Radiation Pressure-Dominated HII Regions
• Herschel Survey of Galactic OH+, H2O+, and H3O+: Probing the Molecular Hydrogen Fraction and Cosmic-Ray Ionization Rate
• HI-to-H2 Transitions and H I Column Densities in Galaxy Star-Forming Regions
• CO-dark gas and molecular filaments in Milky Way type galaxies
• Photodissociation and X-Ray Dominated Regions
• Vertical structure of a supernova-driven turbulent magnetized ISM
• Superbubbles in the Multiphase ISM and the Loading of Galactic Winds
• Magnetohydronamic Evolution of HII Regions in Molecular Clouds: Simulation Methodology, Tests, and Uniform Media
• First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations
• Mutual influence of supernovae and molecular clouds
• Galaxies that Shine: radiation-hydrodynamical simulations of disk galaxies
• The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies
• Star Formation Efficiency and Dispersal of Giant Molecular Clouds with UV Radiation Feedback: Dependence on Gravitational Boundedness and Magnetic Fields
• The SILCC project --- IV. Impact of dissociating and ionising radiation on the interstellar medium and Halpha emission as a tracer of the star formation rate
• The GRIFFIN project -- Formation of star clusters with individual massive stars in a simulated dwarf galaxy starburst
• SILCC VI -- Multi-phase ISM structure, stellar clustering, and outflows with supernovae, stellar winds, ionising radiation and cosmic rays
• The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces
• Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields
• Pressure-Regulated, Feedback-Modulated Star Formation In Disk Galaxies
• Efficiently Cooled Stellar Wind Bubbles in Turbulent Clouds I. Fractal Theory and Application to Star-Forming Clouds
• The ISM in spiral galaxies: can cooling in spiral shocks produce molecular clouds?
• Efficient early stellar feedback can suppress galactic outflows by reducing supernova clustering
• Snap, Crackle, Pop: sub-grid supernova feedback in AMR simulations of disk galaxies
• Efficiently Cooled Stellar Wind Bubbles in Turbulent Clouds II. Validation of Theory with Hydrodynamic Simulations
• The environmental dependence of the X_CO conversion factor
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• Physical Processes in Star Formation
• The Polycyclic Aromatic Hydrocarbon Mass Fraction on a 10 pc scale in the Magellanic Clouds
• Cosmic ray feedback from supernovae in dwarf galaxies
• Thermal and nonthermal dust sputtering in hydrodynamical simulations of the multiphase interstellar medium
• Metallicity dependence of the H/H${}_2$ and C${}^{+}$/C/CO distributions in a resolved self-regulating interstellar medium
• Mapping Luminous Hot Stars in the Galaxy
• The natural emergence of the correlation between H2 and star formation rate surface densities in galaxy simulations
• Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Method and Tests
• The physics of Lyman-alpha escape from disc-like galaxies
• Comparing the pre-SNe feedback and environmental pressures for 6000 HII regions across 19 nearby spiral galaxies
• A Schmidt-Kennicutt law for star formation in the Milky Way disk
• The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies
• Galactic Spiral Shocks with Thermal Instability
• Radiation-MagnetoHydrodynamics simulations of cosmic ray feedback in disc galaxies
• Prestellar Core Formation, Evolution, and Accretion from Gravitational Fragmentation in Turbulent Converging Flows
• Momentum feedback from marginally-resolved HII regions in isolated disc galaxies
• Dependence of $X_{\mathrm{CO}}$ on metallicity, intensity, and spatial scale in a self-regulated interstellar medium
• Diffusive Versus Free-Streaming Cosmic Ray Transport in Molecular Clouds
• RAMSES-RTZ: Non-Equilibrium Metal Chemistry and Cooling Coupled to On-The-Fly Radiation Hydrodynamics
• Cosmic-Ray Transport in Simulations of Star-forming Galactic Disks
• Diffuse Ionized Gas in Simulations of Multiphase, Star-Forming Galactic Disks
• Cosmic-Ray Transport in Varying Galactic Environments
• The global structure of magnetic fields and gas in simulated Milky Way-analogue galaxies
• Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations
• Thermal instability and multi-phase interstellar medium in the first galaxies
• Radiative Supernova Remnants and Supernova Feedback
• Influence of Ion-Neutral Damping on the Cosmic-Ray Streaming Instability: Magnetohydrodynamic Particle-in-cell Simulations
• The Life Cycle of the Central Molecular Zone. II: Distribution of atomic and molecular gas tracers
• Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay
• Rigorous theory for secondary cosmic-ray ionization
• HI-to-H${}_2$ Transitions in Dust-Free Interstellar Gas
• Non-equilibrium ionisation and radiative transport in an evolving supernova remnant
• Smoothed Particle Radiation Hydrodynamics: Two-Moment method with Local Eddington Tensor Closure
• Bimodal Behavior and Convergence Requirement in Macroscopic Properties of the Multiphase Interstellar Medium Formed by Atomic Converging Flows
• A new ionisation network and radiation transport module in PLUTO
• When did the initial mass function become bottom-heavy?

Cited by in corpus (19)

• Introducing TIGRESS-NCR: I. Co-Regulation of the Multiphase Interstellar Medium and Star Formation Rates
• SILCC VII -- Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities
• Disk settling and dynamical heating: histories of Milky Way-mass stellar disks across cosmic time in the FIRE simulations
• RIGEL: Simulating dwarf galaxies at solar mass resolution with radiative transfer and feedback from individual massive stars
• Boosting galactic outflows with enhanced resolution
• Hybrid-chimes: A model for radiative cooling and the abundances of ions and molecules in simulations of galaxy formation
• Radiative and mechanical energies in galaxies I. Contributions of molecular shocks and PDRs in 3C 326 N
• Metallicity Dependence of Molecular Cloud Hierarchical Structure at Early Evolutionary Stages
• Impact of star formation models on the growth of simulated galaxies at high redshifts
• Population III star formation in the presence of turbulence, magnetic fields and ionizing radiation feedback
• The Noctua Suite of Simulations -- The Difficulty of Growing Massive Black Holes in Low-Mass Dwarf Galaxies
• Cloud Formation by Supernova Implosion
• The Milky Way satellite galaxy Leo T: A perturbed cored dwarf
• The origin and evolution of the [CII] deficit in HII regions and star-forming molecular clouds
• Effect of cosmic rays and ionizing radiation on observational ultraviolet plasma diagnostics in the circumgalactic medium
• RIGEL: Feedback-regulated cloud-scale star formation efficiency in a simulated dwarf galaxy merger
• SILCC -- IX. The multi-phase interstellar medium at low metallicity
• Imladris: a detailed and flexible model for galaxy simulations with individual stars
• Stellar associations powering HII regions $\unicode{x2013}$ II. Escape fraction of ionizing photons