[mesa-users] accretion_powered_irradiation

[Chen WC]

Dear All,
I use the “binary_donor_only_implicit_mdot” to calculate the evolution of an intermediate-mass X-ray binary including a 1.4 Msun NS and a 3.0 Msun donor star (with a solar composition),
and the initial orbital period is 1.0 day.

In calculation, I consider the irradiation process by X-ray luminosity of the NS. In the list_project:
****
accretion_powered_irradiation = .true.
col_depth_for_eps_extra = 100.0 ! g/cm^2
****
Other input parameters such as accretor_radius_for_irrd, max_F_irr were adopted to be default.


However, after the mass of the donor star decrease to be ~ 2.3-2.4 Msun, the calculation is very difficult to continue.
The timestep decrease to be <0.1 yr, and the terminal output  for timestep control displays “max increase”.

I would appreciate if someone can help me to solve this problem.


Cheers,


Chen Wen-Cong


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160610/c23bc9b1/attachment.html>