[Mesa-users] error when accreting material similar to earth(high z)

[Francis Timmes]

perhaps take a look at 

https://ui.adsabs.harvard.edu/abs/2023MNRAS.518.5465B/abstract

and its related zenodo entry

https://doi.org/10.5281/zenodo.6946553

for some clues and hints.

also, https://ui.adsabs.harvard.edu/abs/2018A%26A...618A.132K/abstract could be useful.

if you post a minimal working example (inlists, any run_star_extras, etc), then other mesa_user might take to look at it in more detail, and possibly run it to help diagnose. 

fxt






> On Aug 7, 2025, at 7:35 AM, 李宗恒 via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> hello mesa_users:
> following is my issue
> mesa-r23051
> mesasdk-x86_64-linux-23.7.3.tar
> when I try to make the sun-like star to accrete earth-like materials(several earth-masses),some error come out
> retry: adjust_correction failed in eval_equations -- give up in solver 159
> retry: get_T_tau -- L <= 0 226
> the following is the composition I used
> if (s% star_age >1d7 .and. s% mstar < m_limit .and. s% x_logical_ctrl(3))then
>    write(,)"accreting**"
>    s% accrete_same_as_surface = .false.
>    s% accrete_given_mass_fractions = .true.
>    s% num_accretion_species = 4
>    ! s% accretion_species_id(1) = 'li7'
>    ! s% accretion_species_xa(1) = 0.5
>    s% accretion_species_id(1) = 'o16'
>    s% accretion_species_xa(1) = 0.3
>    s% accretion_species_id(2) = 'mg24'
>    s% accretion_species_xa(2) = 0.2
>    s% accretion_species_id(3) = 'si28'
>    s% accretion_species_xa(3) = 0.2
>    s% accretion_species_id(4) = 'fe56'
>    s% accretion_species_xa(4) = 0.3
>    s% mass_change = 1d-11
>    write(,)"mstar=,m_limit=",s% mstar,m_limit
>    write(,)"accreting*"
> else
>    write(,)"mstar=",s% mstar
>    s% mass_change = 0
> endif
> 
> But if I use the abundance below ,everything is fine
> s% num_accretion_species = 4
>       s% accretion_species_id(1) = 'h1'
>       s% accretion_species_xa(1) = 0.75
>       s% accretion_species_id(2) = 'he4'
>       s% accretion_species_xa(2) = 0.2
>       s% accretion_species_id(3) = 'li7'
>       s% accretion_species_xa(3) = 1d-9
>       s% accretion_species_id(4) = 'fe56'
>       s% accretion_species_xa(4) = 0.04999999d0
>       s% mass_change = 1d-11         
> Linux LL 6.6.87.2-microsoft-standard-WSL2 #1 SMP PREEMPT_DYNAMIC Thu Jun  5 18:30:46 UTC 2025 x86_64 x86_64 x86_64 GNU/Linux
> 
> li at LL:~$ gfortran -v
> Using built-in specs.
> COLLECT_GCC=/home/li/mesasdk/bin/gfortran.exec
> COLLECT_LTO_WRAPPER=/home/li/mesasdk/bin/../libexec/gcc/x86_64-pc-linux-gnu/13.1.0/lto-wrapper
> Target: x86_64-pc-linux-gnu
> Configured with: /home/user/sdk2-tmp/build/gcc/configure CC= CXX= --build=x86_64-pc-linux-gnu --host=x86_64-pc-linux-gnu --target=x86_64-pc-linux-gnu --prefix=/home/user/sdk2-tmp/mesasdk --with-gmp=/home/user/sdk2-tmp/mesasdk --with-mpfr=/home/user/sdk2-tmp/mesasdk --with-mpc=/home/user/sdk2-tmp/mesasdk --enable-languages=c,c++,fortran --disable-multilib --disable-nls --disable-libsanitizer --enable-clocale=generic
> Thread model: posix
> Supported LTO compression algorithms: zlibgcc version 13.1.0 (GCC)
> 
> Thanks a lot
> Jason Lee
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>