Abstract
X-ray binaries are known to exhibit different spectral states which are often associated with different black hole accretion modes. The exact geometry and properties of these accretion modes is still uncertain. Recent IXPE measurements of linear polarization of X-ray emission in canonical X-ray binary system Cygnus X-1 allow us to test models for the hard spectral state of accretion in a unique way. We show that general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of accreting stellar black hole in a hard X-ray state may be consistent with the new observational information. In the presented framework, where first-principle models have limited number of free parameters, the polarimetric X-ray observations put constraints on the viewing angle of the inner hot accretion flow.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The direct images of supermassive black holes in M87* and Sgr A* released by the Event Horizon Telescope revealed ring structures which are interpreted as advection dominated relativistically hot accretion flows in an immediate vicinity of the black hole event horizons (Event Horizon Telescope Collaboration et al. 2019, 2022). These images allow us to test complex accretion flow models with unprecedented precision. It is reasonable to ask whether a black hole accretion processes are universal and whether models constructed mostly for the supermassive black holes apply to accreting stellar-mass holes found in our Galaxy.
Stellar-mass black holes accreting matter from a stellar companion (hereafter XRB) are known to display several distinct X-ray spectral states characterized by different luminosities, X-ray spectral indices, multiwavelength components and timing properties (Fender and Belloni 2012; Belloni and Motta 2016). The two main XRB spectral states are a low-hard and a high-soft states. In the low-hard state, the luminosity is lower, the X-ray spectrum has a power-law shape and the source has a persistent radio component. During the high-soft state, the luminosity rises closer the Eddington limit, the spectrum is dominated by a thermal component and the radio counterpart is quenched. Sources can be also found in an intermediate state with thermal and power-law components both visible in the X-ray spectrum. Since many years it is generally understood that the X-ray soft component is produced by a thermal disk and the power-law component is emitted by the so called “corona” which has some, not yet well understood, connection to a relativistic jet emitting photons in the radio (Done et al. 2007). One of the long-standing questions is what is the geometry and properties of this corona and whether the corona is actually a compact jet or rather a part of the disk and, if so, then how is the corona connected to the jet. All these questions are difficult to answer solely observationally because accretion flows in XRBs are spatially unresolved.
In this work we construct a more detailed model of the accretion flow in XRB in the low-hard state using state-of-the-art general relativistic radiation-magnetohydrodynamic (GRRMHD) simulations. Recent Imaging X-ray Polarimetry Explorer (IXPE) measurements of polarization of X-ray emission in XRBs in a hard state (e.g., Krawczynski et al. 2022) opened up a unique opportunity to test accretion flow scenarios in these systems as suggested already by Lightman and Shapiro (1976) and Rees (1975). In this work we show that the GRRMHD simulations indicate consistency with these new observational constraints.
The paper is structured as follows. In Sect. 2 we describe our model for hard-state accretion. In Sect. 3 we show the polarimetric properties of the model. We compare our modeling results to observations and to other models in Sect. 4. We discuss prospects and directions for future studies in Sect. 5.
2 Accretion model for a hard spectral state
We carry out GRRMHD simulations of low-hard state accretion flow onto stellar mass black hole using first-principles GRRMHD code ebhlight (Ryan et al. 2017). The state-of-the-art code couples evolution of plasma and magnetic fields to the evolution of radiation field accounting for radiative looses and radiative forces. The radiation model includes continuum processes such as: synchrotron emission, bremsstrahlung emission, self-absorptions and Compton scatterings of photons. Additionally, the code follows ions and electron temperatures separately allowing for development of two-temperature plasma flows (Ressler et al. 2015). The simulations include evolution of radiation field and hence they are not scale free. Inspired by recently published results on X-ray polarization in XRB Cygnus X-1 (Krawczynski et al. 2022), we scale our model exactly to this system. In our simulation we set black hole mass to \(M_{\mathrm{BH}}=20\mathrm{M_{\odot }}\) (Miller-Jones et al. 2021) and dimensionless spin of the black hole to \(a_{*}=0.9375\) (Fabian et al. 2012; Gou et al. 2014). The GRRMHD simulations require specifying the gas density and magnetic field strength via scaling parameter \(\mathcal{M}_{\mathrm{unit}}\). We run two simulations with \(\mathcal{M}_{\mathrm{unit}} = 1.5 \times 10^{11}\) and \(2.5 \times 10^{11}\). Two values of \(\mathcal{M}_{\mathrm{unit}}\) are chosen to model emission for observers with low (disk face-on) and high (disk edge-on) viewing angles. The mass accretion rates in both models are at the level of \(\dot{m}\equiv \dot{M}/\dot{M}_{\mathrm{Edd}}\sim {\mathrm{few}} \times 10^{-4}\). In X-rays both models are as bright as \(L_{X} \approx 2 \times 10^{36}\) ergs/s consistent with luminosity of Cygnus X-1 during IXPE detection reported by Krawczynski et al. (2022).Footnote 1 The Bolometric luminosities of our models are \(L_{\mathrm{Bol}} \approx 10^{38} \,{\mathrm{ergs\,s^{-1}}}\) (\(L_{\mathrm{Bol}}/L_{ \mathrm{Edd}}\approx 0.05\)).
The simulation initial conditions assume torus in hydrodynamical equilibrium in the equatorial plane of the black hole. The torus is seeded with a loop of poloidal magnetic field. Further details of the simulation setup are provided in the Appendix A.1. It is important to note here that our models initially do not contain a truncated geometrically thin, optically thick diskFootnote 2 component (Shakura and Sunyaev 1973). This means that in our models there is no reflection of the hot accretion flow from the thermal thin disk and hence there is no thermal photons scattering off the hot corona. Our simulations self-consistently evolve magnetic fields. The magnetorotational instability leads to the development of a turbulent accretion flow. The accretion flow is accompanied by wind and a jet with organized magnetic fields. The accretion flow models studied here are magnetically arrested (MADs, Tchekhovskoy et al. 2011, McKinney et al. 2012, and see Appendix A.1). The powers of a jets produced by the MAD simulations are consistent with the estimated power of the relativistic jet observed in Cyg X-1 (\(P_{ \mathrm{jet}}=10^{36}-10^{37}\) ergs/s, Gallo et al. 2005, Russell et al. 2007).
Figure 1 presents maps of plasma density, electron temperature, plasma \(\beta \) (gas to magnetic pressure) and radiation \(\beta _{r}\) (gas to radiation pressure) of the very inner accretion flow at a later time moment of evolution. The very inner accretion flow is geometrically thick, it has two-temperature structure and it is rather strongly magnetised. The optical thickness of the simulations is at most \(\tau _{es} \approx n_{e} \kappa _{es} r_{g} \approx 0.03\) hence the radiation emitted from the gas may have impact on the gas via radiative back reaction rather than Compton forces. Interestingly, in the disk body the regions of colder electrons are spatially coincidental with low plasma \(\beta \) regions suggesting that electrons are cooling rapidly in regions of stronger magnetic fields. These sub-relativistic electrons may enhance Faraday rotation of polarized synchrotron photons (seed photons for Compton process) which may impact the polarization properties of the X-ray emission.
Geometry of the inner accretion flow in the hard spectral state in GRRMHD simulation with \({\mathcal {M}}_{\mathrm{unit}}=1.5 \times 10^{11}\) at time moment \(t=4000\) M. Panels from left to right show the plasma density (in code units), the dimensionless electron temperature, the plasma \(\beta \) parameter indicating dynamically important magnetic fields, and the radiation \(\beta _{r}\) parameter (ratio of gas pressure to radiation pressure). The black contours follow the magnetic field lines
3 Mock spectral energy distribution and polarization
Although the GRRMHD simulations do produce self-consistent emission spectra those spectra do not have polarization information. Hence, we post-process our GRRMHD models with radpol, recently developed, fully relativistic, Monte Carlo code for polarized radiative transfer (Mościbrodzka 2020, 2022). The details of the calculation are provided in the Appendix A.2.
Figure 2 shows polarized broadband SED of the GRRMHD models for a couple of viewing angles around times \(t=4000~{\mathrm{M}}\). The lower viewing angle model has slightly higher \({\mathcal {M}}_{\mathrm{unit}}=2.5\times 10^{11}\) to match the observed X-ray flux and spectral slope. At the shown time the simulation is already relaxed and approximately steady-state with variations due to magnetic turbulence. Each SED has three components: synchrotron hump, X-ray scattered emission and bremsstrahlung of which synchrotron hump and Compton humps are the most prominent. The bottom panels in Fig. 2 show polarimetric information. The fractional linear polarization \(|m|=\sqrt{Q^{2}+U^{2}}/I\) in X-ray band is small, less than 10% for both viewing angles. The fractional linear polarization of the synchrotron emission is naturally higher. The electric vector position angle (\(EVPA \equiv 0.5arg(Q+iU)\)) is zero for close to edge-on inclination angles, meaning that EVPA is aligned with the jet for nearly edge-on view of the model. For lower viewing angles, the EVPA significantly deviates from zero and displays some frequency dependency.
SEDs of the GRRMHD models at a viewing angle for which the SED matches observational constraints (\(i=90\) deg, left panel) and at a viewing angle equal to the estimated inclination angle of the binary system orbital plane (\(i=30\) deg, right panel). Top panels show total intensity spectra, middle panels show the amplitude of the fractional linear polarization \(|m|\%\) and lower panels display the linear polarization angle, EVPA. The gray shaded region marks bandwidth of the X-ray observations (2-8 keV)
Unfortunately, it is not straightforward to measure how much of the X-ray polarization observed in the model is caused by the scatterings or the seed photon polarization because one cannot disentangle different effects easily in these complex simulations (e.g., the scattering process does depend on the polarization of incident photons). Due to the complexity of the model it is even more difficult to explain the behaviour of the EVPA as a function of viewing angle and frequency.
In addition to the spectrum our radiative transfer simulations track the origin of emission detected in different spectral windows. In Fig. 3 we show a snapshot of the GRRMHD simulation together with emission maps produced by our radpol scheme. Notice that in comparison to Fig. 1, Fig. 3 shows the model on larger scales (within 50 M) and plasma density and magnetic field strength are shown in c.g.s. units. Here it is evident that the electrons in the GRRMHD model are mildly relativistic with temperature increasing with latitude. Hence the seed synchrotron emission is produced in the inner parts of the accretion disk but mostly along the jet funnel wall. The inverse-Compton scatterings take place along the jet wall but also in the disk region.
Example maps of plasma density, magnetic field strength (both scaled to Cygnus X-1 system and given in c.g.s. units), dimensionless electron temperature and maps of origin of photons at different energy bands in one snapshot of the GRRMHD simulation. The maps of the photon origin are \(\phi \) angle averaged to improve signal-to-noise ratio in the maps. The color scale in the emission maps (two right-most panels) is proportional to the number of photons emitted. Notice that the size of the maps shown here is larger in comparison to Fig. 1
4 Comparison with observations and other models
The IXPE detected X-ray fractional polarization of Cygnus X-1 in hard state to be on average \(|m|=\)4% (slightly rising with frequency) and the EVPA to be aligned with the radio jet position angle (in our reference system this translates to EVPA=0). Given the simplifications in the model (e.g., the simulations are axisymmetric) our GRRMHD simulations of the hard state are quite consistent with these observational constraints. Low polarization fraction and EVPA aligned with the black hole spin or jet axis can be simultaneously recovered by GRRMHD simulations assuming nearly edge-on viewing angles. In our model the X-ray emission is emitted predominantly by the jet wall although some disk emission is also contributing. For viewing angle of \(i=30\) deg (estimated viewing angle of the Cyg X-1 orbital plane, Miller-Jones et al. 2021) the EVPA can become nearly perpendicular (EVPA\(\approx -100\) deg) compared to the observed value. Overall our results indicate a possible misalignment between the jet axis (black hole spin axis) and the angular momentum of the binary or the outer accretion disk.
In our model the X-ray emission and its polarization are produced by a partially outflowing matter (jet wall). This is consistent with ideas recently put forward by others (Poutanen et al. 2023; Dexter and Begelman 2023) however the first-principles simulations presented here prefer nearly edge-on viewing angles. Our results are in contradiction with the IXPE data interpretation presented in Krawczynski et al. (2022) where authors argue that the measurement made in Cygnus X-1 can be interpreted as a “sandwich”-like corona above the thin disk or “corona” alone (with either external Comptonization or self-synchrotron Comptonization) observed at intermediate inclinations angles (\(i=45-60\) deg). For such high viewing angles our model are depolarized in X-rays and EVPA deviates from zero. All these discrepancies can be attributed to different assumptions made when using different modeling techniques. For example, here we do not include the thin disk component, in contrast to the previous work. On the other hand, we include the previously neglected Faraday rotation which very likely plays a role in depolarization process of synchrotron emission in the disk where sub-relativistic electrons exist. Our models also include self-consistently evolved magnetic field geometry which is not included in the aforementioned geometrical models. Alternatively, the current model may be incorrect and another GRMHD scenario should be considered to explain the IXPE measurements.
5 Future prospects
Overall, the accurate detection of the X-ray polarimetry is a powerful new tool for constraining GRMHD simulations of XRBs and accreting supermassive black holes. In this work, we have made a first attempt to model polarimetric characteristics of XRB, utilizing two state-of-the-art numerical techniques. GRRMHD simulations, although already complex, still have several uncertainties that may impact the interpretation of the observations. At least two types of developments could be pursued.
The sub-grid electron thermodynamics in the GRMHD simulations is uncertain. One could consider alternative to currently used electron heating model which assumes purely Alfvénic cascade for energy dissipation. However it is not clear which model for electron heating is realistic for a given GRRMHD simulation (e.g., Howes 2010; Kawazura et al. 2018; Satapathy et al. 2023). Another missing component in the GRRMHD models are the non-thermal electrons which are likely accelerated along the jet. Including the effect of acceleration using first-principles would require developing new modeling approaches.
A more feasible goal would be to calculate broadband polarization of a fully three dimensional simulations such as the ones presented in Dexter et al. (2021). Our current model in an axisymmetric version of fully three dimensional models M3/M4 in Dexter et al. (2021) but notice that the our model differs in the black hole mass, evolution time and the size of the gas-radiation interaction zone (assumed to be 40 M by Dexter et al. (2021)). Three dimensional simulations allow to study emission from models with accretion disk tilted with respect to the equatorial plane so they may be helpful to explain the discrepancy between the viewing angle of the inner accretion flow and the orbital plane of the system. Finally, only three dimensional GRRMHD models are suitable for time variability studies.
Data Availability
No datasets were generated or analysed during the current study.
Notes
The IXPE X-ray flux reported in Krawczynski et al. 2022 is \(F_{X} = 3.5-4.5 \times 10^{-9} {\mathrm{ergs\, cm^{-2}\, s^{-1}}}\). Given the distance to Cyg X-1, \(d=2.22\) kpc, the source X-ray luminosity \(L_{X} = 2-2.5\times 10^{36} {\mathrm{ergs \, s^{-1}}}\). The reported photon index is \(\Gamma \approx 1.62\) meaning that the X-ray part of the spectrum in \(\nu L_{\nu}\) vs \(\nu \) plots should have negative slope of −0.38.
Formally a thin disk could be self-consistently formed in the GRRMHD model inner region but our accretion rates are still too low for this to happen.
References
Belloni, T.M., Motta, S.E.: Transient black hole binaries. In: Bambi, C. (ed.) Astrophysics of Black Holes: From Fundamental Aspects to Latest Developments, pp. 61 (2016). arXiv:1603.07872. https://doi.org/10.1007/978-3-662-52859-4_2
Dexter, J., Begelman, M.C.: A relativistic outflow model of the X-ray polarization in Cyg X-1. arXiv e-prints https://doi.org/10.48550/arXiv.2308.01963 [astro-ph.HE] (2023)
Dexter, J., Scepi, N., Begelman, M.C.: Radiation GRMHD simulations of the hard state of black hole X-ray binaries and the collapse of a hot accretion flow. Astrophys. J. Lett. 919(2), L20 (2021). arXiv:2109.06239 [astro-ph.HE]. https://doi.org/10.3847/2041-8213/ac2608
Done, C., Gierliński, M., Kubota, A.: Modelling the behaviour of accretion flows in X-ray binaries. Everything you always wanted to know about accretion but were afraid to ask. Astron. Astrophys. Rev. 15(1), 1–66 (2007). arXiv:0708.0148 [astro-ph]. https://doi.org/10.1007/s00159-007-0006-1
Event Horizon Telescope Collaboration et al.: First M87 event horizon telescope results. I. The shadow of the supermassive black hole. Astrophys. J. Lett. 875(1), L1 (2019). arXiv:1906.11238 [astro-ph.GA]. https://doi.org/10.3847/2041-8213/ab0ec7
Event Horizon Telescope Collaboration et al.: First Sagittarius A* event horizon telescope results. I. The shadow of the supermassive black hole in the center of the Milky Way. Astrophys. J. Lett. 930(2), L12 (2022). https://doi.org/10.3847/2041-8213/ac6674
Fabian, A.C., Wilkins, D.R., Miller, J.M., et al.: On the determination of the spin of the black hole in Cyg X-1 from X-ray reflection spectra. Mon. Not. R. Astron. Soc. 424(1), 217–223 (2012). arXiv:1204.5854 [astro-ph.HE]. https://doi.org/10.1111/j.1365-2966.2012.21185.x
Fender, R., Belloni, T.: Stellar-mass black holes and ultraluminous X-ray sources. Science 337(6094), 540 (2012). arXiv:1208.1138 [astro-ph.HE]. https://doi.org/10.1126/science.1221790
Fishbone, L.G., Moncrief, V.: Relativistic fluid disks in orbit around Kerr black holes. Astrophys. J. 207, 962–976 (1976). https://doi.org/10.1086/154565
Gallo, E., Fender, R., Kaiser, C., et al.: A dark jet dominates the power output of the stellar black hole Cygnus X-1. Nature 436(7052), 819–821 (2005). arXiv:astro-ph/0508228 [astro-ph]. https://doi.org/10.1038/nature03879.
Gou, L., McClintock, J.E., Remillard, R.A., et al.: Confirmation via the continuum-fitting method that the spin of the black hole in Cygnus X-1 is extreme. Astrophys. J. 790(1), 29 (2014). arXiv:1308.4760 [astro-ph.HE]. https://doi.org/10.1088/0004-637X/790/1/29
Howes, G.G.: A prescription for the turbulent heating of astrophysical plasmas. Mon. Not. R. Astron. Soc. 409, L104–L108 (2010). arXiv:1009.4212 [astro-ph.HE]. https://doi.org/10.1111/j.1745-3933.2010.00958.x
Kawazura, Y., Barnes, M., Schekochihin, A.A.: Thermal disequilibration of ions and electrons by collisionless plasma turbulence. arXiv e-prints arXiv:1807.07702 [physics.plasm-ph] (2018)
Krawczynski, H., Muleri, F., Dovčiak, M., et al.: Polarized X-rays constrain the disk-jet geometry in the black hole X-ray binary Cygnus X-1. Science 378(6620), 650–654 (2022). arXiv:2206.09972 [astro-ph.HE]. https://doi.org/10.1126/science.add5399
Lightman, A.P., Shapiro, S.L.: Polarization of X-rays from Cygnus X-1: a test of the accretion disk model. Astrophys. J. 203, 701–703 (1976). https://doi.org/10.1086/154131
McKinney, J.C., Tchekhovskoy, A., Blandford, R.D.: General relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around black holes. Mon. Not. R. Astron. Soc. 423(4), 3083–3117 (2012). arXiv:1201.4163 [astro-ph.HE]. https://doi.org/10.1111/j.1365-2966.2012.21074.x
Miller-Jones, J.C.A., Bahramian, A., Orosz, J.A., et al.: Cygnus X-1 contains a 21-solar mass black hole—implications for massive star winds. Science 371(6533), 1046–1049 (2021). arXiv:2102.09091 [astro-ph.HE]. https://doi.org/10.1126/science.abb3363
Mościbrodzka, M.: General relativistic polarized radiative transfer with inverse-Compton scatterings. Mon. Not. R. Astron. Soc. 491(4), 4807–4815 (2020). arXiv:1910.02413 [astro-ph.HE]. https://doi.org/10.1093/mnras/stz3329
Mościbrodzka, M.A.: Polarization-sensitive Compton scattering by accelerated electrons. Astrophys. J. Suppl. Ser. 263(1), 6 (2022). arXiv:2207.10487 [astro-ph.HE]. https://doi.org/10.3847/1538-4365/ac972c
Poutanen, J., Veledina, A., Beloborodov, A.M.: Polarized X-rays from windy accretion in Cygnus X-1. Astrophys. J. Lett. 949(1), L10 (2023). arXiv:2302.11674 [astro-ph.HE]. https://doi.org/10.3847/2041-8213/acd33e
Rees, M.J.: Expected polarization properties of binary X-ray sources. Mon. Not. R. Astron. Soc. 171, 457–465 (1975). https://doi.org/10.1093/mnras/171.3.457
Ressler, S.M., Tchekhovskoy, A., Quataert, E., et al.: Electron thermodynamics in GRMHD simulations of low-luminosity black hole accretion. Mon. Not. R. Astron. Soc. 454, 1848–1870 (2015). arXiv:1509.04717 [astro-ph.HE]. https://doi.org/10.1093/mnras/stv2084
Russell, D.M., Fender, R.P., Gallo, E., et al.: The jet-powered optical nebula of Cygnus X-1. Mon. Not. R. Astron. Soc. 376(3), 1341–1349 (2007). arXiv:astro-ph/0701645 [astro-ph]. https://doi.org/10.1111/j.1365-2966.2007.11539.x
Ryan, B.R., Ressler, S.M., Dolence, J.C., et al.: The radiative efficiency and spectra of slowly accreting black holes from two-temperature GRRMHD simulations. Astrophys. J. Lett. 844, L24 (2017). arXiv:1707.04238 [astro-ph.HE]. https://doi.org/10.3847/2041-8213/aa8034
Satapathy, K., Psaltis, D., Ozel, F.: Global Electron Thermodynamics in Radiatively Inefficient Accretion Flows. arXiv e-prints https://doi.org/10.48550/arXiv.2304.10684 [astro-ph.HE] (2023)
Shakura, N.I., Sunyaev, R.A.: Black holes in binary systems. Observational appearance. Astron. Astrophys. 24, 337–355 (1973)
Tchekhovskoy, A., Narayan, R., McKinney, J.C.: Efficient generation of jets from magnetically arrested accretion on a rapidly spinning black hole. Mon. Not. R. Astron. Soc. 418(1), L79–L83 (2011). arXiv:1108.0412 [astro-ph.HE]. https://doi.org/10.1111/j.1745-3933.2011.01147.x
Acknowledgements
The author thanks Jiri Svoboda, Małgosia Sobolewska, George Wong, Maciek Wielgus, Pierre-Olivier Petrucci, Henric Krawczynski for useful discussions and comments on observations of XRBs and simulations of accreting black holes. Author also thanks Jason Dexter for discussion at early stage of the project in 2021 and sharing one time slice of the 3-D model M4.
Funding
We gratefully acknowledge the HPC RIVR consortium (www.hpc-rivr.si) and EuroHPC JU (eurohpc-ju.europa.eu) for funding this research by providing computing resources of the HPC system Vega at the Institute of Information Science (www.izum.si).
Author information
Authors and Affiliations
Contributions
M.M wrote the main manuscript text, prepared all figures and carried out all calculations.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A: Methods
Appendix A: Methods
1.1 A.1 GRRMHD simulations
The GRRMHD simulations of accretion onto a black hole are carried out by means of public code ebhlightFootnote 3 developed by Ryan et al. (2017). ebhlight integrates equations of ideal radiation-magnetohydrodynamics in Kerr space time where multi-frequency radiative transport is integrated out using Monte Carlo techniques. The code additionally tracks total and electron entropy equations which allow to calculate total and electron heating rates as the model evolves (Ressler et al. 2015). The latter requires a model for distribution of dissipated energies between electrons and protons. We adopt a sub-grid model for turbulent heating of electrons and ions (see e.g., Howes 2010) as implemented in ebhlight by the code authors. The electron model also accounts for Coulomb couplings and radiative cooling. The self-consistently evolved radiative transfer includes synchrotron emission and absorption, Compton processes and bremsstrahlung. The ebhlight scheme allows for all possible matter-radiation interactions (four-momentum exchange during emission, absorption and scatterings or, in other words, cooling and radiative forces).
Our initial conditions is the standard Fishbone and Moncrief torus, often used to initialize advection dominated accretion flow models (Fishbone and Moncrief 1976). The initial torus is seeded with single loop of poloidal magnetic field. The initial magnetic field configuration is defined using the vector potential:
where \(r\) is the radius and \(\theta \) is the polar angle in range from 0 deg to 180 deg. This initial magnetic field configuration leads to the development of the Magnetically Arrested Disks (MADs, (Tchekhovskoy et al. 2011; McKinney et al. 2012)). We integrate the MAD solution using ebhlight until \(t_{f}=4500\) M, where \(M=GM/c^{3}\) is the time unit in geometriesed unit system (\(G=c=1\)). Figure 4 shows evolution of mass accretion rate and magnetic flux over time confirming large accumulation of the magnetic flux at the horizon consistent with the MAD solutions obtained by other authors.
Evolution of the mass accretion rate (\(\dot{M}\)) and magnetic field flux (\(\phi _{\mathrm{B}}\equiv \Phi _{\mathrm{B}}/\sqrt{\dot{M}c^{2}}\)) near the event horizon of the black hole in our axisymmetric GRRMHD simulation. The magnetic field flux at the black hole horizon is normalized by the square root of the accretion rate at the horizon. The horizontal dashed line marks saturation level \(\phi _{B}=15\) for which the accretion flows becomes magnetically arrested (MAD)
In the GRRMHD simulation, our grid is two dimensional and has resolution of \(240\times 240\) zones in radial and \(\theta \) direction. We use modified Kerr Schild coordinates which are logarithmic in radius, and mixed-modified in \(\theta \) direction for a better resolution of the inner part of the disk in the equatorial plane of the black hole (for details see the ebhlight code paper). We assume a ceiling, \(\Theta _{e,max}=10^{3}\), and a floor, \(\Theta _{e,min}=10^{-3}\), for electron temperatures through out the computational domain. The emissivity formulas for synchrotron processes used in ebhlight are inaccurate for \(\Theta _{e}<0.3\) therefore all synchrotron emission from plasma with sub-relativistic electron temperature is suppressed. However, the synchrotron emission from plasma with \(\Theta _{e}<0.3\) is negligible. Compton scatterings, on the other hand, are allowed for all possible electron temperatures. In the GRRMHD simulation, we allow for interactions between gas and radiation field only within 100 M of the simulation domain. In the outer regions these interactions are negligible when it comes to the gas dynamics. The outer radius of the simulation extends to \(R_{\mathrm{out}}=1000\) M.
1.2 A.2 Polarized radiative transfer simulations
General relativistic polarized radiative transfer is carried out using postprocessing Monte Carlo scheme radpol (Mościbrodzka 2020, 2022). radpol includes polarized synchrotron emission, self-absorption and polarization sensitive Compton scatterings. In the postprocessing we make the following simplifications. We assume that synchrotron emissivity for Stokes V and that Faraday conversion are null. These processes are weak and not important for X-ray polarization predictions. (Notice that the Faraday rotation is included and it affects the linear polarization.) To be self-consistent with GRRMHD simulation, we assume that the electron distribution function is thermal. Synchrotron emission from regions with \(\Theta _{e}<0.3\) is suppressed similar to GRRMHD model, but Compton scatterings on relativistically cold electrons are allowed, as well as Faraday rotation. All postprocessing radiative transfer also makes a cut to exclude uncertain emission from highly magnetized regions (defined with \(\sigma \equiv B^{2}/ 4\pi \rho c^{2} >1\)) where ebhlight (similarly to all other GRMHD and GRRMHD codes) fails to accurately solve for internal gas energies). In the post-processing radiative transfer step we calculate emission including entire GRRMHD computational domain also outside of 100 M. All presented spectra are \(\phi \) angle averaged but are not folded around equator.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Moscibrodzka, M. What is the hard spectral state in X-ray binaries? Insights from GRRMHD accretion flows simulations and polarization of their X-ray emission. Astrophys Space Sci 369, 68 (2024). https://doi.org/10.1007/s10509-024-04333-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10509-024-04333-3