[Mesa-users] PISN models for high metallicity

Tue Jul 2 10:10:10 UTC 2024

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240702/4ec28185/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ppisn2.zip
Type: application/x-zip-compressed
Size: 16658978 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240702/4ec28185/attachment-0001.bin>