Anisotropic Linear Maxwell System

• The anisotropic linear Maxwell system is defined by tensorial constitutive relations that account for directional permittivity, permeability, and magnetoelectric coupling.
• It is formulated in a variational H(curl) framework and discretized using Nédélec edge elements for reliable finite element approximations.
• Numerical experiments validate convergence and well-posedness, making the approach applicable to complex metamaterials and advanced photonic structures.

Anisotropic Linear Maxwell System

The anisotropic linear Maxwell system generalizes the classical Maxwell equations to account for the directional dependence (anisotropy) of material parameters such as permittivity and permeability, and permits magnetoelectric coupling in bianisotropic media. This framework underpins the analysis of electromagnetic phenomena in engineered materials, metamaterials, and complex natural media. The formulation encompasses fully tensorial constitutive relations, nontrivial polarization and magnetization fields, possible magnetic charge and current densities, and is foundational for robust finite element and other numerical methods (Fernando et al., 2022).

1. Fundamental Equations and Constitutive Structure

The time-harmonic anisotropic Maxwell equations in phasor notation specify the electric field $E(r)$, magnetic field $H(r)$, electric displacement $D(r)$, and magnetic flux density $B(r)$ as follows: \begin{align} \nabla\cdot D &= \rho_E, \tag{1a}\ \nabla\cdot B &= \rho_M, \tag{1b}\ \nabla\times H &= i\,\omega\,D + J_E, \tag{1c}\ -\,\nabla\times E &= i\,\omega\,B + J_M. \tag{1d} \end{align} Here, $J_E$ and $J_M$ denote impressed electric and magnetic current densities, $\rho_E$ and $\rho_M$ the corresponding charge densities, and $\omega$ is the angular frequency.

General linear constitutive relations for bianisotropic media are \begin{align} D &= \epsilon\,E + \xi\,H + P, \tag{3a}\ B &= \zeta\,E + \mu\,H + M, \tag{3b} \end{align} where $\epsilon(r)$, $\mu(r)$ are symmetric positive-definite tensor fields (permittivity and permeability), $\xi(r), \zeta(r)$ are magnetoelectric coupling tensors, and $P(r), M(r)$ are prescribed polarization and magnetization fields (Fernando et al., 2022).

These equations subsume standard and generalized media, including:

• Anisotropic dielectrics and magnetics ($\epsilon, \mu$ nonproportional to the identity)
• Magnetoelectric media with nonzero off-diagonal $\xi, \zeta$
• Media permitting polarization and magnetization densities, as well as magnetic charge and current

The continuity (charge conservation) constraints inherent in the system read

$$i\,\omega\,\rho_E + \nabla\cdot J_E = 0, \qquad i\,\omega\,\rho_M + \nabla\cdot J_M = 0. \quad\tag{2}$$

2. Variational Formulation and Functional Analysis

To enable rigorous analysis and numerical approximation, the system is recast in a variational (weak) form. Eliminating $D, B$ via constitutive laws, the time-harmonic system becomes a coupled system for the primary fields $(E, H)$: \begin{align} i\,\omega\,\epsilon\,E + i\,\omega\,\xi\,H - \nabla\times H &= -\left(J_E + i\,\omega\,P\right), \tag{4a}\ i\,\omega\,\zeta\,E + i\,\omega\,\mu\,H + \nabla\times E &= -\left(J_M + i\,\omega\,M\right). \tag{4b} \end{align}

The natural energy space for these fields is

$$H(\mathrm{curl}; \Omega) = \left\{ v \in L^2(\Omega)^3 \mid \nabla\times v \in L^2(\Omega)^3 \right\},$$

with appropriate tangential-trace subspaces for the imposed boundary conditions.

The weak form: for auxiliary fields implementing nonhomogeneous boundary conditions, find $E \in H(\mathrm{curl};\Omega, \Sigma^E)$, $H \in H(\mathrm{curl};\Omega, \Sigma^M)$ such that for all admissible test fields, \begin{align} \langle i\omega\epsilon E + i\omega\xi H - \nabla\times H, V\rangle &= -\langle S_E, V\rangle, \tag{5a}\ \langle i\omega\zeta E + i\omega\mu H + \nabla\times E, W\rangle &= -\langle S_M, W\rangle, \tag{5b} \end{align} where $S_E$, $S_M$ collect all source and boundary liftings (Fernando et al., 2022).

Functional analysis in the $H(\mathrm{curl})$ framework is essential for addressing well-posedness, regularity, and compatibility of both interior and boundary traces, as detailed in trace and extension theory for Lipschitz domains (Stachura et al., 2020).

