\volnopage

Vol.0 (20xx) No.0, 000–000

11institutetext: Physics Department, King Fahd University of Petroleum and Minerals, Dhahran 31261, Saudi Arabia; [email protected] ; [email protected]
\vs\noReceived 20xx month day; accepted 20xx month day

Charged Particles Capture Cross-Section by a Weakly Charged Schwarzschild Black Hole

A. M. Al Zahrani    11    A. Al-Jama 221122
Abstract

We study the capture cross-section of charged particles by a weakly charged Schwarzschild black hole. The dependence of the maximum impact parameter for capture on the particle’s energy is investigated numerically for different values of the electromagnetic coupling strength between the particle and the black hole. The capture cross-section is then calculated. We show that the capture cross-section is independent of the electromagnetic coupling for ultra-relativistic particles. The astrophysical implications of our results are discussed.

keywords:
stars: black holes — Stars, accretion, accretion disks — Physical Data and Processes, black hole physics — Physical Data and Processes

1 Introduction

Studying the capture cross-section of black holes is central to understand how matter interacts with them. It helps us understand the process of matter accretion by a black hole which in turn determines how its mass, angular momentum and charge evolve. It can also help us understand the environment near black holes. Moreover, scrutinizing capture cross-section can be used to test theories of gravity in strong gravitational fields.

Astrophysicists generally assume black holes are electrically neutral. This is because they would quickly attract oppositely charged matter to balance out any access charge. However, there are compelling reasons why weakly charged black holes might exist as discussed in Al Zahrani (2021, 2022); Zajaček & Tursunov (2019); Zajaček et al. (2019, 2018); Carter (1973) and the references therein. The differences in how a black hole accretes electrons and protons within its plasma environment, influenced by radiation, could render it charged. Also, the spin of a black hole in the presence of a magnetic field can induce the accretion of charged particles. In fact, using the EHT observations, it was inferred that Sgr A* and M87 can be charged Ghosh & Afrin (2023); Kocherlakota et al. (2021). The black hole’s charge is weak in the sense that it has no tangible effect on spacetime, but its effect on charged particle dynamics is prominent.

There are numerous astrophysical scenarios wherein charged particles are drawn into black holes. Stars within the Roche limits near black holes often contribute matter through tidal interactions. Additionally, stars emit streams of charged particles as stellar winds. Highly energetic charged particles, resulting from supernovae, gamma-ray bursts, and bipolar jets from compact objects, frequently find their way into the vicinity of black holes. These processes collectively enrich the environment around black holes with a significant population of charged particles.

The concept of capture cross-sections has been explored extensively for various black hole types. Foundational treatment which examine photon and neutral particle capture by Schwarzschild black holes was given in several monographs, such as Frolov & Zelnikov (2011). Further work addressed capture cross-sections of charged and neutral particles by Kerr-Newman black holes, including the implications for black hole spin and charge evolution Young (1976). Capture by Reissner-Nordström black holes was also investigated Zakharov (1994). In the context of higher-dimensional black holes, studies have focused on calculating photon critical impact parameters for Schwarzschild-Tangherlini black holes Connell & Frolov (2008); Tsukamoto (2014); Singh & Ghosh (2018); Bugden (2020). The capture cross-section for massive particles was determined in Ahmedov et al. (2021). Additionally, research extends to particle capture in Myers-Perry rotating spacetime which describes rotating black holes in five-dimensions Gooding & Frolov (2008). Moreover, wave capture cross-sections have been studied for various black hole configurations (see Anacleto, et al. (2023) and the references within).

In this research, we examine the capture cross-section of charged particles by a weakly charged Schwarzschild black hole and discuss the astrophysical consequences of our findings. The paper is organized as follows: In Sec. 2, we review the dynamics of charged particles in the background of a weakly charged black hole. We then review the capture cross-section of neutral particles in Sec. 3. The capture cross-section of charged particles is calculated for different coupling strengths and particle energies in Secs. 4. Finally, we summarize our main findings and discuss their astrophysical consequences in Sec 5. We use the sign conventions adopted in Misner et al. (1973) and geometrized units where c=G=k=1𝑐𝐺𝑘1c=G=k=1italic_c = italic_G = italic_k = 1, where k𝑘kitalic_k is the electrostatic constant.

2 Charged Particles near a Weakly Charged Schwarzschild Black Hole

Here, we review the dynamics of charged particles near a weakly charged black hole. The spacetime geometry around a black hole of mass M𝑀Mitalic_M and charge Q𝑄Qitalic_Q is described by the Schwarzschild Reissner-Nordströn metric which reads Misner et al. (1973)

ds2=hdt2+h1dr2+r2dθ2+r2sin2θdϕ2,𝑑superscript𝑠2𝑑superscript𝑡2superscript1𝑑superscript𝑟2superscript𝑟2𝑑superscript𝜃2superscript𝑟2superscript2𝜃𝑑superscriptitalic-ϕ2ds^{2}=-hdt^{2}+h^{-1}dr^{2}+r^{2}d\theta^{2}+r^{2}\sin^{2}{\theta}d\phi^{2},italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - italic_h italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_h start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_d italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (1)

where h=1rS/r+Q2/r21subscript𝑟𝑆𝑟superscript𝑄2superscript𝑟2h=1-r_{S}/r+Q^{2}/r^{2}italic_h = 1 - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT / italic_r + italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT / italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and rS=2Msubscript𝑟𝑆2𝑀r_{S}=2Mitalic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 2 italic_M is the Schwarzschild radius. The electromagnetic 4-potential is

