[Mesa-users] Query about relax_core in MESA
Abha Vishwakarma
abha.vishwakarma at niser.ac.in
Sat Mar 2 07:50:20 UTC 2024
Dear Dr. Earl,
Thanks a lot for your reply. I didn't realise gyre_in_mesa_envelope is
another test suite. That definitely helped clear a lot of confusion. Thanks
again!
Best,
Abha
On Tue, Feb 27, 2024 at 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
> 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. changing the Roche-lobe size from an inlist
> (Temmink, K.D. (Karel))
> 2. Re: Query about relax_core in MESA (Abha Vishwakarma)
> 3. Re: changing the Roche-lobe size from an inlist (Pablo Marchant)
> 4. Visiting Graduate Student Fellowships at the Flatiron
> Institute (Jared Goldberg)
> 5. Re: Query about relax_core in MESA (Earl Bellinger)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 26 Feb 2024 12:37:47 +0000
> From: "Temmink, K.D. (Karel)" <karel.temmink at ru.nl>
> To: "mesa-users at lists.mesastar.org" <mesa-users at lists.mesastar.org>
> Subject: [Mesa-users] changing the Roche-lobe size from an inlist
> Message-ID:
> <
> DB4PR10MB7422DF62C08417CDDB0BFB23F35A2 at DB4PR10MB7422.EURPRD10.PROD.OUTLOOK.COM
> >
>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear MESA users,
>
> I'm struggling a bit with making basic changes to the MESA code (version
> r12115).
> I want to multiply the Roche lobe radius by a number that could be
> provided in the inlists (via a control option such as: x_ctrl(1) = 1.1).
>
> To me, it seemed that the best way would be to make changes to the
> eval_rlobe function in $MESA_DIR/binary/private/binary_utils.f90 .
> However, this function does not have access to the necessary pointers (b%)
> to make the above work.
> I saw that this function is called in a number of places across MESA (such
> as binary_utils.f90, binary_ce.f90, binary_evolve.f90 and binary_lib.f90).
> That would mean that making such a change would also require making
> changes there, and possibly in more places too.
>
> I was wondering if someone knows of a better or more efficient way of
> accomplishing what I want?
> Any help, tips or advice is greatly appreciated!
>
> Kind regards,
>
> Karel Temmink
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/81e4af73/attachment-0001.htm
> >
>
> ------------------------------
>
> Message: 2
> Date: Mon, 26 Feb 2024 18:31:47 +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_p51XYX8ry+FmonWn8QvDUmK8P1VT2pzqgprMu6G8bzu9-A at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> 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
> > 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
> > > 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/
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL: <
> > >
> >
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
> > > >
> > >
> > > ------------------------------
> > >
> > > 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
> > > modelling red supergiants with neutron star cores (Thorne??ytkow
> > objects):
> > > https://arxiv.org/abs/2305.07337
> > > and modelling stars with primordial black holes at their center
> (Hawking
> > > stars):
> > > https://github.com/earlbellinger/black-hole-sun
> > > 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>
> > > 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
> > >
> > > 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
> > >
> > >
> > > 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/
> > > >
> > > >
> > > > _______________________________________________
> > > > mesa-users at lists.mesastar.org
> > > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > > >
> > > >
> > > -------------- next part --------------
> > > An HTML attachment was scrubbed...
> > > URL: <
> > >
> >
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
> > > >
> > >
> > > ------------------------------
> > >
> > > 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
> > > >
> > >
> > > 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
> > >
> > -------------- next part --------------
> > finished doing relax core
> >
> > net name basic.net
> > 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
> > >
> >
> > ------------------------------
> >
> > 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
> > >> 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/
> > >> -------------- next part --------------
> > >> An HTML attachment was scrubbed...
> > >> URL: <
> > >>
> >
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
> > >> >
> > >>
> > >> ------------------------------
> > >>
> > >> 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
> > >> modelling red supergiants with neutron star cores (Thorne??ytkow
> > objects):
> > >> https://arxiv.org/abs/2305.07337
> > >> and modelling stars with primordial black holes at their center
> (Hawking
> > >> stars):
> > >> https://github.com/earlbellinger/black-hole-sun
> > >> 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>
> > >> 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
> > >>
> > >> 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
> > >>
> > >>
> > >> 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/
> > >> >
> > >> >
> > >> > _______________________________________________
> > >> > mesa-users at lists.mesastar.org
> > >> > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > >> >
> > >> >
> > >> -------------- next part --------------
> > >> An HTML attachment was scrubbed...
> > >> URL: <
> > >>
> >
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
> > >> >
> > >>
> > >> ------------------------------
> > >>
> > >> 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
> > >> >
> > >>
> > >> End of Mesa-users Digest, Vol 79, Issue 11
> > >> ******************************************
> > >>
> > > _______________________________________________
> > > mesa-users at lists.mesastar.org
> > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> > >
> > >
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: <
> >
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240224/38372676/attachment-0001.htm
> > >
> >
> > End of Mesa-users Digest, Vol 79, Issue 18
> > ******************************************
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/9e9c2073/attachment-0001.htm
> >
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: removed.mod
> Type: audio/x-mod
> Size: 163126 bytes
> Desc: not available
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/9e9c2073/attachment-0001.bin
> >
>
> ------------------------------
>
> Message: 3
> Date: Mon, 26 Feb 2024 14:27:53 +0100
> From: Pablo Marchant <pamarca at gmail.com>
> To: "Temmink, K.D. (Karel)" <karel.temmink at ru.nl>
> Cc: "mesa-users at lists.mesastar.org" <mesa-users at lists.mesastar.org>
> Subject: Re: [Mesa-users] changing the Roche-lobe size from an inlist
> Message-ID:
> <
> CAF_dd5O34eU07P6Ag4TgZHBU9ybsk+sZaHD0ZpDeQwtT71fHZQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Karel,
>
> Don't think there is an efficient way to achieve that right now. You will
> need to modify the source code to achieve this. I would suggest changing
> the eval_rlobe function to receive the binary_info too:
>
> real(dp) function eval_rlobe(binary_id, m1, m2, a) result(rlobe)
>
> and then make fixes till the compiler is happy. Modifying things at this
> lower level will ensure you catch all places where there could be an issue.
> I guess the ideal approach would be to provide a hook for this, if you feel
> comfortable enough with it you can try to give that a go and submit a pull
> request to add that functionality.
>
> Cheers
>
> On Mon, Feb 26, 2024 at 1:38?PM Temmink, K.D. (Karel) via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
> > Dear MESA users,
> >
> > I'm struggling a bit with making basic changes to the MESA code (version
> > r12115).
> > I want to multiply the Roche lobe radius by a number that could be
> > provided in the inlists (via a *control* option such as: x_ctrl(1) =
> 1.1).
> >
> > To me, it seemed that the best way would be to make changes to the
> > eval_rlobe function in $MESA_DIR/binary/private/binary_utils.f90 .
> > However, this function does not have access to the necessary pointers
> (b%)
> > to make the above work.
> > I saw that this function is called in a number of places across MESA
> (such
> > as binary_utils.f90, binary_ce.f90, binary_evolve.f90 and
> binary_lib.f90).
> > That would mean that making such a change would also require making
> > changes there, and possibly in more places too.
> >
> > I was wondering if someone knows of a better or more efficient way of
> > accomplishing what I want?
> > Any help, tips or advice is greatly appreciated!
> >
> > Kind regards,
> >
> > Karel Temmink
> > _______________________________________________
> > mesa-users at lists.mesastar.org
> > https://lists.mesastar.org/mailman/listinfo/mesa-users
> >
> >
>
> --
> Pablo Marchant Campos
> M.Sc on Astrophysics, Universidad Cat?lica de Chile
> PhD on Astrophysics, Argelander-Institut f?r Astronomie, Universit?t Bonn
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/5a4a66de/attachment-0001.htm
> >
>
> ------------------------------
>
> Message: 4
> Date: Mon, 26 Feb 2024 16:27:59 -0500
> From: Jared Goldberg <jgoldberg at flatironinstitute.org>
> To: mesa-users <mesa-users at lists.mesastar.org>
> Subject: [Mesa-users] Visiting Graduate Student Fellowships at the
> Flatiron Institute
> Message-ID:
> <CADjORnANSUfKFojRdMnkm0Q3M9dp_ramTME6gsaQkg28Og=
> kaQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi mesa-users!
>
> For the students currently enrolled in PhD programs who are looking to
> broaden your horizons and work in NYC for half a year, we're happy to
> announce that applications for the 2024-2025 Flatiron Institute Center for
> Computational Astrophysics Pre-Doctoral Program are now open! I imagine
> there are a number of projects with multiple advisors which might be of
> interest to some of you. *The application process is slightly different
> than in years past; applicants must fill out the Form
> <
> https://docs.google.com/forms/d/e/1FAIpQLSfXpC63LtbU0i9_A6w7mszTs9iDDzr7QF4aYIo7CzLUAhxdOg/viewform
> >
> for the stage 1 application by 11:59 pm EST on March 15, 2024.*
>
> A bit more info:
>
> We are pleased to announce that applications for the 2024-2025 Flatiron
> Institute CCA Pre-Doctoral Program are now open. This program enables
> students currently enrolled in a PhD program to spend 4.5 months employed
> at the CCA in Manhattan collaborating with 1 or more CCA mentors on a
> project complementary to their thesis research. Pre-docs receive a stipend
> and other benefits. For more details and to apply, please visit
>
> https://www.simonsfoundation.org/flatiron-institute-center-for-computational-astrophysics-pre-doctoral-program
> .
> Applicants must fill out the Google Form
> <
> https://docs.google.com/forms/d/e/1FAIpQLSfXpC63LtbU0i9_A6w7mszTs9iDDzr7QF4aYIo7CzLUAhxdOg/viewform
> >
> for the stage 1 application by 11:59 pm EST on March 15, 2024.
>
> Students who currently have limited access to computational astrophysics
> resources and mentorship are especially encouraged to apply. The Simons
> Foundation actively seeks a diverse applicant pool and encourages
> candidates of all backgrounds to apply.
>
> Cheers,
> ~Jared
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/e62ced82/attachment-0001.htm
> >
>
> ------------------------------
>
> Message: 5
> Date: Mon, 26 Feb 2024 16:43:51 -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:
> <CAD7H7L3MgNbMARi7BnjpF_RiPrdzWjYK6zaKSXn=oWh_r0yH=
> A at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> 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
> >> 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
> >> > 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/
> >> > -------------- next part --------------
> >> > An HTML attachment was scrubbed...
> >> > URL: <
> >> >
> >>
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
> >> > >
> >> >
> >> > ------------------------------
> >> >
> >> > 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
> >> > modelling red supergiants with neutron star cores (Thorne??ytkow
> >> objects):
> >> > https://arxiv.org/abs/2305.07337
> >> > and modelling stars with primordial black holes at their center
> (Hawking
> >> > stars):
> >> > https://github.com/earlbellinger/black-hole-sun
> >> > 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>
> >> > 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
> >> >
> >> > 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
> >> >
> >> >
> >> > 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/
> >> > >
> >> > >
> >> > > _______________________________________________
> >> > > mesa-users at lists.mesastar.org
> >> > > https://lists.mesastar.org/mailman/listinfo/mesa-users
> >> > >
> >> > >
> >> > -------------- next part --------------
> >> > An HTML attachment was scrubbed...
> >> > URL: <
> >> >
> >>
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
> >> > >
> >> >
> >> > ------------------------------
> >> >
> >> > 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
> >> > >
> >> >
> >> > 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
> >> >
> >> -------------- next part --------------
> >> finished doing relax core
> >>
> >> net name basic.net
> >> 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
> >> >
> >>
> >> ------------------------------
> >>
> >> 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
> >> >> 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/
> >> >> -------------- next part --------------
> >> >> An HTML attachment was scrubbed...
> >> >> URL: <
> >> >>
> >>
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/84864efb/attachment-0001.htm
> >> >> >
> >> >>
> >> >> ------------------------------
> >> >>
> >> >> 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
> >> >> modelling red supergiants with neutron star cores (Thorne??ytkow
> >> objects):
> >> >> https://arxiv.org/abs/2305.07337
> >> >> and modelling stars with primordial black holes at their center
> >> (Hawking
> >> >> stars):
> >> >> https://github.com/earlbellinger/black-hole-sun
> >> >> 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
> >
> >> >> 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
> >> >>
> >> >> 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
> >> >>
> >> >>
> >> >> 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/
> >> >> >
> >> >> >
> >> >> > _______________________________________________
> >> >> > mesa-users at lists.mesastar.org
> >> >> > https://lists.mesastar.org/mailman/listinfo/mesa-users
> >> >> >
> >> >> >
> >> >> -------------- next part --------------
> >> >> An HTML attachment was scrubbed...
> >> >> URL: <
> >> >>
> >>
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240214/2278ce1c/attachment-0001.htm
> >> >> >
> >> >>
> >> >> ------------------------------
> >> >>
> >> >> 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
> >> >> >
> >> >>
> >> >> End of Mesa-users Digest, Vol 79, Issue 11
> >> >> ******************************************
> >> >>
> >> > _______________________________________________
> >> > mesa-users at lists.mesastar.org
> >> > https://lists.mesastar.org/mailman/listinfo/mesa-users
> >> >
> >> >
> >> -------------- next part --------------
> >> An HTML attachment was scrubbed...
> >> URL: <
> >>
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240224/38372676/attachment-0001.htm
> >> >
> >>
> >> End of Mesa-users Digest, Vol 79, Issue 18
> >> ******************************************
> >>
> > _______________________________________________
> > mesa-users at lists.mesastar.org
> > https://lists.mesastar.org/mailman/listinfo/mesa-users
> >
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://lists.mesastar.org/pipermail/mesa-users/attachments/20240226/f33ad941/attachment-0001.htm
> >
>
> End of Mesa-users Digest, Vol 79, Issue 20
> ******************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20240302/faa22cf2/attachment.htm>
More information about the Mesa-users
mailing list