Gauged Skyrme-Maxwell-Einstein Models

• The model defines a 3+1-dimensional field theory where nonlinear chiral Skyrme fields, a U(1) gauge field, and Einstein gravity interact to stabilize topological solitons.
• It employs solution-generating techniques from Einstein-scalar-Maxwell systems to construct analytic, rotating, and magnetically deformed baryonic solutions.
• Applications include modeling electrically charged Skyrmions and non-topological pion stars, offering new insights into compact astrophysical objects and charge quantization.

The gauged Skyrme-Maxwell-Einstein models constitute a class of (3+1)-dimensional field theories in which nonlinear chiral Skyrme fields, a $U(1)$ Maxwell gauge field, and Einstein gravity interact dynamically. These models extend the canonical Skyrme framework, introducing gauge coupling and gravitational self-interaction, enabling the study of strongly gravitating, rotating, and magnetized baryonic matter. Recent developments establish precise correspondences with Einstein-scalar-Maxwell systems, substantially broadening analytic and numeric access to exact solutions, both topological and non-topological. Such configurations—including electrically charged, magnetically deformed Skyrmions and Q-ball-like pion stars—illuminate cosmological and astrophysical scenarios where baryonic charge, gauge fields, and gravity are dominant.

1. Model Definition and Field Content

The four-dimensional gauged Skyrme–Maxwell–Einstein action is constructed as: $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$ where $R$ is the Ricci scalar, $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$ is the Maxwell field strength, $\Sigma_{\mu}=U^{-1}D_{\mu}U$ and $B_{\mu\nu}=[\Sigma_{\mu},\Sigma_{\nu}]$ define the Skyrme current and quartic term, with $U(x)\in SU(2)$ parametrized by three angles $(\Psi,\Theta,\Phi)$. $D_{\mu}U$ is a gauge-covariant derivative and $K = f_\pi^2/4$ with $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$0. The Maxwell gauge coupling minimally couples $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$1 via commutator structure associated with $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$2. The Skyrme sector contains both quadratic and quartic kinetic terms, achieving stabilization of topological solitons.

In related studies, e.g., (Kirichenkov et al., 2023), the chiral field is represented as $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$3, subject to $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$4, and the Maxwell coupling generically distinguishes charged pion modes through a charge matrix $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$5. The Lagrangian further incorporates a pion mass potential term.

2. Ansatz and Sector Reduction

A fundamental simplification arises by imposing a sector-wise ansatz: $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$6 Under this construction, all non-Abelian currents and quartic Skyrme commutators vanish: $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$7, $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$8, rendering the baryonic current $$S = \int d^{4}x\,\sqrt{-g} \left\{\frac{1}{4}R - \frac{1}{4}F_{\mu\nu}F^{\mu\nu} + \frac{K}{4}\,\mathrm{Tr}\left[\Sigma^{\mu}\Sigma_{\mu} + \frac{\lambda}{8}B^{\mu\nu}B_{\mu\nu}\right]\right\}$$9 identically zero. The system reduces to Einstein gravity coupled to a free Maxwell field and a minimally coupled massless scalar, i.e., the Einstein-scalar-Maxwell system: $R$0 with stress-energy: $R$1 The baryonic (topological) charge density is nontrivial via the Callan–Witten term: $R$2 This means any scalar profile with nonvanishing derivative along magnetic lines corresponds to a nonzero baryonic charge in the uplifted Skyrme-Maxwell-Einstein configuration (Canfora et al., 25 Jan 2026).

3. Solution-Generating Techniques and Analytic Construction

The reduction to Einstein-scalar-Maxwell enables the direct application of solution-generating methodology developed for electrovacua—Ernst potentials, Ehlers and Harrison transformations, inverse scattering, and symmetry-based algorithms. Every axisymmetric electrovacuum or scalar-electrovacuum solution with nontrivial scalar gradient can be uplifted via the above field dictionary to a genuine gauged Skyrme–Maxwell–Einstein solution with quantized baryonic charge.

For instance, a Kerr–Newman–like metric with scalar dressing is specified by

$R$3

with $R$4 and $R$5 as specified in detail, and a conformal factor $R$6 encoding back-reaction. The total baryonic charge is given by: $R$7 and quantization of charge enforces discrete allowed values of the rotation parameter $R$8: $R$9 An upper bound for $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$0 is set by the maximal allowed $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$1 for which $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$2 remains real.

