[Mesa-users] Test Suite: 1M_pre_ms_to_wd

Warrick Ball W.H.Ball at bham.ac.uk
Thu Nov 18 21:07:37 UTC 2021


Hi Erica,

What makes you say the star has reached the end of the main sequence?  Your figure isn't completely clear but the HR diagram inset presumably shows the pre-main-sequence evolution, rather than post-MS evolution.

Cheers,
Warrick

___________

Warrick Ball
Postdoc, School of Physics and Astronomy
University of Birmingham, Edgbaston, Birmingham B15 2TT
W.H.Ball at bham.ac.uk
+44 (0)121 414 4552

On Thu, 18 Nov 2021, Erica A Sawczynec wrote:

> Hi Aaron, Even with commenting out the max_model_number in the inlist I am still having the same issue. The star appears to have reached the end of the MS under 1 Gyr. Any other ideas?
> Thanks,
> Erica
>
>       On Nov 18, 2021, at 2:34 PM, Aaron Dotter via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Hi Erica,
> With a max_years_for_timestep = 1d5 and max_model_number = 6500 your simulation will end with an age of ~650 Myr.  
> 
> That does not mean it has reached any specific evolutionary point, just that it has reached the max_model_number and that is one criterion for MESA to terminate the evolution. 
> 
> A 1 Msun model (at solar Z) should have a main sequence lifetime of ~10 Gyr so you would need about 10^5 timesteps to reach the end of the main sequence with a max_years_for_timestep = 1d5.
> 
> Aaron
> 
> On Thu, Nov 18, 2021 at 3:24 PM Erica A Sawczynec <ericasaw at utexas.edu> wrote:
>       Hi Aaron, Will this cause the star to evolve to the MS turn-off in 6500 steps regardless of if it should be there or not? I tried changing this value to be bigger and ran into the same issue. Here is a plot I made using TULIPS of
>       the radius as a function of time. You can see the star is at the MS turn-off at 638 Myr. Trying a run with turning that off now!
> Thanks,
> Erica
> 
> <600Myr_turnoff.png>
>
>       On Nov 18, 2021, at 2:17 PM, Aaron Dotter <aaron.dotter at gmail.com> wrote:
> 
> Hi Erica,
> You have max_model_number = 6500 in the inlist.  Try removing this from the inlist (or set to -1).
> 
> Aaron
> 
> On Thu, Nov 18, 2021 at 3:08 PM Erica A Sawczynec <ericasaw at utexas.edu> wrote:
>       Hi All! I have been tinkering with the 1M_pre_ms_to_wd test suite to try to look at the depth of the convection zone boundary as a function of opacity. I have an altered version that no longer uses the to_end_wd inlist
>       and I set the envelope mass limit so the to_end_agb inlist doesn’t actually go too far past the MS turn-off. I want to finely sample the  MS so I tried turning on the max_years_for_timestep keyword in the debugging
>       section. I set it to something small like 1d5, which feels reasonable for the MS, but I’ve noticed whenever I do this, the test suite evolves a 1M_sol star to the end of the MS in ~600 Myr (which isn’t possible). 
> I thought that the max_years_for_timestep would keep the time to evolve to the MS turn-off the same, but increase the sampling of the MS, is that not the case? I even tried increasing the mass model number (up to 65000),
> thinking maybe that imposed some sort of restriction on the model, but even then the age of the star at the end of the run is way too young. I’ve attached my changed inlist below, note that it will “fail” with the "failed to
> create end_agb.mod when running inlist_to_end_agb_header” error, but there is still a history.data file where you can look at the results of the run.
> I cannot help but feeling I am missing something silly here, but I just can’t figure out what it is. Any help is greatly appreciated!
> Thanks,
> Erica
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 
> 
> 
>


More information about the Mesa-users mailing list