[Mesa-users] Query about relax_core in MESA
Abha Vishwakarma
abha.vishwakarma at niser.ac.in
Sat Mar 2 07:58:56 UTC 2024
Dear Ebraheem,
Thanks a lot for your email! This is exactly what was confusing me. Thanks
again for sharing the directory!
Sincerely,
Abha
On Wed, Feb 28, 2024 at 12:44 AM Ebraheem Farag <ekfarag at asu.edu> wrote:
> Hello Abha,
>
> Not sure what your goal is but, I've attached a compressed
> $MESA_DIR/star/work directory where I:
> - evolved a 12 Msun from pre-ms to zams, saved the model.
> - loaded in the model and excised 6 Msun of the core and then saved the
> resulting model.
>
> If you follow this approach, you should be able to generate whatever
> desired excised model you have in mind and then try loading it into the
> hse_rieman test.
>
> Hopefully this helps get you off the ground...
> (untar with "tar -xzvf work_test.tar.gz")
>
> -EbF
>
> On Mon, Feb 26, 2024 at 2:44 PM Earl Bellinger via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> Dear Abha,
>>
>> I am not sure. I have not used the hse_riemann test case before. Maybe
>> some others can chime in here.
>>
>> You are right that the hse_riemann test case seems to have a similar
>> setup. Its initial model (removed.mod) seems to be generated
>> using gyre_in_mesa_envelope, so perhaps check the sequence of steps there
>> and see if you can use anything from it. Also, you might want to compare
>> your model with removed.mod, to see if there are any glaring differences
>> that may be triggering your roadblock.
>>
>> P. S. thanks for the praise on shmesa! 😊
>>
>> Best wishes,
>> Earl
>>
>>
>> On Mon, Feb 26, 2024 at 8:02 AM Abha Vishwakarma via Mesa-users <
>> mesa-users at lists.mesastar.org> wrote:
>>
>>> Dear Dr. Earl,
>>>
>>> I am trying to reproduce a model of a 12 solar mass star with the core
>>> excised, which is used as an input in the 'hse_riemann' test suite in the
>>> star module of MESA. I have attached the 'removed.mod' file from the
>>> $MESA_DIR/star/test_suite/hse_riemann directory in this email for your
>>> reference.
>>>
>>> I found the details of the core mass and other core properties by going
>>> through this file, and I tried setting the value of the new core properties
>>> accordingly. I am not sure about the Chandrashekhar limit here. Perhaps I
>>> am misunderstanding the function of the relaxation routine related
>>> parameters. I thought that by setting the new_core_mass to the core mass
>>> that I found in the 'removed.mod' file, MESA would define the boundary of
>>> the core accordingly, but I think I might be wrong. Could you please
>>> clarify?
>>>
>>> Another thought process is that maybe I should first evolve a 12 solar
>>> mass star upto TAMS, save the model, and then do a relaxation on this using
>>> another inlist. Sorry, perhaps this is very trivial, but any clues will be
>>> very helpful! Also, shmesa is great, thanks a lot for that too!
>>>
>>> Sincerely,
>>> Abha
>>>
>>> On Sun, Feb 25, 2024 at 5:30 PM <mesa-users-request at lists.mesastar.org>
>>> wrote:
>>>
>>>> Send Mesa-users mailing list submissions to
>>>> mesa-users at lists.mesastar.org
>>>>
>>>> To subscribe or unsubscribe via the World Wide Web, visit
>>>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> or, via email, send a message with subject or body 'help' to
>>>> mesa-users-request at lists.mesastar.org
>>>>
>>>> You can reach the person managing the list at
>>>> mesa-users-owner at lists.mesastar.org
>>>>
>>>> When replying, please edit your Subject line so it is more specific
>>>> than "Re: Contents of Mesa-users digest..."
>>>>
>>>>
>>>> Today's Topics:
>>>>
>>>> 1. Re: Query about relax_core in MESA (Abha Vishwakarma)
>>>> 2. Re: Query about relax_core in MESA (Earl Bellinger)
>>>>
>>>>
>>>> ----------------------------------------------------------------------
>>>>
>>>> Message: 1
>>>> Date: Sun, 25 Feb 2024 07:28:55 +0530
>>>> From: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> To: mesa-users at lists.mesastar.org
>>>> Subject: Re: [Mesa-users] Query about relax_core in MESA
>>>> Message-ID:
>>>> <CAGU_p53irmUyZMhR8xKY8LdGXiY=
>>>> jVR7uwn+RaLqV8dG7zf+qA at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Dear Dr. Earl,
>>>>
>>>> Apologies for the delay in the reply, and thanks a lot for the help! I
>>>> am
>>>> trying to use the relax_core parameter to define a core of a 12 solar
>>>> mass
>>>> star, and if it's possible, it would help me if you could give me some
>>>> clue
>>>> to fix an error that I keep encountering in the process.
>>>>
>>>> The terminal outputs 'finished doing relax core', but proceeds to show
>>>> some
>>>> timestepping error and outputs 'retry: logT <
>>>> hydro_mtx_min_allowed_logT'
>>>> several times and stops. I have attached this error and the
>>>> inlist corresponding to it in the email.
>>>>
>>>> I tried tweaking the 'dlg_core_mass_per_step' parameter to change the
>>>> mass
>>>> stepping, and I again encountered errors regarding the
>>>> 'hydro_mtx_min_allowed_logT' and 'hydro_mtx_max_allowed_logT' limits. It
>>>> might be something trivial, but since I'm new to this I am having
>>>> trouble
>>>> figuring it out.
>>>>
>>>> Thanks,
>>>> Abha
>>>>
>>>> On Thu, Feb 15, 2024, 5:31?PM <mesa-users-request at lists.mesastar.org>
>>>> wrote:
>>>>
>>>> > Send Mesa-users mailing list submissions to
>>>> > mesa-users at lists.mesastar.org
>>>> >
>>>> > To subscribe or unsubscribe via the World Wide Web, visit
>>>> > https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> > or, via email, send a message with subject or body 'help' to
>>>> > mesa-users-request at lists.mesastar.org
>>>> >
>>>> > You can reach the person managing the list at
>>>> > mesa-users-owner at lists.mesastar.org
>>>> >
>>>> > When replying, please edit your Subject line so it is more specific
>>>> > than "Re: Contents of Mesa-users digest..."
>>>> >
>>>> >
>>>> > Today's Topics:
>>>> >
>>>> > 1. Query about relax_core in MESA (Abha Vishwakarma)
>>>> > 2. Re: Query about relax_core in MESA (Earl Bellinger)
>>>> > 3. A question about internal energy in MESA (Yamaguchi, Natsuko)
>>>> >
>>>> >
>>>> > ----------------------------------------------------------------------
>>>> >
>>>> > Message: 1
>>>> > Date: Wed, 14 Feb 2024 21:54:00 +0530
>>>> > From: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> > To: mesa-users at lists.mesastar.org
>>>> > Subject: [Mesa-users] Query about relax_core in MESA
>>>> > Message-ID:
>>>> > <
>>>> > CAGU_p53e7zkHcrmZVKBtoMJ-RRfs-4B2RkP-_-UuF7oJSupmRg at mail.gmail.com>
>>>> > Content-Type: text/plain; charset="utf-8"
>>>> >
>>>> > Dear Mesa Users,
>>>> >
>>>> > I am new to MESA, and I was curious about the parameter 'relax_core'
>>>> in
>>>> > MESA, that I must switch to true if I want a non-zero core
>>>> mass(please let
>>>> > me know if I'm wrong). I looked up the MESA papers, the documentation
>>>> and
>>>> > older MESA Users emails as well, but I still struggle to understand
>>>> what
>>>> > exactly happens once I set this parameter to true. Is this the same
>>>> as the
>>>> > 'relaxation routines' that I find in other sections of the papers as
>>>> well?
>>>> > If so, what exactly do relaxation routines do? I would like some
>>>> insight
>>>> > into it.
>>>> >
>>>> > Sincerely,
>>>> > Abha Vishwakarma
>>>> >
>>>> > -----------------------------------------------------------------
>>>> > *Abha Vishwakarma*
>>>> > Final Year Int. M.Sc.
>>>> > NISER, India
>>>> > Mail: abha.vishwakarma at niser.ac.in
>>>> > Web: sites.google.com/view/abhavishwakarma/
>>>> <https://urldefense.com/v3/__http://sites.google.com/view/abhavishwakarma/__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ykZhl3w$>
>>>> > -------------- next part --------------
>>>> > An HTML attachment was scrubbed...
>>>> > URL: <
>>>> >
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_580XtdlnA$>
>>>> > >
>>>> >
>>>> > ------------------------------
>>>> >
>>>> > Message: 2
>>>> > Date: Wed, 14 Feb 2024 11:53:26 -0500
>>>> > From: Earl Bellinger <earlbellinger at gmail.com>
>>>> > To: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> > Cc: mesa-users at lists.mesastar.org
>>>> > Subject: Re: [Mesa-users] Query about relax_core in MESA
>>>> > Message-ID:
>>>> > <CAD7H7L0tL1P9LD30dRjxmfyLEV5fAsA_auUdhC=
>>>> > 3+a7VNM1EAg at mail.gmail.com>
>>>> > Content-Type: text/plain; charset="utf-8"
>>>> >
>>>> > Dear Abha,
>>>> >
>>>> > Yes, you can use relax_core to change the mass of the inner boundary
>>>> to be
>>>> > nonzero.
>>>> >
>>>> > Applications include modelling planets with rocky interiors:
>>>> > https://docs.mesastar.org/en/latest/test_suite/make_planets.html
>>>> <https://urldefense.com/v3/__https://docs.mesastar.org/en/latest/test_suite/make_planets.html__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_581H76LrA$>
>>>> > modelling red supergiants with neutron star cores (Thorne??ytkow
>>>> objects):
>>>> > https://arxiv.org/abs/2305.07337
>>>> <https://urldefense.com/v3/__https://arxiv.org/abs/2305.07337__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ubyGR3g$>
>>>> > and modelling stars with primordial black holes at their center
>>>> (Hawking
>>>> > stars):
>>>> > https://github.com/earlbellinger/black-hole-sun
>>>> <https://urldefense.com/v3/__https://github.com/earlbellinger/black-hole-sun__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_582p-EthA$>
>>>> > among others.
>>>> >
>>>> > The relaxation routines iteratively change the properties of the
>>>> stellar
>>>> > model, and recompute the model after each iteration. This helps with
>>>> > convergence, as jumping from e.g. a model with zero enclosed mass at
>>>> the
>>>> > center to one with a large enclosed mass may require a large
>>>> adjustment to
>>>> > the surrounding stellar structure, so it's best to do it in a
>>>> sequence of
>>>> > smaller steps.
>>>> >
>>>> > To look a bit under the hood at what is actually happening, we can
>>>> use the
>>>> > new shmesa <https://github.com/MESAHub/mesa/tree/main/scripts/shmesa
>>>> <https://urldefense.com/v3/__https://github.com/MESAHub/mesa/tree/main/scripts/shmesa__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58ts77Kbw$>
>>>> >
>>>> > functionality to find where this code is called:
>>>> >
>>>> > > shmesa grep relax_core
>>>> >
>>>> > mesa/mesa/star/job/run_star_support.f90: if (s% job%
>>>> relax_core
>>>> > .or. &
>>>> > mesa/mesa/star/job/run_star_support.f90: write(*, 1)
>>>> > 'relax_core', s% job% new_core_mass
>>>> > mesa/mesa/star/job/run_star_support.f90: call
>>>> star_relax_core( &
>>>> > mesa/mesa/star/job/run_star_support.f90: s% job%
>>>> > relax_core_years_for_dt, &
>>>> > mesa/mesa/star/job/run_star_support.f90: if
>>>> > (failed('star_relax_core',ierr)) return
>>>> > mesa/mesa/star/private/relax.f90: subroutine do_relax_core( &
>>>> >
>>>> >
>>>> > We see straight away that it is called in private/relax.f90.
>>>> Navigating
>>>> > there on the web:
>>>> >
>>>> >
>>>> >
>>>> https://github.com/MESAHub/mesa/blob/a965ec699b2b4978c43b0b305f40446b1aff05f0/star/private/relax.f90#L1979
>>>> <https://urldefense.com/v3/__https://github.com/MESAHub/mesa/blob/a965ec699b2b4978c43b0b305f40446b1aff05f0/star/private/relax.f90*L1979__;Iw!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_-l0eMZg$>
>>>> >
>>>> > Here you can see how all the relaxation routines are defined.
>>>> >
>>>> > I hope this is helpful.
>>>> >
>>>> > Best wishes,
>>>> > Earl
>>>> >
>>>> > ?
>>>> > *Earl Patrick Bellinger *
>>>> > Assistant Professor
>>>> > Department of Astronomy
>>>> > Yale University
>>>> > https://earlbellinger.com
>>>> <https://urldefense.com/v3/__https://earlbellinger.com__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_59enFhzFg$>
>>>> >
>>>> >
>>>> > On Wed, Feb 14, 2024 at 11:26?AM Abha Vishwakarma via Mesa-users <
>>>> > mesa-users at lists.mesastar.org> wrote:
>>>> >
>>>> > > Dear Mesa Users,
>>>> > >
>>>> > > I am new to MESA, and I was curious about the parameter
>>>> 'relax_core' in
>>>> > > MESA, that I must switch to true if I want a non-zero core
>>>> mass(please
>>>> > let
>>>> > > me know if I'm wrong). I looked up the MESA papers, the
>>>> documentation and
>>>> > > older MESA Users emails as well, but I still struggle to understand
>>>> what
>>>> > > exactly happens once I set this parameter to true. Is this the same
>>>> as
>>>> > the
>>>> > > 'relaxation routines' that I find in other sections of the papers as
>>>> > well?
>>>> > > If so, what exactly do relaxation routines do? I would like some
>>>> insight
>>>> > > into it.
>>>> > >
>>>> > > Sincerely,
>>>> > > Abha Vishwakarma
>>>> > >
>>>> > > -----------------------------------------------------------------
>>>> > > *Abha Vishwakarma*
>>>> > > Final Year Int. M.Sc.
>>>> > > NISER, India
>>>> > > Mail: abha.vishwakarma at niser.ac.in
>>>> > > Web: sites.google.com/view/abhavishwakarma/
>>>> <https://urldefense.com/v3/__http://sites.google.com/view/abhavishwakarma/__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ykZhl3w$>
>>>> > >
>>>> > >
>>>> > > _______________________________________________
>>>> > > 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!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> > >
>>>> > >
>>>> > -------------- next part --------------
>>>> > An HTML attachment was scrubbed...
>>>> > URL: <
>>>> >
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5-sFv_edg$>
>>>> > >
>>>> >
>>>> > ------------------------------
>>>> >
>>>> > Message: 3
>>>> > Date: Thu, 15 Feb 2024 03:01:20 +0000
>>>> > From: "Yamaguchi, Natsuko" <nyamaguc at caltech.edu>
>>>> > To: Evan Bauer via Mesa-users <mesa-users at lists.mesastar.org>
>>>> > Subject: [Mesa-users] A question about internal energy in MESA
>>>> > Message-ID:
>>>> > <
>>>> >
>>>> BY5PR03MB53169F2A0B85C38C06800E43DF4D2 at BY5PR03MB5316.namprd03.prod.outlook.com
>>>> > >
>>>> >
>>>> > Content-Type: text/plain; charset="iso-8859-1"
>>>> >
>>>> > Hello,
>>>> >
>>>> > I would like to clarify what exactly MESA's (specific) 'internal
>>>> energy'
>>>> > is composed of.
>>>> >
>>>> > Is it the (1) gas pressure, (2) radiation pressure, and (3)
>>>> ionization (or
>>>> > recombination) energy combined? And which species is it considering
>>>> for the
>>>> > ionization energy? It's not obvious from Paxton et al. (2018) which
>>>> > primarily describes it in the context of a simple example. If there
>>>> is a
>>>> > better source on this, please let me know.
>>>> >
>>>> > Thank you,
>>>> >
>>>> > Natsuko Yamaguchi
>>>> >
>>>> >
>>>> >
>>>> >
>>>> >
>>>> > -------------- next part --------------
>>>> > An HTML attachment was scrubbed...
>>>> > URL: <
>>>> >
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240215/15fc3794/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240215/15fc3794/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58BjdW9Jw$>
>>>> > >
>>>> >
>>>> > End of Mesa-users Digest, Vol 79, Issue 11
>>>> > ******************************************
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240225/9347d7ad/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240225/9347d7ad/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58JVlGr5Q$>
>>>> >
>>>> -------------- next part --------------
>>>> finished doing relax core
>>>>
>>>> net name basic.net
>>>> <https://urldefense.com/v3/__http://basic.net__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58kRxkehQ$>
>>>> xmstar/mstar
>>>> 3.0091653151900996D-03
>>>> xmstar (g)
>>>> 7.1801448183434444D+31
>>>> M_center (g)
>>>> 2.3789117000193169D+34
>>>> xmstar/Msun
>>>> 3.6109983782281196D-02
>>>> M_center/Msun
>>>> 1.1963890016217718D+01
>>>> R_center (cm)
>>>> 3.8032489518482001D+11
>>>> R_center/Rsun
>>>> 5.4667945261581146D+00
>>>> core density
>>>> 1.0323464061282679D-01
>>>> L_center/Lsun
>>>> 2.1996552195011729D+04
>>>> kap_option gs98
>>>> kap_CO_option gs98_co
>>>> kap_lowT_option lowT_fa05_gs98
>>>> OMP_NUM_THREADS 4
>>>>
>>>>
>>>> retry: logT < hydro_mtx_min_allowed_logT 86 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 86 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: get_T_tau -- L <= 0 1
>>>> retry: get_T_tau -- L <= 0 1
>>>> retry: get_T_tau -- L <= 0 1
>>>> retry: get_T_tau -- L <= 0 1
>>>> retry: get_T_tau -- L <= 0 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 72 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 72 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 76 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 64 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 74 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 78 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 82 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> retry: logT < hydro_mtx_min_allowed_logT 69 1
>>>> dt
>>>> 5.8780610561370854D-07
>>>> min_timestep_limit
>>>> 9.9999999999999995D-07
>>>> stopping because of problems dt < min_timestep_limit
>>>> terminated evolution: cannot find acceptable model
>>>> termination code: min_timestep_limit
>>>> pause_before_terminate: hit RETURN to continue
>>>> -------------- next part --------------
>>>> A non-text attachment was scrubbed...
>>>> Name: inlist_project
>>>> Type: application/octet-stream
>>>> Size: 2160 bytes
>>>> Desc: not available
>>>> URL: <
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240225/9347d7ad/attachment-0001.obj
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240225/9347d7ad/attachment-0001.obj__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5__hVMfsQ$>
>>>> >
>>>>
>>>> ------------------------------
>>>>
>>>> Message: 2
>>>> Date: Sat, 24 Feb 2024 21:24:45 -0500
>>>> From: Earl Bellinger <earlbellinger at gmail.com>
>>>> To: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> Cc: mesa-users at lists.mesastar.org
>>>> Subject: Re: [Mesa-users] Query about relax_core in MESA
>>>> Message-ID:
>>>> <
>>>> CAD7H7L30OoR7aAdEu8HukN6iAdk13ufA-QJMTC7dCiZB-DCTrA at mail.gmail.com>
>>>> Content-Type: text/plain; charset="utf-8"
>>>>
>>>> Dear Abha,
>>>>
>>>> It looks like you are trying to make a model with a core of 11.96 Msun
>>>> and
>>>> a total mass of 12 Msun. I think this would be in violation of the
>>>> Sch?nberg?Chandrasekhar limit. What are you trying to model?
>>>>
>>>> Best wishes,
>>>> Earl
>>>>
>>>>
>>>> On Sat, Feb 24, 2024 at 9:00?PM Abha Vishwakarma via Mesa-users <
>>>> mesa-users at lists.mesastar.org> wrote:
>>>>
>>>> > Dear Dr. Earl,
>>>> >
>>>> > Apologies for the delay in the reply, and thanks a lot for the help!
>>>> I am
>>>> > trying to use the relax_core parameter to define a core of a 12 solar
>>>> mass
>>>> > star, and if it's possible, it would help me if you could give me
>>>> some clue
>>>> > to fix an error that I keep encountering in the process.
>>>> >
>>>> > The terminal outputs 'finished doing relax core', but proceeds to show
>>>> > some timestepping error and outputs 'retry: logT <
>>>> > hydro_mtx_min_allowed_logT' several times and stops. I have attached
>>>> this
>>>> > error and the inlist corresponding to it in the email.
>>>> >
>>>> > I tried tweaking the 'dlg_core_mass_per_step' parameter to change the
>>>> mass
>>>> > stepping, and I again encountered errors regarding the
>>>> > 'hydro_mtx_min_allowed_logT' and 'hydro_mtx_max_allowed_logT' limits.
>>>> It
>>>> > might be something trivial, but since I'm new to this I am having
>>>> trouble
>>>> > figuring it out.
>>>> >
>>>> > Thanks,
>>>> > Abha
>>>> >
>>>> > On Thu, Feb 15, 2024, 5:31?PM <mesa-users-request at lists.mesastar.org>
>>>> > wrote:
>>>> >
>>>> >> Send Mesa-users mailing list submissions to
>>>> >> mesa-users at lists.mesastar.org
>>>> >>
>>>> >> To subscribe or unsubscribe via the World Wide Web, visit
>>>> >> https://lists.mesastar.org/mailman/listinfo/mesa-users
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> >> or, via email, send a message with subject or body 'help' to
>>>> >> mesa-users-request at lists.mesastar.org
>>>> >>
>>>> >> You can reach the person managing the list at
>>>> >> mesa-users-owner at lists.mesastar.org
>>>> >>
>>>> >> When replying, please edit your Subject line so it is more specific
>>>> >> than "Re: Contents of Mesa-users digest..."
>>>> >>
>>>> >>
>>>> >> Today's Topics:
>>>> >>
>>>> >> 1. Query about relax_core in MESA (Abha Vishwakarma)
>>>> >> 2. Re: Query about relax_core in MESA (Earl Bellinger)
>>>> >> 3. A question about internal energy in MESA (Yamaguchi, Natsuko)
>>>> >>
>>>> >>
>>>> >>
>>>> ----------------------------------------------------------------------
>>>> >>
>>>> >> Message: 1
>>>> >> Date: Wed, 14 Feb 2024 21:54:00 +0530
>>>> >> From: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> >> To: mesa-users at lists.mesastar.org
>>>> >> Subject: [Mesa-users] Query about relax_core in MESA
>>>> >> Message-ID:
>>>> >> <
>>>> >> CAGU_p53e7zkHcrmZVKBtoMJ-RRfs-4B2RkP-_-UuF7oJSupmRg at mail.gmail.com>
>>>> >> Content-Type: text/plain; charset="utf-8"
>>>> >>
>>>> >> Dear Mesa Users,
>>>> >>
>>>> >> I am new to MESA, and I was curious about the parameter 'relax_core'
>>>> in
>>>> >> MESA, that I must switch to true if I want a non-zero core
>>>> mass(please let
>>>> >> me know if I'm wrong). I looked up the MESA papers, the
>>>> documentation and
>>>> >> older MESA Users emails as well, but I still struggle to understand
>>>> what
>>>> >> exactly happens once I set this parameter to true. Is this the same
>>>> as the
>>>> >> 'relaxation routines' that I find in other sections of the papers as
>>>> well?
>>>> >> If so, what exactly do relaxation routines do? I would like some
>>>> insight
>>>> >> into it.
>>>> >>
>>>> >> Sincerely,
>>>> >> Abha Vishwakarma
>>>> >>
>>>> >> -----------------------------------------------------------------
>>>> >> *Abha Vishwakarma*
>>>> >> Final Year Int. M.Sc.
>>>> >> NISER, India
>>>> >> Mail: abha.vishwakarma at niser.ac.in
>>>> >> Web: sites.google.com/view/abhavishwakarma/
>>>> <https://urldefense.com/v3/__http://sites.google.com/view/abhavishwakarma/__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ykZhl3w$>
>>>> >> -------------- next part --------------
>>>> >> An HTML attachment was scrubbed...
>>>> >> URL: <
>>>> >>
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_580XtdlnA$>
>>>> >> >
>>>> >>
>>>> >> ------------------------------
>>>> >>
>>>> >> Message: 2
>>>> >> Date: Wed, 14 Feb 2024 11:53:26 -0500
>>>> >> From: Earl Bellinger <earlbellinger at gmail.com>
>>>> >> To: Abha Vishwakarma <abha.vishwakarma at niser.ac.in>
>>>> >> Cc: mesa-users at lists.mesastar.org
>>>> >> Subject: Re: [Mesa-users] Query about relax_core in MESA
>>>> >> Message-ID:
>>>> >> <CAD7H7L0tL1P9LD30dRjxmfyLEV5fAsA_auUdhC=
>>>> >> 3+a7VNM1EAg at mail.gmail.com>
>>>> >> Content-Type: text/plain; charset="utf-8"
>>>> >>
>>>> >> Dear Abha,
>>>> >>
>>>> >> Yes, you can use relax_core to change the mass of the inner boundary
>>>> to be
>>>> >> nonzero.
>>>> >>
>>>> >> Applications include modelling planets with rocky interiors:
>>>> >>
>>>> https://docs.mesastar.org/en/latest/test_suite/make_planets.html
>>>> <https://urldefense.com/v3/__https://docs.mesastar.org/en/latest/test_suite/make_planets.html__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_581H76LrA$>
>>>> >> modelling red supergiants with neutron star cores (Thorne??ytkow
>>>> objects):
>>>> >> https://arxiv.org/abs/2305.07337
>>>> <https://urldefense.com/v3/__https://arxiv.org/abs/2305.07337__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ubyGR3g$>
>>>> >> and modelling stars with primordial black holes at their center
>>>> (Hawking
>>>> >> stars):
>>>> >> https://github.com/earlbellinger/black-hole-sun
>>>> <https://urldefense.com/v3/__https://github.com/earlbellinger/black-hole-sun__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_582p-EthA$>
>>>> >> among others.
>>>> >>
>>>> >> The relaxation routines iteratively change the properties of the
>>>> stellar
>>>> >> model, and recompute the model after each iteration. This helps with
>>>> >> convergence, as jumping from e.g. a model with zero enclosed mass at
>>>> the
>>>> >> center to one with a large enclosed mass may require a large
>>>> adjustment to
>>>> >> the surrounding stellar structure, so it's best to do it in a
>>>> sequence of
>>>> >> smaller steps.
>>>> >>
>>>> >> To look a bit under the hood at what is actually happening, we can
>>>> use the
>>>> >> new shmesa <https://github.com/MESAHub/mesa/tree/main/scripts/shmesa
>>>> <https://urldefense.com/v3/__https://github.com/MESAHub/mesa/tree/main/scripts/shmesa__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58ts77Kbw$>
>>>> >
>>>> >> functionality to find where this code is called:
>>>> >>
>>>> >> > shmesa grep relax_core
>>>> >>
>>>> >> mesa/mesa/star/job/run_star_support.f90: if (s% job%
>>>> relax_core
>>>> >> .or. &
>>>> >> mesa/mesa/star/job/run_star_support.f90: write(*, 1)
>>>> >> 'relax_core', s% job% new_core_mass
>>>> >> mesa/mesa/star/job/run_star_support.f90: call
>>>> star_relax_core(
>>>> >> &
>>>> >> mesa/mesa/star/job/run_star_support.f90: s% job%
>>>> >> relax_core_years_for_dt, &
>>>> >> mesa/mesa/star/job/run_star_support.f90: if
>>>> >> (failed('star_relax_core',ierr)) return
>>>> >> mesa/mesa/star/private/relax.f90: subroutine do_relax_core( &
>>>> >>
>>>> >>
>>>> >> We see straight away that it is called in private/relax.f90.
>>>> Navigating
>>>> >> there on the web:
>>>> >>
>>>> >>
>>>> >>
>>>> https://github.com/MESAHub/mesa/blob/a965ec699b2b4978c43b0b305f40446b1aff05f0/star/private/relax.f90#L1979
>>>> <https://urldefense.com/v3/__https://github.com/MESAHub/mesa/blob/a965ec699b2b4978c43b0b305f40446b1aff05f0/star/private/relax.f90*L1979__;Iw!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_-l0eMZg$>
>>>> >>
>>>> >> Here you can see how all the relaxation routines are defined.
>>>> >>
>>>> >> I hope this is helpful.
>>>> >>
>>>> >> Best wishes,
>>>> >> Earl
>>>> >>
>>>> >> ?
>>>> >> *Earl Patrick Bellinger *
>>>> >> Assistant Professor
>>>> >> Department of Astronomy
>>>> >> Yale University
>>>> >> https://earlbellinger.com
>>>> <https://urldefense.com/v3/__https://earlbellinger.com__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_59enFhzFg$>
>>>> >>
>>>> >>
>>>> >> On Wed, Feb 14, 2024 at 11:26?AM Abha Vishwakarma via Mesa-users <
>>>> >> mesa-users at lists.mesastar.org> wrote:
>>>> >>
>>>> >> > Dear Mesa Users,
>>>> >> >
>>>> >> > I am new to MESA, and I was curious about the parameter
>>>> 'relax_core' in
>>>> >> > MESA, that I must switch to true if I want a non-zero core
>>>> mass(please
>>>> >> let
>>>> >> > me know if I'm wrong). I looked up the MESA papers, the
>>>> documentation
>>>> >> and
>>>> >> > older MESA Users emails as well, but I still struggle to
>>>> understand what
>>>> >> > exactly happens once I set this parameter to true. Is this the
>>>> same as
>>>> >> the
>>>> >> > 'relaxation routines' that I find in other sections of the papers
>>>> as
>>>> >> well?
>>>> >> > If so, what exactly do relaxation routines do? I would like some
>>>> insight
>>>> >> > into it.
>>>> >> >
>>>> >> > Sincerely,
>>>> >> > Abha Vishwakarma
>>>> >> >
>>>> >> > -----------------------------------------------------------------
>>>> >> > *Abha Vishwakarma*
>>>> >> > Final Year Int. M.Sc.
>>>> >> > NISER, India
>>>> >> > Mail: abha.vishwakarma at niser.ac.in
>>>> >> > Web: sites.google.com/view/abhavishwakarma/
>>>> <https://urldefense.com/v3/__http://sites.google.com/view/abhavishwakarma/__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5_ykZhl3w$>
>>>> >> >
>>>> >> >
>>>> >> > _______________________________________________
>>>> >> > 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!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> >> >
>>>> >> >
>>>> >> -------------- next part --------------
>>>> >> An HTML attachment was scrubbed...
>>>> >> URL: <
>>>> >>
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_5-sFv_edg$>
>>>> >> >
>>>> >>
>>>> >> ------------------------------
>>>> >>
>>>> >> Message: 3
>>>> >> Date: Thu, 15 Feb 2024 03:01:20 +0000
>>>> >> From: "Yamaguchi, Natsuko" <nyamaguc at caltech.edu>
>>>> >> To: Evan Bauer via Mesa-users <mesa-users at lists.mesastar.org>
>>>> >> Subject: [Mesa-users] A question about internal energy in MESA
>>>> >> Message-ID:
>>>> >> <
>>>> >>
>>>> BY5PR03MB53169F2A0B85C38C06800E43DF4D2 at BY5PR03MB5316.namprd03.prod.outlook.com
>>>> >> >
>>>> >>
>>>> >> Content-Type: text/plain; charset="iso-8859-1"
>>>> >>
>>>> >> Hello,
>>>> >>
>>>> >> I would like to clarify what exactly MESA's (specific) 'internal
>>>> energy'
>>>> >> is composed of.
>>>> >>
>>>> >> Is it the (1) gas pressure, (2) radiation pressure, and (3)
>>>> ionization
>>>> >> (or recombination) energy combined? And which species is it
>>>> considering for
>>>> >> the ionization energy? It's not obvious from Paxton et al. (2018)
>>>> which
>>>> >> primarily describes it in the context of a simple example. If there
>>>> is a
>>>> >> better source on this, please let me know.
>>>> >>
>>>> >> Thank you,
>>>> >>
>>>> >> Natsuko Yamaguchi
>>>> >>
>>>> >>
>>>> >>
>>>> >>
>>>> >>
>>>> >> -------------- next part --------------
>>>> >> An HTML attachment was scrubbed...
>>>> >> URL: <
>>>> >>
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240215/15fc3794/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240215/15fc3794/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58BjdW9Jw$>
>>>> >> >
>>>> >>
>>>> >> End of Mesa-users Digest, Vol 79, Issue 11
>>>> >> ******************************************
>>>> >>
>>>> > _______________________________________________
>>>> > 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!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>> >
>>>> >
>>>> -------------- next part --------------
>>>> An HTML attachment was scrubbed...
>>>> URL: <
>>>> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240224/38372676/attachment-0001.htm
>>>> <https://urldefense.com/v3/__https://lists.mesastar.org/pipermail/mesa-users/attachments/20240224/38372676/attachment-0001.htm__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_59n56FnvQ$>
>>>> >
>>>>
>>>> End of Mesa-users Digest, Vol 79, Issue 18
>>>> ******************************************
>>>>
>>> _______________________________________________
>>> 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!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$>
>>>
>>> _______________________________________________
>> mesa-users at lists.mesastar.org
>>
>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!IKRxdwAv5BmarQ!ZxNLLzG0LrrQY9bK4tA_TPAqF44zfjNdwa3su3Sc-mMRtXdH3avowhoXx3fa8rmF0FlujnqIHj0fcvxBqgqr_58OMudtyA$
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240302/20784bca/attachment.htm>
More information about the Mesa-users
mailing list