Quasilinear Maxwell System

• Quasilinear Maxwell System is a set of first-order hyperbolic PDEs featuring nonlinear, anisotropic constitutive laws that govern complex electromagnetic interactions.
• It incorporates detailed boundary and interface conditions, such as perfect conductor, absorbing, and transmission types, to model realistic physical scenarios.
• Advanced techniques like energy estimates, microlocal analysis, and dispersive methods ensure local well-posedness, global existence, and exponential decay.

The quasilinear Maxwell system refers to initial-boundary value problems for the macroscopic Maxwell equations in media with instantaneous, nonlinear, and typically anisotropic constitutive relations between the electromagnetic fields. In such systems, the constitutive laws (relating $\mathbf{D}$ to $\mathbf{E}$ and $\mathbf{B}$ to $\mathbf{H}$) are nonlinear functions of the fields, so the resulting PDE system is genuinely quasilinear, first-order, and hyperbolic. The main mathematical and physical distinction from the linear Maxwell system is that material parameters such as permittivity and permeability depend on the magnitude or direction of the underlying fields, leading to complex nonlinear wave interactions, new stability and blow-up phenomena, and challenging issues of well-posedness at both the local and global levels.

1. Mathematical Formulation and Constitutive Laws

The macroscopic quasilinear Maxwell system, in a domain $G\subset\mathbb{R}^3$, is governed by

$$\begin{aligned} \partial_t \mathbf{D} &= \nabla \times \mathbf{H} - \mathbf{J}\ \partial_t \mathbf{B} &= - \nabla \times \mathbf{E}\ \nabla \cdot \mathbf{D} &= \rho,\quad \nabla \cdot \mathbf{B}=0 \end{aligned}$$

with nonlinear, instantaneous constitutive relations

$$\mathbf{D} = \mathbf{P}(x,\mathbf{E},\mathbf{H}),\quad \mathbf{B} = \mathbf{M}(x,\mathbf{E},\mathbf{H}),$$

where $\mathbf{P},\mathbf{M} \in C^m$ in $x$ and $(\mathbf{E},\mathbf{H})$, and the associated susceptibility tensor

$$\chi = \partial_{(\mathbf{E},\mathbf{H})} (\mathbf{P},\mathbf{M})$$

is uniformly positive-definite and symmetric. The current typically takes the form $$\mathbf{J} = \mathbf{J}_0 + \sigma(\mathbf{E},\mathbf{H})\mathbf{E}$$ with $\sigma$ a bounded, potentially field-dependent conductivity. For systematization, these relations can be written as a first-order quasilinear symmetric hyperbolic system,

$$\chi(\mathbf{u})\partial_t \mathbf{u} + \sum_{j=1}^3 A_j^{\text{co}} \partial_j \mathbf{u} + \sigma(\mathbf{u}) \mathbf{u} = f$$

where $\mathbf{u} = (\mathbf{E},\mathbf{H})^T$ and the $A_j^{\text{co}}$ are the standard Maxwell matrices encoding the curl operation (Spitz, 2018, Schnaubelt et al., 2018, Schnaubelt et al., 2018).

2. Boundary and Interface Conditions

Boundary and interface conditions play a decisive role in the PDE theory.

• Perfect conductor: On a perfectly conducting boundary, the tangential electric field vanishes, i.e., $\mathbf{E} \times \nu = 0$, and the normal component of $\mathbf{B}$ vanishes, i.e., $\mathbf{B} \cdot \nu = 0$ (Spitz, 2018).
• Absorbing/Impedance boundary (Silver–Müller): On boundaries with energy absorption, the condition is $$\mathbf{H} \times \nu + [\lambda(\mathbf{E} \times \nu)]\times \nu = 0$$, where $\lambda$ is a symmetric positive-definite tangential tensor, possibly nonlinear in the fields (Pokojovy et al., 2018, Nutt et al., 14 Jan 2026).
• Transmission/interface: At material interfaces, one imposes continuity of the tangential electric field and normal induction: $$[\mathbf{E}_\text{tan}] = 0, \quad [\mathbf{D}_\text{nor}] = -\rho_\Sigma,\quad [\mathbf{B}_\text{nor}] = 0, \quad [\mathbf{H}_\text{tan}] = \mathbf{J}_\Sigma$$ with $\rho_\Sigma$ and $\mathbf{J}_\Sigma$ as surface charge and current densities (Schnaubelt et al., 2018, Dohnal et al., 2021).

