MESA inlists and extra evolution code for "Can the Anomalous Magnetic Braking of Ap/Bp Stars Explain the Orbital Decay of Algol-type Binaries?"

This repository contains the MESA inlists and extra evolution code that are used in the article "Can the Anomalous Magnetic Braking of Ap/Bp Stars Explain the Orbital Decay of Algol-type Binaries?", which has been accepted by ApJ. Algol-type binaries are semi-detached binary systems, where a less massiver Roche lobe-filling K/F-type giant/subgiant transfers material to a more massive hotter A/B-type main-sequence companion through Roche lobe overflow. According to conservative mass transfer model, the orbits of Algol binaries are expected to continuously expand. However, observations have revealed that several Algol binaries undergo rapid orbital decay. In this study, we explore whether anomalous magnetic braking (MB) can account for the observed orbital decay. The uploaded zip file AMB contains the simulation codes and settings related to anomalous MB. The main input files are: inlist_project, which specifies the binary parameters and evolution settings, such as the masses of the donor and accretor stars, the orbital period, and the mass-transfer scheme; inlist1, which specifies the metallicity and stellar-wind scheme of the donor star; inlist2, which specifies the metallicity of the accretor star. The source code modification is located in the src folder as run_binary_extra.f, where we added a subroutine an_jdot_mb to implement the angular momentum loss driven by anomalous MB. All simulations were carried out with MESA version r24.08.1. To reproduce the results presented in the paper, replace the corresponding files in mesa-24.08.1/binary/work with those provided here, adjust the initial conditions to match those described in the article, and then recompile and run MESA.