# Analytical ray-tracing of synchrotron emission around accreting black holes

Alexandra Veledina<fnref n="1" />,<fnref n="2" />,\* and Matthieu Pélissier<fnref n="3" />

D

Received 10 October 2024 / Accepted 12 December 2024

## ABSTRACT

Polarimetric images of accreting black holes encode important information about laws of strong gravity and relativistic motions of matter. Recent advancements in instrumentation have enabled such studies of two objects: the supermassive black holes M87* and Sagittarius A*. Light coming from these sources is produced by a synchrotron mechanism whose polarization is directly linked to magnetic field lines, and propagates toward the observer in a curved spacetime. We studied the distortions of the gas image by employing the analytical ray-tracing technique for polarized light ARTPOL, which is adapted for the case of synchrotron emission. We derived analytical expressions for fast conversion of the intensity or flux, polarization degree, and polarization angle from the local coordinates to those of the observer. We placed an emphasis on the nonzero matter elevation above the equatorial plane and noncircular matter motions. Applications of the developed formalism include static polarimetric imaging of the black hole vicinity and dynamic polarimetric signatures of matter close to the compact object.

**Key words.** accretion, accretion disks – gravitational lensing: strong – polarization – methods: analytical – stars: black holes – galaxies: active

## 1. Introduction

Accretion onto a compact object, a black hole (BH) or a neutron star (NS), leads to a substantial energy release and excites interest in studying the details of energy liberation. Thanks to recent technological advances in radio, submillimeter, and near-infrared (NIR) interferometry, we are able to probe the emission regions with micro-arc-second angular resolution. The highest resolution achieved in such studies, achieved by the Event Horizon Telescope (EHT), allows one to produce images that are a few Schwarzschild radii ($R_S$) in size. Two supermassive BHs – M87* (Event Horizon Telescope Collaboration 2019) and Sgr A* (Event Horizon Telescope Collaboration 2022) – have now been resolved down to the characteristic scales of energy release. Emission at the EHT wavelengths is thought to be shaped by synchrotron mechanism, indicating the presence of high-energy electrons and magnetic fields of significant magnitude near the event horizon. Additional information coming from polarimetric studies (Event Horizon Telescope Collaboration 2021, 2024a) suggests a rather ordered topology of the field.

Both M87* and Sgr A* are characterized by a low level of accretion, where an interplay between plasma processes and strong gravity jointly affect the observed images and average spectro-polarimetric signatures. Disentangling these effects is a difficult and degenerate problem. A promising avenue through which to understand the role of various effects is to use the time-dependent information.

Variability studies suggest an intermittent nature of accretion onto Sgr A*: the stable quiescent level is interrupted by

incursions into bright flares. During the flare, the observed radio/submillimeter, NIR, and X-ray fluxes increase by about an order of magnitude (Baganoff et al. 2001; Genzel et al. 2003; Eckart et al. 2006a; Ponti et al. 2017; GRAVITY Collaboration 2020a; Wielgus et al. 2022). The flares are characterized by a high linear polarization degree (PD, reaching 40%) and rotating polarization angle (PA) over the course of the flare (Eckart et al. 2006b; Trippe et al. 2007). The flares occur up to a few times per day. Their average duration is about an hour, this being consistent with the characteristic Keplerian timescale near the event horizon, although rich substructure morphology with timescales of minutes has been detected (Dodds-Eden et al. 2009). Individual flares are now being resolved in both astrometric and polarimetric spaces; two flares have good simultaneous coverage in both domains (GRAVITY Collaboration 2018, 2020b; Wielgus et al. 2022; GRAVITY Collaboration 2023).

Several scenarios for the origin of flares have been discussed, either based on a hot spot orbiting the BH close to the innermost stable circular orbit (ISCO, Genzel et al. 2003; Wielgus et al. 2022; Vincent et al. 2024) or an expanding blob of plasma ejected from the BH vicinity (Yusef-Zadeh et al. 2006; Vincent et al. 2014; Aimar et al. 2023). More sophisticated models have also been considered, such as a complex orbiting emission pattern formed, for example, by the local perturbations or by the interplay between the outflow and surrounding disk (Matsumoto et al. 2020). Finally, it remains elusive whether all flares can be attributed to a single mechanism or whether, rather, similar flares can be caused by different physical phenomena. In this context, it is important to note the recent change in the statistical properties of flaring activity in Sgr A*, likely owing to the close approaches by S0-2 stars and the G2 cloud (Do et al. 2009, 2019).

<footnote n="1">¹ Department of Physics and Astronomy, FI-20014 University of Turku, Finland</footnote>

<footnote n="2">² Nordita, Stockholm University and KTH Royal Institute of Technology, Hannes Alfvéns väg 12, SE-10691 Stockholm, Sweden</footnote>

<footnote n="3">³ Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote>\* Corresponding author; alexandra.veledina@gmail.com</footnote>
Intrinsic spectra and polarimetric properties of the hot plasma near the event horizon are altered on the way to the distant observed by the relativistic aberration, light- bending and frame-dragging effects (Connors & Stark 1977; Stark & Connors 1977; Pineault 1977; Pineault & Roeder 1977; Connors et al. 1980). While the PD is an invariant quantity, the polarization signatures of each point or element in the accretion disk is expected to be different because of the PA rotation and alteration of the viewing angle; both of these effects depend on the location of the element in the disk, BH spin, and observer inclination. To understand the intrinsic properties of flares, as well as the geometry of the quiescent-level accretion onto Sgr A* and M87*, it is important to distinguish between the effects of strong gravity, radiative transfer, and plasma processes. Calculations of these effects often involve parallel transport of the Stokes vector along the null geodesics; for example, via the conservation of the Walker & Penrose (1970) constant. The geodesics can in turn be computed using the ray-tracing technique (often used, e.g., Dexter & Agol 2009; Vos et al. 2022). However, this straightforward routine is often nonintuitive with respect to the prediction of parameters from the observed properties; it is also rather computationally expensive.

Recently, a number of semi-analytical methods have been developed to serve the purpose of computationally fast studies of the effects of model parameters on the observed images, polarimetry, and variability. Polarimetric images of the M87* surroundings have been obtained (Narayan et al. 2021) by evaluating the Walker-Penrose constant for an approximate lightbending formula of Beloborodov (2002). An extension to this model for the case of a Kerr BH was also considered (Gelles et al. 2021). Furthermore, the effects of relativistic aberration have been isolated by considering the polarimetric signatures of orbiting hot spots in Minkowski space (Vincent et al. 2024). The aforementioned methods consider motion of matter in the equatorial plane; more recently, images of the vertically extended source have also been computed (Chang et al. 2024). The effects of the nonzero matter elevation on polarimetric properties have not been given detailed consideration.

In this work, we study the importance of effects of the non-equatorial and noncircular motion on the observed imaging and polarimetric characteristics. We utilize the approach of analytical ray tracing for spectro-polarimetric (ARTPOL) calculations and give explicit expressions for the PA rotation. The method is based on calculating the rotation of the PA along the photon trajectory using the higher-order approximation to the lightbending formula (Poutanen 2020). It was proved to be robust in the context of the X-ray polarization formed in the atmospheres of accreting NSs (Loktev et al. 2020), accretion disks around Schwarzschild (Loktev et al. 2022), and Kerr (Loktev et al. 2024) BHs. We extend the previous calculations for (i) the cases of polarized synchrotron emission, (ii) noncircular motion of matter, and (iii) nonzero matter elevation above the equatorial plane and nonzero thickness of the disk. We consider the case of the Schwarzschild metric; however, general expressions derived in this work can also be used in the case of the Kerr metric.

The first extension introduces an extra freedom in terms of magnetic field orientation and adds azimuthal dependence of local quantities. The second extension alters the azimuthal angles where Doppler effects are important. Addition of the third effect plays an important role in comparison to the observations.

The paper is organized as follows. In Sect. 2, we describe the geometry and introduce the coordinate system and auxiliary quantities. We also give an explicit expression for the rotation of PA, as a function of the azimuth, caused by general and spe-

<figure><img src="image_1.png" /><figcaption>Fig. 1. Geometry of the accreting matter around the BH. The fluid element is described by the radius vector, $\hat{r}$, and velocity, $\hat{v}$. Its trajectory is shifted from the midplane of the accretion disk by the angle, $\eta$. The observer is located in the direction of $\hat{d}$, making an angle, $i$, with the disk normal, $\hat{n}$.</figcaption></figure>

cial relativity (GR and SR) as well as by the changes in the local magnetic field orientation. The vector form for calculation of the aforementioned rotation angles, a detailed comparison to the previous methods, and approximations are given in the appendix. In Sect. 3, we show examples of the calculations, with an emphasis on the difference between the resulting quantities with the new features included. In Sect. 4, we consider applications of the developed approach to the BH images and spectro-polarimetric properties of flares. We also detail the difference between our results and previous semi-analytical models. In Sect. 5, we summarize our findings.

## 2. Model setup

In this section, we introduce the formalism describing the accretion geometry and transformation of the basic vectors affected by the strong gravity and fast motions of matter. We note that the formalism partly repeats the description given in Loktev et al. (2022, 2024); however, we give the full description here for completeness. Motivated by the earlier studies of M87* and Sgr A* images and flares showing the secondary effect of the BH spin on the observed characteristics (e.g., Gelles et al. 2021), we considered the case of a non-spinning BH.

### 2.1. Accretion geometry and observed flux

We considered accreting matter around the BH in the Schwarzschild metric. The matter can either constitute a compact blob (hot spot) or form an extended structure (accretion disk or its narrow ring). In both cases, we identified the basic plane related to the equatorial (mid)plane of the accretion disk and introduced the Cartesian coordinates related to this plane (see Fig. 1). We set the z axis to coincide with the normal to the disk plane, $\hat{n}$. The x axis lies along the projection of the line of sight on the disk. The unit vector, $\hat{d}$, points toward the observer and the radius vector of the considered point in the disk (or the hot spot) is $\mathbf{r}$ (with the corresponding unit vector, $\hat{r}$) and the unit azimuthal vector is $\hat{\varphi}$. These vectors have the following coordinates:

$$ \hat{n} = (0, 0, 1), \qquad (1) $$

$$ \hat{d} = (\sin i, 0, \cos i), \qquad (2) $$

$$ \hat{r} = (\cos \varphi \cos \eta, \sin \varphi \cos \eta, \sin \eta), \qquad (3) $$