3. Local Well-Posedness and Regularity Results

The canonical functional setting is $H^m$- or $\mathcal{H}^m$-valued Sobolev spaces with $m\geq3$, in which the maximal regularity required by the nonlinearities and Sobolev embedding holds. The main local well-posedness theorem states:

Let the constitutive tensors $\chi,\,\sigma$ be $C^m$ and satisfy the uniform positivity/symmetry requirements, and let the initial data and inhomogeneity be of class $H^m$ and compatible up to order $m$. Then there exists a unique maximal solution $u \in \bigcap_{j=0}^m C^j((T_-,T_+); H^{m-j}(G))$ to the quasilinear Maxwell system—this solution depends continuously on the data and boundary sources, with lifespan determined either by the approach of $u$ to the boundary of the constitutive domain or by blow-up of the Lipschitz norm (Spitz, 2018, Schnaubelt et al., 2018, Schnaubelt et al., 2018). Explicitly, the blow-up/continuation criterion is: $$\text{if } T_+<\infty \text{ and } \text{dist}(u(t),\partial U)>c>0 \quad \forall t<T_+, \text{ then } \|\nabla_x u(t)\|_{L^\infty(G)} \to \infty \text{ as } t\uparrow T_+\,.$$ Analogous results hold with conservative interface or absorbing boundary conditions, with smallness or regularity conditions adapted to the nonlinear character of those constraints (Schnaubelt et al., 2018, Schnaubelt et al., 2018).

4. Global Existence, Stabilization, and Exponential Decay

Under appropriate dissipativity—either via interior conductivity or absorbing boundary—and smallness of initial data in $H^3$, global existence and exponential decay can be achieved:

• For quasilinear Maxwell with absorbing boundary, if the domain is strictly star-shaped and material tensors are $C^3$ and positively definite, then classical $H^3$-solutions exist globally for sufficiently small initial data. The main mechanism is a stabilized energy identity: $E(t) + \int_0^t D(s)\, ds \leq C E(0)$ where $E$ is a high-order energy and $D$ the boundary dissipation, yielding

$$\max_{0\leq j\leq 3} \left( \|\partial_t^j E(t)\|_{H^{3-j}}^2 + \|\partial_t^j H(t)\|_{H^{3-j}}^2 \right) \leq M e^{-\omega t} \|(E^{(0)},H^{(0)})\|_{H^3}^2$$

(Pokojovy et al., 2018, Nutt et al., 14 Jan 2026).

• For interior damping (positive conductivity) and perfect-conductor boundary, small data likewise yield global existence and exponential decay by a barrier method combining high-order energy and observability-type inequalities: $z(t) + \int_s^t z(\tau)\, d\tau \leq C z(s)$ where $z(t)$ is the complete $H^3$-norm, with exponential decay following by a discrete Grönwall argument (Lasiecka et al., 2018).
• In the quasilinear transmission/interface problem, solutions enforcing div-free initial displacement (for divergence constraints in nonlinear dielectrics) are constructed variationally as corrections to approximate wave-packet ansatzes, yielding asymptotically small corrections in the multiple-scale regime relevant for nonlinear optics (Dohnal et al., 2021).

5. Analysis of Two-Dimensional Quasilinear Maxwell and Dispersive Effects

In two spatial dimensions, Strichartz-type dispersive estimates have been established for the quasilinear Maxwell system with rough (as low as $C^2$) permittivity, exploiting FBI phase-space conjugation and microlocal diagonalization. The system can be reduced to a pair of half-wave equations via diagonalization of the principal symbol. This yields sharp dispersive and nonlinear well-posedness results, most notably local well-posedness in $H^s$ for $s>11/6$ for nonlinear permittivities with Kerr-type response (Schippa et al., 2021). These results directly tie the regularity theory for quasilinear Maxwell to that of nonlinear wave equations with rough coefficients, illustrating the influence of the quasilinear structure on dispersive and regularity properties.

