X-ray binaries (XRBs) display spectral state transitions that are accompanied by substantial changes in the hardness, luminosity, and structure of the accretion flow. We developed a GPU-accelerated cooling toolkit for general relativistic magnetohydrodynamic (GRMHD) simulations of accreting black holes that uses texture memory for fast retrieval of pre-computed values. The toolkit incorporates bremsstrahlung, synchrotron, inverse Compton radiation and Coulomb collision processes. We implemented our toolkit into a GRMHD code and used it to simulate a magnetically arrested disk in the context of the XRB low/hard state around a Kerr black hole. We explored the mass accretion rate in the (10^-6-0. 3) Ṁ ₄₃₃ range, where Ṁ ₄₃₃ is the Eddington accretion rate. Our simulations reveal that for low accretion rates (Ṁ 0. 01 Ṁ ₄₃₃), the flow settles into a geometrically thick, low-density, two-temperature hot accretion flow. At higher accretion rates, the flow turns into a cold single-temperature thin disk at r ₈₍ 50 rg. Inside, the disk breaks up into single-temperature thin filaments embedded into a two-temperature hot thick flow. Our GPU texture memory accelerated cooling prescription is 3-5 times faster than the standard radiation M1 closure methods, and 5 times faster than storing the lookup table in global memory.
Motta et al. (Tue,) studied this question.