computing energy generation rate

[Anonymous]

Dear Brad,

I ran run_single_zone for a rp-prorcess trajectory. Is there a way to get the energy generation rate from the output xml?

Thanks a lot.

X.

[Bradley S. Meyer]

Thanks for your post. Yes, there are two codes that you might find useful. If your output is in nucnet-tools-code/examples/network/my_output.xml, you can do

cd nucnet-tools-code/examples/thermo
./examples_make

You can then compute the energy generation rate per gram by

./compute_energy_generation_rate ../network/my_output.xml

You can select out nuclides, reactions, and/or zones or you can compute the energy generation rate per nucleon. For example, try

./compute_energy_generation_rate --help

for a list of options. As another example, try

./compute_energy_generation_rate ../network/my_output.xml --reac_xpath "[count(reactant) = 2 and reactant = 'h1' and product = 'gamma' and count(product) = 2]" --eg_type nucleon

That will compute the energy generation per nucleon from (p,γ) reactions.

You can also compute the energy generation rate from individual reactions. The code is compute_energy_generation_rate_by_reaction. To get the list of possible options, try

./compute_energy_generation_rate_by_reaction

or

./compute_energy_generation_rate_by_reaction --help

For example, to get the dominant energy generating reactions for time dumps 10 to 20, try

./compute_energy_generation_rate_by_reaction ../network/my_output.xml --zone_xpath "[position() >= 10 and position() <= 20]" --eg_type nucleon --cutoff 1.e-5

The cutoff limits the output. You can also select out nuclides and reactions.

Note that the codes do not include energy lost from neutrinos.

Let me know if this doesn't work for you. I haven't yet blogged about this, but I hope to soon, including using the energy generation within the network calculation. Also, let me know if you see anything amiss or have any ideas for other things to include.

Best wishes.

[Anonymous]

Thank you, Brad. What is the unit for the energy generation rate?

X.

[Bradley S. Meyer]

The default units are ergs per second per gram. If you use the option --eg_type nucleon, the units will be ergs per second per nucleon. I've updated the codes to indicate this more clearly. Best wishes.