A D V E R T |
I S E M E N T |
| Math Sites: lean ages 13-∞ readme referees parents | >>> MAITH GAMES <<< | all 372 compute stand |
|
LEVEL 5 OF 13 · Spacetime Penrose inequalities and rigidity
Electromagnetic tails and the Kerr–Newman Penrose inequality
expertly designed by an internal OpenAI model · released 2026-10-05
· original PDF
IntroductionThe Kerr–Newman mass formula suggests the Kerr–Newman Penrose inequality \[ m^2\ge \frac{A}{16\pi}+\frac{Q^2}{2} +\frac{\pi(Q^4+4J^2)}{A}, \qquad A\ge4\pi\sqrt{Q^4+4J^2}, \tag{1}\] where \(A\) is a suitable enclosing area, \(m\) is mass, \(Q\) is total charge, and \(J\) is angular momentum. Its motivation comes from the gravitational-collapse argument underlying Penrose’s original proposal (Penrose 1973); the area branch in (1) is the branch on which the right side increases with \(A\). The formulation using minimum enclosing area, and the distinction between the corresponding upper and lower area bounds, are discussed in (Dain et al. 2013). In the presence of an electromagnetic field, specifying \(J\) is essential. The gravitational angular-momentum flux through a surface need not be conserved. For the globally defined invariant magnetic potential used below, adding the corresponding electromagnetic boundary term gives a conserved flux under the electrovacuum constraints. Khuri–Weinstein include this term in their definition at infinity and state that its disappearance requires appropriate asymptotic conditions (Khuri and Weinstein 2016, Equation (1.12) and the following paragraph). The identification argument of Dain–Khuri–Weinstein–Yamada uses Coulomb expansions for the electric and magnetic fields (Dain et al. 2013, Equations (2.20) and (2.22)). The weaker bounds \(\mathcal E,\mathcal B=O(r^{-2})\) alone do not impose those leading terms. We give a negative resolution of (1) for the precise formulation below: \(J\) is the gravitational ADM flux at infinity, and the electromagnetic assumptions are only the displayed decay and constraint conditions. Our examples have a connected outermost future marginally outer trapped boundary, minimize area among all enclosing cuts, and lie strictly on the physical area branch. They also disprove the associated nondegenerate Kerr–Newman equality implication. Stronger asymptotic conditions, or a definition of angular momentum retaining its electromagnetic term, lead to different assertions; no conclusion about those assertions follows from these examples. Conventions and the class of exteriorsWe work in three spatial dimensions with \(G=c=1\) and zero cosmological constant. On an oriented Riemannian three-manifold \((\Omega,g)\), the cross product is defined by \[(X\times_gY)^\flat(Z)=\mathop{}\!\mathrm dV_g(X,Y,Z).\] Let \(K\) be a symmetric covariant two-tensor, \(\tau=\mathop{\mathrm{tr}}_gK\), and let \(\mathcal E,\mathcal B\) be electric and magnetic vector fields. We use the source-free electrovacuum constraint equations \[ \begin{split} R_g+\tau^2-\left\lvert K\right\rvert_g^2 &=2\bigl(\left\lvert\mathcal E\right\rvert_g^2+\left\lvert\mathcal B\right\rvert_g^2\bigr),\\ \operatorname{div}_g(K-\tau g) &=2(\mathcal E\times_g\mathcal B)^\flat,\\ \operatorname{div}_g\mathcal E &=\operatorname{div}_g\mathcal B=0. \end{split} \tag{2}\] Thus the total constraint densities are \[\mu=\frac{\left\lvert\mathcal E\right\rvert_g^2+\left\lvert\mathcal B\right\rvert_g^2}{8\pi}, \qquad J_{\mathrm{con}} =\frac{(\mathcal E\times_g\mathcal B)^\flat}{4\pi}.\] The covector \(J_{\mathrm{con}}\) is distinct from the scalar angular momentum \(J\). All data are smooth, including at the compact boundary. The exterior is connected, complete as a metric space with its boundary included, and diffeomorphic to the complement of an open ball in \(\mathbb R^3\). There is one asymptotically flat end and a smooth effective \(U(1)\) action preserving \(g,K,\mathcal E,\mathcal B\) and the boundary \(S\). Its generator \(\eta\) has period \(2\pi\) and agrees on the end with the standard oriented rotation about the third coordinate axis. In asymptotic Euclidean coordinates \(x\), write \(r=\left\lvert x\right\rvert\) and assume \[ g_{ij}-\delta_{ij}=O_4(r^{-1}),\qquad K_{ij}=O_3(r^{-3}),\qquad \mathcal E^i,\mathcal B^i=O_3(r^{-2}). \tag{3}\] Here \(O_k(r^{-a})\) bounds each Cartesian derivative of order \(j\le k\) by \(Cr^{-a-j}\). We require \(\mu\) and \(\left\lvert J_{\mathrm{con}}\right\rvert_g\) to be integrable. On a coordinate sphere \(S_R\), let \(n_\delta\) be its Euclidean outward unit normal and \(\nu_R\) its \(g\)-unit outward normal. The asymptotic quantities are \[\begin{align*} E_{\mathrm{ADM}} &=\frac1{16\pi}\lim_{R\to\infty} \int_{S_R}(\partial_jg_{ij}-\partial_ig_{jj}) n_\delta^i\,\mathop{}\!\mathrm dA_\delta,\tag{4}\\ (P_{\mathrm{ADM}})_i &=\frac1{8\pi}\lim_{R\to\infty} \int_{S_R}(K_{ij}-\tau g_{ij})n_\delta^j\,\mathop{}\!\mathrm dA_\delta, \tag{5}\\ J&=\frac1{8\pi}\lim_{R\to\infty} \int_{S_R}(K_{ij}-\tau g_{ij})\nu_R^i\eta^j\,\mathop{}\!\mathrm dA_g, \tag{6}\\ Q_e&=\frac1{4\pi}\lim_{R\to\infty} \int_{S_R}g(\mathcal E,\nu_R)\,\mathop{}\!\mathrm dA_g, \qquad Q_b=\frac1{4\pi}\lim_{R\to\infty} \int_{S_R}g(\mathcal B,\nu_R)\,\mathop{}\!\mathrm dA_g. \tag{7}\end{align*}\] All these limits are required to exist and be finite. We impose \(P_{\mathrm{ADM}}=0\) and \(E_{\mathrm{ADM}}>0\), and set \(m=E_{\mathrm{ADM}}\), \(Q=\sqrt{Q_e^2+Q_b^2}\). In particular (6) has no electromagnetic surface term. For a two-sided surface with unit normal \(\nu\) pointing toward the end, put \[ H_\Sigma=\operatorname{div}_\Sigma\nu, \qquad \mathop{\mathrm{tr}}_\Sigma K=(g^{ij}-\nu^i\nu^j)K_{ij}, \qquad \theta_+(\Sigma)=H_\Sigma+\mathop{\mathrm{tr}}_\Sigma K. \tag{8}\] Thus outward Euclidean spheres have positive mean curvature. The boundary is a future marginally outer trapped surface, or future MOTS, if \(\theta_+(S)=0\). We require this condition and the following precise outermostness condition: there is no compact smooth embedded enclosing surface entirely in the interior, possibly disconnected, for which \(\theta_+\le0\) everywhere. No condition on the inward expansion is imposed. An enclosing cut is the compact intrinsic boundary \(\Gamma=\partial D\) of a connected smooth codimension-zero submanifold \(D\subset\Omega\), closed as a subset of \(\Omega\), whose interior lies in \(\operatorname{int}\Omega\) and which contains the entire sufficiently distant part of the end. The full intrinsic boundary is counted, including any portion coinciding with \(S\). The cut must be smooth, embedded and two-sided; it may be disconnected and need not be symmetric or trapped. In particular \(D=\Omega\) gives the cut \(S\). The required area condition is \[ \mathop{\mathrm{Area}}_g(\Gamma)\ge A:=\mathop{\mathrm{Area}}_g(S)>0 \qquad\text{for every enclosing cut }\Gamma. \tag{9}\] Thus \(A\) is exactly the minimum enclosing area. Main resultTheorem 1. In the class of smooth axisymmetric source-free electrovacuum exteriors specified in 1.1, both of the following occur.
All these examples are maximal, \(\mathop{\mathrm{tr}}_gK=0\), on \(\Omega=\{x\in\mathbb R^3:\left\lvert x\right\rvert\ge1\}\), with connected boundary \(S=\{r=1\}\). Both electromagnetic fields are nonzero. The first conclusion disproves (1) under (3) and (6), even with every horizon and area condition above. The second disproves the nondegenerate equality implication with a Kerr–Newman realization, including realizations by nonconstant time graphs. The usual spacetime convention is \(K(X,Y)=\overline g(\overline\nabla_X n,Y)\) for the future unit normal \(n\); the induced fields are \(\mathcal E^\flat=(\iota_nF)|_{T\Omega}\) and \(\mathcal B^\flat=\star_g(F|_{T\Omega})\). The obstruction in the second conclusion holds for every spacelike slice, so it also excludes a realization extending smoothly to a cross-section of the future horizon or to the bifurcation sphere. The proof is independent of any neutral or charged Penrose inequality or rigidity theorem. It gives a two-parameter family, with amplitude \(s\in[0,1]\) and a transition radius \(L\), such that \[ J=-\frac{4s^2}{15},\qquad m=2+O(s^2/L),\qquad A=64\pi+O(s^2/L^2). \tag{10}\] The electromagnetic correction to the angular-momentum flux is \(+4s^2/15\) at infinity. Angularly varying radial tails supply this term while their energy tends to zero as \(L\to\infty\). At \(s=1\), the defect in (1) tends to \(-1/225\). For fixed large \(L\), a small positive amplitude gives the opposite sign, yielding the equality example by continuity. Two features make the construction explicit. A mixed radial–azimuthal tensor solves the electromagnetic momentum constraint without a vector elliptic equation. Kelvin reflection then gives a positive Newton-potential operator with the Robin condition that makes the unit sphere minimal. The contraction estimate uses the distance of the source from that sphere, and accommodates a noncompact \(r^{-4}\) energy tail. Finally, a closed two-form obtained by radial projection proves area-minimization against all enclosing cuts. These steps are developed in [sec:construction,sec:elliptic,sec:geometry]; 5 computes the fluxes and proves both conclusions. Explicit electrovacuum seedsWe first solve the momentum and Maxwell divergence constraints explicitly. The electromagnetic fields begin at a large radius; one scalar equation will then determine the metric. Throughout the construction, \[\Omega=\{x\in\mathbb R^3:r=|x|\ge1\},\qquad S=\{r=1\},\] with the Euclidean metric \(\delta\) and standard orientation. Write \((r,\theta,\phi)\) for spherical coordinates, with \(\theta\) the colatitude and \(\eta=\partial_\phi\) the generator of the standard rotation action. Fix a smooth function \(\chi:\mathbb R\to[0,1]\) that is zero on \((-\infty,1]\) and one on \([2,\infty)\). For example, one may take \[\chi(t)=\frac{\zeta(t-1)}{\zeta(t-1)+\zeta(2-t)},\qquad \zeta(t)=\begin{cases}e^{-1/t},&t>0,\\0,&t\le0.\end{cases}\] For \(s\in[0,1]\) and \(L>4\), put \(c(r)=\chi(r/L)\). Constants \(C,C_j\) below depend only on this fixed cutoff and the indicated derivative order, and are independent of \(s\) and \(L\) above a fixed lower threshold. Define the scalar potentials, vector fields, and symmetric covariant tensor by \[ \begin{aligned} f&=sc(r)\sin^2\theta,& h&=sc(r)\sin^2\theta\cos\theta,\\ \iota_{\widetilde E}\mathop{}\!\mathrm dV_\delta&=\mathop{}\!\mathrm d(f\,\mathop{}\!\mathrm d\phi),& \iota_{\widetilde B}\mathop{}\!\mathrm dV_\delta&=\mathop{}\!\mathrm d(h\,\mathop{}\!\mathrm d\phi),\\ \sigma&=-\frac{s^2c(r)^2\sin^4\theta}{r^2} (\mathop{}\!\mathrm dr\otimes\mathop{}\!\mathrm d\phi+\mathop{}\!\mathrm d\phi\otimes\mathop{}\!\mathrm dr). \end{aligned} \tag{11}\] Here contraction with the volume form uniquely determines each vector field. The coordinate formulas in (11) extend smoothly across the rotation axis, as 2 verifies. Lemma 2 (Seed identities and bounds). The fields \(\widetilde E,\widetilde B\) and tensor \(\sigma\) are smooth on \(\Omega\), invariant under the standard effective \(U(1)\) action, and vanish on \(\{1\le r\le L\}\). They satisfy \[ \operatorname{div}_\delta\widetilde E =\operatorname{div}_\delta\widetilde B=0,\qquad \mathop{\mathrm{tr}}_\delta\sigma=0,\qquad \operatorname{div}_\delta\sigma =2\mathop{}\!\mathrm dV_\delta(\widetilde E,\widetilde B,\cdot). \tag{12}\] For every Cartesian multi-index \(\alpha\) of order \(j\ge0\), \[ \begin{aligned} |\partial^\alpha\widetilde E|_\delta +|\partial^\alpha\widetilde B|_\delta &\le C_j s r^{-2-j}\mathbf1_{\{r\ge L\}},\\ |\partial^\alpha\sigma|_\delta &\le C_j s^2 r^{-3-j}\mathbf1_{\{r\ge L\}}. \end{aligned} \tag{13}\] Moreover, \(\operatorname{div}_\delta\sigma\) is supported in \(\{L\le r\le2L\}\), with \[|\partial^\alpha(\operatorname{div}_\delta\sigma)|_\delta \le C_j s^2r^{-4-j}\mathbf1_{\{L\le r\le2L\}}.\] Proof. Let \(X=(x^1,x^2,x^3)\), \(T=(-x^2,x^1,0)\), and \(w=(x^1)^2+(x^2)^2\), with \(\flat\) denoting Euclidean metric dual in the following formulas. The primitives and tensor have the Cartesian expressions \[ \begin{aligned} f\,\mathop{}\!\mathrm d\phi&=\frac{sc}{r^2}(x^1\mathop{}\!\mathrm dx^2-x^2\mathop{}\!\mathrm dx^1),\\ h\,\mathop{}\!\mathrm d\phi&=\frac{scx^3}{r^3}(x^1\mathop{}\!\mathrm dx^2-x^2\mathop{}\!\mathrm dx^1),\\ \sigma&=-\frac{s^2c^2w}{r^7} (X^\flat\otimes T^\flat+T^\flat\otimes X^\flat). \end{aligned} \tag{14}\] These are smooth on \(r\ge1\), with no denominator involving distance to the axis. They also show rotation invariance and, since \(X\cdot T=0\), trace-freeness. The exact flux forms are closed, so \(\mathop{}\!\mathrm d(\iota_{\widetilde E}\mathop{}\!\mathrm dV_\delta)= (\operatorname{div}_\delta\widetilde E)\mathop{}\!\mathrm dV_\delta=0\), and likewise for \(\widetilde B\). For the momentum identity we compute off the axis and then use smoothness. Since \(\mathop{}\!\mathrm dV_\delta=r^2\sin\theta\,\mathop{}\!\mathrm dr\wedge\mathop{}\!\mathrm d\theta\wedge\mathop{}\!\mathrm d\phi\), the spherical coordinate components are \[ \begin{aligned} \widetilde E^r&=\frac{2sc\cos\theta}{r^2},& \widetilde E^\theta&=-\frac{sc'\sin\theta}{r^2},& \widetilde E^\phi&=0,\\ \widetilde B^r&=\frac{sc(3\cos^2\theta-1)}{r^2},& \widetilde B^\theta&=-\frac{sc'\sin\theta\cos\theta}{r^2},& \widetilde B^\phi&=0. \end{aligned} \tag{15}\] In particular an orthonormal polar component is \(r\) times the displayed polar component. The covector \(\mathop{}\!\mathrm dV_\delta(\widetilde E,\widetilde B,\cdot)\) has only an azimuthal component, namely \[ \mathop{}\!\mathrm dV_\delta(\widetilde E,\widetilde B,\partial_\phi) =\frac{f_rh_\theta-f_\theta h_r}{r^2\sin\theta} =-\frac{s^2cc'\sin^4\theta}{r^2}. \tag{16}\] To check all three tensor-divergence components, put \(a=\sigma_{r\phi}=-s^2c^2\sin^4\theta/r^2\) and use the coordinate identity \[(\operatorname{div}_\delta\sigma)_i =\frac1{\sqrt{\det\delta}}\partial_j \bigl(\sqrt{\det\delta}\,\sigma_i{}^j\bigr) -\frac12\sigma^{jk}\partial_i\delta_{jk}.\] The last term vanishes: the metric is diagonal in these coordinates, whereas \(\sigma^{jk}\) has only \(r\phi\) and \(\phi r\) entries. Thus \[ \begin{aligned} (\operatorname{div}_\delta\sigma)_r &=\frac1{r^2\sin\theta}\partial_\phi \left(\frac a{\sin\theta}\right)=0,\\ (\operatorname{div}_\delta\sigma)_\theta&=0,\\ (\operatorname{div}_\delta\sigma)_\phi &=\frac1{r^2\sin\theta}\partial_r(r^2\sin\theta\,a) =a_r+\frac{2a}{r} =-\frac{2s^2cc'\sin^4\theta}{r^2}. \end{aligned} \tag{17}\] This proves (12). The last expression also proves the stated support property; for \(r\ge2L\) both fields are radial and parallel. Finally, for \(j\ge1\) the Cartesian derivatives of \(c\) obey \[|\partial^\alpha c|\le C_j r^{-j} \mathbf1_{\{L\le r\le2L\}}.\] Indeed, all such derivatives are supported where \(r\) is comparable to \(L\), and differentiating \(\chi(r/L)\) gives this bound. The coefficients of the two uncut primitive one-forms in (14) are smooth homogeneous functions of degree \(-1\); exterior differentiation therefore gives the field bounds in (13). The uncut Cartesian tensor coefficients are smooth homogeneous functions of degree \(-3\), which gives the tensor bounds by the same product rule. Smooth cutoff jets give the estimates also at \(r=L\) and \(r=2L\). ◻ Reduction to a scalar equationFor a positive function \(u\), define the physical data by \[ g=u^4\delta,\qquad K=u^{-2}\sigma,\qquad \mathcal E=u^{-6}\widetilde E,\qquad \mathcal B=u^{-6}\widetilde B. \tag{18}\] The next calculation fixes the conformal weights and the normalization of the scalar equation directly. Lemma 3 (Conformal electrovacuum equations). For every smooth positive \(u\), the data (18) satisfy \[\tau=\mathop{\mathrm{tr}}_g K=0,\qquad \operatorname{div}_g\mathcal E=\operatorname{div}_g\mathcal B=0, \qquad \operatorname{div}_gK=2(\mathcal E\times_g\mathcal B)^\flat_g.\] They satisfy the electrovacuum Hamiltonian constraint if and only if \[ -\Delta_\delta u=p u^{-7}+\ell u^{-3}\quad\hbox{on }\Omega, \tag{19}\] where, writing \(t=\cos\theta\), \[ \begin{aligned} p&=\frac18|\sigma|_\delta^2 =\frac{s^4c^4\sin^6\theta}{4r^6},\\ \ell&=\frac14\bigl(|\widetilde E|_\delta^2 +|\widetilde B|_\delta^2\bigr)\\ &=\frac{s^2}{4r^4} \left[c^2(9t^4-2t^2+1)+(rc')^2(1-t^4)\right]. \end{aligned} \tag{20}\] The nonnegative coefficients \(p,\ell\) are smooth and rotation invariant. For every Cartesian multi-index \(\alpha\) of order \(j\ge0\), \[ |\partial^\alpha p|\le C_j s^4r^{-6-j}\mathbf1_{\{r\ge L\}},\qquad |\partial^\alpha\ell|\le C_j s^2r^{-4-j}\mathbf1_{\{r\ge L\}}. \tag{21}\] Proof. Since \(\mathop{}\!\mathrm dV_g=u^6\mathop{}\!\mathrm dV_\delta\), the flux forms are unchanged: \[ \iota_{\mathcal E}\mathop{}\!\mathrm dV_g=\iota_{\widetilde E}\mathop{}\!\mathrm dV_\delta, \qquad \iota_{\mathcal B}\mathop{}\!\mathrm dV_g=\iota_{\widetilde B}\mathop{}\!\mathrm dV_\delta. \tag{22}\] Their closedness proves the two Maxwell divergence constraints. The cross-product convention gives, for any vector \(Z\), \[(\mathcal E\times_g\mathcal B)^\flat_g(Z) =\mathop{}\!\mathrm dV_g(\mathcal E,\mathcal B,Z) =u^{-6}\mathop{}\!\mathrm dV_\delta(\widetilde E,\widetilde B,Z).\] Here is the tensor calculation in Cartesian coordinates. Put \(b=2\log u\), so \(g=e^{2b}\delta\) and \(K=e^{-b}\sigma\). The connection difference is \[D^i_{jk}=\delta^i_j b_k+\delta^i_k b_j-\delta_{jk}b^i, \qquad b_j=\partial_jb,\] with indices on the right raised by \(\delta\). Symmetry and trace-freeness give \[\delta^{jk}D^l_{ki}\sigma_{lj}=0,\qquad \delta^{jk}D^l_{kj}=-b^l.\] Consequently the derivative of \(e^{-b}\) cancels the second connection contraction, yielding \[ \mathop{\mathrm{tr}}_gK=e^{-3b}\mathop{\mathrm{tr}}_\delta\sigma=0,\qquad (\operatorname{div}_gK)_i =e^{-3b}\delta^{jk}\partial_k\sigma_{ij} =u^{-6}(\operatorname{div}_\delta\sigma)_i. \tag{23}\] Together with 2, this is exactly the momentum constraint. For completeness, substituting the same \(D\) into \[R_{ij}=\partial_kD^k_{ij}-\partial_jD^k_{ik} +D^k_{kl}D^l_{ij}-D^k_{jl}D^l_{ik}\] gives derivative terms \(-b_{ij}-\delta_{ij}\Delta_\delta b\) and quadratic terms \(b_i b_j-\delta_{ij}|\nabla_\delta b|^2\). Contracting therefore gives \[ R_g=e^{-2b}\bigl(-4\Delta_\delta b -2|\nabla_\delta b|^2\bigr) =-8u^{-5}\Delta_\delta u. \tag{24}\] The norm identities are \[|K|_g^2=u^{-12}|\sigma|_\delta^2,\qquad |\mathcal E|_g^2=u^{-8}|\widetilde E|_\delta^2,\qquad |\mathcal B|_g^2=u^{-8}|\widetilde B|_\delta^2.\] Thus the Hamiltonian constraint \(R_g-|K|_g^2=2(|\mathcal E|_g^2+|\mathcal B|_g^2)\) becomes (19), with the coefficients in (20). Those explicit expressions follow from (15) and \(|\sigma|_\delta^2=2a^2/(r^2\sin^2\theta)\). Their smoothness, nonnegativity, and derivative bounds follow either from their definitions as squared norms or directly from 2 and the product rule. ◻ We will solve (19) with \(u\to1\) at infinity and \(\partial_ru+u/2=0\) on \(S\). This boundary condition will make \(S\) a future MOTS. A solution and the uniform estimates needed to control all other surfaces are established next. The conformal factorThe sources begin at radius \(L\). This separation makes their effect near \(S\) of order \(L^{-2}\), while their total integral is of order \(L^{-1}\). We prove both estimates by constructing the solution directly. Proposition 4 (Uniform conformal factor). There is \(L_0>4\), depending only on the fixed cutoff, such that for every \(L\ge L_0\) and \(s\in[0,1]\), Equation (19) has a smooth axisymmetric solution \(u=u_0+v\) on \(\Omega\), where \(u_0=1+r^{-1}\), satisfying \[ u_r+\tfrac12u=0\quad\text{on }S, \qquad u\longrightarrow1\quad\text{as }r\longrightarrow\infty. \tag{25}\] It is unique among smooth solutions with \(v\ge0\) and \(v=O(r^{-1})\). One may take a constant \(M=3\), independent of \(s,L\), such that \(1\le u_0\le u\le M\). Uniformly in these parameters, \[ \begin{aligned} |v|+r|v_r|&\le Cs^2L^{-2} &&(r\ge1),\\ |v_{rr}|&\le Cs^2L^{-2} &&(1\le r\le2). \end{aligned} \tag{26}\] For every integer \(k\ge0\), all Cartesian derivatives satisfy \[ |D^kv|\le C_k\frac{s^2}{L}r^{-1-k}, \qquad |D^k(u-1)|\le C_k r^{-1-k}. \tag{27}\] For fixed \(L\), \(s\mapsto v_s\) is continuous in the sup norm, as are \[s\longmapsto\int_{\mathbb R^3}H_s\,\mathop{}\!\mathrm dx, \qquad s\longmapsto N_s(0)=\frac1{4\pi}\int_{\mathbb R^3}\frac{H_s(x)}{|x|}\,\mathop{}\!\mathrm dx,\] where \(H_s=pu_s^{-7}+\ell u_s^{-3}\) is extended by zero inside the unit ball and \(N_s\) is its Newton potential. Moreover, \[ 0\le\int_{\mathbb R^3}H_s\,\mathop{}\!\mathrm dx\le C\frac{s^2}{L}, \qquad 0\le N_s(0)\le C\frac{s^2}{L^2}. \tag{28}\] At \(s=0\) the solution is exactly \(u_0\). Here and below \(D^k\) denotes any Cartesian derivative of order \(k\); constants may depend on this order but are independent of \(s\in[0,1]\) and \(L\ge L_0\). We first record the elementary potential estimates used in the proof. Lemma 5 (Newton potential and Robin reflection). Suppose \(L>4\) and \(H\in C(\mathbb R^3)\) obeys \[|H(y)|\le a|y|^{-4}\mathbf1_{\{|y|\ge L\}}.\] Then \[N_H(x)=\frac1{4\pi}\int_{\mathbb R^3}\frac{H(y)}{|x-y|}\,\mathop{}\!\mathrm dy\] is well defined and \(C^1\), satisfies \(-\Delta N_H=H\) distributionally, and is smooth and harmonic on \(\{|x|<L\}\). For \(j=0,1\), \[ |D^jN_H(x)|\le \begin{cases} C_j aL^{-2-j},& |x|\le L/2,\\ C_j a\bigl(r^{-2-j}+L^{-1}r^{-1-j}\bigr),&r=|x|\ge L/2. \end{cases} \tag{29}\] The first estimate holds for every \(j\ge0\). The operator \[(\mathcal RH)(x)=N_H(x)+r^{-1}N_H(x/r^2),\qquad x\in\Omega,\] preserves nonnegativity, satisfies \(-\Delta\mathcal RH=H\) distributionally and \((\partial_r+\tfrac12)\mathcal RH=0\) on \(S\), and obeys \[ \begin{gathered} \|\mathcal RH\|_\infty\le CaL^{-2},\qquad |\mathcal RH|+r|\partial_r\mathcal RH|\le CaL^{-2},\\ |\partial_r^2\mathcal RH|\le CaL^{-2}\quad(1\le r\le2). \end{gathered} \tag{30}\] The homogeneous harmonic Robin problem on \(\Omega\) has no nonzero solution of order \(O(r^{-1})\). Proof. The bound on \(H\) gives \(\|H\|_{L^1}\le4\pi a/L\). The Newton kernel and its first derivatives are locally integrable. Localizing around a point and separating a small ball about its kernel singularity proves that the potential and its first derivatives are continuous; the singular part of the gradient integral is bounded by a constant times the ball’s radius. The fundamental-solution identity \(-\Delta(4\pi|x|)^{-1}=\delta_0\), obtained by integration by parts around the pole, therefore gives the distributional equation by Fubini’s theorem. For \(|x|\le L/2\), source separation allows differentiation under the integral to every order and gives \[|D^jN_H(x)|\le C_j a\int_L^\infty \rho^{-4}\rho^{-1-j}\rho^2\,\mathop{}\!\mathrm d\rho \le C_j aL^{-2-j}.\] For \(r=|x|\ge L/2\), split the integral at \(|x-y|=r/2\). In the inner part \(|y|\ge r/2\), and for \(j=0,1\) its absolute value is at most \[Ca r^{-4}\int_{|z|\le r/2}|z|^{-1-j}\,\mathop{}\!\mathrm dz \le Ca r^{-2-j}.\] The outer part is bounded by \(Cr^{-1-j}\|H\|_{L^1}\le CaL^{-1}r^{-1-j}\). This proves (29). Write \(I(x)=x/r^2\). The Kelvin identity in dimension three is \[\Delta\bigl(r^{-1}N_H(I(x))\bigr) =r^{-5}(\Delta N_H)(I(x)).\] Its right side vanishes on \(\Omega\) because \(|I(x)|\le1<L\). Along each ray \(x=r\omega\), direct differentiation gives \[(\mathcal RH)(\omega)=2N_H(\omega),\qquad \partial_r(\mathcal RH)(\omega)=-N_H(\omega),\] which proves the Robin condition including its angular dependence. Positivity follows from the two positive kernels. Combining (29) in the two radial regions gives the first two bounds in (30). For the last bound both potential arguments are separated from the source. More explicitly, put \(n(t)=N_H(t\omega)\); then \[\frac{\mathop{}\!\mathrm d^2}{\mathop{}\!\mathrm dr^2}\bigl(r^{-1}n(r^{-1})\bigr) =2r^{-3}n(r^{-1})+4r^{-4}n'(r^{-1})+r^{-5}n''(r^{-1}).\] The separated estimates for \(j=0,1,2\), together with the corresponding bound for \(\partial_r^2N_H(r\omega)\), give the claim. For uniqueness of the homogeneous problem, let \(w\) be harmonic on \(\Omega\), \(C^1\) up to \(S\), with \(w=O(r^{-1})\) and \(w_r+\tfrac12w=0\) on \(S\). Extend it to \(0<r<1\) by \(r^{-1}w(I(x))\). Values and tangential derivatives match at \(S\), and the interior radial derivative there is \(-w-w_r=w_r\). Thus the extension \(\widehat w\) is distributionally harmonic across \(S\) and hence harmonic on \(\mathbb R^3\setminus\{0\}\). It is bounded near zero and is \(O(r^{-1})\) at infinity. For sufficiently small \(\varepsilon\) and large \(R\), the maximum principle applied to \(\pm\widehat w\) on \(\varepsilon<|x|<R\) gives \(|\widehat w(x)|\le C\varepsilon/|x|+C/R\), with one constant \(C\): the harmonic right side bounds the absolute value on both boundary spheres. Letting \(\varepsilon\downarrow0\) and \(R\to\infty\) proves uniqueness. ◻ Lemma 6 (Localized derivative bootstrap). Let \(H\) satisfy the hypotheses of 5. If \(H\in C^j(\mathbb R^3)\), then \(N_H\in C^{j+1}(\mathbb R^3)\). If, in addition, \[|D^iH(y)|\le a_j|y|^{-4-i}\mathbf1_{\{|y|\ge L\}} \quad(0\le i\le j),\] then (29) holds through order \(j+1\), with a constant \(C_j a_j\) in place of \(C_j a\). Consequently the local gain in regularity and the derivative decay use no first moment of \(H\). Proof. Fix a point \(x_0\). Choose a smooth compactly supported cutoff \(\zeta\) equal to one near \(x_0\), and keep this cutoff fixed when differentiating with respect to \(x\). Writing \(G(x)=(4\pi|x|)^{-1}\), for \(|\alpha|=j+1\) choose \(\alpha=\beta+e_i\), \(|\beta|=j\). Distributional integration by parts gives, near \(x_0\), \[ D^\alpha N_H(x) =\int_{\mathbb R^3}\partial_iG(x-y)D_y^\beta(\zeta H)(y)\,\mathop{}\!\mathrm dy +\int_{\mathbb R^3}D^\alpha G(x-y)(1-\zeta(y))H(y)\,\mathop{}\!\mathrm dy. \tag{31}\] There are no boundary terms since \(\zeta H\) has compact support. The first integral is continuous by local integrability of \(\partial_iG\) and the same small-ball argument as above; the second is smooth near \(x_0\) by source separation and integrability at infinity. The formula, also applied at lower orders, proves \(C^{j+1}\) regularity. For the estimate at \(|x_0|=r\ge L/2\), take \(\zeta\) supported in \(B(x_0,r/2)\) and equal to one on \(B(x_0,r/4)\), with \(|D^i\zeta|\le C_i r^{-i}\). On its support every Leibniz term in \(D^\beta(\zeta H)\) is bounded by \(C_j a_jr^{-4-j}\). The first integral in (31) at \(x_0\) is therefore bounded by \[C_j a_jr^{-4-j}\int_{|z|\le r/2}|z|^{-2}\,\mathop{}\!\mathrm dz \le C_j a_jr^{-3-j}.\] The second is bounded by \[C_jr^{-2-j}\|H\|_{L^1} \le C_j a_jL^{-1}r^{-2-j}.\] These are precisely the claimed bounds at order \(j+1\); the same argument gives the lower orders. On \(|x|\le L/2\), differentiation of the separated integral already gives every derivative estimate. ◻ Proof of 4. The coefficient estimates from [lem:seeds,lem:conformal] imply, for every \(j\ge0\), \[ |D^jp|\le C_js^4r^{-6-j}\mathbf1_{\{r\ge L\}},\qquad |D^j\ell|\le C_js^2r^{-4-j}\mathbf1_{\{r\ge L\}}. \tag{32}\] Let \(\mathcal C=\{v\in C_b(\Omega):v\ge0\}\), a complete metric space in the sup norm. For \(v\in\mathcal C\), put \[ \begin{split} H_v&=p(u_0+v)^{-7}+\ell(u_0+v)^{-3},\\ N_v(x)&=\frac1{4\pi}\int_{\mathbb R^3}\frac{H_v(y)}{|x-y|}\,\mathop{}\!\mathrm dy,\\ T_s(v)(x)&=N_v(x)+r^{-1}N_v(x/r^2). \end{split} \tag{33}\] The source is extended by zero to \(\mathbb R^3\); this extension is continuous because it vanishes on \(r<L\). Since \(u_0+v\ge1\), the mean-value theorem for the negative powers gives \[ \begin{split} 0\le H_v&\le Cs^2r^{-4}\mathbf1_{\{r\ge L\}},\\ |H_v-H_w|&\le Cs^2r^{-4}\mathbf1_{\{r\ge L\}} \|v-w\|_\infty. \end{split} \tag{34}\] By 5, \(T_s\) maps the whole cone \(\mathcal C\) into itself, with \[\|T_s(v)\|_\infty\le Cs^2L^{-2},\qquad \|T_s(v)-T_s(w)\|_\infty\le CL^{-2}\|v-w\|_\infty.\] Choose \(L_0>4\) so large that the second constant is at most \(1/2\) and the first, for \(s=1\), is at most one. The contraction theorem now gives an actual fixed point \(v\in\mathcal C\), unique in that cone. The operator commutes with rotation about the prescribed axis, so its unique fixed point is axisymmetric. It also gives \(v=0\) when \(s=0\). 5 yields (26), \(1\le u_0\le u\le3\), the Robin condition, and the scalar equation distributionally. The potential bounds give \(v=O(r^{-1})\), hence \(u\to1\). For completeness we bootstrap both smoothness and the full decay. The fixed point is \(C^1\) by 5; near \(S\) it is already smooth by separation of both potential arguments from the source. Suppose derivatives of \(u-1\) through order \(j\) have the uniform bounds \(C_i r^{-1-i}\). The case \(j=0\) follows from (29) and \(u_0=1+r^{-1}\). Because \(u\ge1\), the chain and product rules with (32) give \[|D^iH_v|\le C_js^2r^{-4-i}\mathbf1_{\{r\ge L\}} \quad(0\le i\le j).\] The zero extension is \(C^j\) since the source vanishes near \(S\). 6 gains one derivative for \(N_v\) with the stated bounds. For the reflected term, the separated estimates and the homogeneity of \(I(x)\) give at every order \(k\) \[\bigl|D^k\bigl(r^{-1}N_v(I(x))\bigr)\bigr| \le C_ks^2L^{-2}r^{-1-k}.\] In the region \(1\le r\le L/2\), the ordinary term is bounded by \(C_ks^2L^{-2-k}\); in the other region it is bounded by \(C_ks^2(r^{-2-k}+L^{-1}r^{-1-k})\). Each bound is at most \(C_ks^2L^{-1}r^{-1-k}\) in its region. This proves the next induction step and (27) at every order, with constants independent of \(L,s\). In particular \(u\) is smooth and solves (19) pointwise. If another solution has \(v\ge0\) and \(v=O(r^{-1})\), subtracting \(\mathcal RH_v\) leaves a decaying homogeneous harmonic Robin solution. 5 makes that difference zero; cone uniqueness then identifies the solutions. Finally, fix \(L\ge L_0\). The coefficients have the form \(p_s=s^4p_1\), \(\ell_s=s^2\ell_1\). Thus uniformly for \(v\in\mathcal C\), \[\|T_s(v)-T_t(v)\|_\infty\le CL^{-2}|s-t|.\] Using the contraction constant \(1/2\) in the fixed point equations gives \[\|v_s-v_t\|_\infty\le 2CL^{-2}|s-t|.\] Combining this with the pointwise source Lipschitz bound yields \[|H_s(y)-H_t(y)|\le C|s-t|\,|y|^{-4} \mathbf1_{\{|y|\ge L\}}.\] The right side is integrable with both weights \(1\) and \(|y|^{-1}\). It proves continuity of the two source integrals, and (34) gives their quantitative bounds (28). Uniform convergence also gives continuity of the restriction \(u_s|_S\), which will imply continuity of the horizon area. No finite first moment or multipole expansion has entered the argument. ◻ The horizon and all enclosing surfacesThe estimates for the conformal factor make the coordinate spheres useful barriers. Their positive expansion outside the boundary also gives a closed two-form that compares the boundary area with every enclosing cut. Proposition 7. For all sufficiently large \(L\), uniformly for \(0\le s\le1\), the data constructed in [sec:construction,sec:elliptic] are smooth source-free electrovacuum data on \[\Omega=\{x\in\mathbb R^3:r\ge1\},\qquad S=\{r=1\},\] with \(\mathop{\mathrm{tr}}_gK=0\). The metric is complete with its boundary included, has one asymptotically flat end, and satisfies \[ g_{ij}-\delta_{ij}=O_4(r^{-1}),\qquad K_{ij}=O_3(r^{-3}),\qquad \mathcal E^i,\mathcal B^i=O_3(r^{-2}). \tag{35}\] The standard effective \(U(1)\) action preserves all the data and \(S\), and \(\mu\) and \(|J_{\mathrm{con}}|_g\) are integrable. With the normal toward the end, \(S\) is a future MOTS, and no compact smooth embedded enclosing surface in \(\operatorname{int}\Omega\), possibly disconnected, has everywhere nonpositive future outward expansion. Every enclosing cut in the sense of 1 has area at least \(A=\mathop{\mathrm{Area}}_g(S)\). Moreover, for a constant \(C_A\) independent of \(s,L\), \[ 64\pi\le A\le64\pi+C_A\frac{s^2}{L^2}, \qquad A=64\pi+O(s^2/L^2), \tag{36}\] and \(A(s,L)\) is continuous in \(s\) for fixed \(L\). Proof. The domain is connected and oriented by the standard orientation of \(\mathbb R^3\); its boundary is a connected two-sphere. Outside each sufficiently large ball it is a single product end. Smoothness up to \(S\) and across the rotation axis follows from [lem:seeds,prop:conformal-factor]; the constraints and maximality follow from 3. The seeds are invariant under the standard rotations, and uniqueness of the fixed point makes \(u\) invariant as well. Thus these rotations preserve \(g,K,\mathcal E,\mathcal B\), with generator \(\eta=\partial_\phi\) of period \(2\pi\). The action is effective because a rotation fixing every point of \(\Omega\) has angle \(0\) modulo \(2\pi\). The bound \(1\le u\le M\) gives \(\delta\le g\le M^4\delta\). The Euclidean intrinsic length metric on the closed exterior is complete: a Cauchy sequence has a Euclidean limit in \(\Omega\), and local smooth boundary charts show convergence also in the intrinsic length metric. Comparison of curve lengths now proves completeness for \(g\) with \(S\) included. The derivative bounds in [lem:seeds,prop:conformal-factor] and the product rule give (35). The electrovacuum constraints identify \[\mu=\frac{|\mathcal E|_g^2+|\mathcal B|_g^2}{8\pi}=O(r^{-4}), \qquad J_{\mathrm{con}}=\frac{(\mathcal E\times_g\mathcal B)^{\flat_g}}{4\pi}.\] The fields vanish for \(r\le L\) and are radial, hence parallel, for \(r\ge2L\). Consequently \(J_{\mathrm{con}}\) is supported in the compact annulus \(L\le r\le2L\). Since \(\mathop{}\!\mathrm dV_g=u^6\mathop{}\!\mathrm dx\le M^6\mathop{}\!\mathrm dx\), both required densities are integrable. The horizon and enclosing-surface assertions are proved in [lem:positive-spheres,lem:area-minimum] below. ◻ Lemma 8. For \(L\) sufficiently large, the coordinate sphere \(S_r=\{r=\mathrm{const}\}\) has future outward expansion zero at \(r=1\) and strictly positive at every \(r>1\). No compact smooth embedded enclosing surface in \(\operatorname{int}\Omega\), with any finite number of components, is weakly future outer trapped. Proof. The normal to \(S_r\) toward infinity is \(\nu_r=u^{-2}\partial_r\). Since \(\sigma\) has only mixed radial–azimuthal components, \(K|_{TS_r\times TS_r}=0\), including on the axis by smoothness. The mean curvature is therefore also the future outward expansion, and direct divergence gives \[ \theta_+(S_r)=H_{S_r} =\frac{1}{u^6r^2}\partial_r(r^2u^4) =\frac{2W}{ru^3},\qquad W=u+2ru_r. \tag{37}\] This formula permits angular dependence of \(u\): the normal has no angular components. The Robin condition (25) gives \(W(1,\omega)=0\). For \(r\ge2\), (26) gives \[W=1-r^{-1}+v+2rv_r\ge\frac12-\frac{Cs^2}{L^2}.\] On \(1\le r\le2\), the same estimates yield \[W_r=r^{-2}+3v_r+2rv_{rr}\ge\frac14-\frac{Cs^2}{L^2}.\] Increasing \(L\) once, uniformly over \(s\in[0,1]\), we obtain \[ \begin{aligned} W_r\ge\frac18,\qquad W\ge\frac{r-1}{8} &\quad(1\le r\le2),\\ W\ge\frac14&\quad(r\ge2). \end{aligned} \tag{38}\] Thus the strict positivity holds even arbitrarily close to \(S\). Let \(\Sigma\subset\operatorname{int}\Omega\) be any compact smooth embedded enclosing surface, without a symmetry or connectedness restriction. To cover also the interpretation allowing components with interior edges, let \(\Sigma_{\mathrm{cl}}\) be the union of its components without boundary. The other components do not affect separation. Indeed, a path avoiding \(\Sigma_{\mathrm{cl}}\) can be perturbed to miss their edges and cross their interiors transversely at finitely many points. For each crossing, choose an embedded arc in that component from the crossing to an edge. A narrow tubular strip along the arc gives a detour on one side of the component, around the edge, and back on the other side, removing the crossing without introducing another. Compactness and disjointness of the components, and their location in \(\operatorname{int}\Omega\), let every detour avoid the other components and \(S\). Thus any path avoiding \(\Sigma_{\mathrm{cl}}\) can avoid all of \(\Sigma\), so \(\Sigma_{\mathrm{cl}}\) still encloses \(S\) and is nonempty. Choose \(p\in\Sigma_{\mathrm{cl}}\) where \(r\) has its global maximum \(R>1\) on \(\Sigma_{\mathrm{cl}}\). The increasing radial ray beyond \(p\) misses \(\Sigma_{\mathrm{cl}}\). Starting just beyond \(p\), the same detours past any edged components connect this side to infinity without meeting \(\Sigma\). Hence the prescribed normal at \(p\) is \(\nu_\Sigma(p)=u(p)^{-2}\partial_r\), and \(T_p\Sigma=T_pS_R\). With the convention \(H_\Sigma=\operatorname{div}_\Sigma \nu_\Sigma\), restriction of the ambient Hessian gives, at \(p\), \[0\ge\Delta_\Sigma r =\mathop{\mathrm{tr}}_{T_p\Sigma}\nabla_g^2r-H_\Sigma\,\nu_\Sigma(r) =|\nabla_g r|_g\bigl(H_{S_R}-H_\Sigma\bigr).\] For the last identity, the tangential Hessian of a defining function for a level surface equals its second fundamental form times the length of its gradient. It follows that \(H_\Sigma(p)\ge H_{S_R}(p)>0\). Tangency also gives \(\mathop{\mathrm{tr}}_\Sigma K(p)=0\), so \(\theta_+(\Sigma)(p)>0\). This excludes precisely the condition \(\theta_+\le0\) everywhere. ◻ Lemma 9. For the same choices of \(L\), every enclosing cut \(\Gamma=\partial D\) from the specified class has \(\mathop{\mathrm{Area}}_g(\Gamma)\ge\mathop{\mathrm{Area}}_g(S)\). Thus \(S\) attains the minimum enclosing area, and (36) holds. Proof. Let \(\pi:\Omega\longrightarrow\mathbb S^2\) be the radial projection \(\pi(r,\omega)=\omega\), and let \(\mathop{}\!\mathrm dA_{\mathbb S^2}\) denote the unit round area form. Define the smooth closed two-form \[\beta=\pi^*\bigl(u(1,\omega)^4\mathop{}\!\mathrm dA_{\mathbb S^2}\bigr).\] It is closed because every two-form on \(\mathbb S^2\) is closed. Its comass, namely the supremum of \(|\beta(e_1,e_2)|\) over \(g\)-orthonormal pairs, is \[ \|\beta\|_{*,g}(r,\omega) =\frac{u(1,\omega)^4}{r^2u(r,\omega)^4}\le1. \tag{39}\] Indeed, \(\beta\) annihilates radial vectors and its value on an oriented orthonormal angular pair is the displayed ratio. Furthermore, \[\partial_r(ru^2)=u(u+2ru_r)=uW>0\qquad(r>1),\] by 8; hence \(ru(r,\omega)^2\ge u(1,\omega)^2\). The comass estimate controls every tangent plane, so it does not require \(\Gamma\) to be a radial graph. For completeness, take any connected smooth codimension-zero submanifold \(D\subset\Omega\) that is closed in \(\Omega\), has interior contained in \(\operatorname{int}\Omega\), contains the entire sufficiently distant part of the end, and has compact smooth embedded two-sided intrinsic boundary \(\Gamma\). Choose \(R\) beyond that boundary and beyond the radius where \(D\) contains the entire end. Then \[D_R=D\cap\{r\le R\} \quad\text{is compact and smooth, with}\quad \partial D_R=\Gamma\sqcup S_R.\] Near \(\Gamma\) the truncation leaves \(D\) unchanged, including at contacts with \(S\); near \(S_R\) it is a regular level-set truncation in the interior of \(D\). These boundary pieces are disjoint, so \(D_R\) has no corners. Here \(D_R\) inherits the orientation of \(\Omega\), and its boundary is the full intrinsic boundary. In particular, any portion coinciding with \(S\) already belongs to \(\Gamma\) and contributes to this boundary integral. Stokes’ theorem with the induced boundary orientations gives \[\int_\Gamma\beta=-\int_{S_R}\beta =-\int_{\mathbb S^2}u(1,\omega)^4\mathop{}\!\mathrm dA_{\mathbb S^2}=-A.\] Using (39), including all components, we conclude \[A=\left|\int_\Gamma\beta\right| \le\int_\Gamma\|\beta\|_{*,g}\,\mathop{}\!\mathrm dA_g \le\mathop{\mathrm{Area}}_g(\Gamma).\] The choice \(D=\Omega\) is allowed and gives \(\Gamma=S\), so the infimum equals \(A\) and is attained. No symmetry or finite-perimeter approximation has entered this comparison. Finally, \(u(1,\omega)=2+v(1,\omega)\) with \(0\le v(1,\omega)\le Cs^2/L^2\). Since \(u\le M\), \[0\le A-64\pi =\int_{\mathbb S^2}\bigl((2+v)^4-16\bigr)\mathop{}\!\mathrm dA_{\mathbb S^2} \le4M^3\int_{\mathbb S^2}v\,\mathop{}\!\mathrm dA_{\mathbb S^2} \le C_A\frac{s^2}{L^2}.\] This proves (36). Supremum-norm continuity of \(v\) in \(s\) from 4 also proves continuity of \(A(s,L)\). ◻ Asymptotic quantities and the counterexamplesFix the cutoff of 2. Throughout this section \(L\) is above a single threshold for the preceding construction, uniformly for \(s\in[0,1]\). Write \(H\) and \(N\) for the source and Newton potential at the fixed point. All constants below are independent of these two parameters once that threshold is fixed. Mass, momentum, and chargesProposition 10. The constructed data have finite asymptotic quantities \[ P_{\mathrm{ADM}}=0,\qquad Q_e=Q_b=0, \tag{40}\] and \[ m=2+2N(0)+\frac1{2\pi}\int_{\mathbb R^3}H(x)\,\mathop{}\!\mathrm dx =2+O(s^2/L). \tag{41}\] In particular \(m\ge2\). If \(M\) is a uniform upper bound for \(u\), then \[ m-2\ge\frac{s^2}{3M^3L}. \tag{42}\] The angular momentum in (6) is \[ J=-\frac{4s^2}{15}. \tag{43}\] For \(s>0\), neither electromagnetic field is identically zero. Proof. The Cartesian decay \(K=O(r^{-3})\), with \(\tau=0\), makes the integrals in (5) of order \(O(r^{-1})\). This proves the first part of (40). The conformal weights leave both Maxwell flux forms unchanged. On every coordinate sphere, the electric flux is \[\int_{S_r}\iota_{\mathcal E}\mathop{}\!\mathrm dV_g =\int_{S_r}\mathop{}\!\mathrm d(f\,\mathop{}\!\mathrm d\phi) =2\pi\bigl(f(r,\pi)-f(r,0)\bigr)=0.\] The identical calculation for \(h\) gives zero magnetic flux. Equivalently, these are integrals of globally smooth exact two-forms over closed surfaces. The charge limits are therefore zero. Substituting \(g_{ij}=u^4\delta_{ij}\) in (4) gives \[ E_{\mathrm{ADM}} =-\frac1{2\pi}\lim_{R\to\infty} \int_{S_R}u^3u_r\,\mathop{}\!\mathrm dA_\delta. \tag{44}\] Since \(u-1=O(R^{-1})\) and \(u_r=O(R^{-2})\), replacing \(u^3\) by one has an integrated error \(O(R^{-1})\). The source \(H\) is integrable and the ordinary Newton potential satisfies the exact flux identity \[\int_{S_R}N_r\,\mathop{}\!\mathrm dA_\delta=-\int_{B_R}H\,\mathop{}\!\mathrm dx.\] For the reflected term, smoothness of \(N\) at the origin gives \[R^{-1}N(\omega/R)=\frac{N(0)}{R}+O(R^{-2}), \qquad \partial_R\bigl(R^{-1}N(\omega/R)\bigr) =-\frac{N(0)}{R^2}+O(R^{-3}).\] Together with \(u_0=1+1/r\), these identities prove the existence of the limit and the exact formula (41). In particular no first-moment assumption on \(H\) or multipole expansion of the ordinary potential is needed. The estimates of 4 give \(N(0)=O(s^2/L^2)\), \(\int H=O(s^2/L)\), and all contributions in (41) are nonnegative. Since the linear momentum is zero, \(m=E_{\mathrm{ADM}}\). On \(r\ge2L\) the cutoff is one, and \[ \widetilde E=\frac{2s\cos\theta}{r^2}\partial_r, \qquad \widetilde B=\frac{s(3\cos^2\theta-1)}{r^2}\partial_r. \tag{45}\] Using just the electric field and \(u\le M\), \[H\ge\ell u^{-3}\ge\frac{s^2\cos^2\theta}{M^3r^4} \qquad(r\ge2L).\] As \(\int_{\mathbb S^2}\cos^2\theta\,\mathop{}\!\mathrm dA=4\pi/3\), this gives \[\int H\,\mathop{}\!\mathrm dx\ge\frac{2\pi s^2}{3M^3L}.\] The mass formula proves (42). Finally \(\nu_r=u^{-2}\partial_r\), \(\mathop{}\!\mathrm dA_g=u^4\mathop{}\!\mathrm dA_\delta\), and \(K=u^{-2}\sigma\), so all conformal weights in the angular-momentum integrand cancel. The gravitational flux at every coordinate sphere is \[\begin{align*} J_{\mathrm{grav}}(r) &:=\frac1{8\pi}\int_{S_r}K(\nu_r,\eta)\,\mathop{}\!\mathrm dA_g\\ &=\frac1{8\pi}\int_{S_r}\sigma(\partial_r,\partial_\phi) \,\mathop{}\!\mathrm dA_\delta =-\frac{s^2c(r)^2}{4}\int_0^\pi\sin^5\theta\,\mathop{}\!\mathrm d\theta =-\frac{4s^2c(r)^2}{15}. \tag{46}\end{align*}\] It is constant for \(r\ge2L\), proving (43) and existence of the prescribed limit. The two nonzero angular polynomials in (45), multiplied by the positive factor \(u^{-6}\), also prove the final assertion. ◻ The conserved electromagnetic correctionThe flux (46) vanishes near the horizon but not at infinity. The following calculation identifies the boundary term responsible for the difference, using precisely our constraint signs. Proposition 11. Let \(A_B=h\,\mathop{}\!\mathrm d\phi\), the smooth magnetic potential one-form in 2. Then \[ \operatorname{div}_g\left(K(\cdot,\eta)^{\sharp_g} +2h\mathcal E\right)=0. \tag{47}\] The electromagnetic boundary contribution on every coordinate sphere is \[ J_{\mathrm{EM}}(r) :=\frac1{4\pi}\int_{S_r}h\,\iota_{\mathcal E}\mathop{}\!\mathrm dV_g =\frac{4s^2c(r)^2}{15}. \tag{48}\] Consequently \(\mathcal J(r):=J_{\mathrm{grav}}(r)+J_{\mathrm{EM}}(r)=0\) for all \(r\ge1\), whereas \(J_{\mathrm{EM}}(r)\to4s^2/15\) at infinity. Proof. Since \(A_B\) is invariant and \(A_B(\eta)=h\), Cartan’s identity gives \[\iota_\eta\mathop{}\!\mathrm dA_B=-\mathop{}\!\mathrm dh.\] Using \(\mathop{}\!\mathrm dA_B=\iota_{\mathcal B}\mathop{}\!\mathrm dV_g\) and evaluating on \(\mathcal E\) gives \[\mathcal E(h)=-\mathop{}\!\mathrm dV_g(\mathcal E,\mathcal B,\eta).\] The symmetry of \(K\) and the Killing equation imply \[\operatorname{div}_g\bigl(K(\cdot,\eta)^{\sharp_g}\bigr) =(\operatorname{div}_gK)(\eta) =2\mathop{}\!\mathrm dV_g(\mathcal E,\mathcal B,\eta).\] Combining this with \(\operatorname{div}_g\mathcal E=0\) proves (47). Directly from the flux form, \[J_{\mathrm{EM}}(r) =\frac12\int_0^\pi h f_\theta\,\mathop{}\!\mathrm d\theta =s^2c(r)^2\int_0^\pi\sin^3\theta\cos^2\theta\,\mathop{}\!\mathrm d\theta =\frac{4s^2c(r)^2}{15}.\] Now apply (46). ◻ There is no singular gauge choice in this computation: \(A_B\) is a globally smooth one-form. The scalar \(h=A_B(\eta)\) vanishes on the axis, where \(\eta=0\). An additive constant \(C\) in a scalar potential would change the boundary contribution only by \(C Q_e=0\). The norm of \(A_B\) is \(O(r^{-1})\), while \(A_B(\eta)=h\) need not tend to zero because \(\left\lvert\eta\right\rvert=O(r)\). Moreover, for \(s>0\), the zero-charge fields in (45) retain angularly varying \(r^{-2}\) terms. Multiplication by \(u^{-6}=1+O(r^{-1})\) preserves these leading terms. They satisfy (3), but they are not \(O(r^{-3})\). The divergence constraint does not prohibit this: in Euclidean space, \[\operatorname{div}_\delta\bigl(r^{-2}a(\omega)\partial_r\bigr)=0\] for every smooth angular function \(a\). No curl-free or stationarity condition is imposed. Thus 11 is consistent with electromagnetic angular-momentum conservation; it is the identification of that conserved flux with (6) which fails here. Strict violation and nondegenerate equalityThe area estimate and angular momentum already give, uniformly in the family, \[ A\ge64\pi>\frac{32\pi}{15}\ge8\pi\left\lvert J\right\rvert. \tag{49}\] Since \(Q=0\), this is the strict physical area branch. Set \[ D(s,L)=m(s,L)^2-\frac{A(s,L)}{16\pi} -\frac{4\pi J(s,L)^2}{A(s,L)}. \tag{50}\] By (36), (41), and (43), \[ \lim_{L\to\infty}D(1,L)=-\frac1{225}. \tag{51}\] Hence \(D(1,L)<0\) for all sufficiently large finite \(L\). Proposition 12. For every sufficiently large fixed \(L\), there exists \(s_*\in(0,1)\) such that \(D(s_*,L)=0\). These equality data satisfy the strict branch (49) and have no Kerr–Newman realization with the specified asymptotic parameters and induced fields. Proof. Write the uniform area estimate as \[64\pi\le A\le64\pi+C_A s^2/L^2.\] The lower mass bound (42) and \(A\ge64\pi\) imply \[ D(s,L)\ge \frac{4s^2}{3M^3L}-\frac{C_A s^2}{16\pi L^2}-\frac{s^4}{225}. \tag{52}\] Choose \(L\) above the thresholds for the previous propositions, (51), and positivity of the quadratic coefficient on the right of (52). For this fixed \(L\), the right side is positive at sufficiently small nonzero \(s\). The fixed point depends continuously on \(s\) in the supremum norm. Restriction to \(S\) therefore gives continuity of \(A\), and the integrable source bounds give continuity of \(N(0)\) and \(\int H\). The mass formula and (43) show that \(D(s,L)\) is continuous. Its positive value at a nonzero small parameter and its negative value at \(s=1\) yield a root \(s_*\in(0,1)\) by the intermediate value theorem. At this root, \(Q_e=Q_b=0\), \(J\ne0\), and both electromagnetic fields are nonzero. The strict branch and the equality equation give \[ m^2-\left\lvert J\right\rvert=\frac{(A-8\pi\left\lvert J\right\rvert)^2}{16\pi A}>0. \tag{53}\] Thus the Kerr–Newman member with these parameters would be a subextremal zero-charge member, namely Kerr with zero Maxwell tensor. Indeed the dyonic Kerr–Newman potential is linear in its electric and magnetic charges; see (Chen et al. 2018, Equation (3)). When both vanish, \(F=0\). Every spacelike slice then has \[(\iota_nF)|_{T\Omega}=0,\qquad \star_g(F|_{T\Omega})=0,\] independently of the slice or its future unit normal. This contradicts the nonzero constructed fields, and proves the asserted obstruction. ◻ Proof of 1. [lem:seeds,lem:conformal,prop:conformal-factor,prop:geometry] provide smooth maximal electrovacuum data with every topological, asymptotic, boundary, outermostness, and area condition in 1.1. 10 supplies the finite ADM and charge limits, positive mass and zero linear momentum. The physical branch is strict by (49). For \(s=1\) and sufficiently large finite \(L\), (51) gives the first conclusion. For the same large \(L\), 12 gives the second. ◻
Chen, Chiang-Mei, Sang Pyo Kim, Jia-Rui Sun, and Fu-Yi Tang. 2018. “Pair Production of Scalar Dyons in Kerr–Newman Black Holes.” Physics Letters B 781: 129–38. https://doi.org/10.1016/j.physletb.2018.03.078.
Dain, Sergio, Marcus Khuri, Gilbert Weinstein, and Sumio Yamada. 2013. “Lower Bounds for the Area of Black Holes in Terms of Mass, Charge, and Angular Momentum.” Physical Review D 88: 024048. https://doi.org/10.1103/PhysRevD.88.024048.
Khuri, Marcus, and Gilbert Weinstein. 2016. “The Positive Mass Theorem for Multiple Rotating Charged Black Holes.” Calculus of Variations and Partial Differential Equations 55 (2): 33. https://doi.org/10.1007/s00526-016-0969-8.
Penrose, Roger. 1973. “Naked Singularities.” Annals of the New York Academy of Sciences 224 (1): 125–34. https://doi.org/10.1111/j.1749-6632.1973.tb41447.x.
|
| ||||||||
|