3. Boundary Conditions and Well-Posedness

The anisotropic Maxwell system admits a variety of boundary conditions:

• Perfect electric conductor (PEC): $n\times E = K_M$ on $\Sigma^E$
• Perfect magnetic conductor (PMC): $n\times H = K_E$ on $\Sigma^M$
• Surface charge/current: $n\cdot D = \sigma_E$, $n\cdot B = \sigma_M$

Boundary lifting is implemented by decomposing the total fields as sums of homogeneous and auxiliary trace-carrying components, enabling all test and trial fields to satisfy homogeneous traces and ensuring coercivity and well-posedness in the usual Galerkin framework (Fernando et al., 2022).

On Lipschitz domains, trace operators for $H(\mathrm{curl})$ fields map continuously into fractional Sobolev spaces of tangential traces, with operator norms quantified explicitly in terms of the domain's Lipschitz constant and geometric parameters (Stachura et al., 2020). The analysis guarantees the coercivity of sesquilinear forms associated to the Maxwell system provided the coefficients satisfy uniform ellipticity and loss (for time-harmonic problems).

Quantitative a priori estimates show that the solution norm in $H(\mathrm{curl}; \Omega)$ grows at most quadratically in the domain's Lipschitz constant, and is controlled by data and norms of the boundary incident fields.

4. Finite Element Approximation

For numerical solution, the variational form is discretized on tetrahedral meshes using Nédélec first-kind edge elements of order one, which are curl-conforming and preserve tangential continuity (Fernando et al., 2022). The discrete trial spaces

$$X_h^E \subset H(\mathrm{curl}; \Omega, \Sigma^E),\quad X_h^M \subset H(\mathrm{curl}; \Omega, \Sigma^M)$$

yield the Galerkin problem: find $(E_h, H_h) \in X_h^E \times X_h^M$ such that, for all $(V_h, W_h)$ in the product space,

$$\mathcal{A}((E_h, H_h), (V_h, W_h)) = \mathcal{L}(V_h, W_h).$$

Numerical tests (unit cube, manufactured solutions) exhibit first-order convergence in the $L^2$ norm for $E$ and $H$ as mesh size $h \to 0$, in line with theory for lowest-order edge elements.

5. Generalized Material Responses: Bianisotropy and Source Terms

The inclusion of nonzero polarization ($P$) and magnetic current ($J_M$), as well as bianisotropic magnetoelectric coupling ($\xi, \zeta$), augments only the right-hand side of the main equations. Compatibility (continuity) conditions adjust to include divergence terms of $P, M$ in Gauss's laws, but do not alter the structure of the variational form (Fernando et al., 2022). For given data, the variational system remains coercive under standard assumptions on $\epsilon, \mu$.

A plausible implication is that polarization and magnetic sources can be incorporated in both analysis and computation without changing the core finite element implementation, provided data consistency for compatibility conditions is verified.

6. Analytical and Numerical Robustness

The formulation supports robust analytical and computational treatment of the anisotropic Maxwell system. Coercivity, well-posedness, and convergence are maintained under minimal regularity on coefficients (piecewise regularity), anisotropic tensor structure, and general boundary conditions. The framework is directly extensible to inhomogeneous and complex materials appearing in advanced photonic, metamaterial, and magnetoelectric applications.

Numerical schemes based on this system are validated by quantitative convergence studies and are particularly suited for simulating electromagnetic fields in media with pronounced anisotropy and magnetoelectric coupling (Fernando et al., 2022).

7. Summary Table of Key Structures

Aspect	Mathematical Structure	Reference [arXiv id]
Constitutive Relations	$D=\epsilon E+\xi H+P,\ B=\zeta E+\mu H+M$	(Fernando et al., 2022)
Functional Spaces	$H(\mathrm{curl};\Omega)$ with trace/sobolev subspaces	(Fernando et al., 2022, Stachura et al., 2020)
Weak Form	$$(E,H) \to \mathcal{A}((E, H), (V, W))=\mathcal{L}(V, W)$$	(Fernando et al., 2022)
Finite Element Basis	Nédélec first-kind edge elements	(Fernando et al., 2022)
Compatibility Conditions	$i\omega\rho_E+\nabla\cdot J_E=0$	(Fernando et al., 2022)

The anisotropic linear Maxwell system, as rigorously formulated in the curl-conforming variational setting and discretized by edge elements, constitutes a mathematically robust and computationally viable framework for the analysis and simulation of electromagnetic phenomena in general bianisotropic media (Fernando et al., 2022).