Chemulator: Fast, accurate thermochemistry for dynamical models through emulation

Abstract

Chemical modelling serves two purposes in dynamical models: accounting for the effect of microphysics on the dynamics and providing observable signatures. Ideally, the former must be done as part of the hydrodynamic simulation but this comes with a prohibitive computational cost which leads to many simplifications being used in practice. To produce a statistical emulator that replicates a full chemical model capable of solving the temperature and abundances of a gas through time. This emulator should suffer only a minor loss of accuracy over including a full chemical solver in a dynamical model but would have a fraction of the computational cost. The gas-grain chemical code UCLCHEM was updated to include heating and cooling processes and a large dataset of model outputs from possible starting conditions was produced. A neural network was then trained to map directly from inputs to outputs. Chemulator replicates the outputs of UCLCHEM with an overall mean squared error (MSE) of 0.0002 for a single time step of 1000 yr and is shown to be stable over 1000 iterations with an MSE of 0.003 the log scaled temperature after one time step and 0.006 after 1000 time steps. Chemulator was found to be approximately 50,000 times faster than the time dependent model it emulates but can introduce a significant error to some models.

arXiv comment16 pages, 12 figures, accepted for publication in A&A

References in corpus (29)

• The chemical composition of the Sun
• Simulating Galaxy Formation with the IllustrisTNG Model
• The UMIST database for astrochemistry 2012
• Photoelectric Emission from Interstellar Dust: Grain Charging and Gas Heating
• A PDR-Code Comparison Study
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• A chemical model for the atmosphere of hot Jupiters
• KROME - a package to embed chemistry in astrophysical simulations
• [CII] 158$\mu$m Emission and Metallicity in PDRs
• On the convergence of the critical cooling timescale for the fragmentation of self-gravitating discs
• CO Depletion in Protoplanetary Disks: A Unified Picture Combining Physical Sequestration and Chemical Processing
• UCLCHEM: A Gas-Grain Chemical Code
• Non-convergence of the critical cooling timescale for fragmentation of self-gravitating discs
• Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests
• Radiative cooling rates, ion fractions, molecule abundances and line emissivities including self-shielding and both local and metagalactic radiation fields
• On the gas temperature in circumstellar disks around A stars
• Convergence of the critical cooling rate for protoplanetary disk fragmentation achieved; the key role of numerical dissipation of angular momentum
• Observable signatures of wind--driven chemistry with a fully consistent three dimensional radiative hydrodynamics model of HD 209458b
• The first multi-dimensional view of mass loss from externally FUV irradiated protoplanetary discs
• The chemistry of protoplanetary fragments formed via gravitational instabilities
• Grand challenges in protoplanetary disc modelling
• RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations
• Modelling massive-star feedback with Monte Carlo radiation hydrodynamics: photoionization and radiation pressure in a turbulent cloud
• The growth of H II regions around massive stars: the role of metallicity and dust
• Radiative cooling implementations in simulations of primordial star formation
• Modelling the ArH${}^{+}$ emission from the Crab Nebula
• Incorporating Astrochemistry into Molecular Line Modelling via Emulation
• Complexity reduction of astrochemical networks
• Modelling thermochemical processes in protoplanetary disks I: numerical methods

Cited by in corpus (16)

• Neural networks: solving the chemistry of the interstellar medium
• The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry
• 3D simulations of AGB stellar winds -- II. Ray-tracer implementation and impact of radiation on the outflow morphology
• Machine learning-accelerated chemistry modeling of protoplanetary disks
• De-noising of galaxy optical spectra with autoencoders
• Neural network-based emulation of interstellar medium models
• Quantifying the informativity of emission lines to infer physical conditions in giant molecular clouds. I. Application to model predictions
• MF-Box: Multi-fidelity and multi-scale emulation for the matter power spectrum
• Mapping Synthetic Observations to Prestellar Core Models: An Interpretable Machine Learning Approach
• Using a neural network approach to accelerate disequilibrium chemistry calculations in exoplanet atmospheres
• Radiative Transfer as a Bayesian Linear Regression problem
• NeuralPDR: Neural Differential Equations as surrogate models for Photodissociation Regions
• Neural-Network Chemical Emulator for First-Star Formation: Robust Iterative Predictions over a Wide Density Range
• Beetroots: spatially-regularized Bayesian inference of physical parameter maps -- Application to Orion
• Systematic selection of surrogate models for nonequilibrium chemistry
• Using Statistical Emulation and Knowledge of Grain-Surface Diffusion for Bayesian Inference of Reaction Rate Parameters: An Application to a Glycine Network