[Mesa-users] Blue and red edges of the instability strip

[Bill Paxton]

You might also want to look at the rsp_check_2nd_crossing test case in star/test_suite.

Here’s the comment you’ll find at the start of inlist_ rsp_check_2nd_crossing

! given a history log constructed by a mesa/star run
! such as in the 5M_cepheid_blue_loop test case
! (e.g., history_logs/history_7M.data),
! uses the rsp model building and linear analysis
! to find the instability strip edges,
! then interpolates along the evolutionary track
! to find points at selected Teff offsets from the blue edge.
! the results are written to a file such as 7B_2nd_crossing.data

Dig into the details in run_star_extras to see exactly how it works.

As it stands it is only looking for the blue edge, but the approach might be extended for red edge as well. 

-Bill



> On Jun 15, 2020, at 11:00 PM, Radek Smolec via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Hi Felipe,
> 
> RSP will not produce the edges of the instability strip at one shot, but for every RSP model you run, it returns linear growth rates. As an example, please run test_suite/rsp_Cepheid/
> in the LOGS directory you will find LINA_period_growth.data in which linear periods (second column) and growth rates (third) are located. In the first column you have the mode id (0 - fundamental, 1 - first ov, ...)
> 
> Positive growth rate means that the mode is linearly unstable (inside linear IS). Linear instability edges are lines along which the growth rate is zero.
> If you run a sequence or a grid of models you can interpolate then. Note that linear growth rates are calculated at the very beginning of rsp run (for static model). There is no need to converge the model to full amplitude pulsation if you are interested in growth rates only.
> Best wishes,
> Radek
> 
> 
> wt., 16 cze 2020 o 04:52 Felipe Espinoza Arancibia via Mesa-users <mesa-users at lists.mesastar.org <mailto:mesa-users at lists.mesastar.org>> napisał(a):
> Hi all, 
> I'm a Master's student and I'm using MESA in a project with Cepheids. I wanted to ask, in figure 13 of instrumental paper 5, how can I calculate the edges of the instability strips with RSP?
> I found the input files of the paper in Zenodo but I didn't find anything about that figure. 
> 
> Regards, 
> Felipe.
> _______________________________________________
> mesa-users at lists.mesastar.org <mailto:mesa-users at lists.mesastar.org>
> https://lists.mesastar.org/mailman/listinfo/mesa-users <https://lists.mesastar.org/mailman/listinfo/mesa-users>
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200616/a57e5aa6/attachment.htm>