[Mesa-users] [EXT] Fwd: Binary evolution residual_norm warning for high masses

[Gökhan Yücel]

Hey Mathieu,

Thanks for your reply. I apologize for not being clear. The evolution stops
at the start, right after stars enter the zams. I edited the image for
better understanding.

The red eclipse shows the timestep, which is extremely small as you can see
and the yellow eclipse the age of the system, which is also very small.

[image: Screenshot from 2025-01-30 12-48-30.png]


I have also tried to do without the "zams patch" but the result was the
same.

It seems the problem with rotation. When I turned it off, the evolution
started as it should be. But that's not the realistic case. So, I added,

relax_surface_rotation_v = .true.
new_surface_rotation_v = 30 !km/s

these commands and it seems the evolution proceeds as intended so far.
Obviously I have not done the testing extensively, but it's a good starting
point.

Much appreciated,
Gökhan





On Thu, Feb 6, 2025 at 8:40 PM Mathieu Renzo via Mesa-users <
mesa-users at lists.mesastar.org> wrote:

> Hi Gökhan,
>
> I am unable to run your inlists because they appear to require "patched"
> ZAMS models in
> `$MESA_DIR/data/star_data/zams_models/zams_z2m2_y28_patched.data`.
>
> Is your binary going through case A mass transfer? That would explain why
> it may crash before TAMS if you're running into issues during mass transfer
> during the main sequence of the donor (and consequently certainly also the
> lower mass accretor).
>
> You may want to check which star limits your timesteps (donor or accretor,
> since you run both), of if it is something related to the orbit. To do this
> check the `dt_limit` in the terminal output (refer to
> $MESA_DIR/binary/defaults/control.defaults and
> $MESA_DIR/star/defaults/controls.defaults for what those strings mean).
>
> You also seem to have rotation on for both stars (great!), in my
> experience running binaries with accretors below ~15Msun, the rotationally
> enhanced wind is too weak to get rid of the accreted angular momentum and
> MESA struggles (sometimes producing omega>omega_crit, which is physically
> not allowed). This could be the cause of your residuals. As an experiment
> to test this you could try cranking up artificially the accretor's mass
> loss rate from inlist2 (maybe only during RLOF), and see if this helps.
> Note that this second problem occurs regardless of whether you have case A
> or case B RLOF (actually, sufficiently strong tides in short periods for
> case A may help preventing this). You can try to play with the "implicit
> wind enhancement" to see if you manage to get through.
>
> Another simple experiment you can try is turn off rotation and tides: do
> your binaries work then?
>
> Hope this helps you a bit tracing the root of the issue...
> Cheers,
> Mathieu
>
>
>
>
>
> On 2/6/25 00:40, Gökhan Yücel via Mesa-users wrote:
>
>
> * External Email *
> ------------------------------
> Hey everyone,
>
> Has anyone any idea about what might cause this?
>
> Gökhan.
>
> ---------- Forwarded message ---------
> From: Gökhan Yücel <gokhannyucel at gmail.com>
> Date: Thu, Jan 30, 2025 at 12:52 PM
> Subject: Binary evolution residual_norm warning for high masses
> To: mesa-users <mesa-users at lists.mesastar.org>
>
>
> Hey everyone,
>
> I have been trying to a basic binary evolution with both stars evolve
> simultaneously. But I get warning and evolution slows to an extreme degree,
> even it stops at the start of the evolution if I try to increase the mass
> of any component more than 6 Msun.
>
> [image: Screenshot from 2025-01-30 12-48-30.png]
>
> I tried to increase tol_resuidual_norm and tol_max_residual parameters,
> but it didn't work.
> I know my inlist is basic for high mass evolutions, but I believe it
> shouldn't give me this warning, considering my calculations don't even
> reach to TAMS. I have attached the necessary files to re-produce the
> warning.
>
> Any help would be highly appreciated.
>
> I'm using MESA version r23.05.1
>
> Cheers,
> Gökhan
>
>
> _______________________________________________mesa-users at lists.mesastar.orghttps://lists.mesastar.org/mailman/listinfo/mesa-users
>
>
> --
> Mathieu Renzo <https://users.flatironinstitute.org/~mrenzo/>
> Assistant Professor
> University of Arizona & Steward Observatory
> _______________________________________________
> 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/20250207/9453c722/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2025-01-30 12-48-30.png
Type: image/png
Size: 55782 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20250207/9453c722/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2025-01-30 12-48-30.png
Type: image/png
Size: 56472 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20250207/9453c722/attachment-0001.png>