[Mesa-users] relax_initial_entropy error

[Charles Gibson]

Hello,

I am currently working on evolving a star involved in a collision from an
SPH simulation. I was able to export the composition, entropy, and angular
momentum data from the star so I could use those files to create a MESA
model. I made sure to follow the formatting of the three files based on
what I found in this &star_job link
<https://docs.mesastar.org/en/latest/reference/star_job.html#relax-entropy-filename>.
The composition and angular momentum files were both read in, but MESA
keeps breaking when it hits the entropy file.

The entropy file provides a density, specific internal energy pair. The
values for the density of the outer layers are pretty low at ~10^-8, and
the specific internal energy values are at ~10^10, and none of the values
are negative.

The related settings in my inlist look like this:

&star_job

    get_entropy_for_relax_from_eos = 'eosDE'

    relax_initial_composition = .true.
    relax_composition_filename = 'composition.dat'
    relax_initial_angular_momentum = .true.
    relax_angular_momentum_filename = 'angular_momentum.dat'
    relax_initial_entropy = .true.
    relax_entropy_filename = 'entropy.dat'

When I run MESA, it gives me the terminal output:

                                  relax_initial_entropy
 failed in eosDT_get_T (as eosDE)
 do_relax_initial_entropy ierr          -1
 do_star_job_controls_after ierr          -1
 do_before_evolve_loop ierr          -1
 do_before_evolve_loop ierr          -1

The version of MESA that I am running is 22.05.1.

How might I locate the error within the model/file? Is there an inlist
parameter that I am missing, or are the values in the entropy file too
small to be accurately used to determine entropy? In the past, I have
successfully used MESA to evolve stellar models from this SPH code that
haven't experienced any collisions, so I know that, at least in those
cases, the inlist settings are the correct. Any help would be much
appreciated!

Thanks!
Charlie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230531/6d4c9b29/attachment.htm>