Aμ=Qrδμ0.subscript𝐴𝜇𝑄𝑟superscriptsubscript𝛿𝜇0A_{\mu}=-\frac{Q}{r}\delta_{\mu}^{0}.italic_A start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = - divide start_ARG italic_Q end_ARG start_ARG italic_r end_ARG italic_δ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT . (2)

However, when the charge is weak we can ignore the curvature due to it and use the Schwarzschild metric, which reads Misner et al. (1973)

ds2=fdt2+f1dr2+r2dθ2+r2sin2θdϕ2,𝑑superscript𝑠2𝑓𝑑superscript𝑡2superscript𝑓1𝑑superscript𝑟2superscript𝑟2𝑑superscript𝜃2superscript𝑟2superscript2𝜃𝑑superscriptitalic-ϕ2ds^{2}=-fdt^{2}+f^{-1}dr^{2}+r^{2}d\theta^{2}+r^{2}\sin^{2}{\theta}d\phi^{2},italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - italic_f italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_f start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_θ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ italic_d italic_ϕ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (3)

where f=1rS/r𝑓1subscript𝑟𝑆𝑟f=1-r_{S}/ritalic_f = 1 - italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT / italic_r and rS=2Msubscript𝑟𝑆2𝑀r_{S}=2Mitalic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT = 2 italic_M is the Schwarzschild radius. This weak charge approximation is valid unless the charge creates curvature comparable to that due to the black hole’s mass. This happens when

Q2M2.similar-tosuperscript𝑄2superscript𝑀2Q^{2}\sim M^{2}.italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∼ italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (4)

In conventional units, the weak charge approximation fails when

QG1/2Mk1/21020MMcoloumbs.similar-to𝑄superscript𝐺12𝑀superscript𝑘12similar-tosuperscript1020𝑀subscript𝑀direct-productcoloumbsQ\sim\frac{G^{1/2}M}{k^{1/2}}\sim 10^{20}\>\frac{M}{\>M_{\odot}}\>\text{% coloumbs}.italic_Q ∼ divide start_ARG italic_G start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT italic_M end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT end_ARG ∼ 10 start_POSTSUPERSCRIPT 20 end_POSTSUPERSCRIPT divide start_ARG italic_M end_ARG start_ARG italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT end_ARG coloumbs . (5)

This charge is way greater than the greatest estimated change on any black hole. Although the black hole charge is tiny, its effect on charged particles dynamics is profound because it is multiplied by the charge-to-mass ratio of these particles (1021m1similar-toabsentsuperscript1021superscriptm1\sim 10^{21}\>{\text{m}}^{-1}∼ 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT m start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT for electrons and 1018m1similar-toabsentsuperscript1018superscriptm1\sim 10^{18}\>{\text{m}}^{-1}∼ 10 start_POSTSUPERSCRIPT 18 end_POSTSUPERSCRIPT m start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT for protons).

The Lagrangian describing a charged particle of charge q𝑞qitalic_q and mass m𝑚mitalic_m in a spacetime described by a metric gμνsubscript𝑔𝜇𝜈g_{\mu\nu}italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT and an electromagnetic field produced by a 4-potential Aμsuperscript𝐴𝜇A^{\mu}italic_A start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT reads Chandrasekhar (1983)

L=12mgμνuμuν+quμAμ,𝐿12𝑚subscript𝑔𝜇𝜈superscript𝑢𝜇superscript𝑢𝜈𝑞superscript𝑢𝜇subscript𝐴𝜇L=\frac{1}{2}mg_{\mu\nu}u^{\mu}u^{\nu}+qu^{\mu}A_{\mu},italic_L = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_m italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT + italic_q italic_u start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_A start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT , (6)

where uμdxμ/dτsuperscript𝑢𝜇𝑑superscript𝑥𝜇𝑑𝜏u^{\mu}\equiv dx^{\mu}/d\tauitalic_u start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT ≡ italic_d italic_x start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT / italic_d italic_τ is the particle’s 4-velocity and τ𝜏\tauitalic_τ is its proper time. In our case, the Lagrangian becomes

