Galactic chemical evolution in python
kimmy contains simple tools to study chemical evolution in galaxies.
Jo Bovy (University of Toronto): bovy - at - astro - dot - utoronto - dot - ca
Install the latest release with pip
pip install kimmy
or install the latest version by cloning/forking/downloading the repository and installing using
sudo python setup.py install
or locally using
python setup.py install --user
For an example of usage, see the example notebook. You can also launch a Binder instance and directly play around with this notebook. Or you can use this pyodide-compatible version of the same notebook in JupyterLite.
Currently, the only implemented feature is a simple one-zone chemical model with two elements O (for oxygen) and Fe (for iron). Initialize this model as
import kimmy
oz= kimmy.OneZone()
then for example compute the evolution of the default model and plot the [O/Fe] vs. [Fe/H] sequence
ts= numpy.linspace(0.001,10.,1001)*u.Gyr
plot(oz.Fe_H(ts),oz.O_Fe(ts))
To compute the distribution of [Fe/H], do for example,
FeHs= numpy.linspace(-1.525,1.225,56)
FeH_dist= [oz.Fe_H_DF(f) for f in FeHs]
and similar for the distribution of [O/H] and [O/Fe]. You can directly update the main parameters of the model and the model will be re-computed. For example, to set the outflow mass-loading parameter to one and plot the [O/Fe] vs. [Fe/H] sequence, do
ts= numpy.linspace(0.001,10.,1001)*u.Gyr
oz.eta= 1.
plot(oz.Fe_H(ts),oz.O_Fe(ts))
Keep in mind that once you change a parameter, it remains changed in the model. If you want to go back to the initial set of parameters that you used to initialize the instance, use oz.initial(); if you want to go back to the default set of parameters, use oz.default(). If you want to print the model you are using at any time, do
print(oz)
which prints a nicely formatted list of all of the model parameters.