Skip to content

[WIP] update the nova networks #1962

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
zingale wants to merge 3 commits into
base: development
Choose a base branch
from
Open

[WIP] update the nova networks #1962

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
ea3bd11
update the nova networks
zingale Mar 11, 2026
c1685a2
Merge branch 'development' into update_nova
zingale Mar 26, 2026
905986c
Merge branch 'development' into update_nova
zingale Apr 12, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
33 changes: 22 additions & 11 deletions Docs/source/networks.rst
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -185,31 +185,37 @@ This network is composed of 17 nuclei: $\isotm{H}{1,2}$,
$\isotm{He}{3,4}$, $\isotm{Be}{7}$, $\isotm{B}{8}$,
$\isotm{C}{12,13}$, $\isotm{N}{13-15}$, $\isotm{O}{14-17}$,
$\isotm{F}{17,18}$ and is used to model the onset of a classical novae
thermonuclear runaway. The first set of nuclei, $\isotm{H}{1,2}$,
$\isotm{He}{3,4}$ represent the pp-chain sector of the reaction
network, while the second set, of $\isotm{Be}{7}$, and $\isotm{B}{8}$,
describe the involvement of the x-process. Finally, all the remaining
nuclei are active participants of the CNO cycle with endpoints at
$\isotm{F}{17}$ and $\isotm{F}{18}$. The triple-$\alpha$ reaction
thermonuclear runaway. These are linked together by 55 rates (including reverse rates).

The first set of nuclei, $\isotm{H}{1,2}$, $\isotm{He}{3,4}$ represent
the pp-chain sector of the reaction network, while the second set, of
$\isotm{Be}{7}$, and $\isotm{B}{8}$, describe the involvement of the
x-process. Finally, all the remaining nuclei are active participants
of the CNO cycle with endpoints at $\isotm{F}{17}$ and
$\isotm{F}{18}$. The triple-$\alpha$ reaction
$\alpha(\alpha\alpha,\gamma)\isotm{C}{12}$, serves as bridge between
the nuclei of first and the last set.

The the cold-CNO chain of reactions of the CN-branch are:
The cold-CNO chain of reactions of the CN-branch are:

* :math:`\isotm{C}{12}(p,\gamma)\isotm{N}{13}(\beta^{+}\nu_e)\isotm{C}{13}(p,\gamma)`

while the NO-branch chain of reactions is:

* :math:`\isotm{N}{14}(p,\gamma)\isotm{O}{15}(\beta^{+})\isotm{N}{15}(p,\gamma)\isotm{O}{16}(p,\gamma)\isotm{F}{17}(\beta^{+}\nu_e)\isotm{O}{17}`

where the isotopes $\isotm{N}{15}$ and $\isotm{O}{17}$ may decay back into $\isotm{C}{12}$ and $\isotm{N}{14}$ through
$\isotm{N}{15}(p,\alpha)\isotm{C}{12}$ and $\isotm{O}{17}(p,\alpha)\isotm{N}{14}$ respectively.
where the isotopes $\isotm{N}{15}$ and $\isotm{O}{17}$ may decay back
into $\isotm{C}{12}$ and $\isotm{N}{14}$ through
$\isotm{N}{15}(p,\alpha)\isotm{C}{12}$ and
$\isotm{O}{17}(p,\alpha)\isotm{N}{14}$ respectively.

.. figure:: ../../networks/nova/nova.png
:align: center

Once the temperature reaches a threshold of $\gtrsim 10^8\,\mathrm{K}$, the fast $p$-captures, for example,
$\isotm{N}{13}(p,\gamma)\isotm{O}{14}$, are more likely than the $\beta^{+}$-decays $\isotm{N}{13}(\beta^{+}\nu_e)\isotm{C}{13}$
Once the temperature reaches a threshold of $\gtrsim
10^8\,\mathrm{K}$, the fast $p$-captures, for example,
$\isotm{N}{13}(p,\gamma)\isotm{O}{14}$, are more likely than the
$\beta^{+}$-decays $\isotm{N}{13}(\beta^{+}\nu_e)\isotm{C}{13}$
reactions. These rates are also included in this network.

``nova-li``
Expand All @@ -220,6 +226,11 @@ nuclei beyond fluorine. It should give a more accurate energy in late
stages of the burst, and can also be used to explore lithium
production.

To keep the network size down, a filtering process is done to remove
rates that are not expected to be important under the thermodynamic
conditions encountered in a nova. The result is that the network has
25 nuclei and 54 rates.

