AstroChemistry.jl is a package that solves a user-defined chemistry network for astrophysical applications (in particular for the interstellar medium). It can be used to post-process snapshots of hydrodynamical simulations with a built-in Octree that calculates the column densities for shielding against the FUV radiation. The chemistry coefficients are based on the UMIST database and the CO shielding is based on Visser et al. 2009. The code is multithreading parallel using Julia's built-in functionality.
AstroChemistry.jl makes use of the following Julia packages:
Open the REPL (type julia in the terminal) and install the relevant packages:
import Pkg
Pkg.add("StaticArrays")
Pkg.add("Parameters")
Pkg.add("DifferentialEquations")
Pkg.add("Healpix")
Pkg.add("AstroChemistry")
using AstroChemistry
#specify the chemical species in the network
all_species = ["H", "H-", "H2", "H+", "H2+", "e-"]
#initialize the network
net, dict = initialize_chemistry_network(all_species, grRec=false)
#specify the physical parameters
nH, temp, ξ, IUV, Zp = 100., 20., 1e-16, 0., 1.
par = Par{1,0,Float64}(nH=nH, temp=temp, ξ=ξ, IUV=IUV, Zp=Zp)
#initialize the total (gas-phase) elemental abundance
abtot = AbundTotal() * Zp
#initialize the chemical abundance array
N_spec = length(all_species)
abund = zeros(N_spec)
init_abund(abund, [], abtot, net)
#solve the network
dt = 1e10 #yr
retcode, rr = solve_equilibrium_abundances(abund, dt, par, abtot, net)
#print the final abundances of H+, H and H2
@show retcode, abund[dict["H+"]], abund[dict["H"]], abund[dict["H2"]]
-
1D PDR calculation: see examples/PDRtest.jl for reproducing the "F1" test problem described in the PDR code comparison project Röllig et al. 2017.
-
post-processing: see examples/chem_octree.jl for how to read a Gadget snapshot (format 3) and calculate the chemical abundances in post-processing.
-
The output image of examples/PDRtest.jl.
-
Post-processing results from a hydrodynamical simulation of the interstellar medium in the solar-neighborhood (Hu et al. 2021).
Chia-Yu Hu @ Max Planck Institute for Extraterrestrial Physics (cyhu.astro@gmail.com)