[Mesa-users] Re2: Calculate the epsilon of nuclear reaction

[Francis Timmes]

hi wenyu,

for versions around 11701 - 12778 one can extract the desired information 
by incorporating some calls into your run_star_extras routines.

https://lists.mesastar.org/pipermail/mesa-users/2022-April/013609.html

http://doi.org/10.5281/zenodo.3372837

this was how it was done in the past, before rob & earl did a nice job 
making the flow more general and a default option in 22.11.1.

fxt





> On Dec 29, 2022, at 12:19 AM, 辛文宇 <wyxin at nao.cas.cn> wrote:
> 
> Hi, Prof. Timmes,
> 
> Thanks very much for the infomation of newest version of MESA. These update may be useful for my
> following work. However, I have calculate a series of models with MESA 12778 and it is time-consumed
> to update those models with newest version of MESA.
> Now I want to ensure the ignition temperature of C+C reaction when X(C)=1, 0.1, 1E-3 and so on, 
> according to the formula of eps_nuc = eps_neu. I have find your code from https://cococubed.com/code_pages/nuloss.shtml
> which can calculate eps_neu and its derivatives. But I need eps_nuc and its derivatives now.
> Thus, it is better for me to find the subrotine in MESA to calculate these value independently.
> 
> I also checked the book of kippenhahn, in wich a formula to calculte eps_nuc of C+C is provided,
> but this formula used the C+C rate of CF88. In our models, the C+C rate from Tumino 2018 is also used in some
> models. I am not sure how to replace CF88 to Tumino 2018 in the fomula.
> 
> Thanks,
> Wenyu
> 
> 
> > -----原始邮件-----
> > 发件人: "Francis Timmes" <fxt44 at mac.com>
> > 发送时间: 2022-12-28 23:56:40 (星期三)
> > 收件人: "辛文宇" <wyxin at nao.cas.cn>, "xinwenyu16 at mails.ucas.ac.cn" <xinwenyu16 at mails.ucas.ac.cn>
> > 抄送: "Frank Timmes" <fxt44 at mac.com>, mesa-users <mesa-users at lists.mesastar.org>
> > 主题: Re: [Mesa-users] Calculate the epsilon of nuclear reaction
> > 
> > the current public release 22.11.1 includes the capability of
> > including individual nuclear reaction rates and related quantities 
> > in the history or profile data files - thank you rob & earl! 
> > 
> > see section 10.3 of the arxiv version of the mesa VI instrument paper
> > https://ui.adsabs.harvard.edu/abs/2022arXiv220803651J/abstract , 
> > and visit $MESA_DIR/r22.11.1/star/defaults.
> > 
> > fxt
> > 
> > 
> > 
> > 
> > 
> > > On Dec 28, 2022, at 8:06 AM, 辛文宇 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> > > 
> > > I know the subotine of  $MESA_DIR/net/test/src/sample_net.f90
> > > However, this subrotine give the epsilon of a centain network, 
> > > instead of one reaction. 
> > > 
> > > Thanks, 
> > > Wenyu 
> > > 
> > > 
> > > -----原始邮件-----
> > > 发件人:"辛文宇 via Mesa-users" <mesa-users at lists.mesastar.org>
> > > 发送时间:2022-12-28 22:57:14 (星期三)
> > > 收件人: mesa-users at lists.mesastar.org
> > > 抄送: "xinwenyu16 at mails.ucas.ac.cn" <xinwenyu16 at mails.ucas.ac.cn>
> > > 主题: [Mesa-users] Calculate the epsilon of nuclear reaction
> > > 
> > > Dear everyone: 
> > > 
> > > I want to calculate the energy generation rate(epsilon) for a certain reaction rate. 
> > > Such as 12C+12C reaction rate, I want to calculate the epsilon of 12C+12C reaction with rate from different references under different logT, logRho and mass fraction of 
> > > 12C. Where this subrotine is in MESA_dir. 
> > > 
> > > Thanks, 
> > > Wenyu 
> > > _______________________________________________
> > > mesa-users at lists.mesastar.org
> > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > > 
> </xinwenyu16 at mails.ucas.ac.cn></mesa-users at lists.mesastar.org></mesa-users at lists.mesastar.org></mesa-users at lists.mesastar.org></fxt44 at mac.com></xinwenyu16 at mails.ucas.ac.cn></wyxin at nao.cas.cn></fxt44 at mac.com>