.. figure:: ../../networks/nova-li/nova-li.png
:align: center

Expand Down
206 changes: 103 additions & 103 deletions networks/nova-li/actual_network.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -243,57 +243,57 @@ namespace Rates
k_N13_to_C13_reaclib = 2,
k_O14_to_N14_reaclib = 3,
k_O15_to_N15_reaclib = 4,
k_B8_to_p_Be7_reaclib = 5,
k_B8_to_He4_He4_reaclib = 6,
k_F17_to_p_O16_reaclib = 7,
k_p_p_to_d_reaclib_bet_pos = 8,
k_p_p_to_d_reaclib_electron_capture = 9,
k_p_d_to_He3_reaclib = 10,
k_d_d_to_He4_reaclib = 11,
k_p_He3_to_He4_reaclib = 12,
k_He4_He3_to_Be7_reaclib = 13,
k_He4_Li7_to_B11_reaclib = 14,
k_p_Be7_to_B8_reaclib = 15,
k_p_B11_to_C12_reaclib = 16,
k_p_C12_to_N13_reaclib = 17,
k_He4_C12_to_O16_reaclib = 18,
k_p_C13_to_N14_reaclib = 19,
k_p_N13_to_O14_reaclib = 20,
k_p_N14_to_O15_reaclib = 21,
k_He4_N14_to_F18_reaclib = 22,
k_p_N15_to_O16_reaclib = 23,
k_He4_N15_to_F19_reaclib = 24,
k_He4_O14_to_Ne18_reaclib = 25,
k_He4_O15_to_Ne19_reaclib = 26,
k_p_O16_to_F17_reaclib = 27,
k_He4_O16_to_Ne20_reaclib = 28,
k_p_O17_to_F18_reaclib = 29,
k_He4_O17_to_Ne21_reaclib = 30,
k_p_O18_to_F19_reaclib = 31,
k_p_F17_to_Ne18_reaclib = 32,
k_p_F18_to_Ne19_reaclib = 33,
k_p_F19_to_Ne20_reaclib = 34,
k_d_He3_to_p_He4_reaclib = 35,
k_p_Li7_to_He4_He4_reaclib = 36,
k_He4_N13_to_p_O16_reaclib = 37,
k_p_N15_to_He4_C12_reaclib = 38,
k_He4_O14_to_p_F17_reaclib = 39,
k_p_O17_to_He4_N14_reaclib = 40,
k_p_O18_to_He4_N15_reaclib = 41,
k_He4_F17_to_p_Ne20_reaclib = 42,
k_p_F18_to_He4_O15_reaclib = 43,
k_He4_F18_to_p_Ne21_reaclib = 44,
k_p_F19_to_He4_O16_reaclib = 45,
k_He3_He3_to_p_p_He4_reaclib = 46,
k_d_Be7_to_p_He4_He4_reaclib = 47,
k_p_B11_to_He4_He4_He4_reaclib = 48,
k_He3_Be7_to_p_p_He4_He4_reaclib = 49,
k_He4_He4_He4_to_C12_reaclib = 50,
k_F17_to_O17_weaktab = 51,
k_F18_to_O18_weaktab = 52,
k_Ne18_to_F18_weaktab = 53,
k_Ne19_to_F19_weaktab = 54,
NumRates = k_Ne19_to_F19_weaktab
k_B8_to_He4_He4_reaclib = 5,
k_p_p_to_d_reaclib_beta_pos = 6,
k_p_p_to_d_reaclib_electron_capture = 7,
k_p_d_to_He3_reaclib = 8,
k_d_d_to_He4_reaclib = 9,
k_p_He3_to_He4_reaclib = 10,
k_He4_He3_to_Be7_reaclib = 11,
k_He4_Li7_to_B11_reaclib = 12,
k_p_Be7_to_B8_reaclib = 13,
k_p_B11_to_C12_reaclib = 14,
k_p_C12_to_N13_reaclib = 15,
k_He4_C12_to_O16_reaclib = 16,
k_p_C13_to_N14_reaclib = 17,
k_p_N13_to_O14_reaclib = 18,
k_p_N14_to_O15_reaclib = 19,
k_He4_N14_to_F18_reaclib = 20,
k_p_N15_to_O16_reaclib = 21,
k_He4_N15_to_F19_reaclib = 22,
k_He4_O14_to_Ne18_reaclib = 23,
k_He4_O15_to_Ne19_reaclib = 24,
k_p_O16_to_F17_reaclib = 25,
k_He4_O16_to_Ne20_reaclib = 26,
k_p_O17_to_F18_reaclib = 27,
k_He4_O17_to_Ne21_reaclib = 28,
k_p_O18_to_F19_reaclib = 29,
k_p_F17_to_Ne18_reaclib = 30,
k_p_F18_to_Ne19_reaclib = 31,
k_p_F19_to_Ne20_reaclib = 32,
k_d_He3_to_p_He4_reaclib = 33,
k_p_Li7_to_He4_He4_reaclib = 34,
k_He4_N13_to_p_O16_reaclib = 35,
k_p_N15_to_He4_C12_reaclib = 36,
k_He4_O14_to_p_F17_reaclib = 37,
k_p_O17_to_He4_N14_reaclib = 38,
k_p_O18_to_He4_N15_reaclib = 39,
k_p_F18_to_He4_O15_reaclib = 40,
k_He4_F18_to_p_Ne21_reaclib = 41,
k_p_F19_to_He4_O16_reaclib = 42,
k_He3_He3_to_p_p_He4_reaclib = 43,
k_d_Be7_to_p_He4_He4_reaclib = 44,
k_p_B11_to_He4_He4_He4_reaclib = 45,
k_He3_Be7_to_p_p_He4_He4_reaclib = 46,
k_He4_He4_He4_to_C12_reaclib = 47,
k_F17_to_O17_weaktab = 48,
k_F18_to_O18_weaktab = 49,
k_Ne18_to_F18_weaktab = 50,
k_Ne19_to_F19_weaktab = 51,
k_B8_to_p_Be7_derived = 52,
k_F17_to_p_O16_derived = 53,
k_He4_F17_to_p_Ne20_derived = 54,
NumRates = k_He4_F17_to_p_Ne20_derived
};

