Unified Surface Energy Matrix

• Unified Surface Energy Matrix is a dimension-agnostic construct that encodes anisotropic surface energy and facilitates structure-preserving numerical methods.
• It integrates the Cahn–Hoffman vector with a stabilization term to guarantee unconditional energy stability and volume conservation across 2D and 3D geometries.
• The framework supports both strong and weak variational formulations in geometric PDEs, enabling robust discretizations for anisotropic diffusion and curvature flows.

The unified surface energy matrix is a mathematical construct central to the variational and computational modeling of anisotropic surface phenomena, particularly in geometric partial differential equations (PDEs) and structure-preserving numerical methods. This object encodes the anisotropic surface energy density and its geometric consequences in a single, dimension-agnostic matrix form, delivering a unified framework for both theoretical analysis and robust discretizations in two and three dimensions. It couples the orientation-dependent surface energy, the Cahn–Hoffman vector, and a stabilizing term to guarantee unconditional energy stability and volume conservation, regardless of the specific anisotropy or geometry in question (Bao et al., 2023).

1. Construction: Anisotropic Surface Energy, Cahn–Hoffman Vector, and Stabilization

Let $\Gamma \subset \mathbb{R}^d$ ($d=2$ or $3$) denote a smooth closed orientable hypersurface, with outward unit normal $\mathbf{n} \in \mathbb{S}^{d-1}$. The anisotropic surface energy density is a strictly positive $C^2$ function

$\gamma : \mathbb{S}^{d-1} \to \mathbb{R}_{>0},$

extended one-homogeneously to $\mathbb{R}^d \setminus \{0\}$ by

$$\gamma(\mathbf{p}) = |\mathbf{p}|\, \gamma(\mathbf{p}/|\mathbf{p}|), \qquad \gamma(0) = 0.$$

The Cahn–Hoffman vector

$$\boldsymbol{\xi}(\mathbf{n}) = \nabla_{\mathbf{p}} \gamma(\mathbf{p})\Big|_{\mathbf{p} = \mathbf{n}}$$

satisfies the crucial identity $$\gamma(\mathbf{n}) = \boldsymbol{\xi}(\mathbf{n}) \cdot \mathbf{n}$$ due to homogeneity.

Unconditional energy stability necessitates the inclusion of a stabilizing function $k : \mathbb{S}^{d-1} \to [0, \infty)$, chosen minimally (via local matrix positivity) or larger, depending on the regularity of the anisotropy. The minimal stabilizer $k_0(\mathbf{n})$ is defined so that a local energy estimate holds on each discrete element (Bao et al., 2023).

2. Definition and Algebraic Structure of the Unified Surface Energy Matrix

With these ingredients, the unified surface energy matrix is defined as

$$\mathbf{G}_k(\mathbf{n}) = \gamma(\mathbf{n})\, I_d - \mathbf{n}\, \boldsymbol{\xi}^\top(\mathbf{n}) + \boldsymbol{\xi}(\mathbf{n})\,\mathbf{n}^\top + k(\mathbf{n})\, \mathbf{n} \mathbf{n}^\top.$$

This decomposition admits a symmetric part

$$\mathbf{G}_k^{(s)}(\mathbf{n}) = \gamma(\mathbf{n})\, I_d + k(\mathbf{n})\, \mathbf{n}\mathbf{n}^\top,$$

and a skew (anti-symmetric) part

$$\mathbf{G}_k^{(a)}(\mathbf{n}) = -\mathbf{n}\, \boldsymbol{\xi}^\top + \boldsymbol{\xi}\,\mathbf{n}^\top.$$

On the tangent plane to $\Gamma$, $\mathbf{G}_k^{(s)}$ reduces to a multiple of $\gamma(\mathbf{n}) I$, as $\mathbf{n}$ is normal to the tangent bundle. The skew part is identically divergence-free in all energy balances, and thus does not contribute to energy dissipation.

3. Variational and Strong Formulations in Geometric PDEs

The matrix $\mathbf{G}_k(\mathbf{n})$ governs both the strong and weak forms of anisotropic surface diffusion and geometric evolution:

• Strong form: The chemical potential is $\mu = \nabla_\Gamma \cdot \boldsymbol{\xi}$. The governing equation for the embedding $X: \Gamma \to \mathbb{R}^d$ is

$$\nabla_\Gamma \cdot [\mathbf{G}_k(\mathbf{n})\, \nabla_\Gamma X] = -\mu\,\mathbf{n}.$$

• Weak (variational) form: For all test functions $\omega \in [H^1(\Gamma)]^d$,

$$\int_\Gamma \mu\,\mathbf{n} \cdot \omega\,dA = \int_\Gamma [\mathbf{G}_k(\mathbf{n})\,\nabla_\Gamma X] : \nabla_\Gamma \omega\,dA,$$

