pith. machine review for the scientific record. sign in

arxiv: 2604.09912 · v1 · submitted 2026-04-10 · 🌀 gr-qc

Recognition: unknown

Relativistic figures of equilibrium in the Wald magnetosphere

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:28 UTC · model grok-4.3

classification 🌀 gr-qc
keywords rotating fluidsWald magnetosphereEinstein-Euler equationscharged perfect fluidgeneral relativitynumerical relativitypolytropic equation of state
0
0 comments X

The pith

A rigidly rotating charged perfect fluid with vanishing conductivity can remain compatible with Wald's magnetosphere in non-vacuum spacetimes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that Wald's electromagnetic solution, built from Killing vectors, stays consistent with the current carried by a rotating charged fluid provided electric conductivity is zero. For two common equations of state—constant energy density and polytropic—the conservation law for the energy-momentum tensor integrates, reducing the problem to the usual Einstein-Euler system except for one change in the Euler-Bernoulli equation. Existing pseudospectral codes for rotating fluid stars can therefore be adapted with small modifications to produce numerical models. A reader would care because the result supplies concrete, computable examples of self-gravitating objects that carry both rotation and magnetic fields without invoking full magnetohydrodynamics. The authors demonstrate solutions for both equations of state.

Core claim

In non-vacuum spacetimes, Wald's vacuum magnetic solution remains compatible with the electric current associated with a rotating charged perfect fluid of vanishing electric conductivity. For rigidly rotating fluids with constant energy density or a polytropic equation of state, the conservation of the energy-momentum tensor can be integrated, so the system is governed by nearly standard Einstein-Euler equations with a modification only in the Euler-Bernoulli equation.

What carries the argument

The Wald magnetosphere, an electromagnetic field constructed from two linearly independent Killing vectors in a stationary axisymmetric spacetime, which supplies the background magnetic field that matches the fluid current when conductivity vanishes.

If this is right

  • The Einstein-Euler system for rigidly rotating fluids acquires only a single extra term in the Bernoulli integral when the Wald field is present.
  • Numerical equilibria for constant-density and polytropic cases can be obtained by minor changes to an existing pseudospectral code.
  • Self-gravitating, rigidly rotating charged fluids can be modeled without solving the full coupled Maxwell-Einstein system.
  • The construction works for any stationary axisymmetric spacetime that admits the required Killing vectors.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same integration trick may apply to other simple equations of state, allowing quick extension of the numerical models.
  • Stability calculations for these magnetized configurations could reuse techniques already developed for unmagnetized rotating stars.
  • The models supply initial data for dynamical simulations of rotating compact objects that carry a prescribed magnetic field.

Load-bearing premise

Wald's vacuum magnetic field remains consistent with the electric current of the rotating fluid once the spacetime is allowed to be non-vacuum.

What would settle it

A direct integration of Maxwell's equations with the four-current derived from the rotating charged fluid that yields an electromagnetic field differing from the Wald solution.

Figures

Figures reproduced from arXiv: 2604.09912 by Audrey Trova, Bakhtinur Juraev, Patryk Mach, Pawe{\l} Doruchowski.