$$ \hat{\varphi} = (-\sin \varphi, \cos \varphi, 0), \qquad (4) $$

where $\eta$ is the angle between the radius-vector, $\hat{r}$, and the equatorial plane, $\varphi$ is the azimuthal angle of the projection of $\hat{r}$ to the
x-y plane, measured from the x axis, and i is the inclination of the disk plane to the line of sight. The vectors, $\hat{\mathbf{r}}$ and $\hat{\mathbf{o}}$, form an angle, $\psi$:

$$ \cos \psi \equiv \hat{\mathbf{r}} \cdot \hat{\mathbf{o}} = \sin i \cos \varphi \cos \eta + \cos i \sin \eta. \quad (5) $$

Due to the effects of the light bending, the photons reaching the observer along $\hat{\mathbf{o}}$ are leaving the disk in the other direction, which we denote as $\hat{\mathbf{k}}_0$. Because the photon trajectories lie in a plane in the Schwarzschild metric, the wave vector, $\hat{\mathbf{k}}_0$, can be described by a linear combination of the radius-vector and the observer direction:

$$ \hat{\mathbf{k}}_0 = [\sin \alpha \hat{\mathbf{o}} + \sin(\psi - \alpha)\hat{\mathbf{r}}]/\sin \psi, \quad (6) $$

where $\alpha$ is the angle between the radius vector and the direction toward the observer,

$$ \cos \alpha \equiv \hat{\mathbf{r}} \cdot \hat{\mathbf{k}}_0. \quad (7) $$

The angle, $\alpha$, is related to $\psi$ through the light-bending integral or an approximate light-bending formula (Pechenick et al. 1983; Beloborodov 2002; Poutanen 2020), introduced explicitly below. The wave vector makes an angle, $\zeta$, with the disk normal,

$$ \cos \zeta \equiv \hat{\mathbf{n}} \cdot \hat{\mathbf{k}}_0 = \frac{\sin \alpha \cos i + \sin(\psi - \alpha) \sin \eta}{\sin \psi}. \quad (8) $$

Plasma orbits around the BH with Keplerian velocity and additionally gradually moves toward the compact object. The unit velocity vector is thus not purely parallel to the azimuthal vector. This deviation is described by the angle, $\sigma$, with respect to the equatorial radius vector, $\hat{\mathbf{r}}_e = (\cos \varphi, \sin \varphi, 0)$:

$$ \cos \sigma \equiv \hat{\mathbf{r}}_e \cdot \hat{\mathbf{v}}. \quad (9) $$

We note that we did not consider vertical motion of matter (along the z axis). The unit velocity vector is described by

$$ \hat{\mathbf{v}} = (\cos(\varphi + \sigma), \sin(\varphi + \sigma), 0). \quad (10) $$

Purely Keplerian rotation can be recovered by choosing $\sigma = 90^\circ$ or $270^\circ$; matter moves toward the BH for $90^\circ < \sigma < 270^\circ$. For the magnitude of the dimensionless velocity, $\beta = v/c$, relative to the distant observer, we adopted the relation of Luminet (1979). In this work, we have considered the case in which the matter above the disk at a height of $h = r \sin \eta$, where $r = R/R_S$ is the circumferential radius measured in units of Schwarzschild radii, $R_S = 2GM/c^2$ ($M$ is the mass of the compact object), and $\eta \neq 0$, moves with Keplerian velocity corresponding to its equatorial radius, $r_e = \sqrt{r^2 - h^2}$. We have

$$ \beta = \frac{1}{\sqrt{2(r_e - 1)}} \quad (11) $$

and the Lorentz factor is then

$$ \gamma \equiv \frac{1}{\sqrt{1 - \beta^2}} = \sqrt{\frac{r_e - 1}{r_e - 1.5}}. \quad (12) $$

The angle between the photon momentum and velocity vector is expressed as

$$ \cos \xi \equiv \hat{\mathbf{v}} \cdot \hat{\mathbf{k}}_0 = \frac{\sin \alpha}{\sin \psi} \sin i \cos(\varphi + \sigma) + \frac{\sin(\psi - \alpha)}{\sin \psi} \cos \eta \cos \sigma. \quad (13) $$

Fast motions of matter close to the BH cause the photon trajectory to be altered by relativistic aberration. To account for this effect in the observed polarization properties, we considered quantities in a local fluid comoving frame. We denote quantities measured in this frame with primes. The photon vector in the comoving frame is

$$ \hat{\mathbf{k}}'_0 = \delta[\hat{\mathbf{k}}_0 - \gamma\beta\hat{\mathbf{v}} + (\gamma - 1)\hat{\mathbf{v}}(\hat{\mathbf{v}} \cdot \hat{\mathbf{k}}_0)], \quad (14) $$

where $\delta$ is the Doppler factor,

$$ \delta \equiv \frac{1}{\gamma(1 - \beta\hat{\mathbf{v}} \cdot \hat{\mathbf{k}}_0)} = \frac{1}{\gamma(1 - \beta \cos \xi)}. \quad (15) $$

In the fluid frame, the photon makes an angle, $\xi'$, with the velocity vector,

$$ \cos \xi' \equiv \hat{\mathbf{v}} \cdot \hat{\mathbf{k}}'_0 = \frac{\cos \xi - \beta}{1 - \beta \cos \xi} = \gamma \delta(\cos \xi - \beta). \quad (16) $$

The angle between the photon momentum and the local normal in the disk is

$$ \cos \zeta' \equiv \hat{\mathbf{n}} \cdot \hat{\mathbf{k}}'_0 = \delta \cos \zeta. \quad (17) $$

To account for the light-bending effects and to transform the local emission and polarization to the observer frame, we need a relation between the angle, $\alpha$, of the outgoing photon and the observer direction, $\psi$. We used the approximate relation (Poutanen 2020)

$$ \cos \alpha \approx 1 - (1 - 1/r)y \left\{ 1 + \frac{y^2}{112r^2} - \frac{ey}{100r} \left[ \ln \left( 1 - \frac{y}{2} \right) + \frac{y}{2} \right] \right\}, \quad (18) $$

where $y = 1 - \cos \psi$. This relation gives accurate results given that the photon trajectories lie in the plane, which is strictly satisfied in the Schwarzschild metric. The transition from the Schwarzschild metric to the Kerr metric can be made by accounting for two additional effects: (i) a change in the orbital velocity profile and (ii) an additional twist in the light trajectories caused by the frame-dragging effects. While the former step is straightforward as the proper velocity profile is known, the latter effect is not accounted for by the present light-bending approximation. However, the second effect introduces deviations only for the trajectories that lie very close to the spinning BH, $<3R_S$. Significant errors can be expected for the trajectories originating from the parts of the disk behind the BH ($\varphi \sim 180^\circ$) and a combination of a high spin parameter ($a \gtrsim 0.8$) and a high inclination ($i \gtrsim 60^\circ$). At the same time, the integrated PD and PA in these cases show modest deviations (0.5% and $7^\circ$), becoming comparable to the observational uncertainties (Loktev et al. 2024). For the cases of prograde spins, $a < 0.8$, retrograde spins up to $a = -1$, and/or orbits at $\sim 3R_S$ or farther (comparable to the scales of the BH shadow probed by EHT), the current light-bending approximation leads to a precision that is higher than the observational errors.

The accretion flow is assumed to be optically thin to its own synchrotron emission and its spectral energy distribution follows the power-law dependence $I_E \propto E^{-\alpha_E}$. This means that the observer sees the whole volume of the flow, rather than its surface emission. To treat this, we split the accretion flow in surfaces in cylindrical coordinates, with the z axis being perpendicular to the equatorial plane of the accretion disk, and computed the total flux from each surface parallel to the equatorial plane by further splitting it into segments in the radial
and azimuthal direction. The observed flux of each such segment depends on the solid angle it occupies on the sky (see more details in Poutanen & Beloborodov 2006; Loktev et al. 2022),

$$d\Omega = \frac{dS \cos \zeta}{D^2} \mathcal{L} = \frac{R_S^2}{D^2} \frac{r dr d\varphi}{\sqrt{1-1/r}} \mathcal{L} \cos \zeta, \quad (19)$$

where the factor $dS \cos \zeta$ comes from the projection of the segment, D is the distance, and $\mathcal{L}$ is the lensing factor (Beloborodov 2002; Poutanen 2020),

$$\mathcal{L} \equiv \frac{1}{1-1/r} \frac{d \cos \alpha}{d \cos \psi} \approx \quad (20)$$

$$\approx 1 + \frac{3y^2}{112r^2} - \frac{ey}{100r} \left[ 2 \ln \left( 1 - \frac{y}{2} \right) + y \frac{1 - 3y/4}{1 - y/2} \right]. \quad (21)$$

The segment area projected onto the plane orthogonal to the photon propagation is the Lorentz invariant: $dS' \cos \zeta' = dS \cos \zeta$. Hence, the observed solid angle can be expressed through the area of the disk segment as

$$d\Omega = \frac{\delta dS' \cos \zeta}{D^2} \mathcal{L}. \quad (22)$$

The latter expression will be used for computing the fluxes from the moving spots (blobs), for which the intrinsic area $dS'$ remains constant during the orbital motion. This area can be expressed through the comoving coordinates $(r', \varphi')$ as (see also Nättilä & Pihajoki 2018; Bogdanov et al. 2019; Suleimanov et al. 2020)

$$dS' = \frac{\gamma R_S^2 r_e' dr_e' d\varphi'}{\sqrt{1-1/r}}. \quad (23)$$

## 2.2. Transformation of polarized synchrotron emission

The magnetic field is described in the comoving frame, where it has radial, azimuthal, and vertical components. It is convenient to define the field through the components lying on the equatorial plane, $B_{eq}$, and along the disk normal, $B_z$. The equatorial field is further divided into the component lying along the gas motion and orthogonal to that, described by the angle, $\phi'$, of the magnetic field direction relative to the gas motion. In the Cartesian coordinates of the comoving frame, the magnetic field vector is given by

$$\mathbf{B}' = B_{\text{eq}} [\cos \phi' \hat{\mathbf{v}} + \sin \phi' (\hat{\mathbf{n}} \times \hat{\mathbf{v}})] + B_z \hat{\mathbf{n}}. \quad (24)$$

The unit vector, $\hat{\mathbf{B}}' \equiv \mathbf{B}'/|\mathbf{B}'|$, has components $b_{eq}$ and $b_z$. The intensity and polarization of the photon depends on the angle between the magnetic field and photon direction,

