Deterministic PML Envelope

• Deterministic PML envelope is a mathematically rigorous framework that defines explicit error bounds for perfectly matched layer absorbing boundaries in wave propagation and privacy contexts.
• It employs complex coordinate stretching and smooth polynomial or exponential damping profiles to achieve exponential decay of errors, ensuring energy stability and optimal convergence rates.
• Explicit parameter selection formulas enable deterministic calibration for physical wave equations and privacy-preserving mechanisms, reducing reflection and leakage in simulations.

A deterministic PML envelope is an explicit, mathematically rigorous bound describing the error or leakage rate resulting from the use of perfectly matched layer (PML) absorbing boundaries in computational wave propagation and related contexts. The envelope quantifies, for a user-prescribed tolerance, the required thickness and damping strength profile to guarantee the desired suppression of reflections or leakage, and it extends to discrete and continuous settings with provable energy stability and optimal convergence rates. This notion appears both in physical wave equations—such as acoustics, electromagnetics, and elastodynamics—and in information-theoretic settings such as the pointwise maximal leakage (PML) metric for privacy-preserving mechanisms.

1. Mathematical Foundations: Complex Coordinate Stretching

The deterministic PML envelope is constructed via complex coordinate stretching in the Laplace or Fourier–Laplace domain. For linear hyperbolic wave equations, the solution $U(x, t)$ obeys

$$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$

with truncation by PML implemented through complex-valued metrics $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$. The spatial derivative is replaced by $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$. Upon inversion to the time domain, auxiliary fields are introduced to preserve causality and energy structure.

For the first-order acoustic system, the un-split modal PML is derived via Laplace stretching $S_j(x_j, s) = 1 + d_j(x_j)/s$, leading to a system with auxiliary damping variables $\sigma$ and $\psi$ (Duru et al., 2018). Similarly, in Maxwell's equations and elastodynamics, a stretched coordinate yields a modified equation where the absorption profile $\sigma(r)$ or $d_j(x_j)$ is nonzero within a deterministic layer, producing exponential decay of propagating and evanescent modes (Wei et al., 2019, Duru et al., 2022).

2. Deterministic Envelope Profiles: Polynomial and Exponential Damping

The envelope function—$\sigma_{\xi}(\xi)$ or $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$0—is chosen to transition smoothly from zero in the interior to a maximum at the boundary over a layer thickness $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$1. A standard choice is a cubic or higher-order polynomial ramp:

$$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$2

with $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$3 (typically $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$4 for spectral $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$5-refinement) (Duru et al., 2018, Duru et al., 2022). The coefficient $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$6 is tuned via

$$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$7

to ensure a target reflection $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$8 across the layer. Exponential or mixed profiles are also feasible. For envelope wave equations in laser–plasma modeling, Smilei PIC employs $$P^{-1}U_{t} = \sum_{\xi=x,y,z}A_{\xi}\,\partial_{\xi}U,$$9, with empirical choices $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$0 for reflection $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$1 (Bouchard et al., 2024).

3. Energy Estimates and Stability Guarantees

Energy stability is established by deriving continuous Laplace-domain energy estimates and showing that discretizations—e.g., DGSEM or SBP–SAT schemes—recover analogous discrete estimates. For the PML system,

$$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$2

ensures asymptotic stability for all nonnegative, monotone ramp profiles (Duru et al., 2018). Discretely, penalty-based upwind fluxes and stabilizing terms are required in the auxiliary equations. For discrete elements of size $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$3 and degree $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$4, the PML reflection error converges as

$$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$5

and energy-stable time integration is achieved with order-$$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$6 ADER or Taylor schemes (Duru et al., 2018, Duru et al., 2022).

4. Explicit Error Envelopes and Convergence Rates

Key deterministic error envelopes take the following explicit forms:

