Gamma decay energy documentation

docs/physics/tardisgamma/decayenergy.ipynb

@@ -0,0 +1,108 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Radioactive Decay Energy\n",
+    "\n",
+    "Within the ejecta of a supernova, the gammarays largely come from the decay of $^{56}Ni$ into $^{56}Co$.\n",
+    "This releases a large amount of energy. To understand how much energy is produced from this decay we look at the amount of energy 1 g of $^{56}Ni$ produces in 10 days"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {
+    "vscode": {
+     "languageId": "plaintext"
+    }
+   },
+   "source": [
+    "The equation for radioactive decay is \n",
+    "$$N(t) = N_\\mathrm{0} \\exp\\left(-\\lambda t \\right)$$\n",
+    "$N(t)$ is the number of atoms after time t, $N_\\mathrm{0}$ is the initial number of atoms, and $\\lambda$ is the half life."
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Using the molar mass, we find that 1g of $^{56}Ni$ has $1.0767e22$ atoms. The half life of $^{56}Ni$ is 6.10 days. Plugging these numbers into the equation we can find the number of atoms that have decayed after 10 days."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "from astropy import units as u\n",
+    "from astropy import constants as const\n",
+    "\n",
+    "#initial number of atoms\n",
+    "n_0 = 1.0767e22\n",
+    "#half life\n",
+    "t_half = 6.10\n",
+    "#time in days\n",
+    "time = 10\n",
+    "\n",
+    "#number of atoms keft after 10 days\n",
+    "n_final = n_0 * np.exp(-t_half * time)\n",
+    "\n",
+    "#number of atoms that have decayed\n",
+    "n_decays = n_0 - n_final"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "The gammarays that are released during this radiocative decay can be released from a number of energy level transitions. Each transition has an associated energy and a probability out of 100 decays. From the combination of all of these tranisitions and their probabilities, the total energy released per decay is $1.75 MeV$ <strong data-cite=\"1994Nadyozhin\">[]</strong>."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "The total energy released from radioactive decay of 1g of 56Ni after 10 days is: 1.884225e+28 eV\n"
+     ]
+    }
+   ],
+   "source": [
+    "#energy per decay in eV\n",
+    "energy_per_decay = 1.75e6\n",
+    "\n",
+    "total_energy = n_decays * energy_per_decay\n",
+    "\n",
+    "print(\"The total energy released from radioactive decay of 1g of 56Ni after 10 days is:\", total_energy, \"eV\")"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "tardis",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.12.4"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}

docs/physics/tardisgamma/index.rst

@@ -5,3 +5,4 @@ TARDIS gamma :math:`\gamma`
 
 .. toctree::
     packetinitialization
+    decayenergy

docs/tardis.bib

@@ -348,3 +348,18 @@ @ARTICLE{Boyle2017
        adsurl = {https://ui.adsabs.harvard.edu/abs/2017A&A...599A..46B},
       adsnote = {Provided by the SAO/NASA Astrophysics Data System}
 }
+
+@ARTICLE{1994Nadyozhin,
+       author = {{Nadyozhin}, D.~K.},
+        title = "{The Properties of NI CO Fe Decay}",
+      journal = {\apjs},
+     keywords = {Cobalt Isotopes, Decay, Electron Capture, Gamma Rays, Nickel Isotopes, Nuclear Astrophysics, Nuclear Fusion, Half Life, Iron Isotopes, Neutrinos, Astronomy, ATOMIC DATA, NUCLEAR REACTIONS, NUCLEOSYNTHESIS, ABUNDANCES},
+         year = 1994,
+        month = jun,
+       volume = {92},
+        pages = {527},
+          doi = {10.1086/192008},
+       adsurl = {https://ui.adsabs.harvard.edu/abs/1994ApJS...92..527N},
+      adsnote = {Provided by the SAO/NASA Astrophysics Data System}
+}
+