[Mesa-users] Mesa-users Digest, Vol 52, Issue 4

Tue Nov 9 13:13:48 UTC 2021

Please keep your replies on mesa-users.

> Thank you very much for answering my question!
> I found there is flag  'dump_missing_metals_into_heaviest' in star_job and
> 'accretion_dump_missing_metals_into_heaviest' in controls.The explanations
> of these two flags is same as follow:
>
> this controls the treatment metals that are not included in the current
> net. if this flag is true, then the mass fractions of missing metals are
> added to the mass fraction of the most massive metal included in the net.
> if this flag is false, then the mass fractions of the metals in the net are
> renormalized to make up for the total mass fraction of missing metals.
>
> What difference between these two flags?
>

One is for when you're dealing with accretion (the composition of the
accreted material) and the other is not.

Aaron

>
>
>
> 在 2021-11-08 00:52:41，"Aaron Dotter" <aaron.dotter at gmail.com> 写道：
>
> Hi Liguo,
>
> what is the conversion formula between  number density and mass fraction?
>>
>
> number fraction of isotope = (mass fraction of isotope) / (mass of isotope)
>
> and you have to (re)normalize so that sum(number fractions) = 1.0 and
> sum(mass fractions) = 1.0.
>
>
>
>
>>
>>
>>
>> At 2021-11-05 22:19:17, "Aaron Dotter" <aaron.dotter at gmail.com> wrote:
>>
>>
>> So I have some questions as follow:
>>> 1.In different metallicity Z,is the ratio among all isotopes fixed?what
>>> is the ratio? Which paper do this ratio refer?
>>>
>>
>> Yes, it is the scaled the same for all Z.  The ratio depends on your
>> choice of the heavy element distribution, which can be set in the star_job
>> section of the inlist.
>>
>> https://docs.mesastar.org/en/r15140/reference/star_job.html#initial-zfracs
>>
>>
>>
>>> 2.Why didn't the flag 'dump_missing_metals_into_heaviest = .true.' work?
>>>
>>
>> As I noted in my previous email, the default value for this flag is true,
>> so by setting it to true in your inlist you have not changed anything --
>> and therefore should not expect to see any difference.
>>
>>
>> Aaron
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>>
>>>
>>>
>>>
>>> At 2021-11-04 20:00:03, 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: [Fe/H] Calculation with MESA Outputs (Aaron Dotter)
>>> >   2. Re: c13_pocket stop running during TACHeB (Rob Farmer)
>>> >   3. Re: c13_pocket stop running during TACHeB (Liguo)
>>> >   4. Fw: University of Sydney Postdoctoral Research Associate in
>>> >      Stellar Astrophysics (women only) (Dennis Stello)
>>> >
>>> >
>>> >----------------------------------------------------------------------
>>> >
>>> >Message: 1
>>> >Date: Wed, 3 Nov 2021 08:11:14 -0400
>>> >From: Aaron Dotter <aaron.dotter at gmail.com>
>>> >To: Rob Farmer <robert.j.farmer37 at gmail.com>
>>> >Cc: Mitchell Dennis <mtde226 at hawaii.edu>, mesa-users
>>> >	<mesa-users at lists.mesastar.org>
>>> >Subject: Re: [Mesa-users] [Fe/H] Calculation with MESA Outputs
>>> >Message-ID:
>>> >	<CAN5--jNBGjtMpbLr1RpLnZRJpOA7TOGZbZFC5vnXfc+pqoJruQ at mail.gmail.com>
>>> >Content-Type: text/plain; charset="utf-8"
>>> >
>>> >Hi Mitchell,
>>> >
>>> >A couple of additions to Rob's good suggestion.
>>> >
>>> >- You can add isotopes to a network without adding the nuclear reactions.
>>> >For example, if you're interested in tracking the surface evolution of Fe56
>>> >in solar type stars, which won't create or destroy Fe56, then you can still
>>> >use a simple net like basic for energy generation but also track the mixing
>>> >history of other isotopes.
>>> >
>>> >- Whether you do the above or not, beware of a little known feature in
>>> >MESA: In the star_job inlist there is a flag called
>>> >
>>> >      ! dump_missing_metals_into_heaviest
>>> >      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>>> >
>>> >      ! this controls the treatment metals that are not included in the
>>> >current net.
>>> >      ! if this flag is true, then the mass fractions of missing metals are
>>> >added
>>> >      ! to the mass fraction of the most massive metal included in the net.
>>> >      ! if this flag is false, then the mass fractions of the metals in the
>>> >net
>>> >      ! are renormalized to make up for the total mass fraction of missing
>>> >metals.
>>> >
>>> >      ! ::
>>> >
>>> >    dump_missing_metals_into_heaviest = .true.
>>> >
>>> >What this means is that any elements in the solar mixture that are not
>>> >explicitly included in the network are going to be added to the heaviest
>>> >element in the network.  Depending on the network you choose, this may be a
>>> >small fraction.  If you use an approx21 net, or follow my suggestion above,
>>> >make sure that the element you want to track is not the heaviest in the net
>>> >(or turn the above flag to false but then you have a different problem . .
>>> >.).
>>> >
>>> >See attached example.  I use Zn64 as a placeholder for all the extras so
>>> >that Fe56 will just be Fe56.
>>> >
>>> >Aaron
>>> >-------------- next part --------------
>>> >An HTML attachment was scrubbed...
>>> >URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211103/3906833e/attachment-0001.htm>
>>> >-------------- next part --------------
>>> >A non-text attachment was scrubbed...
>>> >Name: basic_extra.net
>>> >Type: application/octet-stream
>>> >Size: 78 bytes
>>> >Desc: not available
>>> >URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211103/3906833e/attachment-0001.obj>
>>> >
>>> >------------------------------
>>> >
>>> >Message: 2
>>> >Date: Wed, 3 Nov 2021 14:03:13 +0100
>>> >From: Rob Farmer <robert.j.farmer37 at gmail.com>
>>> >To: Liguo <liguo070923 at 163.com>
>>> >Cc: mesa-users <mesa-users at lists.mesastar.org>
>>> >Subject: Re: [Mesa-users] c13_pocket stop running during TACHeB
>>> >Message-ID:
>>> >	<CAPvTSya47h2RX054jJ11oh_=09+Py3gxXvfaimpsZMGs1dQE1w at mail.gmail.com>
>>> >Content-Type: text/plain; charset="utf-8"
>>> >
>>> >Hi,
>>> >That's not enough information to diagnose a problem.
>>> >
>>> >You need to show your inlists (and your nuclear network).
>>> >
>>> >What is your star doing when it stops?
>>> >Have you made any plots to try to diagnose the problem yourself?
>>> >Are there regions of the star where it seems like mesa is having
>>> >difficulties?
>>> >What have you tried to do to fix the problem yourself?
>>> >
>>> >Rob
>>> >
>>> >On Tue, 2 Nov 2021 at 09:47, Liguo via Mesa-users <
>>> >mesa-users at lists.mesastar.org> wrote:
>>> >
>>> >> Dear MESA Mailing List,
>>> >>   I can run c13_pocket successfully with it's default setting.But when I
>>> >> changed it's net to bigger net than agb.net,It will stop running in
>>> >> TACHeB phase with error as follow:
>>> >>   stopping because of problems dt < min_timestep_limit
>>> >>   My net have about 250 isotopes.
>>> >>   Any suggestions or help would be very appreciated!
>>> >>
>>> >>   sincerely!
>>> >>   Liguo
>>> >>
>>> >>
>>> >>
>>> >> _______________________________________________
>>> >> 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/20211103/68b4531b/attachment-0001.htm>
>>> >
>>> >------------------------------
>>> >
>>> >Message: 3
>>> >Date: Thu, 4 Nov 2021 10:49:58 +0800 (CST)
>>> >From: Liguo  <liguo070923 at 163.com>
>>> >To: "Rob Farmer" <robert.j.farmer37 at gmail.com>,  mesa-users
>>> >	<mesa-users at lists.mesastar.org>
>>> >Subject: Re: [Mesa-users] c13_pocket stop running during TACHeB
>>> >Message-ID: <668e3033.2111.17ce8da62d2.Coremail.liguo070923 at 163.com>
>>> >Content-Type: text/plain; charset="gbk"
>>> >
>>> >Dear Rob
>>> >Thank you for your reply!
>>> >My star stopped running after h burning and before he burning.
>>> >I have tried to set min_timestep_limit = 1d-20.Star then stopped running with error segmentation fault-invalidmemory reference.
>>> >I attach my inlist and net file in the end.
>>> >I only change net to my net in inlist_common of c13_pocket.
>>> >
>>> >
>>> >sincerely!
>>> >Liguo
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >? 2021-11-03 21:03:13?"Rob Farmer" <robert.j.farmer37 at gmail.com> ???
>>> >
>>> >Hi,
>>> >
>>> >That's not enough information to diagnose a problem.
>>> >
>>> >
>>> >You need to show your inlists (and your nuclear network).
>>> >
>>> >
>>> >
>>> >What is your star doing when it stops?
>>> >
>>> >Have you made any plots to try to diagnose the problem yourself?
>>> >
>>> >Are there regions of the star where it seems like mesa is having difficulties?
>>> >
>>> >What have you tried to do to fix the problem yourself?
>>> >
>>> >
>>> >Rob
>>> >
>>> >
>>> >
>>> >On Tue, 2 Nov 2021 at 09:47, Liguo via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>>> >
>>> >Dear MESA Mailing List,
>>> >  I can run c13_pocket successfully with it's default setting.But when I changed it's net to bigger net than agb.net,It will stop running in TACHeB phase with error as follow:
>>> >  stopping because of problems dt < min_timestep_limit
>>> >  My net have about 250 isotopes.
>>> >  Any suggestions or help would be very appreciated!
>>> >
>>> >
>>> >  sincerely!
>>> >  Liguo
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >
>>> >_______________________________________________
>>> >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/20211104/df109bb2/attachment-0001.htm>
>>> >-------------- next part --------------
>>> >A non-text attachment was scrubbed...
>>> >Name: inlist_common
>>> >Type: application/octet-stream
>>> >Size: 2038 bytes
>>> >Desc: not available
>>> >URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211104/df109bb2/attachment-0002.obj>
>>> >-------------- next part --------------
>>> >A non-text attachment was scrubbed...
>>> >Name: SprocessWang_agb.net
>>> >Type: application/octet-stream
>>> >Size: 1677 bytes
>>> >Desc: not available
>>> >URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211104/df109bb2/attachment-0003.obj>
>>> >
>>> >------------------------------
>>> >
>>> >Message: 4
>>> >Date: Thu, 4 Nov 2021 15:24:55 +1100
>>> >From: Dennis Stello <dennis.stello at sydney.edu.au>
>>> >To: Mesa-users <mesa-users at lists.mesastar.org>
>>> >Subject: [Mesa-users] Fw: University of Sydney Postdoctoral Research
>>> >	Associate in Stellar Astrophysics (women only)
>>> >Message-ID: <88bd20ef-ec99-2671-51a4-b0b8dd9e134e at sydney.edu.au>
>>> >Content-Type: text/plain; charset="utf-8"; Format="flowed"
>>> >
>>> >Dear mesa users,
>>> >
>>> >A potential opportunity to work on stars in Sydney.
>>> >
>>> >Best
>>> >Dennis/
>>> >
>>> >
>>> >-------- Forwarded Message --------
>>> >Subject: 	Fw: University of Sydney Postdoctoral Research Associate in
>>> >Stellar Astrophysics (women only)
>>> >Date: 	Thu, 4 Nov 2021 04:20:24 +0000
>>> >From: 	Tim Bedding <tim.bedding at sydney.edu.au>
>>> >To: 	Dennis Stello <d.stello at unsw.edu.au>
>>> >
>>> >
>>> >
>>> >Hi Dennis,
>>> >
>>> >Can you please forward to the MESA list?
>>> >
>>> >thanks
>>> >
>>> > ? ?Tim
>>> >
>>> >------------------------------------------------------------------------
>>> >
>>> >The University of Sydney is welcoming applications for a three-year
>>> >*Postdoctoral Research Associate in Stellar Astrophysics (women only)*.
>>> > ? The successful applicant will be working with Tim Bedding
>>> ><mailto:tim.bedding at sydney.edu.au>?and Daniel Huber?in
>>> >*asteroseismology, exoplanets or other related areas of stellar
>>> >astrophysics*.??More details can be found here
>>> ><https://protect-au.mimecast.com/s/dqdXCjZ1N7in8k2JXfW0iOZ?domain=usyd.wd3.myworkdayjobs.com>?or
>>> >on the AAS jobs register. Applications close 28-Nov-2021.
>>> >-------------- next part --------------
>>> >An HTML attachment was scrubbed...
>>> >URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211104/0f01a90f/attachment-0001.htm>
>>> >-------------- next part --------------
>>> >_______________________________________________
>>> >ASA mailing list
>>> >ASA at mailman.sydney.edu.au
>>> >Change membership status or contact information via ASA's Edit Membership page
>>> >https://protect-au.mimecast.com/s/j2HgClx1Nji2xDyK8uyhTIH?domain=asa.wildapricot.org
>>> >
>>> >End of Mesa-users Digest, Vol 52, Issue 4
>>> >*****************************************
>>>
>>>
>>>
>>>
>>>
>>
>>
>>
>>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20211109/b7dbcb17/attachment.htm>