Problem Type	Envelope Bound	Parameters
2D acoustic wave (disk)	$$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$7	Absorption $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$8, layer thickness $$S_{\xi}(s, \xi) = 1 + \sigma_{\xi}(\xi)/(s + \alpha_{\xi}(\xi))$$9, source harmonics (Bryan et al., 9 Mar 2025)
3D electromagnetics	$S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$0	Layer thickness $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$1, absorbing strength $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$2, time $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$3, permeability $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$4 (Wei et al., 2019)
Polynomial PML Profile	$S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$5	$S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$6: profile order, $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$7: max wave speed, $S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$8: thickness, tol (Duru et al., 2022)
DG spectral element	$S_{\xi}^{-1}(s, \xi) \partial/\partial\xi$9	Element size $S_j(x_j, s) = 1 + d_j(x_j)/s$0, PML width $S_j(x_j, s) = 1 + d_j(x_j)/s$1, polynomial degree $S_j(x_j, s) = 1 + d_j(x_j)/s$2 (Duru et al., 2018)

All deterministic envelopes imply exponential decay in error with increasing damping strength or layer thickness, with explicit formulas guiding parameter selection.

5. Deterministic PML Envelope in Information-Theoretic Contexts

An analogous notion arises in privacy analysis as the deterministic pointwise maximal leakage (PML) envelope. For the Gaussian mechanism $S_j(x_j, s) = 1 + d_j(x_j)/s$3, with $S_j(x_j, s) = 1 + d_j(x_j)/s$4 Gaussian and $S_j(x_j, s) = 1 + d_j(x_j)/s$5, the closed form

$S_j(x_j, s) = 1 + d_j(x_j)/s$6

quantifies the minimal information leakage over all deterministic post-processings for failure probability $S_j(x_j, s) = 1 + d_j(x_j)/s$7 (Saeidian, 13 Jan 2026). Extensions hold for general priors with posterior-variance control or strong log-concavity. In contrast to $S_j(x_j, s) = 1 + d_j(x_j)/s$8-differential privacy (DP), which yields $S_j(x_j, s) = 1 + d_j(x_j)/s$9, the deterministic PML envelope gives $\sigma$0 and is closed under arbitrary post-processing.

6. Practical Guidelines and Implementation Considerations

For physical wave problems:

• Choose smooth polynomial (typically cubic or higher-order) damping profiles; layer thicknesses of $\sigma$1–$\sigma$2 wavelengths are generally sufficient for tolerances $\sigma$3.
• Upwind penalty fluxes and auxiliary equation stabilization are essential for discrete energy stability.
• In complex envelope equations, auxiliary differential equations are required to preserve explicit FDTD solvers (Bouchard et al., 2024).
• Explicit parameter selection formulas allow for deterministic calibration, avoiding trial-and-error tuning.

For privacy mechanisms:

• The PML envelope provides a fully explicit privacy–failure probability trade-off, supporting post-processing robustness, with optimal calibration given by $\sigma$4 for target leakage $\sigma$5.

Numerical benchmarks confirm that observed error decay closely follows the deterministic envelope predictions, with exponential reduction to round-off or discretization floors for sufficient thickness and strength (Bryan et al., 9 Mar 2025, Bouchard et al., 2024).

7. Limitations, Open Problems, and Extensions

Certain scenarios present challenges for deterministic envelope guarantees:

• Anisotropic media or slowness-surface branches violating the geometric stability condition can yield instabilities unless modified PML metrics are used (Duru et al., 2022).
• Interface/guided wave modes (e.g., Scholte, Stoneley waves) lack a comprehensive stability analysis.
• Second-order (displacement-only) PML formulations in general heterogeneous media remain analytically challenging (Duru et al., 2022).
• Implementation on unstructured or curvilinear meshes and adaptive layers requires further investigation to ensure absence of discrete reflections.

A plausible implication is that ongoing research may extend deterministic envelope results to broader classes of equations, mesh structures, and interface phenomena. Practical implementation of PML layers in high-performance codes (e.g., Smilei PIC, WaveQLab3D, ExaHyPE) demonstrates resource usage reduction as high as $\sigma$6 for appropriate envelope configuration in multidimensional problems (Duru et al., 2022, Bouchard et al., 2024).

(Duru et al., 2018, Saeidian, 13 Jan 2026, Wei et al., 2019, Duru et al., 2022, Bouchard et al., 2024, Bryan et al., 9 Mar 2025)