Figure 1
Figure 1. Figure 1: FIG. 1: Variation of the surface of the star with respect to [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Distribution of the normalized pressure ¯p [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Magnetic (left) and electric (right) field lines, as seen by ZAMO observers, for the case [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The relative deviation of the mass corresponding to Tab. [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Left: Variation of the surface of the star with respect to [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
read the original abstract

We consider a self-gravitating, rigidly rotating charged perfect fluid immersed in the Wald magnetosphere, constructed out of two linearly independent Killing vectors present in stationary and axially-symmetric spacetimes. We show that in non-vacuum spacetimes, Wald's solution can be compatible with the electric current associated with a rotating charged perfect fluid characterized by the vanishing electric conductivity. We prove that for rigidly rotating fluids with a constant energy density or described by the polytropic equation of state, the resulting equations expressing the conservation of the energy-momentum tensor can be integrated. Consequently, the system can be described by nearly standard Einstein-Euler equations known from the theory of general-relativistic rotating fluids, with modifications introduced in the Euler-Bernoulli equation. Numerical solutions of the Einstein-Euler equations are provided for these two cases by introducing suitable modifications in the pseudospectral code by Ansorg, Kleinw\"{a}chter, and Meinel.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper considers self-gravitating, rigidly rotating charged perfect fluids immersed in the Wald magnetosphere constructed from the two Killing vectors of a stationary axisymmetric spacetime. It shows that Wald's vacuum electromagnetic field remains compatible with the fluid's electric current when electric conductivity vanishes, and proves that for constant energy density or polytropic equations of state the conservation of the energy-momentum tensor integrates to a modified Euler-Bernoulli relation. This reduces the system to the standard Einstein-Euler equations with only that relation altered, and the authors supply numerical equilibria by inserting the modified relation into an existing pseudospectral code.

Significance. If the integrability and numerical implementation hold, the work supplies a concrete, reproducible extension of the theory of relativistic rotating stars to include a Wald-type magnetic field while preserving the structure of the Einstein-Euler system for two standard equations of state. The reuse and modification of the established Ansorg-Kleinwächter-Meinel pseudospectral solver is a strength that supports direct comparison with non-magnetized equilibria.

major comments (2)
  1. [Numerical solutions] Numerical solutions section: the abstract and text assert that solutions exist for the two equations of state, yet no information is given on grid resolution, convergence tests, residual norms, or comparison against the non-magnetized limit (e.g., the known constant-density or polytropic sequences without the Wald field). This omission prevents assessment of whether the reported equilibria are accurate to the claimed precision.
  2. [Derivation of the integrated Euler-Bernoulli equation] The compatibility argument relies on the current four-vector being exactly parallel to the fluid four-velocity when conductivity vanishes; however, the manuscript does not explicitly verify that the resulting Maxwell equations remain satisfied once the fluid back-reacts on the metric (i.e., in the non-vacuum spacetime). A concrete check against the sourced Maxwell equations in the numerical solutions would strengthen the central claim.
minor comments (2)
  1. [Integration of the conservation laws] The notation for the modified Bernoulli integral should be introduced with an explicit equation number immediately after the integration step, rather than being referred to only descriptively.
  2. [Numerical implementation] A brief statement of the boundary conditions imposed at the stellar surface and at spatial infinity would clarify how the Wald field is matched to the interior solution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and insightful comments on our manuscript. We address each major comment below and have prepared revisions to strengthen the presentation of our results.

read point-by-point responses
  1. Referee: Numerical solutions section: the abstract and text assert that solutions exist for the two equations of state, yet no information is given on grid resolution, convergence tests, residual norms, or comparison against the non-magnetized limit (e.g., the known constant-density or polytropic sequences without the Wald field). This omission prevents assessment of whether the reported equilibria are accurate to the claimed precision.

    Authors: We acknowledge that the numerical section would benefit from more detailed information on the computational parameters and validation. In the revised manuscript, we will include the grid resolutions employed in the pseudospectral code, results from convergence tests demonstrating the expected spectral accuracy, the norms of the residuals for the Einstein-Euler equations, and direct comparisons of key quantities (such as mass, angular momentum, and magnetic field strength) with the corresponding non-magnetized equilibria from the literature. These additions will allow readers to assess the accuracy of the reported solutions. revision: yes

  2. Referee: The compatibility argument relies on the current four-vector being exactly parallel to the fluid four-velocity when conductivity vanishes; however, the manuscript does not explicitly verify that the resulting Maxwell equations remain satisfied once the fluid back-reacts on the metric (i.e., in the non-vacuum spacetime). A concrete check against the sourced Maxwell equations in the numerical solutions would strengthen the central claim.

    Authors: The analytical derivation demonstrates that the Wald electromagnetic field, defined using the two Killing vectors, remains compatible with the fluid current (which is parallel to the four-velocity for vanishing conductivity) even in the presence of the fluid's energy-momentum tensor. This compatibility follows from the properties of the Killing vectors and the form of the Maxwell tensor, independent of the specific metric solution. Nevertheless, to address the referee's concern, we will add an explicit verification in the numerical solutions section, showing that the sourced Maxwell equations are satisfied to the level of the numerical truncation error in our computed equilibria. This will be done by computing the divergence of the electromagnetic field tensor and comparing it to the current four-vector. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from standard GR plus Wald construction

full rationale

The paper begins with the Einstein equations, the Wald electromagnetic field built from the two Killing vectors of a stationary axisymmetric spacetime, and the assumption of vanishing conductivity. It then demonstrates that the resulting current is parallel to the fluid four-velocity, shows that the electromagnetic term in the Euler equation becomes a gradient for constant-density and polytropic EOS, integrates to a modified Bernoulli relation, and solves the resulting system numerically by adapting an existing pseudospectral code. None of these steps reduces to a self-definition, a fitted parameter renamed as a prediction, or a load-bearing self-citation. The central compatibility and integrability results are direct mathematical consequences of the stated assumptions and do not presuppose the final equations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard general-relativistic assumptions for stationary axisymmetric spacetimes plus the specific choice of vanishing conductivity to keep the magnetic field solution valid.

axioms (2)
  • domain assumption Spacetime is stationary and axially symmetric, admitting two linearly independent Killing vectors
    Required to construct the Wald magnetosphere from the Killing vectors.
  • ad hoc to paper Electric conductivity of the fluid vanishes
    Invoked to ensure compatibility between Wald's vacuum solution and the fluid's electric current in non-vacuum regions.

pith-pipeline@v0.9.0 · 5471 in / 1215 out tokens · 31964 ms · 2026-05-10T16:28:08.506768+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 35 canonical work pages · 1 internal anchor

  1. [1]

    They can be computed directly from Eq

    = 1 (this corresponds top= 0). They can be computed directly from Eq. (44), or from the explicit solution (45). We set C= 1. a= 0,u t 0 = 1a= 0.1,u t 0 = 0.965166a= 0.2,u t 0 = 0.833049 H(0)(ut

  2. [2]

    1 0.991602 0.923646 H(2)(ut

  3. [3]

    0−0.0401569−0.102436 H(3)(ut

  4. [4]

    0 0.308907 1.4975 H(4)(ut

  5. [5]

    0−2.08438−13.8529 H(5)(ut

  6. [6]

    0 15.2849 136.129 H(6)(ut

  7. [7]

    0−125.617−1506.29 H(7)(ut

  8. [8]

    0 1158.14 18845.7 H(8)(ut

  9. [9]

    0−11906.4−264557. H(9)(ut

  10. [10]

    0 135467.4.12572×10 6 H(10)(ut

  11. [11]

    0−1.6931×10 6 −7.07924×10 7 and p ρ0 = p ε =h−1.(43) Inserting this into Eq. (41), we obtain the equation H′ − H 1 ut + 4a2(3H −2) 1 (ut)3 = 0,(44) which, upon integration, yields H(ut) =− 1 18a " √ 6πut erf √ 6a ut ! exp 6a2 (ut)2 −12a # +Cu t exp 6a2 (ut)2 ,(45) whereCis the integration constant. Using the property lim x→0(erf(x)/x) = 2/√π, it can be ch...

  12. [12]

    :=∂ nH(ut 0)/∂(ut)n, whereu t 0 satisfiesH(u t

  13. [13]

    This choice corresponds to the vanishing pressure,p= 0

    = 1. This choice corresponds to the vanishing pressure,p= 0. In the examples collected in Tab. I, we setC= 1 and choose the values ofa= 0.1 anda= 0.2. In contrast to this behavior, fora= 0 one has the limitH(u t) =Cu t (Eq. (46)). In this caseH ′(ut) =C, and all higher derivatives vanish. B. Polytropic fluids For the polytropic equation of state,p=Kρ Γ 0 ...

  14. [14]

    0.98495 0.905376 0.983643 0.901805 H(2)(ut

  15. [15]

    0.00378691 0.116285 0.0122771 0.153011 H(3)(ut

  16. [16]

    0.0722495−0.0887735 0.0273865−0.326636 H(4)(ut 0)−0.648803−0.977188−0.382514 0.735717 H(5)(ut

  17. [17]

    5.20721 13.9595 3.38265 0.388689 H(6)(ut 0)−44.2599−161.221−29.9133−21.9982 H(7)(ut

  18. [18]

    411.595 1921.32 283.591 376.375 H(8)(ut 0)−4210.37−24781.9−2928.97−5436.9 H(9)(ut

  19. [19]

    Support for Talented Doctoral Students at the Silesian University in Opava 2022

    47268.1 349591.33021.2 79576.8 H(10)(ut 0)−579707−5.39059×10 6 −405334−1.23344×10 6 and consequently p ρ0 = Γ−1 Γ (h−1).(49) Therefore, condition (41) becomes H′ − H 1 ut + 4a2 H+ 2 Γ−1 Γ (H −1) 1 (ut)3 = 0.(50) The above equation can be integrated as H(ut) =− (Γ−1) √ 2πΓut erf a √ 2(3Γ−2)√ Γut exp 2a2(3Γ−2) Γ(ut)2 −4a √ 3Γ−2 2a(3Γ−2) 3/2 +Cu t exp 2a2(3Γ...

  20. [20]

    R. M. Wald, Black hole in a uniform magnetic field, Phys. Rev. D 10, 1680 (1974). doi: 10.1103/PhysRevD.10.1680. 19

  21. [21]

    Ben Achour, I

    J. Ben Achour, I. El Mellah, and E. Gourgoulhon, Jet launching from the Kerr black hole magnetosphere: An electrogeodesic approach, arXiv:2601.05048 (2026). doi: 10.48550/arXiv.2601.05048

  22. [22]

    Meinel, M

    R. Meinel, M. Ansorg, A. Kleinw¨ achter, G. Neugebauer, and D. Petroff, Relativistic Figures of Equilibrium, Cambridge University Press, Cambridge 2008. doi: 10.1017/CBO9780511535154

  23. [23]

    Ansorg, A

    M. Ansorg, A. Kleinw¨ achter, and R. Meinel, Highly accurate calculation of rotating neutron stars: Detailed description of the numerical methods, Astron. Astrophys. 405, 711 (2003). doi: 10.1051/0004-6361:20030618

  24. [24]

    J. D. Bekenstein and E. Oron, New conservation laws in general-relativistic magnetohydrodynamics, Phys. Rev. D 18, 1809 (1978). doi: 10.1103/PhysRevD.18.1809

  25. [25]

    J. D. Bekenstein and E. Oron, Interior magnetohydrodynamic structure of a rotating relativistic star, Phys. Rev. D 19, 2827 (1979). doi: 10.1103/PhysRevD.19.2827

  26. [26]

    ROTATING NEUTRON STAR MODELS WITH MAGNETIC FIELD

    M. Bocquet, S. Bonazzola, E. Gourgoulhon, and J. Novak, Rotating neutron star models with a magnetic field, Astron. Astrophys. 301, 757 (1995). doi: 10.48550/arXiv.gr-qc/9503044

  27. [27]

    Konno, T

    K. Konno, T. Obata, and Y. Kojima, Deformation of relativistic magnetized stars, Astron. Astrophys. 352, 211 (1999). doi: 10.48550/arXiv.gr-qc/9910038

  28. [28]

    Ch. Y. Cardall, M. Prakash, and J. M. Lattimer, Effects of strong magnetic fields on neutron star structure, Astrophys. J. 554, 322 (2001). doi: 10.1086/321370

  29. [29]

    Bonazzola, E

    S. Bonazzola, E. Gourgoulhon, M. Salgado, and J. A. Marck, Axisymmetric rotating relativistic bodies: A new numerical approach for ‘exact’ solutions, Astron. Astrophys. 278, 421 (1993)

  30. [30]

    2012, MNRAS, 427, 127, doi: 10.1111/j.1365-2966.2012.21948.x

    J. Frieben and L. Rezzolla, Equilibrium models of relativistic magnetized neutron stars: Axisymmetric solutions with poloidal and toroidal magnetic fields, Mon. Not. R. Astron. Soc. 427, 3406 (2012). doi: 10.1111/j.1365-2966.2012.22027.x

  31. [31]

    Kiuchi and S

    K. Kiuchi and S. Yoshida, Relativistic stars with purely toroidal magnetic fields, Phys. Rev. D 78, 044045 (2008). doi: 10.1103/PhysRevD.78.044045

  32. [32]

    2010, MNRAS, 401, 1670, doi: 10.1111/j.1365-2966.2009.15794.x

    R. Ciolfi, V. Ferrari, L. Gualtieri, and J. A. Pons, Relativistic models of magnetars: The twisted torus magnetic field configuration, Mon. Not. R. Astron. Soc. 397, 913 (2009). doi: 10.1111/j.1365-2966.2009.14990.x

  33. [33]

    Ciolfi and L

    R. Ciolfi and L. Rezzolla, Twisted-torus configurations with large toroidal magnetic fields in relativistic stars, Mon. Not. R. Astron. Soc. Lett. 435, L43 (2013). doi: 10.1093/mnrasl/slt092

  34. [34]

    Ioka and M

    K. Ioka and M. Sasaki, Relativistic stars with poloidal and toroidal magnetic fields and meridional flow, Astrophys. J. 600, 296 (2004). doi: 10.1086/379650

  35. [35]

    , keywords =

    O. Zanotti and L. Rezzolla, General relativistic electromagnetic fields of a slowly rotating magnetized neutron star — II. Solution of the induction equations, Mon. Not. R. Astron. Soc. 331, 376 (2002). doi: 10.1046/j.1365-8711.2002.05196.x

  36. [36]

    Ioka and M

    K. Ioka and M. Sasaki, Grad–Shafranov equation in non-circular stationary axisymmetric spacetimes, Phys. Rev. D 67, 124026 (2003). doi: 10.1103/PhysRevD.67.124026

  37. [37]

    Gourgoulhon, Ch

    E. Gourgoulhon, Ch. Markakis, K. Uryu, and Y. Eriguchi, Magnetohydrodynamics in stationary and axisymmetric space- times: A fully covariant approach, Phys. Rev. D 83, 104007 (2011). doi: 10.1103/PhysRevD.83.104007

  38. [38]

    K. Uryu, E. Gourgoulhon, Ch. M. Markakis, K. Fujisawa, A. Tsokaros, and Y. Eriguchi, Equilibrium solutions of relativis- tic rotating stars with mixed poloidal and toroidal magnetic fields, Phys. Rev. D 90, 101501 (2014). doi: 10.1103/Phys- RevD.90.101501

  39. [39]

    2010, MNRAS, 401, 1670, doi: 10.1111/j.1365-2966.2009.15794.x

    J. Otani, R. Takahashi, and Y. Eriguchi, Equilibrium states of magnetized toroid–central compact object systems, Mon. Not. R. Astron. Soc. 396, 2152 (2009). doi: 10.1111/j.1365-2966.2009.14842.x

  40. [40]

    Zanotti and D

    O. Zanotti and D. Pugliese, Von Zeipel’s theorem for a magnetized circular flow around a compact object, Gen. Relativ. Gravit. 47, 44 (2015). doi: 10.1007/s10714-015-1886-4

  41. [41]

    P. Mach, S. Gimeno-Soler, J. A. Font, A. Odrzywo lek, and M. Pir´ og, Self-gravitating magnetized tori around black holes in general relativity, Phys. Rev. D 99, 104063 (2019). doi: 10.1103/PhysRevD.99.104063

  42. [42]

    W. Dyba, P. Mach, and M. Pietrzy´ nski, Toroidal magnetic fields in self-gravitating disks around black holes, Phys. Rev. D 104, 044058 (2021). doi: 10.1103/PhysRevD.104.044058

  43. [43]

    Shibata, Rotating black hole surrounded by self-gravitating torus in the puncture framework, Phys

    M. Shibata, Rotating black hole surrounded by self-gravitating torus in the puncture framework, Phys. Rev. D 76, 064035 (2007). doi: 10.1103/PhysRevD.76.064035

  44. [44]

    Karkowski, W

    J. Karkowski, W. Kulczycki, P. Mach, E. Malec, A. Odrzywo lek, and M. Pir´ og, General-relativistic rotation: Self-gravitating fluid tori in Keplerian motion around black holes, Phys. Rev. D 97, 104034 (2018). doi: 10.1103/PhysRevD.97.104034

  45. [45]

    Karkowski, W

    J. Karkowski, W. Kulczycki, P. Mach, E. Malec, A. Odrzywo lek, and M. Pir´ og, Self-gravitating axially symmetric disks in general-relativistic rotation, Phys. Rev. D 97, 104017 (2018). doi: 10.1103/PhysRevD.97.104017

  46. [46]

    Ch. G. Tsagas, Electromagnetic fields in curved spacetimes, Class. Quantum Grav. 22, 393 (2005). doi: 10.1088/0264- 9381/22/2/011

  47. [47]

    Kov´ aˇ r, P

    J. Kov´ aˇ r, P. Slan´ y, Z. Stuchl´ ık, V. Karas, C. Cremaschini, and J. C. Miller, Role of electric charge in shaping equilibrium configurations of fluid tori encircling black holes, Phys. Rev. D 84, 084002 (2011). doi: 10.1103/PhysRevD.84.084002

  48. [48]

    Trova, Equilibrium of charged fluid around a Kerr black hole immersed in a magnetic field: Variation of angular momentum, Int

    A. Trova, Equilibrium of charged fluid around a Kerr black hole immersed in a magnetic field: Variation of angular momentum, Int. J. Mod. Phys. D 33, 2350112 (2024). doi: 10.1142/S0218271823501122

  49. [49]

    Trova, E

    A. Trova, E. Hackmann, V. Karas, K. Schroven, J. Kov´ aˇ r, and P. Slan´ y, Influence of test charge and uniform magnetic field on charged fluid equilibrium structures, Phys. Rev. D 101, 083027 (2020). doi: 10.1103/PhysRevD.101.083027

  50. [50]

    L. I. Petrich, S. L. Shapiro, and S. A. Teukolsky, Accretion onto a moving black hole: An exact solution, Phys. Rev. Lett. 60, 1781 (1988). doi: 10.1103/PhysRevLett.60.1781

  51. [51]

    Straumann, General Relativity (2nd edition), Springer, Dordrecht 2013

    N. Straumann, General Relativity (2nd edition), Springer, Dordrecht 2013. doi: 10.1007/978-94-007-5410-2

  52. [52]

    R. M. Wald, General Relativity, The University of Chicago Press, Chicago 1984. doi: 10.7208/chicago/9780226870373.001.0001. 20

  53. [53]

    Ansorg and D

    M. Ansorg and D. Petroff, Black holes surrounded by uniformly rotating rings, Phys. Rev. D 72, 024019 (2005). doi: 10.1103/PhysRevD.72.024019

  54. [54]

    J. M. Bardeen and R. V. Wagoner, Relativistic disks. I. Uniform rotation, Astrophys. J. 167, 359 (1971). doi: 10.1086/151039