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

辛文宇 wyxin at nao.cas.cn
Thu Dec 29 07:19:34 UTC 2022

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.


> -----原始邮件-----
> 发件人: "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>

More information about the Mesa-users mailing list