Born–Infeld Models: A Nonlinear Framework

• Born–Infeld models are nonlinear deformations characterized by square-root or determinant structures that impose finite field scales while recovering linear behavior in weak-field limits.
• They extend classical theories to diverse settings including electrodynamics, gravity, cosmology, and reaction–diffusion, offering regularized phenomena such as bounded point-charge fields and modified wave propagation.
• Formulations using auxiliary-field linearization and symplectic embedding reveal deep symmetry properties—like self-duality and electric-magnetic duality—that distinguish these models from generic nonlinear deformations.

Searching arXiv for recent Born–Infeld-related papers to supplement the provided corpus. arxiv_search(query="Born-Infeld models review electrodynamics gravity DBI", max_results=10) Born–Infeld models are nonlinear square-root or determinant deformations of otherwise linear theories, introduced in electrodynamics to soften the large-field regime and later generalized to scalar, multifield, supersymmetric, gravitational, kinetic, cosmological, and reaction–diffusion settings. Across these realizations, the recurring structural themes are a finite field or slope scale, nonlinear constitutive relations, weak-field limits that recover the linear parent theory, and symmetry constraints—especially self-duality or special symplectic organization—that sharply distinguish Born–Infeld theories from generic nonlinear deformations (Cerchiai et al., 2016, 0812.1981, Garrione, 2023).

1. Defining structures

In nonlinear electrodynamics, a standard starting point is a Lagrangian $L=K(Y,Z)+A_\mu J^\mu$ built from the Lorentz invariants

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$

with equations of motion

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$

Within this class, the one-field Born and Born–Infeld theories differ by the presence or absence of the pseudoscalar invariant $F_{\mu\nu}{}^{*}F^{\mu\nu}$. In the conventions used for multifield constructions,

$$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$

whereas

$$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$

Both are self-dual under Legendre transform, but only Born–Infeld has the enhancement to continuous $U(1)$ electric-magnetic duality (Cerchiai et al., 2016).

A second, more general defining viewpoint is auxiliary-field linearization. For $n$ Abelian field strengths, the quadratic parent Lagrangian

$$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$

packages the nonlinear theory into a scalar-dependent symplectic matrix

$$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$

Eliminating the nondynamical matrices $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$0 and $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$1 reproduces nonlinear Born or Born–Infeld models, while their symmetry content becomes a statement about the homogeneous scalar manifold $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$2, its embedding into $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$3, and the chosen symplectic frame (Cerchiai et al., 2016).

This suggests that “Born–Infeld model” names less a single equation than a construction principle: a nonlinear completion with a finite field scale, typically organized so that weak fields reproduce the original theory while strong fields are controlled by square-root or determinant structure.

2. Nonlinear electrodynamics

In electrodynamics, Born–Infeld theory regularizes the self-field of point charges by replacing Maxwell’s linear constitutive law with a bounded nonlinear one. For a static point electric charge $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$4, the Born–Infeld field is

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$5

so $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$6 for $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$7, while near the origin the field does not diverge and saturates at order $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$8. The same theory is exceptional among generic nonlinear electrodynamics because the two effective optical metrics seen by perturbations coincide up to conformal rescaling, so the eikonal theory is no-birefringent. On the point-charge background the resulting optical geometry is

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$9

which has a throat at $$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$0. The paper emphasizing this result is explicit that the geometry is an analogue optical wormhole for perturbative photons, not a spacetime wormhole of the underlying Minkowski background, and that the same optical metric arises for a magnetic monopole and for a dyon because of Born–Infeld duality invariance (Jiménez et al., 2024).

Born–Infeld theory also appears as a distinguished point inside wider Born–Infeld-type families. A three-parameter nonlinear electrodynamics model with

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$1

contains standard Born–Infeld electrodynamics at

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$2

and exponential electrodynamics in the limit $$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$3. In this family, finite point-charge fields and finite electrostatic self-energy persist for $$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$4, but exact no-birefringence and exact electric-magnetic duality survive only at the Born–Infeld point. The exact all-orders no-birefringence condition reduces to the unique solution

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$5

so Born–Infeld is singled out not merely by regularized electrostatics but by its symmetry structure (Kruglov, 2016).