6. Quasilinear Maxwell Systems in Geometric and General Relativistic Settings

The Maxwell–Born–Infeld (MBI) system is a prototype genuinely quasilinear electromagnetic model arising in nonlinear electrodynamics and string-theoretic contexts. On curved backgrounds such as Schwarzschild spacetime, the system is governed by a quasilinear Lagrangian density and the field equations

$$\nabla_{[\alpha}F_{\beta\gamma]}=0,\qquad \nabla^\mu \left( (1+F-G^2)^{-1/2}[F_{\mu\nu}-G\,{}^*F_{\mu\nu}] \right)=0,$$

where $F,\,G$ are field invariants (Pasqualotto, 2017). The quasilinear nature manifests in the principal part being exactly Maxwellian but with crucial nonlinear coupling. The MBI system admits global existence and quantitative decay for small initial data via a combination of canonical stress energies, bootstrap assumptions, and the $r^p$ method adapted from hyperbolic PDE theory on black hole backgrounds, with the nonlinear error terms treated perturbatively using the structure of the system and null decompositions. The methods generalize to broader classes of quasilinear, gauge-invariant Maxwell-type systems provided the principal part is hyperbolic and appropriate energy structure is present (Pasqualotto, 2017).

7. Techniques and Structural Phenomena

Across the different settings, several methodological advances characterize the contemporary analysis of the quasilinear Maxwell system:

• Energy and a priori estimates: Fundamental at linear and nonlinear levels, employing the symmetric hyperbolic structure, tangential/normal regularity splitting, and Moser-type calculus (Spitz, 2018, Pokojovy et al., 2018, Lasiecka et al., 2018).
• Barrier and bootstrapping methods: Employed for global existence and exponential decay, combining energy-dissipation and observability-type inequalities, and regularity bootstrapping to handle the quasilinear error terms (Pokojovy et al., 2018, Lasiecka et al., 2018, Nutt et al., 14 Jan 2026).
• Normal trace estimates: New microlocal and functional analytic approaches for control of normal traces are crucial for closing energy-type and observability inequalities, especially in the presence of (nonlinear) absorbing boundary conditions (Nutt et al., 14 Jan 2026).
• Compatibility conditions: For preservation of constraints and full $H^m$-wellposedness, compatibility up to order $m$ in initial and boundary data is necessary (Spitz, 2018, Schnaubelt et al., 2018).
• Dispersive/FBI techniques in low dimensions: The application of FBI phase-space conjugation, dyadic parametrix reduction, and sharp Strichartz estimates extends to quasilinear Maxwell systems, proving regularity thresholds consistent with scalar rough-coefficient wave theory (Schippa et al., 2021).

These features highlight the deep interplay between hyperbolic PDE theory, nonlinear functional analysis, geometric insights, and material science in the study of quasilinear Maxwell systems.

---

References:

• (Spitz, 2018) – Local wellposedness of nonlinear Maxwell equations with perfectly conducting boundary conditions
• (Pokojovy et al., 2018) – Boundary Stabilization of Quasilinear Maxwell Equations
• (Lasiecka et al., 2018) – Exponential Decay of Quasilinear Maxwell Equations with Interior Conductivity
• (Schnaubelt et al., 2018) – Local wellposedness of quasilinear Maxwell equations with conservative interface conditions
• (Schnaubelt et al., 2018) – Local wellposedness of quasilinear Maxwell equations with absorbing boundary conditions
• (Dohnal et al., 2021) – A quasilinear transmission problem with application to Maxwell equations with a divergence-free $\mathcal D$-field
• (Nutt et al., 14 Jan 2026) – Normal trace inequalities and decay of solutions to the nonlinear Maxwell system with absorbing boundary
• (Schippa et al., 2021) – On quasilinear Maxwell equations in two dimensions
• (Pasqualotto, 2017) – Nonlinear stability for the Maxwell-Born-Infeld system on a Schwarzschild background