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

Sun Nov 7 16:52:41 UTC 2021

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/20211107/811cb49a/attachment.htm>