Exact electrostatic solutions in simple geometries make the regularization mechanism concrete. In SI units, Abelian Born–Infeld electrostatics with parameter $$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$6 obeys

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$7

with energy density

$$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$8

For an infinite line of charge density $$\partial_\nu\!\left(K_Y F^{\mu\nu}+K_Z\widetilde F^{\mu\nu}\right)=J^\mu, \qquad \partial_\mu \widetilde F^{\mu\nu}=0.$$9,

$F_{\mu\nu}{}^{*}F^{\mu\nu}$0

so the Maxwell divergence is replaced by the finite limit $F_{\mu\nu}{}^{*}F^{\mu\nu}$1 as $F_{\mu\nu}{}^{*}F^{\mu\nu}$2. For an infinitely long uniformly charged cylinder, the field is likewise bounded and reduces to the Maxwell expression only in the large-$F_{\mu\nu}{}^{*}F^{\mu\nu}$3 limit (Moayedi et al., 2017).

A complementary two-dimensional formulation uses complex analysis. In the complex plane $F_{\mu\nu}{}^{*}F^{\mu\nu}$4, the electrostatic equations admit a complex potential $F_{\mu\nu}{}^{*}F^{\mu\nu}$5 with

$F_{\mu\nu}{}^{*}F^{\mu\nu}$6

and a holomorphic seed $F_{\mu\nu}{}^{*}F^{\mu\nu}$7 or $F_{\mu\nu}{}^{*}F^{\mu\nu}$8 such that

$F_{\mu\nu}{}^{*}F^{\mu\nu}$9

This reproduces the Coulombian complex potential in the weak-field limit while accommodating the Born–Infeld bound $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$0. The construction yields explicit monopolar and multipolar solutions, and for two equal but opposite charges it leads to an intrinsically Born–Infeld effect: the attractive force is lower than its Coulombian value and decreases to zero when the charges approach each other below a distance controlled by the Born–Infeld constant (Ferraro, 2010).

3. Scalar, multifield, and supersymmetric realizations

Dirac–Born–Infeld scalar theories replace the canonical kinetic term by a square root while preserving a first-order structure for appropriate choices of the potential. In the $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$1-dimensional models studied through

$$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$2

the first-order relation

$$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$3

implies topological energy

$$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$4

and the fluctuation operator factorizes as

$$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$5

so $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$6. In the explicit $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$7-, $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$8-, sine-Gordon-, double-sine-Gordon-, and multi-sine-Gordon-like examples, the kink profiles often coincide with the canonical ones, while the DBI parameter $$\mathcal L_{\mathrm{Born}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu}}\right),$$9 alters the potential, the energy density, and the stability potential (Bazeia et al., 2017).

A cosmological generalized DBI model promotes the constant $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$0 of generalized Chaplygin gas Born–Infeld matter to a field-dependent quantity $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$1, giving

$$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$2

$$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$3

This unifies the rolling tachyon limit $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$4 and the generalized Chaplygin gas limit $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$5. Because the varying potential induces $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$6-variation,

$$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$7

the potential-driven deviation from $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$8 is constrained much more strongly than the Chaplygin component. The reported $$\mathcal L_{\mathrm{BI}}=\mu^2\left(1-\sqrt{1+\frac{1}{2\mu^2}F_{\mu\nu}F^{\mu\nu} -\frac{1}{16\mu^4}(F_{\mu\nu}{}^{*}F^{\mu\nu})^2}\right).$$9 confidence-level bounds are

$U(1)$0

and, alternatively,

$U(1)$1

The paper’s conclusion is that the potential must be extremely flat (Tavares et al., 2021).

Multifield Born and Born–Infeld theories emerge naturally from the auxiliary symplectic framework. Choosing the scalar manifold $U(1)$2 yields the $U(1)$3-covariant multifield Born–Infeld action

$U(1)$4

while the diagonal embedding $U(1)$5 gives a new $U(1)$6-field Born theory,

$U(1)$7

with manifest $U(1)$8 symmetry and Legendre self-duality (Cerchiai et al., 2016).

Supersymmetric generalizations place the nonlinear constraints under the control of $U(1)$9 special geometry. In the $n$0 models built from a cubic prepotential

$n$1

partial $n$2 breaking leads, in the nonlinear limit, to the tensorial constraint

$n$3

