[Mesa-users] How to calculate the neutrinos flux in different radii

[peter zhang]

Greetings mesa-users,
I am using MESA  to calculate the flux of different types of neutrinos in
different radii. There do not exist similar options in
`profile_columns.list`.  I try to calculate the number of reactions of
different types to represent the flux of different types of neutrinos.

I have extracted the screened rates in `run_star_extras.f90`, which is
attached in* the *attachment. The mass fractions of different isotopes are
also acquired from the model file. The unit of screened rates is related to
the number of reactants which has referred in the `net_def.f90`:
```
! so for a 2-body reaction, there are 2 Y factors, each with units
[moles/gram]
! and the rate_screened units for such a reaction are [grams/(mole-sec)],
! which when multiplied by [moles/gram]^2 gives the same units as dydt.
```

For example, the number (/cm3/s) of  reaction `pep`(p + e^- + p -> H2
+nu_e), which is a 3-body reaction, is the product of screened rate, the
molar density of p, and* the molar density of electrons*. The number
density of electrons is needed but *I cannot find an option or variable to
output the number density of electrons. Maybe there exist simpler methods
or I misunderstand something.* Any help or clarification you could provide
would be greatly appreciated.

I am looking forward to your reply.

Thank you for your time,

Peter
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221109/047a1a43/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f90
Type: application/octet-stream
Size: 21894 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221109/047a1a43/attachment.obj>