$$\cos \theta' \equiv \hat{\mathbf{B}}' \cdot \hat{\mathbf{k}}_0 = b_{\text{eq}}(\cos \phi' \cos \xi' - \delta \tilde{b} \sin \phi') + \delta b_z \cos \zeta, \quad (25)$$

where

$$\tilde{b} = \frac{\sin \alpha}{\sin \psi} \sin i \sin(\varphi + \sigma) + \frac{\sin(\psi - \alpha)}{\sin \psi} \cos \eta \sin \sigma. \quad (26)$$

The radiation field in the comoving frame can be generally described by the Stokes vector,

$$\mathbf{I}'_{E'}(\theta') \equiv \begin{bmatrix} I'_{E'} \\ Q'_{E'} \\ U'_{E'} \end{bmatrix} = I'_{E'} \begin{bmatrix} 1 \\ P_s(E') \cos 2\chi_0 \\ P_s(E') \sin 2\chi_0 \end{bmatrix}, \quad (27)$$

where $I'_{E'}$ is the local rest-frame specific intensity, $P(E')$ corresponds to the PD in this frame, and $\chi_0$ gives the PA relative to the chosen axis. A natural choice for synchrotron emission is the axis going along the magnetic field vector. The basis vectors are

$$\hat{\mathbf{e}}_1 = \frac{\hat{\mathbf{B}}' - \cos \theta' \hat{\mathbf{k}}_0'}{\sin \theta'}, \quad \hat{\mathbf{e}}_2 = \frac{\hat{\mathbf{k}}_0' \times \hat{\mathbf{B}}'}{\sin \theta'}. \quad (28)$$

In this basis, $\chi_0 = 90^\circ$ and synchrotron emission and polarization can be expressed as

$$\mathbf{I}'_{E'}(\theta') = c_0 E'^{-\alpha_E} I'_{\text{ph}} (B \sin \theta')^{1+\alpha_E} \begin{bmatrix} 1 \\ -P_s \\ 0 \end{bmatrix}, \quad (29)$$

where $c_0$ is a constant (e.g., Eq. (6.36) in Rybicki & Lightman 1979), $I'_{\text{ph}}$ is the photon path length through the synchrotron-emitting medium (measured in the comoving frame), and $P_s = (p+1)/(p+7/3)$ gives the PD of the power-law electrons with index $p = 2\alpha_E + 1$ (Ginzburg & Syrovatskii 1965). The typical value considered for the optically thin synchrotron emission is $\alpha_E = 0.7$, while the observations of bright flares in Sgr A* tend to give a somewhat steeper slope, $\alpha_E \approx 1$ (Gillessen et al. 2006). The latter value was used in the modeling of the Sgr A* and M87* (Narayan et al. 2021; Vincent et al. 2024), and we adopted this value for the calculations presented in this work. The advantage of such a representation of the observed Stokes parameters is that it explicitly distinguishes the PA from the PD. The PD is a Lorentz invariant; hence, the transformation of polarization signatures reduces to the tracing of the changes in the photon electric field orientation (as a result of aberration and light bending) along the photon trajectory.

The intensity of the steady matter flow in the comoving frame is related to the observed intensity as

$$I_E = g^3 I'_{E'} \propto g^{3+\alpha_E} E^{-\alpha_E}, \quad (30)$$

where $g \equiv E/E' = \delta \sqrt{1-1/r}$ is the total redshift factor (Luminet 1979; Chen & Shaham 1989) and the latter expression corresponds to the power-law spectrum $I'_{E'} \propto E'^{-\alpha_E}$. The photon path can be expressed as

$$I'_{\text{ph}} = g l_{\text{ph}} = \frac{h_0}{\hat{\mathbf{k}}_0^{(z)}} = \frac{h_0}{\cos \zeta'}, \quad (31)$$

where $h_0$ is the characteristic vertical thickness of each synchrotron-emitting layer in the disk, seen in the observer's frame, and $\hat{\mathbf{k}}_0^{(z)}$ is the z component of the photon momentum vector in the comoving frame.

The observed flux from a segment of the disk is related to the observed intensity of radiation as

$$dF_E = I_E d\Omega = \delta I_E \mathcal{L} \cos \zeta \frac{dS'}{D^2}. \quad (32)$$

Stokes parameters (in flux units) of each layer in the disk were obtained by integrating over the equatorial radii and azimuthal coordinates.

$$\begin{align} \mathbf{F}_E &\equiv \begin{bmatrix} F_I(E) \\ F_Q(E) \\ F_U(E) \end{bmatrix} = \int d\Omega \begin{bmatrix} I_E \\ Q_E \\ U_E \end{bmatrix} = \tag{33} \\ &= \int d\Omega g^3 \mathbf{M}(r, \varphi) \mathbf{I}'_{E'}(\theta') = \nonumber \\ &= \int \frac{dS'}{D^2} g^3 \delta \mathcal{L} \cos \zeta \mathbf{M}(r, \varphi) \mathbf{I}'_{E'}(\theta'), \nonumber \end{align}$$
where the transformation from the vector, $I'_{E'}(\theta')$, to the observed Stokes vector was done through the rotation Mueller matrix,

$$ \mathbf{M}(r, \varphi) = \begin{bmatrix} 1 & 0 & 0 \\ 0 & \cos 2\chi^{\text{tot}} & -\sin 2\chi^{\text{tot}} \\ 0 & \sin 2\chi^{\text{tot}} & \cos 2\chi^{\text{tot}} \end{bmatrix}. \quad (34) $$

Here, $\chi^{\text{tot}}$ corresponds to the total rotation of the polarization (Stokes) vector from the comoving frame to the observer's frame. For the synchrotron emission of the power-law electrons, the expression in Eq. (33) can be simplified to

$$ \begin{aligned} F_E &= f_{E,0} \int d\Omega \frac{g^{3+\alpha_E}}{\delta \cos \zeta} (\sin \theta')^{1+\alpha_E} \begin{bmatrix} 1 \\ -P_s \cos 2\chi^{\text{tot}} \\ -P_s \sin 2\chi^{\text{tot}} \end{bmatrix} = \\ &= f_{E,0} \iint \frac{\gamma R_s^2 r'_c dr'_c d\varphi'}{\sqrt{1-1/r}} g^{3+\alpha_E} \mathcal{L}(\sin \theta')^{1+\alpha_E} \begin{bmatrix} 1 \\ -P_s \cos 2\chi^{\text{tot}} \\ -P_s \sin 2\chi^{\text{tot}} \end{bmatrix}, \end{aligned} \quad (35) $$

where $f_{E,0} = c_0 h_0 E^{-\alpha_E} B^{1+\alpha_E}$. In the latter representation, we have preserved the expression of $dS'$, as it is most useful for moving spots, while the former expression is useful for calculations of the synchrotron-emitting disk image. Below, we have also used the relative Stokes parameters, obtained as $q = F_Q(E)/F_I(E)$ and $u = F_U(E)/F_I(E)$. The developed formalism can be used to construct the spectro-polarimetric images (see also Loktev et al. 2022, 2024). For this, we introduced the Cartesian coordinates, X and Y, in the plane of the sky, scaled to units of $R_s$, as

$$ \begin{bmatrix} X \\ Y \end{bmatrix} = b \begin{bmatrix} -\sin \Phi \\ \cos \Phi \end{bmatrix} = \frac{r}{\sqrt{1-1/r}} \frac{\sin \alpha}{\sin \psi} \begin{bmatrix} \sin \varphi \cos \eta \\ \sin i \sin \eta - \cos i \cos \varphi \cos \eta \end{bmatrix}. \quad (36) $$

Here, $\Phi$ corresponds to the position angle of the point in the plane of the sky, measured counterclockwise from the projection of the disk axis on the sky, and $b = r \sin \alpha / \sqrt{1-1/r}$ corresponds to the impact parameter.

## 2.3. Rotation of polarization angle

The method described above allowed us to compute the observed Stokes parameters once the total rotation of the polarization plane along the photon trajectory ($\chi^{\text{tot}}$) was known. The PA seen by the observer (with respect to the direction of the disk axis) was computed as a sum of the PA of emitted photon and its total rotation:

$$ \chi = \chi_0 + \chi^{\text{tot}}, \quad (37) $$

where $\chi_0 = 90^\circ$ for synchrotron emission. Under the assumption of photon trajectories lying in plane (corresponding to the Schwarzschild metric), the total rotation can be presented as the sum of rotations caused by the light-bending and relativistic aberration effects (the aberration effects need to be computed in the curved space, i.e., accounting for the deflection of the photon path, which is different from the case of the flat space Pineault 1977; Loktev et al. 2022). We followed the methodology described in Loktev et al. (2022, 2024), in which the rotation angles were computed relative to the normal of the disk (in our case, the equatorial plane). As the PA of synchrotron radiation is related to the magnetic field direction, we first need to link this direction to the disk normal via the rotation angle, $\chi^B$. Then, the rotation of polarization plane resulting from light bending is

described by the angle $\chi^{\text{GR}}$ and the special relativity (aberration) effects are accounted for by the angle $\chi^{\text{SR}}$. This leads to the total rotation angle through the sum

$$ \chi^{\text{tot}} = \chi^B + \chi^{\text{GR}} + \chi^{\text{SR}}. \quad (38) $$

The explicit expression for the rotation caused by the effects of special relativity is given in Appendix A:

$$ \tan \chi^{\text{SR}} = -\beta \frac{\cos \alpha \cos \zeta}{\sin^2 \zeta - \beta \cos \xi}. \quad (39) $$

The form of this expression is the same as in Eq. (31) of Loktev et al. (2022); however, the terms are not the same. Namely, in this work we explicitly considered the cases of a nonzero accretion velocity ($\sigma \neq 90^\circ$ or $270^\circ$) and an out-of-equatorial-plane matter location ($\eta \neq 0$); hence, the expressions for the terms $\cos \zeta$, $\sin \zeta$, and $\cos \xi$ are generally different from the ones given in the former work (see Eqs. (8) and (13)). In the case of the flat (Minkowski) space, the expression reduces to

$$ \tan \chi_{\text{flat}}^{\text{SR}} = -\beta \tan i \frac{\sin i \cos \varphi \cos \eta + \cos i \sin \eta}{\sin i - \beta \cos(\varphi + \sigma)}. \quad (40) $$

For $\eta = 0$ and $\sigma = 90^\circ$, it reduces to Eq. (32) of Loktev et al. (2022).

The expression for $\chi^{\text{GR}}$ becomes lengthy once we consider photons originating from the plane that is not equatorial. In coordinate form, it reads as (see Appendix A)

$$ \tan \chi^{\text{GR}} = \cos \eta \frac{\tilde{a}_1 \sin \psi - \tilde{a}_2 \sin i \sin \varphi}{\tilde{a}_1 \tilde{a}_2 - \cos^2 \eta \sin^2 i \sin^2 \varphi \sin \alpha}, \quad (41) $$

where

$$ \tilde{a}_1 = \sin \eta + \cot \psi \sin i \cos \alpha \sin \varphi, \quad (42) $$

$$ \tilde{a}_2 = \sin \alpha [\cos \alpha \cos i - \sin \eta \cos(\psi - \alpha)]. \quad (43) $$

For $\eta = 0$, this expression reduces to Eq. (29) of Loktev et al. (2022).

The rotation angle transforming the polarization direction from the local magnetic field direction relative to the disk normal can be expressed as (see Appendix A)

$$ \chi^B = \arctan \left( \frac{\hat{k}'_0 \cdot (\hat{n} \times \hat{B}')}{\hat{n} \cdot \hat{B}' - (\hat{k}'_0 \cdot \hat{n})(\hat{k}'_0 \cdot \hat{B}')} \right). \quad (44) $$

The expression for $\chi^B$ takes into account both the light-bending effects ($\psi \neq \alpha$) and the Lorentz transformation of the magnetic field (primed quantities). This angle is related to both general and special relativity effects, and depends on $\alpha$ and $\delta$. The explicit expression for the angle $\chi^B$ reads as

$$ \tan \chi^B = \frac{[\tilde{b} \cos \phi' + \gamma \sin \phi' (\cos \xi - \beta)]}{\delta \cos \zeta [\gamma \cos \phi' (\cos \xi - \beta) + \tilde{b} \sin \phi'] + \frac{b_c}{b_{\text{eq}}} [\delta \cos^2 \zeta - \frac{1}{\delta}]}, \quad (45) $$

where $\tilde{b}$ is taken from Eq. (26).

The expression for $\chi^B$ for purely circular rotation ($\sigma = 90^\circ$) in the equatorial plane ($\eta = 0$), far from the compact object ($\gamma \to 1$, $\beta \to 0$) and small, yet nonzero inclinations ($i \to 0$), reduces to

$$ \tan \chi^B \approx \tan \varphi, \quad \text{if } \phi' = 90^\circ, \quad (46) $$

$$ \tan \chi^B \approx -\cot \varphi, \quad \text{if } \phi' = 0. $$
Hence, we obtain

$$ \chi^B \approx \varphi \quad \text{for the radial field,} \qquad (47) $$

$$ \chi^B \approx \varphi - 90^{\circ} \quad \text{for the toroidal field.} $$

The observed PA is then $\chi = \varphi + 90^{\circ}$ for the radial field and $\chi = \varphi$ for the toroidal one, and in both cases the result matches the intuitive expectations.

## 3. Results

The developed formalism allowed us to analytically compute the polarimetric signatures of matter close to the event horizon. The power of such a formalism is that we can separately understand each effect contributing to the total modification of polarization, and hence develop an intuition about the range of the model parameters for the given observed set of data. Furthermore, using the developed formalism, we can compare the predictions of the polarization rotation to the synchrotron source in an orbit around the BH.

Many parameters affect the polarimetric properties of the synchrotron source: distance to the BH, $r$, B-field topology, observer inclination, $i$, the location of the source with respect to the disk plane, $\eta$, and the direction of velocity vector determined by $\sigma$. Generally, with increasing distance to the BH, $r$, all relativistic effects fade and only the topology of the magnetic field affects the PA for different parts of the disk. Other effects are discussed below in detail. To produce the figures, we have considered the point source at given coordinates, unless stated otherwise.

### 3.1. Rotation of polarization angle due to relativistic effects

To disentangle the rotation of the polarization plane caused by relativistic effects from the effect of the B-field orientation changes, we first considered the case of a vertical field, $B = (0, 0, B_z)$. In the absence of the relativistic effects, this case is trivial: the PA is constant and equal to $90^{\circ}$ (for all inclinations except for the degenerate $i = 0$). In the absence of relativistic effects, which can alter the direction of the magnetic field, the angle is the same across the entire disk. We verified that our calculations reproduce this value with high accuracy.

We then included the special and general relativistic effects; we show the contours of constant PA as a function of the disk coordinates $(r, \varphi)$ in Fig. 2. We assumed circular counterclockwise rotation, $\sigma = 90^{\circ}$, and considered different inclinations, $i = 30^{\circ}, 60^{\circ},$ and $80^{\circ}$, and various elevations above the disk plane, $\eta = 0, \pm 20^{\circ}$. We note that the cases $\eta \neq 0$ correspond to conical surfaces above or below the disk plane, but very similar contours can be obtained by fixing the height to $h = 3 \tan \eta$ for all radii. All of the parameters are listed in Table A.1. The azimuthal angle in the disk grows from the direction $(X, -Y)$. For the case of rotation in the equatorial plane ($\eta = 0$) the contours of $\chi^{GR} + \chi^{SR}$ are essentially the same as the ones shown in Fig. 4 of Loktev et al. (2022). We note the presence of a critical point at $i = 30^{\circ}$ near coordinates $(-3, 2)$, where all lines of constant PA rotation intersect. This point corresponds to the place in the disk that an observer sees face-on as a result of the light-bending and SR effects; hence, PD=0 at this point, leading to a degeneracy of PA.

Movement around the origin in Fig. 2 corresponds to a rotation of a point (infinitely small spot) in a disk at a certain radius. The radius containing the critical point separates the region where the rotation of the PA crosses the $\pm 90^{\circ}$ cycle point from

the region where the PA does not make a cycle. These two regions correspond to either one or two loops made by Stokes parameters per one orbital cycle of a spot, and the existence of such critical point has previously been noticed in a number of works (e.g., Dovčiak et al. 2008; GRAVITY Collaboration 2018, 2020b; Gelles et al. 2021; Vos et al. 2022; Loktev et al. 2022; Vincent et al. 2024).

The location of a critical point has been discussed in detail in Vincent et al. (2024) under the assumption of Minkowski space. However, its appearance is not entirely caused by aberration effects, and there are essential differences with respect to the contour morphology in the Schwarzschild metric. First, the SR rotation itself is different in the Minkowski and Schwarzschild metric (compare Eqs. (39) and (40)), as the $\hat{k}_0$ and $\hat{\delta}$ directions are not aligned. Inclusion of the light bending within the calculation of the PA rotation caused by aberration mainly leads to a strong up-down contour asymmetry (around the horizontal line, $y = 0$, in Fig. 2) and a shift in the position of the critical point from the azimuth, $\varphi = 270^{\circ}$, for flat space toward $\varphi \approx 230^{\circ}$ for the Schwarzschild space.

Interestingly, pure aberration effects in the Schwarzschild metric can lead to a formation of two critical points at different radii (Fig. 6 in Loktev et al. 2022); however, the additional rotation due to light-bending effects alters their location and may lead to a shift of one of them below ISCO. This effect is seen in Fig. 2: only one critical point is observed in the panels (a), (d), and (g), corresponding to $i = 30^{\circ}$. At higher inclinations, the critical point disappears because the joint action of aberration and light-bending effects are not strong enough to direct the photons escaping vertically from the disk into the observer's direction, and hence the PA rotation is always smaller than one complete cycle.

The radius and azimuthal coordinate of the critical point, along with the contour morphology, is also altered if the matter is located above or below the equatorial plane (middle and lower panels in Fig. 2). Since in our assumptions we considered the same speed of matter at different elevations, the effect purely scales with the magnitude of light bending.

### 3.2. Effects of magnetic field orientation

Additional rotation of PA can be gained because of the varying orientation of magnetic field with the azimuth; for example, in the cases of purely toroidal and radial magnetic fields. In Fig. 3, we show the results of calculations of $\chi^B$, as a function of the disk azimuth, $\varphi$, for these two cases and three different inclinations of the observer: $i = 30^{\circ}, 60^{\circ},$ and $80^{\circ}$. Circular counterclockwise ($\sigma = 90^{\circ}$) rotation of matter at the equatorial plane of the disk ($\eta = 0$) is assumed. All of the parameters are listed in Table A.1. Different lines correspond to matter located at different equatorial radii: $r_e = 3$ (solid black), 5 (dashed red), 15 (dash-dotted green) and 50 (dotted blue).

For radii of $r_e \ge 50$, the angle makes two full loops per azimuthal cycle for any inclination. In the simplest case of the disk at a relatively low inclination, $i = 30^{\circ}$, the rotation, $\chi^B$, increases with the azimuth nearly linearly. This replicates the behavior for small inclinations and large distances to the compact object described in Eq. (47).

For radii of $r_e \le 5$, there are noticeable deviations from the behavior at large radii thanks to the pronounced effects of relativistic aberration, modified by the light bending. The difference is particularly great at an azimuth of $\varphi \sim 270^{\circ}$; that is, for the approaching side of the disk. For $i = 30^{\circ}$, the changes are severe, as rotation reduces to one loop over the azimuthal cycle,
<figure><img src="image_2.png" /></figure>

**Fig. 2.** Contours of PA rotation ($\chi^{\text{SR}} + \chi^{\text{GR}}$) for different $i = 30^{\circ}, 60^{\circ}$, and $80^{\circ}$ and $\eta = 0, \pm 20^{\circ}$ in the disk coordinates $(r, \varphi)$, where $r$ is measured in units of Schwarzschild radii. Azimuth increases in the counterclockwise direction, as is shown by the arrow in panel (a). Numbers on the contours correspond to the rotation angle (in degrees). The thick red line without a mark corresponds to $\pm 90^{\circ}$ rotation. Matter is assumed to have a purely circular Keplerian velocity in the counterclockwise direction ($\sigma = 90^{\circ}$). The dashed line indicates the location of the ISCO. Contours corresponding to $\eta = 0$ (upper row) match those shown in Loktev et al. (2022). The effects of matter elevation above (or below) the orbital plane are clearly visible when comparing to the case of in-plane rotation. If matter is located in front of the BH (middle row), for the same radius from the BH the GR and SR effects appear to be less pronounced.

and hence the whole topology of rotation changes. The behavior of the PA with the azimuth is now different for the radial and toroidal fields, making it potentially possible to distinguish these fields via orbital (azimuthal) phase-resolved polarimetric data.

For both the toroidal and radial fields, the total PA rotation ($\chi^{\text{tot}}$) is dominated by the rotation caused by the different magnetic field orientation ($\chi^B$), rather than the rotation caused by relativistic effects ($\chi^{\text{GR}} + \chi^{\text{SR}}$), and the contours of constant rotation do not contain the critical point (examples of contours for the case of the toroidal field and different $i$ are shown in Fig. 4). However, the lines of PA rotation have nearly the same morphology for pure radial and toroidal fields, albeit shifted in azimuth. Hence, their linear combination may lead to a substantial reduction of the PA rotation with the azimuth, essentially mimicking the vertical field, and may replicate the behavior near the critical point (Fig. 2). We find that the critical point appears, for example, for the combination of radial and toroidal field components ($\phi' = 45^{\circ}$) at high inclinations ($i = 80^{\circ}$) and matter below the equatorial plane ($\eta = -20^{\circ}$).

### 3.3. Effects of non-equatorial motion

The GR and SR effects on rotation of polarization plane have been studied, using the aforementioned formalism, in the case of razor-thin matter in an equatorial plane of BHs (Loktev et al. 2022, 2024). Emitting matter may, however, be elevated above the equatorial plane in the case of geometrically thick disks or if emission arises at the interface between the inflow and outflow (edge-brightened jet found, e.g. in high-resolution images of Cen A, Janssen et al. 2021). Our formalism enables a direct comparison between the polarization rotation for the matter rotating in the plane containing the BH itself (equatorial plane, $\eta = 0$) and that of the matter elevated by the height, $h$ (set by $\eta \neq 0$) above or below this plane.

In Fig. 5, we show the resulting observed PA ($\chi$) for the specific case of $\eta = 20^{\circ}$. The topology of the lines is the same as for the case $\eta = 0$, but the absolute values differ. In the cases of the equatorial field (two upper rows of panels), the PA rotation is dominated by the $\chi^B$ term and GR and SR effects are less pronounced, as lines corresponding to different radii nearly over-
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Rotation of PA due to magnetic field orientation, as a function of azimuth (in disk coordinates) for different equatorial radii (measured in units of Schwarzschild radii): $r_e = 3$ (solid black), 5 (dashed red), 15 (dash-dotted green), and 50 (dotted blue). The matter is assumed to undergo circular, counterclockwise rotation ($\sigma = 90^\circ$) with velocity following Luminet's law. The upper panels correspond to a purely radial comoving magnetic field ($\phi' = 90^\circ$) and the lower panels correspond to a purely toroidal field ($\phi' = 0$). The columns (left to right) correspond to the observer inclination, $i = 30^\circ, 60^\circ$, and $80^\circ$.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 4.** Contours of constant total rotation angle ($\chi^{\text{tot}}$) for the case of a toroidal magnetic field ($\sigma = 90^\circ, \phi' = 0$) and the equatorial plane matter movement, in disk coordinates. We note the disappearance of the critical point, compared to the case of a vertical field, as the rotation is dominated by the changing magnetic field orientation ($\chi^B$). The observed PA, $\chi$, was obtained by adding $\chi_0 = 90^\circ$.</figcaption></figure>

lap, especially for smaller inclinations. Figure 2 (d-f) shows this effect in more detail. The contours of the same $\chi^{\text{GR}} + \chi^{\text{SR}}$ are generally closer to the event horizon for $\eta = 20^\circ$, compared to $\eta = 0$. This means that for the same distance from the BH, the polarization plane rotation is more severe for the case of equatorial rotation than for the matter located between this plane and the observer. The main reason for this is that for higher $\eta$, a smaller part of photon path lies close to the strong gravitational field of the BH; hence, bending effects are less pronounced. Likewise, if the matter is located below the equatorial plane – that is, behind the BH – relativistic effects become most pronounced and contours of higher PA rotation appear (Fig. 2g-i).

To quantify the effects of matter elevation at different radii, in Fig. 6 we show the difference between the total rotation ($\chi^{\text{tot}}$) for the case of $\eta = 20^\circ$ with respect to the case $\eta = 0$ for various inclinations and magnetic field configurations, observer inclinations, and radii to the BH. For large distances, $r$, the effect of elevated matter is minor, as was expected. However, for distances

of $r \lesssim 10$, the difference can reach $30^\circ$, becoming comparable to the total rotation, and hence cannot be neglected. The difference can reach $90^\circ$ for certain parameters and azimuths. It is generally asymmetric with respect to $\varphi$ and has a complex dependence on the magnetic field orientation, making it difficult to a priori predict the effect of elevated matter on the PA rotation. However, it is generally clear that the effect is highest at azimuths of $180^\circ \le \varphi \le 270^\circ$, where the GR and SR effects are most pronounced. Since these effects strongly depend on the distance to the compact object and matter velocity, any changes to the distance caused by nonzero elevation result in substantial changes to the total PA rotation.

### 3.4. Effects of noncircular motion

Accretion flows around BHs have a nonzero radial velocity component – the accretion speed. In our formalism, this is parameterized by the angle, $\sigma$, which may potentially affect the result-
<figure><img src="image_5.png" /><figcaption>**Fig. 5.** PA seen by the observer ($\chi$) obtained for the angle $\eta = 20^\circ$, for different radii, observer inclinations, and magnetic fields, as a function of the azimuth in the disk. The upper panels show the case of the radial field, the middle panels correspond to the purely toroidal field, and the lower panels correspond to the vertical field. The columns correspond, from left to right, to an increasing observer inclination: $i = 30^\circ, 60^\circ$, and $80^\circ$. The matter is assumed to rotate counterclockwise ($\sigma = 90^\circ$), with a velocity following Luminet's law.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 6.** Difference between the PA obtained for the angle $\eta = 20^\circ$ and $\eta = 0$ ($\chi^{\text{tot}}(\eta = 20^\circ) - \chi^{\text{tot}}(\eta = 0)$), for different radii, observer inclinations, and magnetic fields, as a function of the azimuth in the disk coordinates. The upper panels show the case of the radial field, the middle panels correspond to the purely toroidal field, and the lower panels correspond to the vertical field. The columns, from left to right, correspond to an increasing observer inclination: $i = 30^\circ, 60^\circ$, and $80^\circ$. The matter is assumed to rotate counterclockwise ($\sigma = 90^\circ$), with a velocity following Luminet's law.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 7.** Effects of noncircular motion ($\sigma$) on the observed PA ($\chi^{\text{rot}}$) for different inclinations (left to right: $i = 30^{\circ}$, $60^{\circ}$, and $80^{\circ}$) and magnetic fields (upper panels: radial field, middle panels: toroidal field, lower panels: vertical field), as a function of the azimuth in the disk coordinates. Different directions of motion correspond to $\sigma = 90^{\circ}$ (solid black lines), $110^{\circ}$ (dashed red), $150^{\circ}$ (dash-dotted green), and $180^{\circ}$ (dotted blue). The emitting matter is assumed to be located in the equatorial plane ($\eta = 0$) at a distance of $r = 5$, with a dimensionless velocity obeying Luminet's law.</figcaption></figure>

ing rotation of polarization, as the highest aberration may be achieved at a different angle. The trivial case of changing direction is the clockwise rotation, $\sigma = 270^{\circ}$, which can be visualized by simply changing the azimuthal angle, $\varphi_{\text{cl-w}} = 360^{\circ} - \varphi$. Another important example includes spherical (Bondi-Hoyle) accretion, $\sigma = 180^{\circ}$. We consider the effects of noncircular motions of matter for different magnetic field geometries and varying inclinations in Fig. 7.

In equatorial field cases (panels a–f), the field direction is always fixed either parallel or orthogonal to the matter motion, and hence when the velocity direction changes, the field direction follows, introducing constant shifts to the lines. Apart from the constant shifts, there are minor deviations of the PA dependence on the azimuthal angle owing to the different combined SR and GR effects. Again, these are mostly pronounced for the locations of matter behind the BH ($\varphi \sim 180^{\circ}$).

In the case of the vertical field, the changes are generally likewise minor (Fig. 7h,i). The highest difference with respect to the default counterclockwise rotation occurs for $i = 30^{\circ}$ (Fig. 7g), close to the location of the critical point. By changing the velocity vector, it is possible to shift this point, changing the character of the PA rotation from the limited range (solid black line in Fig. 7g) to two complete loops per orbit (dashed

red, dash-dotted green, and dotted blue). Hence, the $\sigma$ parameter may be degenerate with inclination.

# 4. Applications

In this section, we compare the topology of static and dynamic signatures of BHs. The developed formalism can be applied to reproduce the horizon-scale images and polarimetric properties of M87* and Sgr A* (Event Horizon Telescope Collaboration 2019, 2021, 2022, 2024a), as well as dynamic polarimetric signatures of Sgr A* (Wielgus et al. 2022; GRAVITY Collaboration 2023). We start with the expected polarimetric signatures of the orbiting compact and extended spots. We then discuss polarimetric signatures of the narrow disk ring. After that, we show examples of computed images of the geometrically thick accretion disk.

## 4.1. Hot spot rotation

The quiescent level of Sgr A* accretion is known to be interrupted by the bright flares, when the observed radio (submillimeter) and NIR fluxes can suddenly increase by 1–2 orders
<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Contours of constant GR and SR rotation angles ($\chi^{\text{GR}} + \chi^{\text{SR}}$), in disk coordinates. The position of the critical point is shown as a function of inclination for rotation caused solely by GR and SR effects (equivalent to $B = B_z$). The disk rotates counterclockwise ($\sigma = 90^\circ$) in the equatorial plane, but at inclinations of $i = 165^\circ$ (a), $155^\circ$ (b), and $140^\circ$ (c) it is seen as rotating in the clockwise direction.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Flux Stokes parameters, $F_Q$ and $F_U$, obtained using $\eta = 0$, for different radii, $B$ and $i$. Keplerian velocity and $\sigma = 90^\circ$ are assumed (note that for the considered inclinations, the observed matter movement is clockwise on the sky). Increasing saturation of color corresponds to progression of the spot along the azimuth. The time when the line intersects zero PD corresponds to the passage through the critical point. Note different scales in figures.</figcaption></figure>

of magnitude (e.g. Baganoff et al. 2001; Genzel et al. 2003). The characteristic timescales of the flares are roughly consistent with Keplerian time close to the event horizon, and the hot spot rotation scenario was proposed to explain the flares (Broderick & Loeb 2006; Trippe et al. 2007). Alternatively, a pattern movement could be responsible for the production of similar characteristics (Matsumoto et al. 2020; Aimar et al. 2023). Recently, observed movement of the brightness centroid

and loops of polarized fluxes have given additional support to these scenarios (GRAVITY Collaboration 2018; Wielgus et al. 2022; GRAVITY Collaboration 2023). Perhaps the most intriguing property found in those data is the one observed polarimetric loop in Stokes parameters per one orbital period, with a zero-polarization point (where absolute Stokes parameters $Q = U = 0$) being located within this loop (in the NIR, where interstellar polarization is expected to be low). Moreover, the PA was found
<figure><img src="image_10.png" /><figcaption>**Fig. 10.** Stokes parameters, *U* and *Q*, for matter flying toward the observer ($\eta_{\text{max}} = 45^{\circ}$), in the direction orthogonal to rotation plane, at different launching azimuthal angles: $\varphi_0 = 90^{\circ}$ (upper panels), $180^{\circ}$ (middle panels), and $270^{\circ}$ (lower panels). Note different scales in panels.</figcaption></figure>

to change at a nearly constant rate as a function of time, presum-
ably corresponding to the azimuthal angle in the accretion disk
(GRAVITY Collaboration 2023).

One complete rotation of PA per orbit is not expected to
appear for a large fraction of the parameter space and is tightly
linked to the location of the critical point (where the observed
polarization turns zero). If the orbit does not intersect the loca-
tion of a critical point ($r_{cr}$, $\varphi_{cr}$), either two complete loops of
Stokes parameters (for orbits within $r_{cr}$) or one loop that does not
enclose the origin $Q = U = 0$ (orbits beyond the radius of the
zero-polarization point) are generally expected (see Sect. 3.1).
The appearance of the critical point is also tightly related to the
magnetic field direction (Sect. 3.2): if the spot is moving in the
equatorial plane, the critical point is only observed in the case
of the vertical field $B = B_z$. In order to reproduce the observed
topology, the radius and inclination of the hot spot should clearly
coincide with the location of the critical point, so that the sec-
ond loop, which occupies a small fraction of the orbital motion,
is small or completely absent, and the rotation around the sec-
ondary loop is fast enough to not affect the constant rate of PA
changes. In other words, for the data that are limited in time and
angular resolution, it should not be feasible to resolve the move-
ment along the second polarimetric loop.

The limited range of radii and inclinations for which the
critical point exists, as well as the strong dependence on the
magnetic field topology, then translates into highly constrain-

ing estimates on these parameters that come from polarimetry
(e.g., GRAVITY Collaboration 2018, 2023). In Fig. 8, we show
the contours of PA rotation for the case of the vertical field,
with inclinations over 90° that are believed to be relevant to
the accreting matter around Sgr A*. For the inclinations outside
140° < *i* < 165°, the location of the critical point is either out-
side of the ~10*R<sub>S</sub*, corresponding to scales that the instruments
are sensitive to, or within the ISCO.

In Fig. 9, we show the patterns traced by the spot rotating
in the equatorial plane. Increasing time (or azimuth) is shown as
saturating color: the pale part of the curves corresponds to the
beginning of the rotation from $\varphi = 0$. To produce this figure,
we used the latter representation of (Stokes) fluxes in Eq. (35),
and kept $dS' = \text{const}$ throughout the rotation. While Eq. (35)
was derived assuming the spot emits through the upper surface
of the disk, the same expression is also valid for the synchrotron
emission from a spherical, optically thin spot of radius $h_0$, as in
this case the factors, $\cos \zeta$, entering the area, $dS'$, and the photon
path, $l'_{\text{ph}}$, are omitted.

The loop topology becomes altered if the matter gradu-
ally rises above or below the equatorial plane. This scenario is
expected; for example, if the spots correspond to ejecting plas-
moids (Aimar et al. 2023). We illustrate this case by assum-
ing a gradual increase in matter elevation above the orbital
plane, keeping the equatorial radius of matter the same (helical
<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Same as in Fig. 10, but for matter flying away from the observer ($\eta_{max} = -45^{\circ}$) in the direction orthogonal to the rotation plane. Note different scales in panels.</figcaption></figure>

motion). The results of calculations for the case of the vertical magnetic field and different inclinations are shown in Figs. 10 and 11. The topology of the loops strongly depends on the azimuthal angle, $\varphi_0$, at which the spot starts rising, as is evident from the comparison of the top, middle, and bottom panels.

Next, we considered the behavior of relative Stokes parameters in the regime of strong gravity and fast motions. The PD is Lorentz-invariant, and hence compact spots located in any part of the disk are expected to have the PD that corresponds to the PD of intrinsic emission from the corresponding disk segment viewed at an angle, $\zeta'$. Angular dependence of the intrinsic PD, for example in the case of the optically thick disk atmospheres, leads to a modulation of the PD with orbital phase, since the viewing angle, $\zeta'$, for a given observer's inclination and radius, changes with the azimuth (e.g., Fig. 7 of Loktev et al. 2022). However, in the case of the synchrotron emission, the PD does not depend on the viewing angle, but only on the index of the power-law electrons (Ginzburg & Syrovatskii 1965). Hence, in the simple scenario, we do not expect changes in the PD with the orbital phase (azimuth), which translates into the topology of two complete cycles or a part of the cycle pattern in the *qu* plane. In the lower panels of Fig. 12, we show the direct calculation of the relative Stokes parameters in the case of an infinitely small spot at different radii and inclinations (black, red, and green

curves) that confirm the two-cycle or part-cycle pattern. A more complex pattern may be obtained if the angular dependence of the PD is caused by the angular-dependent (with respect to the magnetic field orientation) electron spectral index, *p*.

## 4.2. Polarimetry of an extended spot and disk ring

An important property of the observed signal is the considerably lower PD detected in the flares, ~25%, as compared to the theoretical expectations, ~70%. The reduction in the PD may be caused by the beam depolarization (Narayan et al. 2021), or if the spot is extended in the azimuthal, vertical, and/or radial direction (see examples of simulations in Yfantis et al. 2024). In Fig. 12, we compare the absolute ($F_Q, F_U$) and relative ($q, u$) Stokes parameters for the infinitely compact spot (black, red, and green curves) to those that are averaged in the azimuthal direction (infinitely thin disk segment), assuming the spot size of a quarter circle (blue, violet, and yellow curves). The spot is assumed to move and each point in the plot corresponds to the azimuth of the middle of the spot. The depolarization essentially occurs because of the averaging over Stokes vectors with different angles, but a large spot size (comparable to the orbital scale) is needed. The single loop and distorted shapes of the *qu* parameters can be recovered.
<figure><img src="image_12.png" /><figcaption>**Fig. 12.** Stokes parameters of the infinitely compact spots and disk segments. Averaging was made over $\Delta\varphi = 45^\circ$ (i.e., a quarter of the circle). The upper and lower panels correspond to the absolute and relative Stokes parameters, respectively. Black, red, and green lines denote calculations of the infinitely small spot, while blue, violet, and yellow lines correspond to the parameters averaged over the segment. Note that several lines overlap in the lower panels.</figcaption></figure>

It is noteworthy that the behavior of NIR polarization during the flares does not always follow the one-loop pattern over the ~70 min timescale (thought to correspond to the orbital period of the spots). Namely, the pattern detected in flares from August 18, 2019, and May 19, 2022 (Fig. 1 of GRAVITY Collaboration 2023), resembles a line rather than a loop. Furthermore, simultaneous high-precision astrometric measurements available for the flare on May 19, 2022, do not resemble a circle, expected for the circular motion around the compact object, but likewise are more consistent with a line. Similar behavior is observed for the astrometric flare from July 28, 2018. The patterns detected in polarization are consistent with those shown in the right columns of Figs. 10 and 11; that is, hot spots at high inclinations and flying away from the orbital plane. Likewise, the path of hot spots at inclinations closer to 90° is more squeezed, compared to the nearly circular one at low (~0° or 180°) inclinations (e.g., images of accreting BHs in Fig. 9 of Loktev et al. 2022). We suggest that the flares showing elongated paths (or elongated loops) in astrometry and/or polarimetry can correspond to hot spots passing the BH at high inclinations.

If the extent of the hot spot becomes large, the PD drops, but it does not generally become zero. In the case of the disk ring, the net polarization arises owing to the broken symmetry as a result of the combined GR and SR effects. The net polarization of ~5% was observed in Sgr A* (Event Horizon Telescope Collaboration 2024a). In Fig. 13, we show the net polarization (Stokes parameters) of the narrow disk ring at various radii, as a function of inclination. All points intersect at the origin, $(q, u) = (0, 0)$, for $i = 0$ and the polarization generally increases with increasing $i$. The closer the disk ring is located to the BH, the higher the net polarization is, and the more the PA deviates from the direction along the disk normal. Hence, the $u$ parameter of the net polarization can be used to obtain information about the proximity of matter to the compact object.

The polarimetric signatures of the disk ring depend on the inclination; hence, they are expected to change if the ring under-

<figure><img src="image_13.png" /><figcaption>**Fig. 13.** Relative Stokes parameters (in fractional units) of the disk ring (in the equatorial plane, $\eta = 0$, and rotating in the counterclockwise direction, $\sigma = 90^\circ$) with vertical field $B = B_z$, as a function of inclination for different ring radii. Zero polarization corresponds to the case of $i = 0$ and increases, along the line, to $i = 90^\circ$. Inclinations of $i = 10^\circ, 20^\circ$, and $30^\circ...80^\circ$ are marked with crosses. For inclinations of $i > 90^\circ$, the curve is reflected around the $u = 0$ axis.</figcaption></figure>

goes Lense-Thirring precession. The precessing disk rings scenario has been considered in the context of quasi-periodic eruptions (Raj & Nixon 2021; Arcodia et al. 2021). In the course of the Lense-Thirring precession, the inclination to the disk changes in a periodic way within limits set by the observer's orientation relative to the BH spin and the misalignment angle between the BH spin and the plane of matter inflow. Intrinsic dependence of the polarimetric signatures on the inclination in Fig. 13 can then be converted to the observed profiles simply by
<figure><img src="image_14.png" /><figcaption>**Fig. 14.** PD as a function of inclination, scaled to its maximal value, $P_s$, in the case of toroidal magnetic field ($\phi' = 90^\circ$), equatorial counterclockwise rotation ($\eta = 0, \sigma = 90^\circ$). The case $r_e = 10^4$ shows the dependence in the absence of relativistic effects and can be approximated well by $\sin^\pi i$.</figcaption></figure>

computing the inclination of the disk at each precession phase (Veledina et al. 2013).

In Fig. 14, we show the dependence of the PD on the inclination angle in the case of the toroidal magnetic field, $\phi' = 90^\circ$, and different radii. The case of $r = 10^4$ might be relevant to the outer parts of the jet, as the ratio of the poloidal and toroidal components of the field decreases linearly with distance in the standard Blandford & Königl (1979) jet model. At small viewing angles, the PD is small because positive contributions to polarization coming from azimuthal angles $\varphi$ and $\varphi + 180^\circ$ are (nearly) canceled by the negative contributions coming from the azimuthal angles $\varphi \pm 90^\circ$, for any $\varphi$. Light bending and aberration alter the PD at high inclinations, where otherwise (for slow motions and no light bending) the contributions from different parts of the disk have nearly the same PA; hence, the PD reaches a maximal value, $P_s$. For small distances from the compact object, we observe the effect of depolarization, caused by the different viewing angles of different parts of the disk; it is mostly pronounced at large inclinations. A good approximation for the case of a high radius is given by the $\sin^\pi i$ function (solid green line in Fig. 14).

For toroidal fields, the net polarization is aligned with the disk (or jet) axis for matter located at large distances. Close to the BH, the net PA deviates from this direction by up to $12^\circ$, achieved at $r_e = 3$. The PA deviation weakly depends on the inclination.

### 4.3. Image of the geometrically thick disk

Polarimetric images of two accreting supermassive BH environments recently obtained by the Event Horizon Telescope (Event Horizon Telescope Collaboration 2021, 2024a) show a consistent picture of spiral polarization pattern that is dominated by the azimuthally symmetric structure. The images correspond to the quiescent level of Sgr A* and M 87* emission and both show the rotationally symmetric pattern of PA rotation.

In Fig. 15, we plot the intensity (Eq. (35)) in the upper panel and the PD in the lower panel. The disk is assumed to be geometrically thick,

<figure><img src="image_15.png" /><figcaption>**Fig. 15.** Color-coded image of the thick accretion disk around a BH. In the upper panel, the color refers to the monochromatic intensity; in the lower panel, it refers to PD. The drop in PD close to the critical point can be seen in the upper left part of the image; otherwise, some depolarization is caused by the nonzero disk thickness.</figcaption></figure>

with $\eta = 30^\circ$ at the ISCO, which translates into a constant height of $h = 1.73R_s$. We split the disk into 11 horizontal layers (one layer in the disk midplane and five layers below and above it) and assumed the elevation of matter above the midplane to be constant within each layer. We verified that this resolution is sufficient for convergence purposes. The observer inclination, $i = 30^\circ$, and vertical magnetic field, $B = B_z$, were assumed. We obtained the characteristic image of an asymmetric ring-like structure, similar to the observed ones in M 87* and Sgr A*, as well as to the ones obtained in modeling of earlier works (Narayan et al. 2021; Vincent et al. 2024; Event Horizon Telescope Collaboration 2024b). However, the area of the region with high intensity is larger in the case of the thick disk, compared to the flat disk, and the shape of the BH shadow is squashed in the case of the thick disk. Both effects are caused by the presence of matter above the equatorial plane, at azimuthal angles close to zero (i.e., between the observer and the BH).

The image of the disk is depolarized compared to the results for the narrow, flat ring. The depolarization is caused by the substantial vertical thickness: emission coming from different layers of the disk has different rotation angles of the polarization plane
(see contours of PA rotation for different η in Fig. 2); hence, the net polarization is reduced. The dark area around coordinates (−4, 3) corresponds to the place where the disk is viewed face-on as a result of light-bending and aberration effects. For the vertical field, both intensity and PD decrease to zero here, and nonzero values appear as a result of the different position of this point for different disk layers. For the direct comparison to the data, additional beam depolarization and internal Faraday rotation (that is known to affect the polarization signatures in the submillimeter range, Event Horizon Telescope Collaboration 2021, 2024a; Wielgus et al. 2024) need to be considered. The latter effect can be taken into account by assuming that each layer between the emission point and the observer acts as a Faraday screen.

## 5. Summary

In this study, we have derived explicit analytical expressions for the transformation of polarization signatures for a fast-moving matter in the Schwarzschild metric in the case of the synchrotron emission. The rotation of the PA is caused by the joint effects of general and special relativity, as well as magnetic field orientation, which were considered separately. We considered different cases of moving matter direction, its elevation above and below the equatorial plane and magnetic field orientation, and computed the resulting PAs as a function of the azimuthal coordinate in the disk. We have further discussed potential applications of the developed formalism for the dynamic polarimetric signatures of spots orbiting around the compact object, as well as static images of geometrically thick accretion flows. We find that a number of the observed signatures can be explained by the nonzero thickness of the disk, which was generally omitted in previous analytical results. We highlight an important difference of the calculations presented in this work, compared to the polarization signatures in Minkowski space. The simple form of the derived expressions makes them attractive for the acceleration of the post-processing calculations, which can be done for a direct comparison with observations.

*Acknowledgements.* AV acknowledges support the Academy of Finland grant 355672. Nordita is supported in part by NordForsk. Authors thank V. Loktev, J. Poutanen and M. Wielgus for the useful discussions and comments on the paper.

## References

Aimar, N., Dmytriiev, A., Vincent, F. H., et al. 2023, A&A, 672, A62

Arcodia, R., Merloni, A., Nandra, K., et al. 2021, Nature, 592, 704

Baganoff, F. K., Bautz, M. W., Brandt, W. N., et al. 2001, Nature, 413, 45

Beloborodov, A. M. 2002, ApJ, 566, L85

Blandford, R. D., & Königl, A. 1979, ApJ, 232, 34

Bogdanov, S., Lamb, F. K., Mahmoodifar, S., et al. 2019, ApJ, 887, L26

Broderick, A. E., & Loeb, A. 2006, MNRAS, 367, 905

Chang, D. O., Johnson, M. D., Tiede, P., & Palumbo, D. C. M. 2024, ApJ, 974, 143

Chen, K., & Shaham, J. 1989, ApJ, 339, 279

Connors, P. A., & Stark, R. F. 1977, Nature, 269, 128

Connors, P. A., Piran, T., & Stark, R. F. 1980, ApJ, 235, 224

Dexter, J., & Agol, E. 2009, ApJ, 696, 1616

Do, T., Ghez, A. M., Morris, M. R., et al. 2009, ApJ, 691, 1021

Do, T., Witzel, G., Gautam, A. K., et al. 2019, ApJ, 882, L27

Dodds-Eden, K., Porquet, D., Trap, G., et al. 2009, ApJ, 698, 676

Dovčiak, M., Muleri, F., Goosmann, R. W., Karas, V., & Matt, G. 2008, MNRAS, 391, 32

Eckart, A., Baganoff, F. K., Schödel, R., et al. 2006a, A&A, 450, 535

Eckart, A., Schödel, R., Meyer, L., et al. 2006b, A&A, 455, 1

Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2019, ApJ, 875, L1

Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2021, ApJ, 910, L12

Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2022, ApJ, 930, L12

Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2024a, ApJ, 964, L25

Event Horizon Telescope Collaboration (Akiyama, K., et al.) 2024b, ApJ, 964, L26

Gelles, Z., Himwich, E., Johnson, M. D., & Palumbo, D. C. M. 2021, Phys. Rev. D, 104, 044060

Genzel, R., Schödel, R., Ott, T., et al. 2003, Nature, 425, 934

Gillessen, S., Eisenhauer, F., Quataert, E., et al. 2006, ApJ, 640, L163

Ginzburg, V. L., & Syrovatskii, S. I. 1965, ARA&A, 3, 297

GRAVITY Collaboration (Abuter, R., et al.) 2018, A&A, 618, L10

GRAVITY Collaboration (Abuter, R., et al.) 2020a, A&A, 638, A2

GRAVITY Collaboration (Jiménez-Rosales, A., et al.) 2020b, A&A, 643, A56

GRAVITY Collaboration (Abuter, R., et al.) 2023, A&A, 677, L10

Janssen, M., Falcke, H., Kadler, M., et al. 2021, Nat. Astron., 5, 1017

Loktev, V., Salmi, T., Nättilä, J., & Poutanen, J. 2020, A&A, 643, A84

Loktev, V., Veledina, A., & Poutanen, J. 2022, A&A, 660, A25

Loktev, V., Veledina, A., Poutanen, J., Nättilä, J., & Suleimanov, V. F. 2024, A&A, 685, A84

Luminet, J. P. 1979, A&A, 75, 228

Matsumoto, T., Chan, C.-H., & Piran, T. 2020, MNRAS, 497, 2385

Narayan, R., Palumbo, D. C. M., Johnson, M. D., et al. 2021, ApJ, 912, 35

Nättilä, J., & Pihajoki, P. 2018, A&A, 615, A50

Pechenick, K. R., Ftaclas, C., & Cohen, J. M. 1983, ApJ, 274, 846

Pineault, S. 1977, MNRAS, 179, 691

Pineault, S., & Roeder, R. C. 1977, ApJ, 212, 541

Ponti, G., George, E., Scaringi, S., et al. 2017, MNRAS, 468, 2447

Poutanen, J. 2020, A&A, 640, A24

Poutanen, J., & Beloborodov, A. M. 2006, MNRAS, 373, 836

Raj, A., & Nixon, C. J. 2021, ApJ, 909, 82

Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: Wiley-Interscience)

Stark, R. F., & Connors, P. A. 1977, Nature, 266, 429

Suleimanov, V. F., Poutanen, J., & Werner, K. 2020, A&A, 639, A33

Trippe, S., Paumard, T., Ott, T., et al. 2007, MNRAS, 375, 764

Veledina, A., Poutanen, J., & Ingram, A. 2013, ApJ, 778, 165

Vincent, F. H., Paumard, T., Perrin, G., et al. 2014, MNRAS, 441, 3477

Vincent, F. H., Wielgus, M., Aimar, N., Paumard, T., & Perrin, G. 2024, A&A, 684, A194

Vos, J., Mościbrodzka, M. A., & Wielgus, M. 2022, A&A, 668, A185

Walker, M., & Penrose, R. 1970, Commun. Math. Phys., 18, 265

Wielgus, M., Moscibrodzka, M., Vos, J., et al. 2022, A&A, 665, L6

Wielgus, M., Issaoun, S., Martí-Vidal, I., et al. 2024, A&A, 682, A97

Yfantis, A. I., Wielgus, M., & Mościbrodzka, M. 2024, A&A, 691, A327

Yusef-Zadeh, F., Roberts, D., Wardle, M., Heinke, C. O., & Bower, G. C. 2006, ApJ, 650, 189
Appendix A: Rotation of polarization angle

Light bending and relativistic aberration alter the direction of the photon path and hence change the direction of its electric field orientation, as seen by the observer. To determine the rotation of PA, we follow the method described in Loktev et al. (2022). For an arbitrary vector $l$, not collinear with the photon propagation direction $\hat{k}$, we introduce the polarization basis, formed by two vectors orthogonal to each other and to the vector $\hat{k}$ as

$$ \mathcal{B}(\hat{k}, l) = \{\hat{e}_1^{k,l}, \hat{e}_2^{k,l}\} = \left\{ \frac{l - (\hat{k} \cdot l)\hat{k}}{|\hat{k} \times l|}, \frac{\hat{k} \times l}{|\hat{k} \times l|} \right\}. \quad (A.1) $$

The PA is given by the angle between the polarization plane of the photon and $\hat{e}_1^{k,l}$, in the counterclockwise direction, as viewed by the observer intercepting the photon. For any other vector $l_2$, the new basis is $\mathcal{B}(\hat{k}, l_2)$, and the PA is given by the angle with respect to the new axis $\hat{e}_1^{k,l_2}$. The new basis is built around the same photon direction, hence the PAs are different by the angle of rotation between the bases,

$$ \chi(\hat{k}, l_2 \rightarrow l) = \mathcal{B}(\hat{k}'_0, \hat{B}') \rightarrow \mathcal{B}(\hat{k}'_0, \hat{n}), \quad (A.2) $$

that can be computed as

$$ \tan \chi(\hat{k}, l_2 \rightarrow l) = \frac{\hat{k} \cdot (l \times l_2)}{l \cdot l_2 - (\hat{k} \cdot l)(\hat{k} \cdot l_2)}. \quad (A.3) $$

For the synchrotron radiation of matter close to the BH, the photon is initially emitted along the direction $\hat{k}'_0$ with polarization orthogonal to the magnetic field lines $\hat{B}'$, the latter defines the basis $\mathcal{B}(\hat{k}'_0, \hat{B}')$. Next, we shift to the description of the photon polarization plane orientation with respect to the normal to the plane of accreting matter $\hat{n}$. Further, the relativistic aberration lead to the transition from $\hat{k}'_0$ to $\hat{k}_0$. Finally, the light bending changes the photon vector $\hat{k}_0$ to $\hat{o}$. Hence, we make the following steps between the bases:

$$ \mathcal{B}(\hat{k}'_0, \hat{B}') \rightarrow \mathcal{B}(\hat{k}'_0, \hat{n}) \rightarrow \mathcal{B}(\hat{k}_0, \hat{n}) \rightarrow \mathcal{B}(\hat{o}, \hat{n}). \quad (A.4) $$

Each step is associated with the rotation of PA by the angle between the bases and is associated with different effects: magnetic field orientation $\chi^B$, aberration $\chi^{SR}$ and light bending $\chi^{GR}$. This can be explicitly written as

$$ \chi = \chi_0 + \chi^{\text{tot}} = \chi_0 + \chi^B + \chi^{\text{SR}} + \chi^{\text{GR}}, \quad (A.5) $$

where $\chi_0 = 90^\circ$.

Rotation associated with the magnetic field orientation is computed as the shift $\mathcal{B}(\hat{k}'_0, \hat{B}) \rightarrow \mathcal{B}(\hat{k}'_0, \hat{n})$. Using Eq. A.3, $\chi^B$ is computed as

$$ \chi^B = \arctan \left( \frac{\hat{k}'_0 \cdot (\hat{n} \times \hat{B}')}{\hat{n} \cdot \hat{B}' - (\hat{k}'_0 \cdot \hat{n})(\hat{k}'_0 \cdot \hat{B}')} \right). \quad (A.6) $$

Explicit analytical expression for the considered geometry and matter motions is given in Eq. (45).

The second and the third transformations have been described in Loktev et al. (2022) and we briefly mention them below for the sake of completeness. The transformation $\mathcal{B}(\hat{k}'_0, \hat{n}) \rightarrow \mathcal{B}(\hat{k}_0, \hat{n})$ is done via three steps:

$$ \mathcal{B}(\hat{k}'_0, \hat{n}) \rightarrow \mathcal{B}(\hat{k}'_0, \hat{v}) \rightarrow \mathcal{B}(\hat{k}_0, \hat{v}) \rightarrow \mathcal{B}(\hat{k}_0, \hat{n}), \quad (A.7) $$

where the second step introduces zero rotation, as the vectors $\hat{k}'_0$, $\hat{v}$ and $\hat{k}_0$ lie in the same plane. For the considered case $\hat{v} \perp \hat{n}$, the rotation angle is obtained as

$$ \tan \chi^{\text{SR}} = - \frac{\beta(\hat{k}_0 \cdot (\hat{v} \times \hat{n}))(\hat{k}_0 \cdot \hat{n})}{1 - (\hat{k}_0 \cdot \hat{n})^2 - \beta(\hat{k}_0 \cdot \hat{v})}. \quad (A.8) $$

The transformation $\mathcal{B}(\hat{k}_0, \hat{n}) \rightarrow \mathcal{B}(\hat{o}, \hat{n})$ is done via an additional step involving the radial vector direction $\hat{r}$:

$$ \mathcal{B}(\hat{k}_0, \hat{n}) \rightarrow \mathcal{B}(\hat{k}_0, \hat{r}) \rightarrow \mathcal{B}(\hat{o}, \hat{r}) \rightarrow \mathcal{B}(\hat{o}, \hat{n}), \quad (A.9) $$

where, again, the second step introduces zero rotation, as the vectors $\hat{k}_0$, $\hat{r}$ and $\hat{o}$ belong to the same plane in Schwarzschild metric (photon trajectories lie in plane). The resulting vector expression for the rotation angle is

$$ \tan \chi^{\text{GR}} = \cos \eta \frac{(\hat{o} \cdot \hat{\varphi})\tilde{g}_1 - (\hat{k}_0 \cdot \hat{\varphi})\tilde{g}_2}{\tilde{g}_1\tilde{g}_2 + \cos^2 \eta (\hat{o} \cdot \hat{\varphi})(\hat{k}_0 \cdot \hat{\varphi})}, \quad (A.10) $$

where

$$ \tilde{g}_1 = (\hat{r} \cdot \hat{n}) - (\hat{k}_0 \cdot \hat{n})(\hat{k}_0 \cdot \hat{r}), \quad (A.11) $$

$$ \tilde{g}_2 = (\hat{r} \cdot \hat{n}) - (\hat{o} \cdot \hat{n})(\hat{o} \cdot \hat{r}). \quad (A.12) $$
<table><caption>**Table A.1.** Parameters of the models: inclination (i), velocity direction (σ), magnetic field ($b_{eq}$, $\phi'$, $b_z$), elevation (η) and equatorial radius ($r_e$).</caption><thead><tr><th>Figure</th><th><i>i</i><br>(deg)</th><th>σ<br>(deg)</th><th><i>b</i><sub>eq</sub></th><th>φ'<br>(deg)</th><th><i>b</i><sub>z</sub></th><th>η<br>(deg)</th><th><i>r</i><sub>e</sub></th></tr></thead><tbody><tr><td>Fig. 3a</td><td>30</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 3b</td><td>60</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 3c</td><td>80</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 3d</td><td>30</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 3e</td><td>60</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 3f</td><td>80</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5a</td><td>30</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5b</td><td>60</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5c</td><td>80</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5d</td><td>30</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5e</td><td>60</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5f</td><td>80</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5g</td><td>30</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5h</td><td>60</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 5i</td><td>80</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6a</td><td>30</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6b</td><td>60</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6c</td><td>80</td><td>90</td><td>1</td><td>90</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6d</td><td>30</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6e</td><td>60</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6f</td><td>80</td><td>90</td><td>1</td><td>0</td><td>0</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6g</td><td>30</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6h</td><td>60</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 6i</td><td>80</td><td>90</td><td>0</td><td>0</td><td>1</td><td>20</td><td>3, 5, 15, 50</td></tr><tr><td>Fig. 7a</td><td>30</td><td>90, 120, 150, 180</td><td>1</td><td>90</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7b</td><td>60</td><td>90, 120, 150, 180</td><td>1</td><td>90</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7c</td><td>80</td><td>90, 120, 150, 180</td><td>1</td><td>90</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7d</td><td>30</td><td>90, 120, 150, 180</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7e</td><td>60</td><td>90, 120, 150, 180</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7f</td><td>80</td><td>90, 120, 150, 180</td><td>1</td><td>0</td><td>0</td><td>0</td><td>5</td></tr><tr><td>Fig. 7g</td><td>30</td><td>90, 120, 150, 180</td><td>0</td><td>0</td><td>1</td><td>0</td><td>5</td></tr><tr><td>Fig. 7h</td><td>60</td><td>90, 120, 150, 180</td><td>0</td><td>0</td><td>1</td><td>0</td><td>5</td></tr><tr><td>Fig. 7i</td><td>80</td><td>90, 120, 150, 180</td><td>0</td><td>0</td><td>1</td><td>0</td><td>5</td></tr><tr><td>Fig. 9a</td><td>165</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9b</td><td>150</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9c</td><td>140</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9d</td><td>165</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9e</td><td>150</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9f</td><td>140</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9g</td><td>165</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9h</td><td>150</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 9i</td><td>140</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)a</td><td>165</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)b</td><td>150</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)c</td><td>140</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)d</td><td>165</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)e</td><td>150</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)f</td><td>140</td><td>90</td><td>1</td><td>0</td><td>0</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)g</td><td>165</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)h</td><td>150</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 10(11)i</td><td>140</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0...45(-45)</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 12a,d</td><td>165</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 12b,e</td><td>150</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 12c,f</td><td>130</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 13</td><td>0...90</td><td>90</td><td>0</td><td>0</td><td>1</td><td>0</td><td>3, 4, 5, 7</td></tr><tr><td>Fig. 14</td><td>0...90</td><td>90</td><td>1</td><td>90</td><td>0</td><td>0</td><td>3, 7, 10<sup>4</sup></td></tr></tbody></table>