A closed set of fluid equations with radiation reaction force (RRF) are constructed from the moments of the appropriate single particle kinetic equation describing a relativistic degenerate (high density) electron plasma. The closure, in analogy with the Maxwellian closure for non-degenerate plasmas, is affected via a parametrized Fermi-Dirac distribution. It is shown that the degeneracy increases RRF just as will be predicted from the so-called “Compton Rocket” effect.

A relativistically hot plasma, exposed to intense radiation, experiences a higher radiation pressure force compared to a cold plasma. The phenomenon known as the “Compton Rocket” was originally discovered by O'Dell.^{1} In O'Dell's work, the force acting on the plasma was derived using a phenomenological, test-particle approach. Energy-momentum conservation is employed to address particle–photon interaction, and the resulting force is integrated over the distribution function. It was shown that incident high-frequency photon flux, entering a hot plasma, undergoes Thomson scattering that redirects it back toward the primary source. The plasma exhibits a tendency to propel itself away from the radiation source, drawing momentum primarily from the anisotropic loss of its internal energy. Applying similar but a relativistic covariant formalism, it was demonstrated in Refs. 2–4 that this process is accompanied by plasma cooling. Note that the force acting on the plasma fluid may not only accelerate the plasma but also induce deceleration at velocities surpassing certain (model-dependent) bulk flow values. This deceleration occurs due to the photon flux aberration, commonly known as the radiative drag. Methods and results obtained in Refs. 1–4 have been employed to study the dynamics of producing relativistic outflows (jets) from highly luminous radiation sources, such as active galactic nuclei (AGNs) or compact galactic objects.^{5}

In this paper, we explore how the radiation reaction force (RRF) would be modified, likely, intensified by the “Compton Rocket” effect. Not considered in the original calculation, the RRF, by increasing the degeneracy and effectively raising the electron degeneracy “temperature” (even in a thermally cold plasma), might be a stimulant for a stronger effect.

The basic model we investigate is contained in the manifestly covariant relativistic hydrodynamic equations that were derived taking into account the radiation reaction force (RRF) for an electron plasma.^{6,7} The model equations are simply the moments of the relativistic kinetic equation with RRF in the Landau–Lifshitz (LL) form. In the limit where the plasma is immersed in the field of incoherent photon fluxes, the derived equations can be readily reduced to those presented in Refs. 1–4.

It is interesting that in Ref. 7, we did demonstrate that RRF, when radiation consists of coherent ultra-strong electromagnetic (EM) pulses, can lead to an effective acceleration of electron (electron–positron) plasma. Since the relativistic hot plasmas have excess internal energy, the radiative thrust of RRF increases just as in “Compton Rocket” effect.

Before getting into detailed calculations, let us review the phenomenology of the compact astrophysical objects—neutron stars, white dwarfs, and AGN—in which the electron number density is believed to be in the range from 10^{26} to $ 10 34 \u2009 cm \u2212 3$. The average inter particle distance can be considerably smaller than the electron thermal de Broglie wavelength. The electron gas in this case is degenerate, must obey Fermi-Dirac statistics, and plasma particles are weakly coupled.^{8,9}

The electrons in a high density degenerate plasma can have relativistic energies, even when they are thermally cold. The Fermi energy of degenerate electrons is $ \u03f5 F = m e c 2 ( \gamma F \u2212 1 )$ with $ \gamma F = 1 + p F 2 / m e 2 c 2$, where *p _{F}* is the Fermi momentum determined by the rest-frame electron density

*n*via $ p F = m e c ( n / n c ) 1 / 3$,

^{10}where $ n c = m e 3 c 3 / 3 \pi 2 \u210f 3 = 5.9 \xd7 10 29 \u2009 cm \u2212 3$ is the normalizing critical number-density. If $ n \u2265 n c$, the electrons acquire non-thermal energy that could become relativistic when

*n*gets larger. When electrons have relativistic speeds, the power radiated due to acceleration increases dramatically. Relativistic electrons emit radiation anisotropically effectively enhancing the radiation reaction force acting on the plasma.

It has been widely accepted that x-ray emission might appear from accreting white dwarfs and pulsars. Details of this process are discussed in the review article of Mukai^{11} and references therein. The accretion matter, falling onto the star's surface, generates high-frequency radiation via Bremsstrahlung. This radiation penetrates the star's interior, which is mostly made up of very dense degenerate plasma. Consequently, the investigation of relativistic degenerate plasma (present in the interiors of compact objects) immersed in strong radiation, is a subject of considerable astrophysical interest—it could, in fact, be an excellent laboratory for the manifestation of the Compton Rocket effect.

*α*th component of the contravariant reduced four-momentum $ u \alpha = [ \gamma , \gamma u / c ] , \u2009 \gamma = ( 1 \u2212 u 2 / c 2 ) \u2212 1 / 2$, where

**u**is the particle velocity, $ d s = cdt / \gamma , \u2009 F \alpha \beta $ is the electromagnetic tensor [ $ u \alpha u \alpha = 1$ since the metric is $ g \alpha \beta ( 1 , \u2212 1 , \u2212 1 , \u2212 1 )$], and $ g \alpha $ is the contravariant component of the radiation reaction force. The following expression for $ g \alpha $ was first derived, perturbatively, by Landau–Lifshitz;