The coefficients $n$4 classify inequivalent multifield Born–Infeld systems, while the vacuum values are fixed by attractor equations

$n$5

This construction shows that coupled multifield supersymmetric Born–Infeld theories are dictated by special geometry rather than by arbitrary nonlinear couplings (Ferrara et al., 2014).

4. Gravitational Born–Infeld models

Born–Infeld ideas enter gravity in several inequivalent ways. In teleparallel Born–Infeld gravity, the starting point is the torsion scalar $n$6 of the Teleparallel Equivalent of General Relativity, which depends only on first derivatives of the vielbein. Replacing the TEGR Lagrangian by a Born–Infeld-type square root preserves second-order field equations and yields explicit modified solutions. In $n$7 dimensions, the BTZ sector is deformed mainly through an effective cosmological constant

$n$8

so that BTZ black holes can exist even when the original $n$9. In spatially flat FRW cosmology with matter, the modified Friedmann equation bounds the Hubble rate,

$$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$0

and the early-time big-bang singularity is replaced by a past-eternal de Sitter-like phase with bounded curvature invariants (0812.1981).

A different nonrelativistic direction is Born–Infeld–Hořava gravity, where determinant-based spatial-curvature potentials are constructed so that their small-curvature expansion reproduces Hořava gravity at quadratic order. The exact actions contain infinitely many higher-spatial-curvature terms and are described as $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$1 extensions, whereas truncations produce finite-$$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$2 models, including half-integer values because Cotton-tensor terms contribute odd numbers of spatial derivatives. The direct $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$3-dimensional action

$$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$4

was chosen to reproduce the Hořava potential through $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$5, but the paper also records a striking exact result: in the minimal model, static spherically symmetric solutions are ruled out (Gullu et al., 2010).

Three-dimensional determinant completions of new massive gravity display yet another pattern. For the Born-Infeld extension of NMG on an $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$6 background, the transverse-traceless linearized equation factorizes as

$$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$7

so the theory propagates a massless and a massive graviton. Unlike TMG, NMG, or GMG, however, the would-be critical point $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$8 requires the singular parameter value $$\mathcal{L} = -\frac14\,F_{\mu\nu}^{\,T}\, g\, F^{\mu\nu} + \frac14\,F_{\mu\nu}^{\,T}\,\theta\,{}^{*}F^{\mu\nu} -\frac{\mu^2}{2}\,\mathrm{Tr}(N\mathcal{M})+\text{const.}$$9, so pure Born–Infeld gravity has no regular critical point and no regular logarithmic bulk modes (Setare et al., 2014). In the related analysis of $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$0-dimensional Born–Infeld gravity and its Chern–Simons extension, pure Born–Infeld gravity again has only a limiting logarithmic solution as $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$1, whereas Born–Infeld–Chern–Simons gravity admits a genuine logarithmic AdS-wave solution along the chiral line

$$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$2

where the left central charge vanishes (Alishahiha et al., 2010).

Eddington-inspired Born–Infeld gravity provides a metric-affine realization in which the physical metric $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$3 and an auxiliary metric $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$4 are related algebraically through matter. Coupling nonlinear $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$5-models to EiBI gravity yields “minimal modifications” of the corresponding GR geometries in the sense that the large-distance form remains close to the GR global-monopole or Reissner–Nordström-with-deficit-angle solution, while the interior geometry changes qualitatively. Wormhole structures always arise, but the paper is explicit that this does not guarantee geodesic completeness. For the quadratic matter model $$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$6, a tuned mass–charge relation produces a subset of solutions that are regular everywhere and geodesically complete (Nascimento et al., 2019).

5. Reaction–diffusion, blow-up, and kinetic Born–Infeld systems

Born–Infeld nonlinearities also appear in non-electromagnetic partial differential equations. A reaction–diffusion model driven by the one-dimensional Born–Infeld, or Minkowski-curvature, operator

$$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$7

inherits the structural gradient bound

$$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$8

and admits a first-order reduction via

$$\mathcal{M}[g,\theta] = \begin{pmatrix} g+\theta g^{-1}\theta & -\theta g^{-1}\ -g^{-1}\theta & g^{-1} \end{pmatrix}\in Sp(2n).$$9

which converts traveling fronts into a scalar two-point problem with bounded reduced flux

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$00