and, paralleling the diffusion equation,

$$\partial_t X \cdot \mathbf{n},\, \varphi + \nabla_\Gamma \mu \cdot \nabla_\Gamma \varphi = 0,\qquad \forall \varphi \in H^1(\Gamma).$$

These forms apply identically for $d=2$ (curves) and $d=3$ (surfaces), achieving a true unification across dimensions (Bao et al., 2023).

4. Structure-Preservation: Volume Conservation and Unconditional Energy Stability

The algebraic construction yields critical structure-preserving features:

• The skew component $\mathbf{G}_k^{(a)}$ exactly cancels in discrete volume balances, ensuring volume conservation.
• The symmetric component, together with the choice $k(\mathbf{n}) \ge k_0(\mathbf{n})$, guarantees that the bilinear form $v \mapsto v^T \mathbf{G}_k(\mathbf{n}) v$ is positive semidefinite on the tangent bundle, essential for the discrete “local energy estimate.”
• Summing over all elements in a PFEM discretization,

$$W^{m+1} - W^m \le -\tau \|\nabla_\Gamma \mu^{m+1}\|^2 \le 0,$$

establishes unconditional energy dissipation without CFL-type time-step restrictions.

• The same matrix handles both strongly and weakly anisotropic surface energies, provided the minimal stability condition

$\gamma(-\mathbf{n}) < (5-d)\, \gamma(\mathbf{n})$

is met—i.e., $\gamma(-\mathbf{n}) < 3\, \gamma(\mathbf{n})$ for $d=2$, $\gamma(-\mathbf{n}) < 2\,\gamma(\mathbf{n})$ for $d=3$.

5. Generalization: Symmetrization, Parametric Families, and Connections to Other Energies

Recent work (Bao et al., 31 Dec 2025) extends the concept by introducing a symmetrization parameter $\alpha$: $$\hat{\mathbf{G}}_k^\alpha(\theta) = \hat\gamma(\theta)\, I_2 - \mathbf{n}(\theta)\,\boldsymbol{\xi}^\top(\theta) + \alpha\, \boldsymbol{\xi}(\theta)\,\mathbf{n}^\top(\theta) + k(\theta)\,\mathbf{n}(\theta)\,\mathbf{n}^\top(\theta).$$ This two-parameter family encompasses all known forms of surface energy matrices. The optimal symmetrization $\alpha = -1$ minimizes the stabilizer $k$ and relaxes the energy stability condition to $3\hat\gamma(\theta) \ge \hat\gamma(\theta-\pi)$, while $\alpha \neq -1$ enforces a strictly stronger condition or additional constraints.

The matrix framework recovers previous structure-preserving finite element discretizations for curvature flow, surface diffusion, and area-preserved flow under both classical and anisotropic settings (Bao et al., 31 Dec 2025, Bao et al., 2022, Bao et al., 2022).

6. Practical Implementation in Structure-Preserving PFEMs

The unified surface energy matrix directly informs both the strong PDE formulation and the parametric finite element discretization (SP-PFEM):

• In spatially and temporally discretized schemes, the matrix $G_k(\mathbf{n})$ couples with mass-lumped inner products and discrete derivatives to achieve energy stability and mesh quality.
• Exact discrete area or volume conservation is enforced via midpoint-type normal approximations.
• Local algebraic estimates, crucial for unconditional stability, rely on the positive definiteness of $G_k^{(s)}$ and the controlled addition of $k(\mathbf{n})$ chosen per element.

Extensive numerical benchmarks validate the theoretical properties—area decay or conservation, monotone energy dissipation, and robust mesh evolution—across a wide range of anisotropic energies and flow regimes (Bao et al., 2023, Bao et al., 31 Dec 2025, Bao et al., 2022).

7. Broader Applications and Unified Operator Formulations

Analogous unified matrix constructions characterize surface free energy and stress in thermodynamics and statistical mechanics (Pasquale et al., 2019), geometric mesh generation and adaptation (Kolasinski et al., 2019), continuum-lattice energy bridging (Rosakis, 2012), and operator-theoretic approaches to Casimir forces (Bimonte et al., 2021). In each context, a single matrix or tensor collects orientation-dependent thermodynamic, atomistic, or operator contributions, providing a concise quadratic or block-matrix form that renders the energy computation, stability analysis, or response prediction both unified and dimension-agnostic. These matrices underpin Shuttleworth's equation in interfacial stress analysis, the alignment denominator in moving mesh PDEs, and the fluctuation-dissipation trace–log formula in quantum field theory.

In summary, the unified surface energy matrix is fundamental for encoding anisotropic energetics and their induced variational structures, enabling robust, structure-preserving schemes and universal operator-theoretic formulations in mathematical physics, computational geometry, and numerical analysis.