[Mesa-users] [EXT] Re: PISN models for high metallicity
Ebraheem Farag
ekfarag at asu.edu
Sun Jul 14 23:44:33 UTC 2024
Hi Wenyu,
I was able to reproduce the issue you encountered. It seems r12778 has an
interesting quark in that if the atmosphere builder returns a negative
luminosity, MESA will draw an error and terminate the model. This is not
the case in r15140, and instead the model will retry the step to achieve a
solution. (I'm not sure if this was or wasn't the case in r11701, but it's
an opportunity for you to investigate).
For this reason, I decided to move to r15140 to demonstrate that a high
metallicity model is possible. I was a little confused by your inlists as
inlist_before is almost identical to inlist_ppisn and I wasn't always sure
when each inlist was being read, so I decided to return to a more
simplified setup for this example. I copied over
"$MESA_DIR/star/test_suite/ppisn", uncommented a couple things in the
run_star extras and added a few lines at the beginning of inlist_ppisn. To
see the changes between this model and the test_suite, you;ll want to do a
diff on inlist_ppisn, and run_star_extras.f90. I turned off mass loss, set
the initial mass to 120Msun, set the initial metallicity to z=0.005,
Y=0.995, (and associated Zbase), then evolved the model until it exploded
as a pair instability supernova as expected with the following termination
code:
" Entire star is expanding above the escape velocity, PISN!"
I've attached a clean model directory and pgstar movie. Note the surface is
quite dense, this is because Pextra_factor = 2 and also because I prevented
remeshing in the surface zones during hydro, to help with convergence. You
can still try switching to zoning in log_tau and see what happens (in which
case you should return Pextra_factor to ''1"). It would be interesting for
you to try running this model with Pextra_factor = 1 or 1.5 so you can see
how it affects the profile of the surface regions (pgstar_flag = .true. is
always recommended).
To compare to your inlists for r12778, you'll want to do a diff between the
inlist_ppisn and run_star_extras.f90, and inlist_extra your inlists. This
is just a test_suite model, so probably many things should be tightened
including the mesh/timesteps and other controls you previously adopted.
Hopefully this setup in combination with yours can help you experiment with
what does and doesn't work, and you can go off and explore some interesting
science :).
p.s. this test_suite example was run with rotation on (w/wcrit = 0.1), but
it should run similarly with rotation off.
best of luck,
EbF
On Sat, Jul 6, 2024 at 11:49 AM 辛文宇 <xinwenyu16 at mails.ucas.ac.cn> wrote:
> I forget to attach the new log files for the output.
>
>
> -----原始邮件-----
> *发件人:*"辛文宇" <xinwenyu16 at mails.ucas.ac.cn>
> *发送时间:*2024-07-07 00:55:28 (星期日)
> *收件人:* "Ebraheem Farag" <ekfarag at asu.edu>
> *抄送:* mesa-users <mesa-users at lists.mesastar.org>, "Mathieu Renzo" <
> mrenzo at arizona.edu>, "Frank Timmes" <fxtimmes at gmail.com>
> *主题:* Re: [EXT] Re: PISN models for high metallicity
>
> Hi EbF,
>
>
> Thanks for your so many suggestions very much.
>
> I testing all of your comments.
>
> 1. turn off the mlt++ in inlist_hydro_on, failure.
>
> 2. set Pextra_factor = 1 in inlist_hydro_on, failure.
>
> 3. try switching from zoning in logR to zoning to logtau, failure.
>
> 4. okay_to_remesh = .false. in inlist_hydro_on.
>
> The Riemann hydro is turn on in step 2205. After some steps of retries,
> log dt decrease to -6.
>
> Then after ~200 steps, the crash takes place again.
>
> I attached the logRho and entropy near the surface. There are small
> changes.
>
> The inlist_hydro_on is also attached. Only this file is changed.
>
>
> Do you have any other comments?
>
>
> Thanks,
>
> Wenyu
>
>
> -----原始邮件-----
> *发件人:*"Ebraheem Farag" <ekfarag at asu.edu>
> *发送时间:*2024-07-06 01:19:07 (星期六)
> *收件人:* "Mathieu Renzo" <mrenzo at arizona.edu>
> *抄送:* mesa-users at lists.mesastar.org, xinwenyu16 at mails.ucas.ac.cn
> *主题:* Re: [Mesa-users] [EXT] Re: PISN models for high metallicity
>
> Hi Wenyu,
>
> Looking at your inlists and your "120M.log" file I have a few comments.
>
> It appears your model is crashing immediately after Riemann hydro turns
> on. Riemann hydro turns on at timestep 2207, and the model crashes at
> timestep 2208. Check that mlt++ is turn off entirely when Riemann hydro
> turns on (perhaps by placing controls in 'inlist_hydro_on') as you are
> still reading inlist_before at logTcen ~ 9d0.
>
> Try returning Pextra_factor to 1 when Riemann hydro turns on if you think
> that will help.
>
> Try turning off remeshing for some few time steps to see if that helps
> "okay_to_remesh = .false.", this will smooth the transition from cell face
> to cell center velocities.
>
> As a last resort, try switching from zoning in logR to zoning
> to logtau --> set :
> split_merge_amr_log_zoning = .false.
> split_merge_amr_logtau_zoning = .true. ! adjust
> 'split_merge_amr_nz_baseline' accordingly.
>
> Happy debugging :)
>
> -EbF
>
>
>
>
>
>
>
> On Fri, Jul 5, 2024 at 2:47 AM Mathieu Renzo via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> Hi Wenyu,
>>
>> I agree with Eb that this smells like a surface issue -- likely opacity
>> related since it appears at large Z but not at lower Z. I suppose (also
>> based on the Kippenhahn you sent) the issue is occurring during the first
>> PPISN pulse.
>>
>> As a test, have you tried running a model with a H-rich envelope? This
>> would change quite a bit the surface properties (unless your wind setup is
>> such that at these mass and metallicity the H-rich envelope is removed by
>> radiatively driven winds before the pulses). If you find that models get
>> farther in their evolution (maybe though the first pulse) with an H rich
>> envelope, it would strengthen the suspicion that it's boundary layers that
>> are the problem.
>>
>> In that case, there are *many* things one could try to play with, from
>> outer boundary conditions, to wind mass loss rates pre-pulse, to doing the
>> re-mapping to a new star earlier. But you would need to carry out extensive
>> tests of any of these, which can impact the evolution significantly,
>> introducing numerical artifacts.
>>
>> Also, you probably are aware that at high Z PPISN is not expected (cf. Langer
>> et al. 2007
>> <https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2007A*26A...475L..19L/abstract__;JQ!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S27tiFJVA$>
>> and more recently Farmer et al. 2019
>> <https://urldefense.com/v3/__https://doi.org/10.3847/1538-4357/ab518b__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S2MagqY9A$>),
>> although one could invoke B-fields to quench the winds at high Z and bypass
>> this limit (e.g., Georgy et al. 2017
>> <https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2017A*26A...599L...5G/abstract__;JQ!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S0yHF-5Og$>).
>>
>>
>> Hope this helps, and good luck!
>> Mathieu
>>
>> On 7/5/24 10:15, 辛文宇 wrote:
>>
>>
>>
>> * External Email *
>> ------------------------------
>> Hi, fxt,
>>
>>
>> Thank you very much for giving me this paper. In fact I have noticed this
>> paper before this email.
>>
>> And I also use the similar resolution settings with this paper.
>>
>> As I have mentioned, for the models with Z<1d-3, the models works well.
>>
>> The crashes happen only for Z>.1d-3 case.
>>
>> > maybe in future postings remove the star executable before assembling a
>> zip archive of a mesa work directory (wink).
>>
>> OK, thanks for your reminder.
>>
>>
>> Hi, EbF,
>>
>>
>> Also thanks for you suggestions. I attached the logRho and entropy as a
>> function of Mr near the surface before crashes,
>>
>> as well as the mass fraction of several isotopes. I also show the
>> kippenhan diagram.
>>
>> There is a sharply increase (decrease) in logRho (entropy). Do the
>> crashes occur because this？
>>
>> I guess this sharply increase in logRho is because I set Pextra_factor =
>> 2, and I need change is as defacult at the end
>>
>> of He burning?
>>
>>
>> Thanks,
>>
>> Wenyu
>>
>>
>>
>> -----原始邮件-----
>> *发件人:*"Ebraheem Farag" <ekfarag at asu.edu> <ekfarag at asu.edu>
>> *发送时间:*2024-07-05 06:26:30 (星期五)
>> *收件人:* "Frank Timmes" <fxtimmes at gmail.com> <fxtimmes at gmail.com>
>> *抄送:* mesa-users <mesa-users at lists.mesastar.org>
>> <mesa-users at lists.mesastar.org>, "辛文宇" <xinwenyu16 at mails.ucas.ac.cn>
>> <xinwenyu16 at mails.ucas.ac.cn>
>> *主题:* Re: [Mesa-users] PISN models for high metallicity
>>
>> I can’t run 11701 or 12778 on my current machine, but i would suggest
>> looking at the surface regions of the trho profile to figure out if it it’s
>> just the surface going haywire, especially at higher metallicity. It’s good
>> to see what’s going on before the model crashes.
>>
>> -EbF
>>
>> On Thu, Jul 4, 2024 at 4:00 PM Frank Timmes via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>>
>>> hi wenyu,
>>>
>>> my first guess is mass or time resolution. there might be some useful
>>> resolution settings to consider in
>>> https://urldefense.com/v3/__https://zenodo.org/records/6930577__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJggObPg5ig$
>>> , with an accompanying article
>>> https://urldefense.com/v3/__https://ui.adsabs.harvard.edu/abs/2022ApJ...937..112F/abstract__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgiZuHIOyg$
>>> . these models were run with r11701, somewhat close to your r12778, but at
>>> a smaller metallicity.
>>>
>>> maybe in future postings remove the star executable before assembling a
>>> zip archive of a mesa work directory (wink).
>>>
>>> fxt
>>>
>>>
>>>
>>>
>>>
>>> > On Jul 2, 2024, at 3:10 AM, 辛文宇 <xinwenyu16 at mails.ucas.ac.cn> wrote:
>>> >
>>> > Hi everyone,
>>> >
>>> > I am calculating the PISN models with ppisn example in test_suite and
>>> in mesa-r12778.
>>> > For Z<1E-3, the models work well. However, when I use Z>1d-3, the
>>> models stop at the end of He burning.
>>> > I set okay_to_reduce_gradT_excess = .true. to make the models go
>>> through pre-MS, which is turn off at
>>> > power_he_burn > 0.9 L_surf. I also use Pextra_factor = 2 to ensure the
>>> model go through He burning.
>>> > However, when T_c reaches 8.988, there are some retries like:
>>> > hydro_mtx: logRho too large 2208 1258 1
>>> 1.2040579673786977D+01 4.2901366099600260D+00
>>> > hydro_mtx: logRho too large 2208
>>> 1335 1 1.2588366322584088D+01 4.3060459950999634D+00
>>> > hydro_mtx: logRho too large 2208
>>> 1258 1 1.2040579673786977D+01 4.2901366099600260D+00
>>> > hydro_mtx: logRho too large 2208
>>> 1335 1 1.2588366322584088D+01 4.3060459950999634D+00
>>> > set_vars_for_solver failed: model, nz 2208
>>> 1533
>>> > eval_equations: set_newton_vars returned ierr -1
>>> > Finally, the error are as follows:
>>> > Error termination. Backtrace:
>>> >
>>> > Could not print backtrace: unrecognized DWARF version in .debug_info
>>> at 6
>>> > #0 0xa90697 in __utils_lib_MOD_mesa_error
>>> > at ../public/utils_lib.f:1021
>>> > #1 0x4f47e2 in get_t_tau
>>> > at ../private/atm_support.f90:310
>>> > #2 0x4f70cc in __atm_support_MOD_get_atm_pt
>>> > at ../private/atm_support.f90:127
>>> > #3 0x50ab78 in __hydro_vars_MOD_get_surf_pt
>>> > at ../private/hydro_vars.f90:1459
>>> > #4 0x50c068 in __hydro_vars_MOD_set_teff_info_for_eqns
>>> > at ../private/hydro_vars.f90:936
>>> > #5 0x8599e5 in pt_eqns_surf
>>> > at ../private/hydro_eqns.f90:2372
>>> > #6 0x8599e5 in eval_equ_for_solver
>>> > at ../private/hydro_eqns.f90:376
>>> > #7 0x85b5bf in __hydro_eqns_MOD_eval_equ
>>> > at ../private/hydro_eqns.f90:64
>>> > #8 0x81f640 in __hydro_newton_procs_MOD_eval_equations
>>> > at ../private/hydro_newton_procs.f90:143
>>> > #9 0x7f82d7 in adjust_correction
>>> > at ../private/star_newton.f90:955
>>> > #10 0x7f82d7 in do_newton
>>> > at ../private/star_newton.f90:631
>>> > #11 0x7f82d7 in __star_newton_MOD_newton
>>> > at ../private/star_newton.f90:183
>>> > #12 0x799795 in newt
>>> > at ../private/solve_hydro.f90:830
>>> > #13 0x799f0d in hydro_newton_step
>>> > at ../private/solve_hydro.f90:766
>>> > #14 0x79bd7e in do_hydro_newton
>>> > at ../private/solve_hydro.f90:163
>>> > #15 0x79bd7e in __solve_hydro_MOD_do_hydro_converge
>>> > at ../private/solve_hydro.f90:163
>>> > #16 0x79eb54 in __struct_burn_mix_MOD_do_struct_burn_mix
>>> > at ../private/struct_burn_mix.f90:207
>>> > #17 0x5d793e in __evolve_MOD_do_evolve_step_part2
>>> > at ../private/evolve.f90:599
>>> > #18 0x4220fd in __star_lib_MOD_star_evolve_step
>>> > at ../public/star_lib.f90:550
>>> > #19 0x43d091 in __run_star_support_MOD_run1_star
>>> > at ../job/run_star_support.f90:130
>>> > #20 0x41b67f in __run_star_MOD_do_run_star
>>> > at
>>> /home/wyxin/work/mesa/12778/mesa-r12778/star/job/run_star.f:26
>>> > #21 0x41b746 in run
>>> > at ../src/run.f:16
>>> > #22 0x41b788 in main
>>> > at ../src/run.f:2
>>> > 2.5691E+05 21.505411 7.756651 7.400769 7.643802 -5.045201
>>> 0.000000 0.050090 0.087833 0.000177 9.254E-01 -0.221E-01
>>> hold
>>> > log_rel_run_E_err -4.7322703837780633
>>> > rel_E_err -8.2532431125648629D-07
>>> >
>>> > check gamma integral 2.6093959149115335E-002
>>> > newton solver rejected trial model
>>> > s% model_number 2210
>>> > s% hydro_call_number 2953
>>> > nz 2396
>>> > s% num_retries 4
>>> > s% num_backups 0
>>> > s% number_of_backups_in_a_row 0
>>> > log dt/secyer
>>> -3.3598229247300733D+00
>>> >
>>> > retry 2210
>>> > retry log10(dt/yr), log10(dt), retry_factor 1 2210
>>> -3.961883 3.537229 0.500000
>>> > get_T_tau: L <= 0 2210
>>> -5.6720755514084073D+42
>>> > File: ../private/atm_support.f90, Line: 310
>>> > DATE: 2024-07-02
>>> > TIME: 08:57:57
>>> >
>>> > I am not sure how to make the models go through the further evolution.
>>> > Is there anyone can give me some suggestions?
>>> > The output file and inlist files are attached.
>>> >
>>> > Thanks,
>>> > Wenyu
>>> > <ppisn2.zip>_______________________________________________
>>> > mesa-users at lists.mesastar.org
>>> >
>>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgjNV7oIZg$
>>> >
>>>
>>> _______________________________________________
>>> mesa-users at lists.mesastar.org
>>>
>>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!bQd8zAP4sa19edBT_5ot5y-ETo-NtUApGFIVwnCdA6mSRTndGfv0e4ehJKYCLY2Tpmgdy2zyO907bq5vgdSWJgjNV7oIZg$
>>>
>>>
>> _______________________________________________ mesa-users at lists.mesastar.org https://lists.mesastar.org/mailman/listinfo/mesa-users <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S1vjR60vg$>
>>
>>
>> --
>> Mathieu Renzo
>> <https://urldefense.com/v3/__https://users.flatironinstitute.org/*mrenzo/__;fg!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S0IR04qPQ$>
>> Assistant Professor
>> University of Arizona & Steward Observatory
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>>
>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!eJty_qjigdgXhcYFpBCRby3xzoKN9OJCcn86gVf0wSVpe_e_Kb0JU25_YIHhQqkaapsl5pq1BAmkWVxes_6o8S1vjR60vg$
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240714/313aad0b/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ppisn_15140_120M.tar.gz
Type: application/x-gzip
Size: 34847 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240714/313aad0b/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 120Msun_0.005Z_ppisn.mp4
Type: video/mp4
Size: 3726030 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240714/313aad0b/attachment-0001.mp4>
More information about the Mesa-users
mailing list