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

[Brian Jackson]

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
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist
Type: application/octet-stream
Size: 458 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_project
Type: application/octet-stream
Size: 888 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist1
Type: application/octet-stream
Size: 939 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist2
Type: application/octet-stream
Size: 914 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: planet_init.mod
Type: application/octet-stream
Size: 189621 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: star_init_20Myr.mod
Type: application/octet-stream
Size: 228381 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170228/fa907156/attachment-0005.obj>