This boundedness produces several asymptotic behaviors not present for linear diffusion. In the varying-field-strength regime $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$01, $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$02, the large-field limit $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$03 rigorously recovers the Maxwell or linear-diffusion front speed and profile in $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$04. In the singular perturbation regime $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$05,

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$06

the critical speed does not vanish: $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$07 and the limiting front becomes one-sided sharp, given by the $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$08-gluing of a piecewise linear branch of slope $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$09 and an inviscid branch solving $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$10. This is presented as a new phenomenon specific to Born–Infeld diffusion (Garrione, 2023).

A scalar hyperbolic Born–Infeld equation,

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$11

arises as the timelike minimal-surface equation for graphs in Lorentz–Minkowski space. In one spatial dimension it admits the explicit self-similar blow-up family

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$12

and the same family also solves the linear wave equation after similarity reduction. The paper proves Lyapunov nonlinear stability of these timelike self-similar blow-up solutions inside a strictly proper subset of the backward light cone, using weighted energy estimates and a Nash–Moser iteration (Yan, 2018).

The low-dimensional Vlasov–Born–Infeld system couples collisionless matter to nonlinear electromagnetic fields in one-and-one-half dimensions. After introducing angular variables

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$13

the Born–Infeld field subsystem becomes the diagonal quasilinear system

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$14

with characteristic speeds

$$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$15

This reveals a strictly hyperbolic, linearly degenerate structure and leads to local existence and uniqueness of $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$16 solutions under compact-support and smallness assumptions. The main obstruction to global theory is that, unlike Maxwell theory, Born–Infeld characteristic speeds depend on the solution and may resonate with particle velocities (Lee, 2015).

6. Recurring mechanisms and limits of the Born–Infeld paradigm

Several mechanisms recur across these disparate models. First, bounded field or slope scales are central: $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$17 in electrostatics, $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$18 for smooth Born–Infeld diffusion fronts, and $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$19 at the 2D electrostatic singular set in the complex-plane construction (Moayedi et al., 2017, Garrione, 2023, Ferraro, 2010). Second, weak-field or low-curvature limits recover the underlying linear theory, as in the Maxwell limit of Born–Infeld electrodynamics, the TEGR limit of teleparallel Born–Infeld gravity, and the canonical limit $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$20 of scalar DBI kink models (Jiménez et al., 2024, 0812.1981, Bazeia et al., 2017).

At the same time, the literature repeatedly stresses that Born–Infeld regularization is selective rather than universal. Born–Infeld electrodynamics regularizes the electric field of a point charge, but the effective optical wormhole metric of perturbations remains singular at the center for the ideal point-source background (Jiménez et al., 2024). EiBI gravity generically produces wormhole structures, but geodesic completeness requires additional conditions and may fail in untuned branches (Nascimento et al., 2019). Three-dimensional Born–Infeld gravity removes the regular critical/logarithmic point familiar from some $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$21 models rather than preserving it (Setare et al., 2014). Born–Infeld reaction–diffusion fronts can retain a strictly positive limiting speed in a singular perturbation regime where linear and saturating diffusions instead slow to zero (Garrione, 2023).

A common misconception is therefore that every Born–Infeld deformation simply “regularizes singularities.” The comparative evidence is more precise. What Born–Infeld structure reliably supplies is a nonlinear scale that constrains constitutive response, modifies characteristic propagation, and often preserves a distinguished symmetry pattern—such as no birefringence, $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$22 or $$Y=-\frac14 F_{\mu\nu}F^{\mu\nu},\qquad Z=-\frac14 F_{\mu\nu}\widetilde F^{\mu\nu},$$23 duality, or symplectic self-duality. Whether this leads to finite self-energy, bounded curvature, a wormhole throat, or a positive selected front speed depends on the sector, the choice of auxiliary geometry or constitutive manifold, and the matter content (Kruglov, 2016, Cerchiai et al., 2016).

Taken together, these results suggest that Born–Infeld models form a coherent but heterogeneous class. Their unity lies in nonlinear completion by square-root or determinant structure and in the finite-field philosophy inherited from electrodynamics; their diversity lies in the different ways that bounded response, duality, and nonlinearity reorganize dynamics in electrodynamics, scalar theory, supersymmetry, gravity, kinetic theory, and nonlinear diffusion.