Implementing Eddington limit for BHs based on spinup/ISCO #570
Conversation
…sts the accretion based on the BH ISCO due to spinup since BH formation
Codecov Report
@@ Coverage Diff @@
## develop #570 +/- ##
===========================================
+ Coverage 87.13% 88.12% +0.99%
===========================================
Files 40 40
Lines 25393 25432 +39
===========================================
+ Hits 22124 22410 +286
+ Misses 3269 3022 -247
Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here. |
| @@ -800,6 +800,11 @@ SUBROUTINE hrdiag(mass,aj,mt,tm,tn,tscls,lums,GB,zpars, | |||
| endif | |||
| endif | |||
| mt = mrem | |||
| * Store the initial BH mass for calculating the ISCO later | |||
There was a problem hiding this comment.
What do you think about initializing Mbh_initial to 0, then checking if Mbh_initial.eq.0 instead of BHbirth_ind?
There was a problem hiding this comment.
Yep makes total sense and is definitely cleaner, I’ll push that change to this PR.
|
@katiebreivik I've done what you suggested. Looks like the unit tests never encountered super-Eddington accretion so they still pass. There's still the issue with the deltam1_bcm for the first timestep of the BCM array sometimes still being super-Eddington, though I don't think this is a huge issue. |
|
Hey @michaelzevin sorry for being slow! -- can you open an issue with an example binary that still goes super Eddington on the first time step? I just want to make sure we don't lose it in the mix. Once that issue's open, feel free to merge! |
|
Good call. I've opened an Issue at #571. I'll now go ahead and merge. |
This PR adjusts the Eddington limit for BHs based on the implementation in Marchant et al. 2017, following Equations 3-5. It introduces a new global parameter
Mbh_initialthat remembers the initial mass of the BH at formation, and calculates theetaparameter that alters the mass transfer rate based on where the ISCO is, assuming the BH is spinning up from stable accretion as in Thorne 1970.As of now, I added an additional global parameter
BHbirth_indthat is 0 before the BH forms and 1 after so thatMbh_initialdoesn't get overwritten inhrdiag.f, though there's probably a more glamorous way to ensure this.This actually causes the Eddington limit for BHs to be ~1 order of magnitude higher for low-spinning BHs than the standard BSE prescription!!! It is unclear where Hurley got the expression for the Eddington limit that is implemented in stock BSE (Equation 67 of Hurley et al. 2002). The two expressions are identical for
eta = 0.5, which is slightly above the maximum value ofeta=0.42which is attained for a maximally-spinning BH (one that accreted more than sqrt(6) times its initial mass).Attached are the mass transfer rates for a few binaries, comparing the current BH accretion limit to the one in this PR. There is still an issue of deltam_1 for the first timestep written into the bcm array not obeying the imposed Eddington limit, which can be seen in a few of these systems.
There are probably other changes that need to be made to the unit test files, so let me know what else needs to be done!