[mesa-users] trying to run RLO calculations for hot Jupiter systems

[Brian Jackson]

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