# binary_c reaction networks:
# to be read by the script found in
# src/perl/scripts/make_reaction_network_code.pl
#
#
############################################################
# pp chains
############################################################

Network pp

# as the full ppI/II/III (below) with H2 (and more?) in eq

SIGMAV_PP H1(H1,beta)H2 # (1.39e-20)
SIGMAV_DP H2(H1,gamma)He3 # (1.43e-3)
SIGMAV_He3_He3 He3(He3,2H1)He4 # (2.16e-13)
SIGMAV_He3_He4 He3(He4,gamma)Be7 # (1.64e-18)
SIGMAV_Be7_e Be7(e,gamma)Li7 # 
SIGMAV_Be7_p Be7(H1,He4)He4 # neglects B8 and Be8 formation and decay
SIGMAV_Li7_p_alpha Li7(H1,He4)He4 # (1.31e-7) 

# assume some reactions are in equilibrium
#H2 equilibrium
#Li7 equilibrium
#Be7 equilibrium
#He3 equilibrium # fails because it requires a quadratic solution

skip



Network coldCNO

# cold CNO with H,He (instantaneous beta decays)

SIGMAV_T12 C12(H1,gamma)C13
SIGMAV_T13 C13(H1,gamma)N14
SIGMAV_T14 N14(H1,gamma)N15
SIGMAV_T15 N15(H1,He4)C12
SIGMAV_T15_BRANCH N15(H1,gamma)O16
SIGMAV_T16 O16(H1,gamma)O17
SIGMAV_T17 O17(H1,He4)N14
SIGMAV_O17G O17(H1,gamma)O18
SIGMAV_O18A O18(H1,He4)N15
SIGMAV_O18G O18(H1,gamma)F19
SIGMAV_F19A F19(H1,He4)O16

exit

# setting any of these in equilibrium fails!

Network hotCNO

# hot CNO with H,He (with some beta decays) 

SIGMAV_T12 C12(H1,gamma)N13
SIGMAV_TBETA13 N13(beta,nu)C13
SIGMAV_T13 C13(H1,gamma)N14

# full reaction
SIGMAV_TPN13 N13(H1,gamma)O14
SIGMAV_TO14BETA O14(beta,nu)N14
# simplified (cold) form
#SIGMAV_TPN13 N13(H1,gamma)N14

SIGMAV_T14 N14(H1,gamma)O15
SIGMAV_TBETA15 O15(beta,nu)N15
SIGMAV_T15 N15(H1,He4)C12
SIGMAV_T15_BRANCH N15(H1,gamma)O16

# full reaction
SIGMAV_T16 O16(H1,gamma)F17
SIGMAV_TBETA17 F17(beta,nu)O17
# simplified (cold) form
#SIGMAV_T16 O16(H1,gamma)O17

# the next reaction is F17(H1,gamma)Ne18 then instant Ne18(beta)F18
#SIGMAV_TF17 F17(H1,gamma)F18

SIGMAV_T17 O17(H1,He4)N14

# full reaction
SIGMAV_O17G O17(H1,gamma)F18
SIGMAV_F18BETA F18(beta,nu)O18
# simplified (cold) form
#SIGMAV_O17G O17(H1,gamma)O18

SIGMAV_O18A O18(H1,He4)N15
SIGMAV_O18G O18(H1,gamma)F19
SIGMAV_F19A F19(H1,He4)O16

Network NeNa

SIGMAV_T20 Ne20(H1,gamma)Ne21
SIGMAV_T21 Ne21(H1,gamma)Ne22
SIGMAV_T22 Ne22(H1,gamma)Na23
SIGMAV_T23 Na23(H1,He4)Ne20

Network MgAl

SIGMAV_T24 Mg24(H1,gamma)Mg25
SIGMAV_T25 Mg25(H1,gamma)Al26
SIGMAV_TBETA26 Al26(beta,nu)Mg26
SIGMAV_T26P Al26(H1,gamma)Al27
SIGMAV_T26 Mg26(H1,gamma)Al27
SIGMAV_T27 Al27(H1,gamma)Si28
SIGMAV_T27a Al27(H1,He4)Mg24

Network NeNaMgAlnoleak

# NeNa/MgAl with no leak
SIGMAV_T20 Ne20(H1,gamma)Ne21
SIGMAV_T21 Ne21(H1,gamma)Ne22
SIGMAV_T22 Ne22(H1,gamma)Na23
SIGMAV_T23 Na23(H1,He4)Ne20
SIGMAV_T24 Mg24(H1,gamma)Mg25
SIGMAV_T25 Mg25(H1,gamma)Al26
SIGMAV_TBETA26 Al26(beta,nu)Mg26
SIGMAV_T26P Al26(H1,gamma)Al27
SIGMAV_T26 Mg26(H1,gamma)Al27
SIGMAV_T27 Al27(H1,gamma)Si28
SIGMAV_T27a Al27(H1,He4)Mg24

Network NeNaMgAl

# full NeNa/MgAl (with leak)
SIGMAV_T20 Ne20(H1,gamma)Ne21
SIGMAV_T21 Ne21(H1,gamma)Ne22
SIGMAV_T22 Ne22(H1,gamma)Na23
SIGMAV_T23 Na23(H1,He4)Ne20
SIGMAV_T23P Na23(H1,gamma)Mg24
SIGMAV_T24 Mg24(H1,gamma)Mg25
SIGMAV_T25 Mg25(H1,gamma)Al26
SIGMAV_TBETA26 Al26(beta,nu)Mg26
SIGMAV_T26P Al26(H1,gamma)Al27
SIGMAV_T26 Mg26(H1,gamma)Al27
SIGMAV_T27 Al27(H1,gamma)Si28
SIGMAV_T27a Al27(H1,He4)Mg24