// rate names -- note: the rates are 1-based, not zero-based, so we pad
Expand All @@ -306,56 +306,56 @@ namespace Rates
"N13_to_C13_reaclib", // 2,
"O14_to_N14_reaclib", // 3,
"O15_to_N15_reaclib", // 4,
"B8_to_p_Be7_reaclib", // 5,
"B8_to_He4_He4_reaclib", // 6,
"F17_to_p_O16_reaclib", // 7,
"p_p_to_d_reaclib_bet_pos", // 8,
"p_p_to_d_reaclib_electron_capture", // 9,
"p_d_to_He3_reaclib", // 10,
"d_d_to_He4_reaclib", // 11,
"p_He3_to_He4_reaclib", // 12,
"He4_He3_to_Be7_reaclib", // 13,
"He4_Li7_to_B11_reaclib", // 14,
"p_Be7_to_B8_reaclib", // 15,
"p_B11_to_C12_reaclib", // 16,
"p_C12_to_N13_reaclib", // 17,
"He4_C12_to_O16_reaclib", // 18,
"p_C13_to_N14_reaclib", // 19,
"p_N13_to_O14_reaclib", // 20,
"p_N14_to_O15_reaclib", // 21,
"He4_N14_to_F18_reaclib", // 22,
"p_N15_to_O16_reaclib", // 23,
"He4_N15_to_F19_reaclib", // 24,
"He4_O14_to_Ne18_reaclib", // 25,
"He4_O15_to_Ne19_reaclib", // 26,
"p_O16_to_F17_reaclib", // 27,
"He4_O16_to_Ne20_reaclib", // 28,
"p_O17_to_F18_reaclib", // 29,
"He4_O17_to_Ne21_reaclib", // 30,
"p_O18_to_F19_reaclib", // 31,
"p_F17_to_Ne18_reaclib", // 32,
"p_F18_to_Ne19_reaclib", // 33,
"p_F19_to_Ne20_reaclib", // 34,
"d_He3_to_p_He4_reaclib", // 35,
"p_Li7_to_He4_He4_reaclib", // 36,
"He4_N13_to_p_O16_reaclib", // 37,
"p_N15_to_He4_C12_reaclib", // 38,
"He4_O14_to_p_F17_reaclib", // 39,
"p_O17_to_He4_N14_reaclib", // 40,
"p_O18_to_He4_N15_reaclib", // 41,
"He4_F17_to_p_Ne20_reaclib", // 42,
"p_F18_to_He4_O15_reaclib", // 43,
"He4_F18_to_p_Ne21_reaclib", // 44,
"p_F19_to_He4_O16_reaclib", // 45,
"He3_He3_to_p_p_He4_reaclib", // 46,
"d_Be7_to_p_He4_He4_reaclib", // 47,
"p_B11_to_He4_He4_He4_reaclib", // 48,
"He3_Be7_to_p_p_He4_He4_reaclib", // 49,
"He4_He4_He4_to_C12_reaclib", // 50,
"F17_to_O17_weaktab", // 51,
"F18_to_O18_weaktab", // 52,
"Ne18_to_F18_weaktab", // 53,
"Ne19_to_F19_weaktab" // 54,
"B8_to_He4_He4_reaclib", // 5,
"p_p_to_d_reaclib_beta_pos", // 6,
"p_p_to_d_reaclib_electron_capture", // 7,
"p_d_to_He3_reaclib", // 8,
"d_d_to_He4_reaclib", // 9,
"p_He3_to_He4_reaclib", // 10,
"He4_He3_to_Be7_reaclib", // 11,
"He4_Li7_to_B11_reaclib", // 12,
"p_Be7_to_B8_reaclib", // 13,
"p_B11_to_C12_reaclib", // 14,
"p_C12_to_N13_reaclib", // 15,
"He4_C12_to_O16_reaclib", // 16,
"p_C13_to_N14_reaclib", // 17,
"p_N13_to_O14_reaclib", // 18,
"p_N14_to_O15_reaclib", // 19,
"He4_N14_to_F18_reaclib", // 20,
"p_N15_to_O16_reaclib", // 21,
"He4_N15_to_F19_reaclib", // 22,
"He4_O14_to_Ne18_reaclib", // 23,
"He4_O15_to_Ne19_reaclib", // 24,
"p_O16_to_F17_reaclib", // 25,
"He4_O16_to_Ne20_reaclib", // 26,
"p_O17_to_F18_reaclib", // 27,
"He4_O17_to_Ne21_reaclib", // 28,
"p_O18_to_F19_reaclib", // 29,
"p_F17_to_Ne18_reaclib", // 30,
"p_F18_to_Ne19_reaclib", // 31,
"p_F19_to_Ne20_reaclib", // 32,
"d_He3_to_p_He4_reaclib", // 33,
"p_Li7_to_He4_He4_reaclib", // 34,
"He4_N13_to_p_O16_reaclib", // 35,
"p_N15_to_He4_C12_reaclib", // 36,
"He4_O14_to_p_F17_reaclib", // 37,
"p_O17_to_He4_N14_reaclib", // 38,
"p_O18_to_He4_N15_reaclib", // 39,
"p_F18_to_He4_O15_reaclib", // 40,
"He4_F18_to_p_Ne21_reaclib", // 41,
"p_F19_to_He4_O16_reaclib", // 42,
"He3_He3_to_p_p_He4_reaclib", // 43,
"d_Be7_to_p_He4_He4_reaclib", // 44,
"p_B11_to_He4_He4_He4_reaclib", // 45,
"He3_Be7_to_p_p_He4_He4_reaclib", // 46,
"He4_He4_He4_to_C12_reaclib", // 47,
"F17_to_O17_weaktab", // 48,
"F18_to_O18_weaktab", // 49,
"Ne18_to_F18_weaktab", // 50,
"Ne19_to_F19_weaktab", // 51,
"B8_to_p_Be7_derived", // 52,
"F17_to_p_O16_derived", // 53,
"He4_F17_to_p_Ne20_derived" // 54,
};

}
Expand All @@ -382,8 +382,8 @@ namespace NSE_INDEX
constexpr int NumNSERatePairs = 2;

inline AMREX_GPU_MANAGED amrex::Array2D<std::int8_t, 1, NumNSERatePairs, 1, 10, amrex::Order::C> rate_pair_data {
H1, Be7, -1, B8, -1, -1, Be7, B8, k_p_Be7_to_B8_reaclib, k_B8_to_p_Be7_reaclib,
H1, O16, -1, F17, -1, -1, O16, F17, k_p_O16_to_F17_reaclib, k_F17_to_p_O16_reaclib
H1, Be7, -1, B8, -1, -1, Be7, B8, k_p_Be7_to_B8_reaclib, k_B8_to_p_Be7_derived,
H1, O16, -1, F17, -1, -1, O16, F17, k_p_O16_to_F17_reaclib, k_F17_to_p_O16_derived
};
}
#endif
Expand Down
Loading
Loading