This route enables a transfer of the full analytic machinery of known scalar-Maxwell backgrounds into the topologically nontrivial sector of Skyrme-Maxwell-Einstein theory, effectively mapping solution branches, quantization phenomena, and rotation-charge relations (Canfora et al., 25 Jan 2026).

4. Numerical Solutions: Skyrmions and Pion Stars

Comprehensive numerical investigations (Kirichenkov et al., 2023) address two broad classes:

• Topological solutions (Skyrmions, $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$3): Electrically charged, magnetically deformed objects characterized by toroidal magnetic flux and two distinct solution branches in gravitational coupling $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$4. The lower branch, connected to flat-space Skyrmions, increases in mass with $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$5; the upper, in strong gravity ($$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$6), connects to Bartnik–McKinnon-type solutions. Gauging explicitly breaks spherical symmetry due to magnetic field topology. Electromagnetic energy lifts the mass branches and alters domain boundaries.
• Non-topological solutions (pion stars, $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$7): Self-gravitating Q-ball-like structures without a flat-space limit, appearing in curved spacetime for $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$8. The mass and electric charge rise as frequency decreases, with solution branches terminating in spiral fashions at critical $$F_{\mu\nu}=\partial_{\mu}A_{\nu}-\partial_{\nu}A_{\mu}$$9. For nonzero gauge coupling, strong gravity leads to singular solutions, and critical values grow with $\Sigma_{\mu}=U^{-1}D_{\mu}U$0.

Both types are constructed via ansätze for metric, Skyrme, and Maxwell fields, yielding ODEs or PDEs solved by Newton–Raphson collocation methods (accuracy $\Sigma_{\mu}=U^{-1}D_{\mu}U$1–$\Sigma_{\mu}=U^{-1}D_{\mu}U$2) across parameter grids.

5. Physical Significance, Charges, and Domains of Existence

The domain structure is controlled by $\Sigma_{\mu}=U^{-1}D_{\mu}U$3:

• Skyrmions ($\Sigma_{\mu}=U^{-1}D_{\mu}U$4) possess both baryonic and electromagnetic charges, plus magnetic dipole moments. Maximal allowed gravitational coupling increases with frequency and gauge strength, with solution branches meeting at $\Sigma_{\mu}=U^{-1}D_{\mu}U$5.
• Pion stars ($\Sigma_{\mu}=U^{-1}D_{\mu}U$6) are analogous to boson stars but lack flat-space Q-ball limits due to the quartic Skyrme term. Only one branch exists, which emerges from vacuum in curved space.
• Both configurations display universal spiraling in $\Sigma_{\mu}=U^{-1}D_{\mu}U$7 space.

Quantization: For uplifted Kerr–Newman-like solutions, baryonic charge quantization implies discrete allowed rotation, and there exists an upper $\Sigma_{\mu}=U^{-1}D_{\mu}U$8 determined by model integration constants.

Boundary conditions and charges conform to regularity at the origin and asymptotic flatness ($\Sigma_{\mu}=U^{-1}D_{\mu}U$9); standard constraints (e.g., $B_{\mu\nu}=[\Sigma_{\mu},\Sigma_{\nu}]$0 on symmetry axes) prevent singularities.

6. Outlook and Research Directions

Gauged Skyrme–Maxwell–Einstein models illuminate astrophysical scenarios involving compact baryonic objects with strong electromagnetic and gravitational fields. Topological Skyrmion stars may serve as baryonic models of electrically and magnetically charged, rotating compact stars. Non-topological pion stars extend the paradigm of boson stars, stabilized by nonlinear self-interactions.

Systematic utilization of solution-generating techniques now available through the scalar-Maxwell correspondence enables exact construction of further analytic and rotating solutions. Prospective research avenues include extensions to multi-Skyrmion systems, hairy black holes, time-dependent Q-clouds, and exploration of holographic or nuclear physics implications. Numerical evidence supports a generic spiraling mechanism in gravitating soliton systems involving both quadratic and quartic kinetic terms under gauge and gravitational interaction (Canfora et al., 25 Jan 2026, Kirichenkov et al., 2023).

A plausible implication is the accessibility of genuinely analytic, highly magnetized, rotating baryonic matter configurations in General Relativity, potentially relevant to strongly gravitating astrophysical objects.