L𝐿\displaystyle Litalic_L =\displaystyle== 12m[f(dtdτ)2+f1(drdτ)2+r2(dθdτ)2+r2sin2θ(dϕdτ)2]12𝑚delimited-[]𝑓superscript𝑑𝑡𝑑𝜏2superscript𝑓1superscript𝑑𝑟𝑑𝜏2superscript𝑟2superscript𝑑𝜃𝑑𝜏2superscript𝑟2superscript2𝜃superscript𝑑italic-ϕ𝑑𝜏2\displaystyle\frac{1}{2}m\left[-f\left(\frac{dt}{d\tau}\right)^{2}+f^{-1}\left% (\frac{dr}{d\tau}\right)^{2}+r^{2}\left(\frac{d\theta}{d\tau}\right)^{2}+r^{2}% \sin^{2}{\theta}\left(\frac{d\phi}{d\tau}\right)^{2}\right]divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_m [ - italic_f ( divide start_ARG italic_d italic_t end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_f start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG italic_d italic_r end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_d italic_θ end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ ( divide start_ARG italic_d italic_ϕ end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] (7)
qQdtdτ.𝑞𝑄𝑑𝑡𝑑𝜏\displaystyle-qQ\frac{dt}{d\tau}.- italic_q italic_Q divide start_ARG italic_d italic_t end_ARG start_ARG italic_d italic_τ end_ARG .

This Lagrangian is cyclic in t𝑡titalic_t and ϕitalic-ϕ\phiitalic_ϕ, which means that the particle’s energy and azimuthal angular momentum are constants of motion. The specific energy and azimuthal angular momentum are, respectively, given by

\displaystyle{\cal E}caligraphic_E =\displaystyle== 1mL(dtdτ)=fdtdτ+qQmr,1𝑚𝐿𝑑𝑡𝑑𝜏𝑓𝑑𝑡𝑑𝜏𝑞𝑄𝑚𝑟\displaystyle-\frac{1}{m}\frac{\partial L}{\partial\left(\frac{dt}{d\tau}% \right)}=f\frac{dt}{d\tau}+\frac{qQ}{mr},- divide start_ARG 1 end_ARG start_ARG italic_m end_ARG divide start_ARG ∂ italic_L end_ARG start_ARG ∂ ( divide start_ARG italic_d italic_t end_ARG start_ARG italic_d italic_τ end_ARG ) end_ARG = italic_f divide start_ARG italic_d italic_t end_ARG start_ARG italic_d italic_τ end_ARG + divide start_ARG italic_q italic_Q end_ARG start_ARG italic_m italic_r end_ARG , (8)
\displaystyle\ellroman_ℓ =\displaystyle== 1mL(dϕdτ)=r2sin2θdϕdτ1𝑚𝐿𝑑italic-ϕ𝑑𝜏superscript𝑟2superscript2𝜃𝑑italic-ϕ𝑑𝜏\displaystyle\frac{1}{m}\frac{\partial L}{\partial\left(\frac{d\phi}{d\tau}% \right)}=r^{2}\sin^{2}{\theta}\frac{d\phi}{d\tau}divide start_ARG 1 end_ARG start_ARG italic_m end_ARG divide start_ARG ∂ italic_L end_ARG start_ARG ∂ ( divide start_ARG italic_d italic_ϕ end_ARG start_ARG italic_d italic_τ end_ARG ) end_ARG = italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ divide start_ARG italic_d italic_ϕ end_ARG start_ARG italic_d italic_τ end_ARG (9)

Combining these equations with the normalization condition gμνuμuν=1subscript𝑔𝜇𝜈superscript𝑢𝜇superscript𝑢𝜈1g_{\mu\nu}u^{\mu}u^{\nu}=-1italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT italic_u start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT italic_u start_POSTSUPERSCRIPT italic_ν end_POSTSUPERSCRIPT = - 1 and solving for dr/dτ𝑑𝑟𝑑𝜏dr/d\tauitalic_d italic_r / italic_d italic_τ give

(drdτ)2=(qQmr)2f[r2(dθdτ)2+2r2sin2θ+1].superscript𝑑𝑟𝑑𝜏2superscript𝑞𝑄𝑚𝑟2𝑓delimited-[]superscript𝑟2superscript𝑑𝜃𝑑𝜏2superscript2superscript𝑟2superscript2𝜃1\left(\frac{dr}{d\tau}\right)^{2}=\left({\cal E}-\frac{qQ}{mr}\right)^{2}-f% \left[r^{2}\left(\frac{d\theta}{d\tau}\right)^{2}+\frac{\ell^{2}}{r^{2}\sin^{2% }{\theta}}+1\right].( divide start_ARG italic_d italic_r end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( caligraphic_E - divide start_ARG italic_q italic_Q end_ARG start_ARG italic_m italic_r end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_f [ italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_d italic_θ end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG roman_ℓ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_θ end_ARG + 1 ] . (10)

In the equatorial plane where θ=π/2𝜃𝜋2\theta=\pi/2italic_θ = italic_π / 2, the equation becomes

(drdτ)2=(qQmr)2f(2r2+1).superscript𝑑𝑟𝑑𝜏2superscript𝑞𝑄𝑚𝑟2𝑓superscript2superscript𝑟21\left(\frac{dr}{d\tau}\right)^{2}=\left({\cal E}-\frac{qQ}{mr}\right)^{2}-f% \left(\frac{\ell^{2}}{r^{2}}+1\right).( divide start_ARG italic_d italic_r end_ARG start_ARG italic_d italic_τ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( caligraphic_E - divide start_ARG italic_q italic_Q end_ARG start_ARG italic_m italic_r end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_f ( divide start_ARG roman_ℓ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 1 ) . (11)

Let us rewrite the last equation in a dimensionless form. We first introduce the following dimensionless quantities:

𝒯=τrS,ρ=rrS,=rS.formulae-sequence𝒯𝜏subscript𝑟𝑆formulae-sequence𝜌𝑟subscript𝑟𝑆subscript𝑟𝑆{\cal T}=\frac{\tau}{r_{S}},\hskip 28.45274pt\rho=\frac{r}{r_{S}},\hskip 28.45% 274pt{\cal L}=\frac{\ell}{r_{S}}.caligraphic_T = divide start_ARG italic_τ end_ARG start_ARG italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG , italic_ρ = divide start_ARG italic_r end_ARG start_ARG italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG , caligraphic_L = divide start_ARG roman_ℓ end_ARG start_ARG italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG . (12)

Equation 11 then becomes

(dρd𝒯)2=(αρ)2f(2ρ2+1),superscript𝑑𝜌𝑑𝒯2superscript𝛼𝜌2𝑓superscript2superscript𝜌21\left(\frac{d\rho}{d{\cal T}}\right)^{2}=\left({\cal E}-\frac{\alpha}{\rho}% \right)^{2}-f\left(\frac{{\cal L}^{2}}{\rho^{2}}+1\right),( divide start_ARG italic_d italic_ρ end_ARG start_ARG italic_d caligraphic_T end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( caligraphic_E - divide start_ARG italic_α end_ARG start_ARG italic_ρ end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_f ( divide start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 1 ) , (13)

where

α=qQmrS.𝛼𝑞𝑄𝑚subscript𝑟𝑆\alpha=\frac{qQ}{mr_{S}}.italic_α = divide start_ARG italic_q italic_Q end_ARG start_ARG italic_m italic_r start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT end_ARG . (14)

The parameter α𝛼\alphaitalic_α represent the relative strength of the electromagnetic force to the Newtonian gravitational force. We can rewrite Eq. 13 as

(dρd𝒯)2=(V+)(V),superscript𝑑𝜌𝑑𝒯2subscript𝑉subscript𝑉\left(\frac{d\rho}{d{\cal T}}\right)^{2}=({\cal E}-V_{+})({\cal E}-V_{-}),( divide start_ARG italic_d italic_ρ end_ARG start_ARG italic_d caligraphic_T end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( caligraphic_E - italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ) ( caligraphic_E - italic_V start_POSTSUBSCRIPT - end_POSTSUBSCRIPT ) , (15)

where

V±=αρ±f(2ρ2+1),subscript𝑉plus-or-minusplus-or-minus𝛼𝜌𝑓superscript2superscript𝜌21V_{\pm}=\frac{\alpha}{\rho}\pm\sqrt{f\left(\frac{{\cal L}^{2}}{\rho^{2}}+1% \right)},italic_V start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = divide start_ARG italic_α end_ARG start_ARG italic_ρ end_ARG ± square-root start_ARG italic_f ( divide start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 1 ) end_ARG , (16)

is an effective potential. It is V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT that corresponds to physical, future-directed motion and hence will be used in all of the analyses below. Without loss of generality, we will consider >00{\cal L}>0caligraphic_L > 0 only.

It was estimated in Ref. Zajaček et al. (2019) that the charge of Sgr A* is 1081015superscript108superscript101510^{8}-10^{15}10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT coulomb. Using the lower limit of charge, the coupling constant for electrons αesubscript𝛼𝑒\alpha_{e}italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT and protons αpsubscript𝛼𝑝\alpha_{p}italic_α start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT near Sgr A*, which has a mass of M=4.3×106M𝑀4.3superscript106subscript𝑀direct-productM=4.3\times 10^{6}M_{\odot}italic_M = 4.3 × 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT according to Ref. GRAVITY Collaboration (2023), are

αe109,similar-tosubscript𝛼𝑒superscript109\displaystyle\alpha_{e}\sim 10^{9},italic_α start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT , (17)
αp106.similar-tosubscript𝛼𝑝superscript106\displaystyle\alpha_{p}\sim 10^{6}.italic_α start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT . (18)

3 Capture Cross-Section of Neutral Particles

Before we tackle the main problem, let us find the capture cross-cross section for neutral particles first. Setting α=0𝛼0\alpha=0italic_α = 0, the effective potential V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT reduces to

V+=f(2ρ2+1).subscript𝑉𝑓superscript2superscript𝜌21V_{+}=\sqrt{f\left(\frac{{\cal L}^{2}}{\rho^{2}}+1\right)}.italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT = square-root start_ARG italic_f ( divide start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 1 ) end_ARG . (19)

Capture occurs whenever the particle’s energy is greater than the maximum of V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. The function V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT is at an extremum when dV+/dρ=0𝑑subscript𝑉𝑑𝜌0dV_{+}/d\rho=0italic_d italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT / italic_d italic_ρ = 0 oder

ρ2+(32ρ)2=0,superscript𝜌232𝜌superscript20\rho^{2}+(3-2\rho){\cal L}^{2}=0,italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 3 - 2 italic_ρ ) caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 , (20)

which gives the position of the extrema in terms of {\cal L}caligraphic_L as

ρ±=2±23,subscript𝜌plus-or-minusplus-or-minussuperscript2superscript23\rho_{\pm}={\cal L}^{2}\pm{\cal L}\sqrt{{\cal L}^{2}-3},italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ± caligraphic_L square-root start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 end_ARG , (21)

where [3,)3{\cal L}\in[\sqrt{3},\infty)caligraphic_L ∈ [ square-root start_ARG 3 end_ARG , ∞ ). When =33{\cal L}=\sqrt{3}caligraphic_L = square-root start_ARG 3 end_ARG (minabsentsubscriptmin\equiv{\cal L}_{\text{min}}≡ caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT), ρ+subscript𝜌\rho_{+}italic_ρ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ρsubscript𝜌\rho_{-}italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT meet at a saddle point. Inspecting d2V+/dρ2superscript𝑑2subscript𝑉𝑑superscript𝜌2d^{2}V_{+}/d\rho^{2}italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT / italic_d italic_ρ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT reveals that ρsubscript𝜌\rho_{-}italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT corresponds to the position of the local maximum of V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT. In terms of {\cal L}caligraphic_L, the escape condition =V+|ρ=ρevaluated-atsubscript𝑉𝜌subscript𝜌{\cal E}=V_{+}|_{\rho=\rho_{-}}caligraphic_E = italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT | start_POSTSUBSCRIPT italic_ρ = italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT end_POSTSUBSCRIPT becomes

=227[(23+)323+9]1/2,227superscriptdelimited-[]superscript233superscript23912{\cal E}=\sqrt{\frac{2}{27}}\left[{\cal L}\left(\sqrt{{\cal L}^{2}-3}+{\cal L}% \right)-\frac{3\sqrt{{\cal L}^{2}-3}}{{\cal L}}+9\right]^{1/2},caligraphic_E = square-root start_ARG divide start_ARG 2 end_ARG start_ARG 27 end_ARG end_ARG [ caligraphic_L ( square-root start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 end_ARG + caligraphic_L ) - divide start_ARG 3 square-root start_ARG caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 end_ARG end_ARG start_ARG caligraphic_L end_ARG + 9 ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , (22)

Inverting this equation gives

=[274362+(928)3/2+88(21)]1/2superscriptdelimited-[]27superscript436superscript2superscript9superscript283288superscript2112\displaystyle{\cal L}=\left[\frac{27{\cal E}^{4}-36{\cal E}^{2}+\left(9{\cal E% }^{2}-8\right)^{3/2}{\cal E}+8}{8\left({\cal E}^{2}-1\right)}\right]^{1/2}caligraphic_L = [ divide start_ARG 27 caligraphic_E start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 36 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 9 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 8 ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT caligraphic_E + 8 end_ARG start_ARG 8 ( caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) end_ARG ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT (23)

The impact parameter b𝑏bitalic_b is defined as the perpendicular distance between the center of force and the incident velocity Goldstein et al. (2001). It can be written as

b=𝒫=21,𝑏𝒫superscript21b=\frac{{\cal L}}{\cal P}=\frac{{\cal L}}{\sqrt{{\cal E}^{2}-1}},italic_b = divide start_ARG caligraphic_L end_ARG start_ARG caligraphic_P end_ARG = divide start_ARG caligraphic_L end_ARG start_ARG square-root start_ARG caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 end_ARG end_ARG , (24)

where 𝒫𝒫\cal Pcaligraphic_P is the specific linear momentum. The maximum impact parameter for capture bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT is given by

bmax=[274362+(928)3/2+8]22(21)1/2.subscript𝑏maxsuperscriptdelimited-[]27superscript436superscript2superscript9superscript2832822superscript2112b_{\text{max}}=\frac{\left[27{\cal E}^{4}-36{\cal E}^{2}+\left(9{\cal E}^{2}-8% \right)^{3/2}{\cal E}+8\right]}{2\sqrt{2}\left({\cal E}^{2}-1\right)}^{1/2}.italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT = divide start_ARG [ 27 caligraphic_E start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 36 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 9 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 8 ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT caligraphic_E + 8 ] end_ARG start_ARG 2 square-root start_ARG 2 end_ARG ( caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) end_ARG start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT . (25)

The capture cross-section σcapsubscript𝜎cap{\sigma}_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT is given by

σcap=πbmax2=π8274362+(928)3/2+8(21)2.subscript𝜎cap𝜋superscriptsubscript𝑏max2𝜋827superscript436superscript2superscript9superscript28328superscriptsuperscript212{\sigma}_{\text{cap}}=\pi b_{\text{max}}^{2}=\frac{\pi}{8}\frac{27{\cal E}^{4}% -36{\cal E}^{2}+\left(9{\cal E}^{2}-8\right)^{3/2}{\cal E}+8}{\left({\cal E}^{% 2}-1\right)^{2}}.italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT = italic_π italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = divide start_ARG italic_π end_ARG start_ARG 8 end_ARG divide start_ARG 27 caligraphic_E start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 36 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( 9 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 8 ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT caligraphic_E + 8 end_ARG start_ARG ( caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (26)
Refer to caption
Figure 1: The maximum impact parameter for capture bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a neutral particle.
Refer to caption
Figure 2: The capture cross-section σcapsubscript𝜎cap{\sigma}_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a neutral particle.

Figures 1 and 2 are plots of bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT and the capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E, respectively. For ultra-relativistic particles (1much-greater-than1{\cal E}\gg 1caligraphic_E ≫ 1),

bmax=332+322+𝒪(13).subscript𝑏max33232superscript2𝒪1superscript3b_{\text{max}}=\frac{3\sqrt{3}}{2}+\frac{\sqrt{3}}{2{\cal E}^{2}}+{\mathcal{O}% }\left(\frac{1}{{\cal E}^{3}}\right).italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT = divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG + divide start_ARG square-root start_ARG 3 end_ARG end_ARG start_ARG 2 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG caligraphic_E start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ) . (27)

The corresponding capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT is therefore

σcap=27π4+9π22+𝒪(13).subscript𝜎cap27𝜋49𝜋2superscript2𝒪1superscript3\sigma_{\text{cap}}=\frac{27\pi}{4}+\frac{9\pi}{2{\cal E}^{2}}+{\mathcal{O}}% \left(\frac{1}{{\cal E}^{3}}\right).italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT = divide start_ARG 27 italic_π end_ARG start_ARG 4 end_ARG + divide start_ARG 9 italic_π end_ARG start_ARG 2 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG caligraphic_E start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ) . (28)

For a slowly moving particle with speed v1much-less-than𝑣1v\ll 1italic_v ≪ 1,

1+v22,1superscript𝑣22{\cal E}\approx 1+\frac{v^{2}}{2},caligraphic_E ≈ 1 + divide start_ARG italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG , (29)

and thus

bmax=21+𝒪(1)=2v+𝒪(v),subscript𝑏max21𝒪12𝑣𝒪𝑣b_{\text{max}}=\frac{\sqrt{2}}{\sqrt{{\cal E}-1}}+{\mathcal{O}}(\sqrt{{\cal E}% -1})=\frac{2}{v}+{\mathcal{O}}(v),italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT = divide start_ARG square-root start_ARG 2 end_ARG end_ARG start_ARG square-root start_ARG caligraphic_E - 1 end_ARG end_ARG + caligraphic_O ( square-root start_ARG caligraphic_E - 1 end_ARG ) = divide start_ARG 2 end_ARG start_ARG italic_v end_ARG + caligraphic_O ( italic_v ) , (30)

and the capture cross-section becomes

σcap=4πv2+𝒪(v0).subscript𝜎cap4𝜋superscript𝑣2𝒪superscript𝑣0\sigma_{\text{cap}}=\frac{4\pi}{v^{2}}+{\mathcal{O}}(v^{0}).italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT = divide start_ARG 4 italic_π end_ARG start_ARG italic_v start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( italic_v start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT ) . (31)

4 Capture Cross-Section of Charged Particles

We will now follow the same procedure we used for the neutral particle. However, analytic expressions are not viable in this case and we will resort to numerical solutions, except in the ultra-relativistic particle case. The structure of the effective potential V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT is generically similar to the neutral particle’s. The effect of α𝛼\alphaitalic_α is to raise (lower) the peak of V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT for positive (negative) α𝛼\alphaitalic_α. The effective potential V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT is at an extremum when

2α(ρ±1)(2+ρ±2)ρ±2(32ρ±)ρ±2=0.2𝛼subscript𝜌plus-or-minus1superscript2superscriptsubscript𝜌plus-or-minus2subscript𝜌plus-or-minussuperscript232subscript𝜌plus-or-minussuperscriptsubscript𝜌plus-or-minus202\alpha\sqrt{(\rho_{\pm}-1)\left({\cal L}^{2}+\rho_{\pm}^{2}\right)\rho_{\pm}}% -{\cal L}^{2}(3-2\rho_{\pm})-\rho_{\pm}^{2}=0.2 italic_α square-root start_ARG ( italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT - 1 ) ( caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT end_ARG - caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 3 - 2 italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ) - italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 . (32)

The extremum is a maximum when

α[2(12ρ±)+(34ρ±)ρ±2](ρ±1)(2+ρ±2)ρ±+2ρ±22<0.𝛼delimited-[]superscript212subscript𝜌plus-or-minus34subscript𝜌plus-or-minussuperscriptsubscript𝜌plus-or-minus2subscript𝜌plus-or-minus1superscript2superscriptsubscript𝜌plus-or-minus2subscript𝜌plus-or-minus2subscript𝜌plus-or-minus2superscript20\frac{\alpha\left[{\cal L}^{2}(1-2\rho_{\pm})+(3-4\rho_{\pm})\rho_{\pm}^{2}% \right]}{\sqrt{(\rho_{\pm}-1)\left({\cal L}^{2}+\rho_{\pm}^{2}\right)\rho_{\pm% }}}+2\rho_{\pm}-2{\cal L}^{2}<0.divide start_ARG italic_α [ caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( 1 - 2 italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ) + ( 3 - 4 italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ) italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] end_ARG start_ARG square-root start_ARG ( italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT - 1 ) ( caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT end_ARG end_ARG + 2 italic_ρ start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT - 2 caligraphic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < 0 . (33)

To be consistent with the notation of the previous section, we let ρ+subscript𝜌\rho_{+}italic_ρ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT correspond to the minimum of V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ρsubscript𝜌\rho_{-}italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT correspond to the maximum. Here, minsubscriptmin{\cal L}_{\text{min}}caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT (the value at which ρsubscript𝜌\rho_{-}italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT and ρ+subscript𝜌\rho_{+}italic_ρ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT meet) depends on the value of α𝛼\alphaitalic_α. The two parameters are related by the relation

α8+6α4min2(min23)8α2min4(min2+9)+3min4(min23)2=0.superscript𝛼86superscript𝛼4superscriptsubscriptmin2superscriptsubscriptmin238superscript𝛼2superscriptsubscriptmin4superscriptsubscriptmin293superscriptsubscriptmin4superscriptsuperscriptsubscriptmin2320-\alpha^{8}+6\alpha^{4}{\cal L}_{\text{min}}^{2}\left({\cal L}_{\text{min}}^{2% }-3\right)-8\alpha^{2}{\cal L}_{\text{min}}^{4}\left({\cal L}_{\text{min}}^{2}% +9\right)+3{\cal L}_{\text{min}}^{4}\left({\cal L}_{\text{min}}^{2}-3\right)^{% 2}=0.- italic_α start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT + 6 italic_α start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 ) - 8 italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 9 ) + 3 caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ( caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 . (34)

Figure 3 is a plot of minsubscriptmin{\cal L}_{\text{min}}caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT vs α𝛼\alphaitalic_α. When α=1/2𝛼12\alpha=1/2italic_α = 1 / 2, minsubscriptmin{\cal L}_{\text{min}}caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT approaches zero. This is because V+subscript𝑉V_{+}italic_V start_POSTSUBSCRIPT + end_POSTSUBSCRIPT ceases to have a local minimum for α1/2𝛼12\alpha\geq 1/2italic_α ≥ 1 / 2. Physically, this limit corresponds to the case when the Coulomb repulsion becomes too strong for stable orbits to exist as discussed in Ref. Al Zahrani (2021).

Refer to caption
Figure 3: The value of {\cal L}caligraphic_L at which ρ+subscript𝜌\rho_{+}italic_ρ start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and ρsubscript𝜌\rho_{-}italic_ρ start_POSTSUBSCRIPT - end_POSTSUBSCRIPT meet (minsubscriptmin{\cal L}_{\text{min}}caligraphic_L start_POSTSUBSCRIPT min end_POSTSUBSCRIPT) vs. the electromagnetic coupling parameter α𝛼\alphaitalic_α.

Figure 4 shows how bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT depends on {\cal E}caligraphic_E for several negative values of the coupling parameter α𝛼\alphaitalic_α. The effect of increasing |α|𝛼|\alpha|| italic_α | is to increase the values of bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT for all energies. This is expected because the Coulombs attraction makes it easier for a charges particle to get captured. In all cases, bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT is a monotonic function of {\cal E}caligraphic_E. In the ultra-relativistic limit, bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT approaches 33/23323\sqrt{3}/23 square-root start_ARG 3 end_ARG / 2, the limit in the neutral particle case, for any finite value of α𝛼\alphaitalic_α, provided that α𝛼\alphaitalic_α is not too large compared to {\cal E}caligraphic_E.

Refer to caption
Figure 4: The maximum impact parameter for capture bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with α=0𝛼0\alpha=0italic_α = 0 (black), α=1𝛼1\alpha=-1italic_α = - 1 (blue), α=2𝛼2\alpha=-2italic_α = - 2 (green), α=3𝛼3\alpha=-3italic_α = - 3 (red).

Figure 5 shows how bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT depends on {\cal E}caligraphic_E for several values of α𝛼\alphaitalic_α between 00 and 0.50.50.50.5. In this range, there is competition between that gravitational ’attraction’ and the Coulomb repulsion. The curves have richer structure. They falls quickly as {\cal E}caligraphic_E goes beyond 1111 and reach a minimum. After that, the curves rise and reach 33/23323\sqrt{3}/23 square-root start_ARG 3 end_ARG / 2 asymptotically.

Refer to caption
Figure 5: The maximum impact parameter for capture bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with α=0𝛼0\alpha=0italic_α = 0 (black), α=0.1𝛼0.1\alpha=0.1italic_α = 0.1 (blue), α=0.3𝛼0.3\alpha=0.3italic_α = 0.3 (green), α=0.5𝛼0.5\alpha=0.5italic_α = 0.5 (red).

Figure 6 shows how bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT depends on {\cal E}caligraphic_E for several positive values of α𝛼\alphaitalic_α greater than 0.50.50.50.5. Generally, bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT becomes smaller as α𝛼\alphaitalic_α increases. This is expected because the greater the Coulomb repulsion the more difficult it is for a charged particle to be captured. In fact, there is a threshold energy thrsubscriptthr{\cal E}_{\text{thr}}caligraphic_E start_POSTSUBSCRIPT thr end_POSTSUBSCRIPT below which capture cannot occur. It is given by

thr=α+14α.subscriptthr𝛼14𝛼{\cal E}_{\text{thr}}=\alpha+\frac{1}{4\alpha}.caligraphic_E start_POSTSUBSCRIPT thr end_POSTSUBSCRIPT = italic_α + divide start_ARG 1 end_ARG start_ARG 4 italic_α end_ARG . (35)

This equation is valid for α0.5𝛼0.5\alpha\geq 0.5italic_α ≥ 0.5 only. Fig. 7 shows how thrsubscriptthr{\cal E}_{\text{thr}}caligraphic_E start_POSTSUBSCRIPT thr end_POSTSUBSCRIPT vary with α𝛼\alphaitalic_α.
 
The capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT corresponding to Figs. 4, 5 and 6 is shown in Figs. 8, 10 and 9, respectively. In all cases, σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E curves inherent the features of the bminsubscript𝑏minb_{\text{min}}italic_b start_POSTSUBSCRIPT min end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E curves.

Refer to caption
Figure 6: The max impact parameter for capture bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with with α=0𝛼0\alpha=0italic_α = 0 (black), α=1𝛼1\alpha=1italic_α = 1 (blue), α=2𝛼2\alpha=2italic_α = 2 (green), α=3𝛼3\alpha=3italic_α = 3 (red).
Refer to caption
Figure 7: The energy threshold for escape thrsubscriptthr{\cal E}_{\text{thr}}caligraphic_E start_POSTSUBSCRIPT thr end_POSTSUBSCRIPT vs. electromagnetic coupling parameter α𝛼\alphaitalic_α.
Refer to caption
Figure 8: The capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with α=0𝛼0\alpha=0italic_α = 0 (black), α=1𝛼1\alpha=-1italic_α = - 1 (blue), α=2𝛼2\alpha=-2italic_α = - 2 (green), α=3𝛼3\alpha=-3italic_α = - 3 (red).
Refer to caption
Figure 9: The capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with α=0𝛼0\alpha=0italic_α = 0 (black), α=0.1𝛼0.1\alpha=0.1italic_α = 0.1 (blue), α=0.3𝛼0.3\alpha=0.3italic_α = 0.3 (green), α=0.5𝛼0.5\alpha=0.5italic_α = 0.5 (red).
Refer to caption
Figure 10: The capture cross-section σcapsubscript𝜎cap\sigma_{\text{cap}}italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT vs. {\cal E}caligraphic_E for a charged particle with α=0𝛼0\alpha=0italic_α = 0 (black), α=1𝛼1\alpha=1italic_α = 1 (blue), α=2𝛼2\alpha=2italic_α = 2 (green), α=3𝛼3\alpha=3italic_α = 3 (red).

For ultra-relativistic particles, we can write bmaxsubscript𝑏maxb_{\text{max}}italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT as

bmax=3323α+92α2632+𝒪(13).subscript𝑏max3323𝛼92superscript𝛼263superscript2𝒪1superscript3b_{\text{max}}=\frac{3\sqrt{3}}{2}-\frac{\sqrt{3}\alpha}{{\cal E}}+\frac{9-2% \alpha^{2}}{6\sqrt{3}{\cal E}^{2}}+{\mathcal{O}}\left(\frac{1}{{\cal E}^{3}}% \right).italic_b start_POSTSUBSCRIPT max end_POSTSUBSCRIPT = divide start_ARG 3 square-root start_ARG 3 end_ARG end_ARG start_ARG 2 end_ARG - divide start_ARG square-root start_ARG 3 end_ARG italic_α end_ARG start_ARG caligraphic_E end_ARG + divide start_ARG 9 - 2 italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 6 square-root start_ARG 3 end_ARG caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG caligraphic_E start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ) . (36)

The corresponding capture cross-section is then

σcap=27π49πα+(4α2+9)π22+𝒪(13).subscript𝜎cap27𝜋49𝜋𝛼4superscript𝛼29𝜋2superscript2𝒪1superscript3\sigma_{\text{cap}}=\frac{27\pi}{4}-\frac{9\pi\alpha}{{\cal E}}+\frac{\left(4% \alpha^{2}+9\right)\pi}{2{\cal E}^{2}}+{\mathcal{O}}\left(\frac{1}{{\cal E}^{3% }}\right).italic_σ start_POSTSUBSCRIPT cap end_POSTSUBSCRIPT = divide start_ARG 27 italic_π end_ARG start_ARG 4 end_ARG - divide start_ARG 9 italic_π italic_α end_ARG start_ARG caligraphic_E end_ARG + divide start_ARG ( 4 italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 9 ) italic_π end_ARG start_ARG 2 caligraphic_E start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + caligraphic_O ( divide start_ARG 1 end_ARG start_ARG caligraphic_E start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ) . (37)

These limiting results are in agreement with our numerical findings.

5 Conclusion

We have studied the capture cross-section of charged particles by a weakly charged Schwarzschild black hole. We have shown that a trace charge on the black hole can have prominent effects.

When the Coulomb force between a charged particle and the black hole is attractive, it enlarges the capture cross-section significantly. This is expected since the Coulomb attraction enhances the capture of charged particles. However, when the Coulomb force between a charged particle and the black hole is repulsive, it shrinks the capture cross-section significantly. When the electromagnetic coupling strength is below a critical value, capture is possible for all values of the particle’s energy. When the electromagnetic coupling strength is above the critical value, there is a minimum value of the particle’s energy below which capture is impossible. This is because the Coulomb repulsion surpasses the gravitational attraction unless the particle’s radial momentum is large enough.

Our results emphasizes the assertion that charged black holes will favorably accretes charges of the opposite sign. However, it is still possible for the black hole charge to grow if the plunging charged particles are energetic enough to the limit that the capture cross-section becomes independent of the sign of the charges. Moreover, the fact that the electromagnetic coupling constant is three orders of magnitudes greater for electrons than protons suggests that it is relatively easier for a black hole to accumulate positive charge than negative charge.

It will be an astrophysically interesting to study the energies of charged particles near an astrophysical black hole to understand better how the black hole’s charge evolves. The problem can be astrophyically more viable when other astrophysical black holes, such as rotating black holes, are studied (in progress).

References

  • Ahmedov et al. (2021) Ahmedov B., Rahimov O., Toshmatov B., 2021, Universe, 7(8) , 307.
  • Al Zahrani (2021) Al Zahrani A., 2021, Phys. Rev. D, 103.
  • Al Zahrani (2022) Al Zahrani A., 2022, ApJ, 937.
  • Anacleto, et al. (2023) Anacleto M., et al., 2023, arXiv:2307.09536v1 [gr-qc].
  • Bugden (2020) Bugden M., 2020, Class. Quantum Gravity, 37, 015001.
  • Carter (1973) Carter B., 1973, Black Hole Equilibrium States, Black Holes, eds. C. DeWitt and B. S. DeWitt (Gordon and Breach Science Publishers, Inc. New York, p. 57.
  • Chandrasekhar (1983) Chandrasekhar S., 1983, The Mathematical Theory of Black Holes, Oxford University Press.
  • Connell & Frolov (2008) Connell P., Frolov V., 2008, Phys. Rev. D, 78, 024032.
  • Frolov & Zelnikov (2011) Frolov V., Zelnikov A., 2011, Introduction to black hole physics, Oxford university press.
  • Ghosh & Afrin (2023) Ghosh S. and Afrin M., 2023, ApJ, 944, 174.
  • Goldstein et al. (2001) Goldstein H., Poole C. and Safko J., 2001, Classical Mechanics, third eddition, Pearson.
  • Gooding & Frolov (2008) Gooding C., Frolov A., 2008, Phys. Rev. D, 77, 104026.
  • GRAVITY Collaboration (2023) GRAVITY Collaboration, 2023, A & A, 677, L10.
  • Kocherlakota et al. (2021) Kocherlakota P. et al., 2021, (EHT Collaboration), Phys. Rev. D, 103, 104047.
  • Misner et al. (1973) Misner C., Thorne K., Wheeler J., 1973, Gravitation, W. H. Freeman and Co., San Francisco.
  • Singh & Ghosh (2018) Singh B. , Ghosh S., 2018, Annals of Physics 395, 127.
  • Tsukamoto (2014) Tsukamoto N., et al., 2014, Phys. Rev. D, 90, 064043.
  • Young (1976) Young P., 1976, Phys. Rev. D, 14, 3281.
  • Zajaček et al. (2018) Zajaček M. et al., 2018, Monthly Notices of the Royal Astronomical Society, 480 4, 4408.
  • Zajaček et al. (2019) Zajaček M. et al., 2019, J. Phys.: Conf. Ser., 1258, 012031.
  • Zajaček & Tursunov (2019) Zajaček M., Tursunov A., arXiv:1904.04654.
  • Zakharov (1994) Zakharov A., 1994, Class. Quantum Grav., 11, 1027.