[Mesa-users] PISN models for high metallicity

[Frank Timmes]

hi wenyu,

my first guess is mass or time resolution. there might be some useful resolution settings to consider in https://zenodo.org/records/6930577 ,  with an accompanying article https://ui.adsabs.harvard.edu/abs/2022ApJ...937..112F/abstract . 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://lists.mesastar.org/mailman/listinfo/mesa-users
>