[Mesa-users] The nuclear burning after log T = 9.85 is stop ?

[Francis Timmes]

please keep all messages on mesa-users.

reaclib rates are fit up to 10e9 k. bad things can (and do) happen when 
these polynomial/rational fits are extrapolated beyond their fitting range.
mesa limits the temperature to 10e9 k to avoid such extrapolations and issues a warning.

fxt



> On Sep 18, 2022, at 8:24 PM, 辛文宇 <wyxin at nao.cas.cn> wrote:
> 
> Hi, Pro. Timmes,
> 
> Thanks very much for your help.
> I tried to uncomment them and re-run the inlist_finish.
> when log T_c = 9.97, the following warning appears:
> 
> WARNING: evaluating rates with lgT=10.001623 which is above 
> lgT=10.000000, rates have been truncated
> 
> I use the reaction rates for jina reaclib. Are the rates in reaclib
> only include the the part for T < 10 GK ?
> 
> Thanks,
> Wenyu 
> 
> > -----Original Messages-----
> > From: "Francis Timmes" <fxt44 at mac.com>
> > Sent Time: 2022-09-18 00:32:48 (Sunday)
> > To: "辛文宇" <wyxin at nao.cas.cn>, xinwenyu16 at mails.ucas.ac.cn
> > Cc: "Frank Timmes" <fxt44 at mac.com>, mesa-users <mesa-users at lists.mesastar.org>
> > Subject: Re: [Mesa-users] The nuclear burning after log T = 9.85 is stop ?
> > 
> > > ... nuclear burning in the outer shell seems to stopped.
> > 
> > inlist_finish turns off nuclear reactions for this dynamical phase of evolution:
> > 
> >       eps_nuc_factor = 0d0
> >       dxdt_nuc_factor = 0d0
> > 
> > https://docs.mesastar.org/en/release-r22.05.1/reference/controls.html#eps-nuc-factor
> > https://docs.mesastar.org/en/release-r22.05.1/reference/controls.html#dxdt-nuc-factor
> > 
> > fxt
> > 
> > 
> > 
> > 
> > > On Sep 13, 2022, at 8:04 PM, 辛文宇 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> > > 
> > > Hi, 
> > > 
> > > 
> > > 
> > > I am using  a example '25M_pre_ms_to_core_collapse' in test_suite of MESA 12778
> > > 
> > > to evolve a star form ZAMS to Fe core collapse.
> > > 
> > > Some of the inlist files are attached. The following three header files are used successively.
> > > 
> > > inlist_to_lgT_9.9_header
> > > 
> > > inlist_convert_header
> > > 
> > > inlist_finish_header
> > > 
> > > When log T_c reaches 9.85, the 'inlist_to_lgT_9.9_header' is finished and 
> > > 
> > > the 'inlist_convert_header' is continue. This file seems evolve the models forwards 100 steps.
> > > 
> > > Then the file 'inlist_finish_header' will finish the folloing evolution until the log Rho_c reaches
> > > 
> > > 10.4 g cm^-3. However, in this part, the logT and logRho increases, but the nuclear burning
> > > 
> > > in the outer shell seems to stopped. I am sure why. 
> > > 
> > > Here also attached a Kippenhahn diagram for M = 28 Msun. 
> > > 
> > > When logT_c =9.85, log t_final - t = -5, at the end of the calculation log t_final - t = -8.
> > > 
> > > It seems that from the center  to M_r=8 Msun, all the nuclear burning is stopped.
> > > 
> > > I also checked the energy generation rate of C, O and Si burning in the profileXX.data
> > > 
> > > they are 0 in this stage.
> > > 
> > > 
> > > 
> > > Thanks,
> > > 
> > > Wenyu Xin
> > > 
> > > <inlist.zip><28M_kipp.pdf>_______________________________________________
> > > mesa-users at lists.mesastar.org
> > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > > 
> </inlist.zip></mesa-users at lists.mesastar.org></mesa-users at lists.mesastar.org></fxt44 at mac.com></wyxin at nao.cas.cn></fxt44 at mac.com>