# C-burning with A=23 URCA rate module generator import pynucastro as pyna from pynucastro.networks import StarKillerNetwork library_file = "20180319default2" mylibrary = pyna.rates.Library(library_file) all_nuclei = ["p", "he4", "c12", "c13", "n13", "n14", "n15", "o14", "o15", "o16", "o17", "f17", "f18"] nova_library = mylibrary.linking_nuclei(all_nuclei, with_reverse=False) print(len(nova_library._rates)) rc = pyna.RateCollection(libraries=[nova_library]) comp = pyna.Composition(rc.get_nuclei()) comp.set_solar_like() rc.plot(outfile="nova.png", rho=1.e4, T=9.e7, comp=comp)
import pynucastro as pyna
import matplotlib.pyplot as plt
files = ["c12-pg-n13-ls09",
"c13-pg-n14-nacr",
"n13--c13-wc12",
"n13-pg-o14-lg06",
"n14-pg-o15-im05",
"n15-pa-c12-nacr",
"o14--n14-wc12",
"o15--n15-wc12",
"o14-ap-f17-Ha96c",
"f17-pg-ne18-cb09",
"ne18--f18-wc12",
"f18-pa-o15-il10"]
rc = pyna.RateCollection(files)
comp = pyna.Composition(rc.get_nuclei())
comp.set_solar_like()
rc.plot(rho=2.e6, T=3.e7, comp=comp, outfile="cno_flow.png")
library_file = "20180319default2"
mylibrary = pyna.rates.Library(library_file)
data_list = mylibrary.get_rates()
all_nuclei = [
"p", "he4", "ne20", "o20", "f20", "mg24", "al27", "o16", "si28", "s32",
"p31"
]
escn_library = mylibrary.linking_nuclei(all_nuclei, with_reverse=True)
escn_tabular = [
"f20--o20-toki", "ne20--f20-toki", "o20--f20-toki", "f20--ne20-toki"
]
rc = pyna.RateCollection(libraries=[escn_library])
comp = pyna.Composition(rc.get_nuclei())
comp.set_nuc("o16", 0.5)
comp.set_nuc("ne20", 0.3)
comp.set_nuc("mg24", 0.1)
comp.set_nuc("o20", 1.e-5)
comp.set_nuc("f20", 1.e-5)
comp.set_nuc("p", 1.e-5)
comp.set_nuc("he4", 1.e-2)
comp.set_nuc("al27", 1.e-2)
comp.set_nuc("si28", 1.e-2)
comp.set_nuc("s32", 1.e-2)
comp.set_nuc("p31", 1.e-2)
comp.normalize()
# C-burning with A=23 URCA rate module generator
import pynucastro as pyna
library_file = "20180319default2"
mylibrary = pyna.rates.Library(library_file)
all_reactants = [
"p", "he4", "c12", "o16", "ne20", "mg24", "si28", "s32", "ar36", "ca40",
"ti44", "cr48", "fe52", "ni56", "al27", "p31", "cl35", "k39", "sc43",
"v47", "mn51", "co55", "c14", "n13", "n14", "o18", "f18", "ne21"
]
subCh = mylibrary.linking_nuclei(all_reactants)
rc = pyna.RateCollection(libraries=[subCh])
comp = pyna.Composition(rc.get_nuclei())
comp.set_all(0.1)
comp.set_nuc("he4", 0.95)
comp.normalize()
rc.plot(outfile="subch2.pdf",
rho=1.e6,
T=1.e9,
comp=comp,
hide_xalpha=True,
size=(1500, 450),
node_size=500,
node_font_size=11,
node_color="#337dff",