^{8}it was, however, later shown by Spohn

^{12}and Rohrlich

^{13}to be exact for a single point particle. Following Ref. 6, LL force is written conveniently as

The LL equations [Eq. (1)] offer several advantages over the Lorentz–Abraham–Dirac formulation (LAD).^{14} Specifically, in contrast to LAD, these equations are of the second rather than third order, and they do not permit a self-force in the absence of an external field. Additionally, as demonstrated in Ref. 15 (see also Refs. 7 and 16), Eq. (1) possesses exact solutions for plane electromagnetic fields. We would like to emphasize that from expression (2), one can easily recover the recoil force acting on electrons during Thomson scattering of photons, provided that the photon energy in the particle rest frame ( $ \u210f \omega $) is much less than the electron's rest mass energy ( $ m c 2$).

^{17}Assuming the validity of the weak coupling condition among plasma particles, the kinetic collisionless equation for a single-particle electron (or positron) distribution function, influenced by RRF as follows: evolves as

*g*= 2 for electrons (positrons)], the chemical potential

*μ*, the temperature

*T*(measured in energy units), and the hydrodynamic four velocity $ U \alpha $: $ U \alpha = ( c \gamma v , \gamma v \u2009 V ) , \u2009 \gamma v = ( 1 \u2212 V 2 / c 2 ) \u2212 1 / 2 ( U \alpha U \alpha = c 2 )$. This distribution yields the flux four vector: $ \Gamma \alpha = n U \alpha $, while for energy momentum tensor of fluid, we have

*p*is the pressure. As shown by Cercignani and Kremer,

^{18}for these parameters, the distribution (11) yields the following general relations:

*A*

_{1}and $ A 2$ are found to be

^{19}In this limit, the integral (16) simplifies to:

^{18},

*A*

_{1}and

*A*

_{2}, we get

*A*

_{1}and

*A*

_{2}into Eq. (21), we obtain

**E**and

**B**are the generated low frequency electric and magnetic fields, respectively. Note that in the right hand side of (33), the index

*α*pertains only to spatial components.

In weakly degenerate case where $ x = p F / m c \u226a 1$, the expression (33) for RRF simplifies to the one obtained for a classical Maxwell distribution in the nonrelativistic regime ( $ T \u226a m c 2$). Since $ p F = h ( 3 n / 8 \pi ) 1 / 3$, the condition of weak degeneracy can be written in terms of electron density $ n \u226a n c = 5.9 \xd7 10 29 \u2009 cm \u2212 3$. However, it is noteworthy that the gas (fluid) approximation for plasma remains valid only if the average potential energy of particles is less than the Fermi energy, implying that the density falls within the solid density range $ 10 23 \u2009 cm \u2212 3 \u226a n \u226a 5.9 \xd7 10 29 \u2009 cm \u2212 3$.

For a relativistic degenerate plasma ( $ n \u2265 5.9 \xd7 10 29 \u2009 cm \u2212 3$), with Fermi momentum greater than the rest mass, i.e., $ x 2 \u2265 1$ (i.e., $ p F \u2265 m c )$, the RRF is strongly boosted up, scaling quadratically with *x*; it is a most spectacular manifestation of the Compton Rocket effect.

Just for perspective, we remind the reader that in a relativistic dilute (non-degenerate) plasma, the enhancement of RRF would occur only when the kinematic temperature is relativistic ( $ T \u2265 m c 2$). However, it would be nowhere near as strong as in the strongly degenerate case.

*U*, becomes

_{rad}We would like to emphasize that the model, described above, is valid in the Thomson limit. The dynamics of the Compton scattering regime, pertaining a cold classical plasma embedded in the x-ray photon fluxes, have been studied for some time.^{20–22} More recently, the importance of photon plasma interaction in the Compton scattering regime was demonstrated^{23,24} for a variety of astrophysical conditions. To extend our work to the Compton regime, one should include quantum effects in radiation reaction force as well as introduce photon particle collisions in the relativistic kinetic equation. Such a work will constitute a much more involved study, which is beyond the scope of the current letter.

By invoking a local Fermi-Dirac (thermodynamic) closure, we derived hydrodynamic equations that include the radiation reaction force, RRF. These equations are valid for arbitrary (physically allowed) temperature and densities of electron (positron) plasmas. It is shown that RRF is strongly amplified in strongly degenerate plasmas—it is, perhaps, one of the strongest expressions of the Compton Rocket effect in a well-defined physical system.

Working out the consequences of this phenomenon for compact astrophysical objects is under way.

This research was supported by the Shota Rustaveli National Science Foundation of Georgia under Grant No. FR22-8273. The work of S.M.M. was supported by USDOE under Contract No. DE-FG 03–96ER-54366.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**V. I. Berezhiani:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). **S. M. Mahajan:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal).

## DATA AVAILABILITY

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

## REFERENCES

*Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects*

*Quantum Plasmas: An Hydrodynamic Approach*

*The Relativistic Boltzmann Equation: Theory and Applications*