[mesa-users] trying to run RLO calculations for hot Jupiter systems
Evan Bauer
ebauer at physics.ucsb.edu
Tue Feb 28 20:13:52 EST 2017
Hi Brian,
Sorry I’m a little late jumping in on this one, but this issue is already partially fixed in the latest public MESA release (9575). It will now throw out a warning about mdot_edd at you and prevent the code from calculating NaN’s for mdot_edd_eta. Perhaps just a little motivation to upgrade to the latest and greatest version ;)
Cheers,
Evan
> On Feb 28, 2017, at 4:52 PM, Brian Jackson <bjackson at boisestate.edu> wrote:
>
> That got it. Thanks!
>
> Brian
>
> On Tue, Feb 28, 2017 at 5:38 PM, Pablo Marchant <pamarca at gmail.com> wrote:
>> Besides from the nuclear network discrepancy, here is an issue that Evan ran
>> into and reported to me a couple of weeks ago. The default implementation
>> for the Eddington mass accretion limit applies to BH accretors, and it
>> breaks down and produces NaNs when evolving both stars simultaneously with
>> limit_retention_by_mdot_edd = .true.. This translates into the obscure error
>> messages you are getting Brian.
>>
>> Just commited a fix that will be up in the next MESA version. For the
>> moment, just set limit_retention_by_mdot_edd = .false. in your inlist and
>> things should work.
>>
>> Cheers!
>>
>> On Tue, Feb 28, 2017 at 5:47 PM, Robert Farmer <rjfarmer at asu.edu> wrote:
>>>
>>> Hi
>>> The nuclear networks (and hence the isotopes each model has defined) in
>>> the two model files are different, can you try in the planet inlist
>>>
>>> new_net_name = 'o18_and_ne22.net'
>>> change_initial_net = .true.
>>>
>>> so the planet model will have the same network as the star
>>>
>>> Rob
>>>
>>>
>>> On Tue, Feb 28, 2017 at 4:30 PM, Brian Jackson <bjackson at boisestate.edu>
>>> wrote:
>>>>
>>>> Dear MESA Users,
>>>> In my group, we're trying to continue our analysis of planetary
>>>> Roche-lobe overflow (RLO).
>>>>
>>>> For that, we generate starting planetary and stellar models
>>>> (attached). The star model is generated using the 1M_pre_ms_to_WD test
>>>> suite with only "max_age = 20e6" added to inlist_1.0 file. The
>>>> planetary model is generated using the make_planets test suite with
>>>> only "max_age = 20d6" changed in inlist_evolve_1.0_MJ_10.0_ME_2.0_RJ
>>>> file.
>>>>
>>>> We make some modifications to inlists for the binary_both_stars test
>>>> suite to load in the saved model files but otherwise use the standard
>>>> test suite.
>>>>
>>>> The following shows the last few lines from MESA before it stops
>>>> running our calculation:
>>>>
>>>> get_xa_for_accretion: accretion species mass fractions do not add to 1.0
>>>> sum(xa(1:species))
>>>> 0.0000000000000000D+00
>>>> h1 1
>>>> 0.0000000000000000D+00
>>>> he3 2
>>>> 0.0000000000000000D+00
>>>> he4 3
>>>> 0.0000000000000000D+00
>>>> c12 4
>>>> 0.0000000000000000D+00
>>>> n14 5
>>>> 0.0000000000000000D+00
>>>> o16 6
>>>> 0.0000000000000000D+00
>>>> ne20 7
>>>> 0.0000000000000000D+00
>>>> mg24 8
>>>> 0.0000000000000000D+00
>>>> 1 1 2.000000000E+07
>>>> 9.244518984E-04 -9.900000000E+01 -3.245769587E-01
>>>> 2 1 2.000000000E+07
>>>> 9.999853565E-01 -9.900000000E+01 -6.799606683E-01
>>>> have reached retry limit so now backup
>>>> have reached retry limit so now backup
>>>>
>>>> dt
>>>> 6.3566365930376634D-07
>>>> min_timestep_limit
>>>> 9.9999999999999995D-07
>>>>
>>>> stopping because of convergence problems dt < min_timestep_limit
>>>>
>>>> terminated evolution: convergence problems
>>>>
>>>> termination code: min_timestep_limit
>>>>
>>>> ---
>>>>
>>>> The species mass fractions look like nonsense, even though they look
>>>> fine in the initial model files. And strangely, they change depending
>>>> on the initial orbital period we choose.
>>>>
>>>> I am running Mac OS X 10.9.5, MESA version 8845, and the MESA SDK
>>>> version from April 4 2016. I have attached the binary inlists and
>>>> initial models we use to run the calculation.
>>>>
>>>> Please let me know if you need additional information and if you have
>>>> any advice.
>>>>
>>>> Many thanks,
>>>> Brian Jackson
>>>>
>>>>
>>>> ------------------------------------------------------------------------------
>>>> Check out the vibrant tech community on one of the world's most
>>>> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
>>>> _______________________________________________
>>>> mesa-users mailing list
>>>> mesa-users at lists.sourceforge.net
>>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>>
>>>
>>>
>>>
>>> ------------------------------------------------------------------------------
>>> Check out the vibrant tech community on one of the world's most
>>> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>
>>
>>
>>
>> --
>> Pablo Marchant Campos
>> M.Sc on Astrophysics, Universidad Católica de Chile
>> PhD student, Argelander-Institut für Astronomie
>
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, SlashDot.org! http://sdm.link/slashdot
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
More information about the Mesa-users
mailing list