Finite Element Conformal Hessian Complexes

• The paper introduces finite element conformal Hessian complexes as exact discrete realizations of BGG-type differential sequences, using an enlarged kernel (CH = P1 ⊕ Span{|x|²}) to enforce symmetry and tracelessness.
• A novel methodology employs local bubble spaces, commuting projection operators, and precise polynomial exact sequences to handle nonstandard continuity and trace constraints on tetrahedral meshes.
• Applications to the linearized Einstein–Bianchi system and related FEEC frameworks demonstrate that the conformal modification yields stable, structure-preserving discretizations for problems in relativity, elasticity, and fluid mechanics.

Searching arXiv for papers on finite element conformal Hessian complexes and closely related BGG/FEEC constructions. First search: exact phrase and key terms. Finite element conformal Hessian complexes are conforming discrete realizations of a BGG-type differential complex in which an $H^2$ scalar field is mapped by the trace-free Hessian into symmetric, traceless tensor fields and then propagated by higher differential operators through spaces with nonstandard continuity and trace constraints. In three dimensions, one formulation on a contractible polyhedral domain is

${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$

where $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$ is the $6$-dimensional kernel of the dev–Hessian (Huang, 2 Aug 2025). A closely related presentation used for the linearized Einstein–Bianchi system writes the third slot as $H(\div,\Omega;\mathbb S\cap\mathbb T)$ and the terminal operator as $\div$ (Guo et al., 6 Aug 2025). Across these formulations, the central objectives are exactness, commuting interpolation, preservation of symmetry and tracelessness, and the construction of stable finite-element discretizations on tetrahedral or simplicial meshes.

1. Continuous conformal Hessian structure

The continuous conformal Hessian complex is distinguished from the ordinary Hessian complex by the replacement of the full Hessian with its deviatoric part. In the three-dimensional formulation summarized in (Huang, 2 Aug 2025), the leftmost kernel is not merely $P_1(\Omega)$ but

$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$

and the sequence is exact on a contractible polyhedral domain. Exactness means that the kernel of $\dev\circ\hess$ is precisely ${\rm CH}$, the image of ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$0 equals the kernel of ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$1, the image of ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$2 equals the kernel of ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$3, and the final operator is onto (Huang, 2 Aug 2025).

The Einstein–Bianchi-oriented formulation in (Guo et al., 6 Aug 2025) uses the spaces ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$4, ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$5, ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$6, and ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$7, again over a bounded, contractible Lipschitz domain. There the kernel of ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$8 is written as

${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$9

and the exactness relations are stated as $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$0, $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$1, and $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$2 (Guo et al., 6 Aug 2025). A persistent misconception is that the conformal complex is obtained from the ordinary Hessian complex by simply “dropping the trace.” The explicit kernel enlargement from $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$3 to $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$4 or $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$5 shows that the conformal modification changes the algebraic and cohomological structure at the left end of the sequence (Huang, 2 Aug 2025).

This conformal viewpoint also has an $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$6-dimensional BGG formulation. In the form-valued-form construction, the continuous conformal-Hessian complex is the special case $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$7 of a BGG sequence that zig-zags between rows $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$8 and $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\}$$9 in the big BGG diagram, with $6$0 and the later differentials given by exterior differentiation in the first index followed by projection (Hu et al., 5 Mar 2025). This places conformal Hessian complexes alongside other BGG-derived complexes for elasticity, divdiv, and related tensorial PDEs.

2. BGG derivation and polynomial exact sequences

The finite-element literature treats conformal Hessian complexes as part of a broader BGG program in which new complexes are produced from de Rham-type ingredients. In three dimensions, the nonconformal Hessian, elasticity, and divdiv complexes were systematically derived from smooth finite element de Rham complexes, the $6$1-$6$2 decomposition, and trace complexes, together with two reduction operations and one augmentation operation (Chen et al., 2022). The conformal constructions in (Huang, 2 Aug 2025) inherit this logic but specialize it to symmetric traceless tensors and the operator $6$3.

At the algebraic level, (Huang, 2 Aug 2025) proves an exact polynomial conformal Hessian sequence

$6$4

for each integer $6$5. The paper also records the relevant dimension data and kernel/image relations. This polynomial exactness is the local algebraic model for the global discrete complex and is the counterpart, in the conformal setting, of the polynomial exactness results used in ordinary Hessian complexes on tetrahedral, cuboidal, or spline-based meshes (Hu et al., 2023, Arf et al., 2021).

The local BGG derivation in (Huang, 2 Aug 2025) is especially notable because it is carried out on bubble spaces. The exact local diagrams are obtained by combining bubble de Rham and bubble div–div complexes with bubble Hessian and bubble elasticity complexes, together with a reduction map $6$6 (Huang, 2 Aug 2025). This suggests a methodological shift from purely global BGG arguments toward local bubble-complex constructions, a perspective that becomes explicit in the statement that the paper introduces “a novel application of the discrete BGG framework, combined with the geometric decomposition of bubble spaces and a reduction operation, to local bubble finite element complexes” (Huang, 2 Aug 2025).

3. Bubble spaces, trace conditions, and extra smoothness

A defining analytic difficulty of conformal Hessian discretization is the coexistence of algebraic constraints and extra smoothness. The trace-structure framework describes this difficulty directly: for nonstandard finite element complexes such as Hessian, Elasticity, divdiv, “the main challenge is the existence of extra smoothness,” and exact bubble complexes in different dimensions imply correct cohomology and permit $6$7 bounded interpolation (Hu et al., 28 Sep 2025). Although (Hu et al., 28 Sep 2025) presents a two-dimensional Hessian model, its conceptual emphasis on bubble exactness and trace structures is closely aligned with the conformal three-dimensional constructions.

In (Huang, 2 Aug 2025), local spaces are indexed by a smoothness vector $6$8. Scalar bubbles are defined by the vanishing of derivatives on vertices, edges, and faces up to the prescribed orders, and tensorial bubble spaces are then refined by trace conditions appropriate to $6$9 or $H(\div,\Omega;\mathbb S\cap\mathbb T)$0. For instance, $H(\div,\Omega;\mathbb S\cap\mathbb T)$1 imposes vanishing tangential-tangential and related boundary traces, whereas $H(\div,\Omega;\mathbb S\cap\mathbb T)$2 imposes vanishing normal-normal data, vanishing second div–div trace, and edge vanishing (Huang, 2 Aug 2025). A geometric decomposition of the symmetric-traceless subspace of $H(\div,\Omega;\mathbb S\cap\mathbb T)$3 into interior bubbles and facewise contributions is stated explicitly in that work.

This bubble analysis is not merely technical bookkeeping. It is the mechanism by which local exactness is proved and by which global degrees of freedom are organized. The same general principle appears in related constructions. In the form-valued-form setting, local exactness and an Euler-characteristic identity imply global exactness of the discrete conformal-Hessian complex on contractible domains (Hu et al., 5 Mar 2025). In the Einstein–Bianchi discretization, the proof of discrete exactness uses local bubble-complexes on each tetrahedron, a BGG-type splitting, and a cell-by-cell dimension count together with a global Euler identity (Guo et al., 6 Aug 2025).

The middle slot also illustrates how nonstandard trace spaces enter the theory. Earlier work constructed conforming finite elements for $H(\div,\Omega;\mathbb S\cap\mathbb T)$4 and $H(\div,\Omega;\mathbb S\cap\mathbb T)$5, proved conformity and unisolvence, and showed that these elements are not $H(\div,\Omega;\mathbb S\cap\mathbb T)$6-conforming (Sander, 2021). Since conformal Hessian complexes use $H(\div,\Omega;\mathbb S\cap\mathbb T)$7-type tensor fields with symmetry and tracelessness constraints, these results form part of the technical background for designing the second space of the conformal sequence.

4. Global finite-element realizations

The three-dimensional conformal finite-element complex in (Huang, 2 Aug 2025) is built on a shape-regular triangulation $H(\div,\Omega;\mathbb S\cap\mathbb T)$8 using four global spaces: $H(\div,\Omega;\mathbb S\cap\mathbb T)$9 The scalar space uses $\div$0 shape functions together with vertex, edge, face, and cell moments. The tensor spaces use $\div$1 and $\div$2, with degrees of freedom distributed across vertices, edges, faces, and cells so as to enforce the trace/symmetry constraints of $\div$3 and $\div$4 conformity (Huang, 2 Aug 2025).

A more explicit tetrahedral realization is given in (Guo et al., 6 Aug 2025). There one fixes an integer $\div$5 and, on each tetrahedron $\div$6, defines an $\div$7-conforming space $\div$8 with shape functions $\div$9 and inter-element continuity described as $P_1(\Omega)$0 at vertices, $P_1(\Omega)$1 on edges, and $P_1(\Omega)$2 on faces; an $P_1(\Omega)$3-conforming space $P_1(\Omega)$4 with shape functions $P_1(\Omega)$5; an $P_1(\Omega)$6-conforming space $P_1(\Omega)$7 with shape functions $P_1(\Omega)$8; and a discontinuous $P_1(\Omega)$9 space (Guo et al., 6 Aug 2025). The local degrees of freedom are written in terms of vertex derivatives, edge moments, face moments, and interior bubble moments, and unisolvence is established for each slot.

Construction	Discrete sequence	Characteristic feature
(Huang, 2 Aug 2025)	$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$0	smoothness vectors, bubble BGG diagrams, exactness theorem
(Guo et al., 6 Aug 2025)	$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$1	explicit tetrahedral DoFs and Einstein–Bianchi discretization
(Hu et al., 5 Mar 2025)	$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$2 with $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$3 or symmetry-reduced variants	general $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$4-dimensional form-valued-form BGG construction

The global spaces in these constructions are assembled by single-valued geometric degrees of freedom. This is the standard FEEC/BGG mechanism for turning local polynomial exactness into a conforming subcomplex, but in the conformal setting the traces are substantially more intricate because symmetry, tracelessness, and extra smoothness are imposed simultaneously (Huang, 2 Aug 2025, Guo et al., 6 Aug 2025).

5. Commuting projections, exactness, and approximation

A central requirement for a finite element conformal Hessian complex is the existence of canonical interpolation or projection operators that commute with the differential operators. In (Huang, 2 Aug 2025), canonical interpolation operators $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$5, $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$6, and analogous maps for the later slots are defined using the same moments and point-values as the degrees of freedom. In (Guo et al., 6 Aug 2025), the corresponding canonical interpolations $${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$7 are obtained by setting each local degree of freedom equal to the corresponding continuous one, and the commutation identities

$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$8

are verified directly from the degrees of freedom (Guo et al., 6 Aug 2025).

Exactness of the discrete sequence is established under explicit smoothness constraints in (Huang, 2 Aug 2025). If

$${\rm CH}=P_1(\Omega)\oplus\operatorname{Span}\{|\mathbf x|^2\},$$9

then the finite-element conformal Hessian sequence is a complex and is exact (Huang, 2 Aug 2025). The exactness theorem identifies the discrete kernel of $\dev\circ\hess$0 with $\dev\circ\hess$1 and the kernel of $\dev\circ\hess$2 with the image of $\dev\circ\hess$3; the paper also states that each pair of successive differential operators has a stable right-inverse on the discrete spaces (Huang, 2 Aug 2025).

Approximation theory in (Huang, 2 Aug 2025) is stated in local Sobolev norms. For $\dev\circ\hess$4,

$\dev\circ\hess$5

and analogous bounds hold on the tensor slots in the graph norms of $\dev\circ\hess$6, $\dev\circ\hess$7, $\dev\circ\hess$8, and $\dev\circ\hess$9 (Huang, 2 Aug 2025). The paper therefore states that one obtains optimal-order convergence of Galerkin methods that respect the differential complex. In the form-valued-form construction, exactness can also be derived from local exactness plus the Euler-characteristic identity

${\rm CH}$0

and the interpolation operators satisfy ${\rm CH}$1 (Hu et al., 5 Mar 2025).

This emphasis on commuting projectors links the conformal theory to earlier FEEC discretizations of ordinary Hessian complexes. Spline-based structure-preserving discretization on affine images of the unit cube proved closed exactness, uniform inf–sup stability for the associated Hodge–Laplacians, and optimal ${\rm CH}$2-type convergence in graph norms (Arf et al., 2021). The conformal theory extends that structure-preserving agenda to symmetric and traceless tensor variables.

6. Applications and broader research landscape

The most explicit application presently attached to a conformal Hessian finite-element complex is the discretization of the linearized Einstein–Bianchi system. In (Guo et al., 6 Aug 2025), the system is reformulated as the Hodge wave equation associated with the conformal Hessian complex, with unknowns ${\rm CH}$3, ${\rm CH}$4, and ${\rm CH}$5. The semidiscrete scheme uses the discrete spaces ${\rm CH}$6, ${\rm CH}$7, and ${\rm CH}$8, and standard energy arguments are stated to yield stability and an ${\rm CH}$9 error estimate in the ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$00-energy norm under the regularity assumptions ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$01, ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$02, and ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$03 (Guo et al., 6 Aug 2025).

More broadly, (Huang, 2 Aug 2025) states that the resulting conformal complexes support stable and structure-preserving numerical methods for applications in relativity, Cosserat elasticity, and fluid mechanics. The form-valued-form framework extends the construction to ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$04 for general ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$05, ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$06, and ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$07, and lists applications to strain and stress tensors in continuum mechanics and metric and curvature tensors in differential geometry in any dimension (Hu et al., 5 Mar 2025). The same work also states that the entire finite-element complex is invariant under conformal scalings of the ambient metric (Hu et al., 5 Mar 2025).

The conformal constructions sit within a wider program of Hessian-complex discretization. Nonconformal but structurally adjacent developments include conforming ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$08 and ${\rm CH}\xrightarrow{\subset}H^2(\Omega)\xrightarrow{\dev\circ\hess}H(\symcurl;\mathbb S\cap\mathbb T)\xrightarrow{\symcurl}H(\div\div;\mathbb S\cap\mathbb T)\xrightarrow{\div\circ\div}L^2(\Omega)\to0,$09 elements (Sander, 2021), three-dimensional finite element Hessian complexes derived from BGG machinery (Chen et al., 2022), exact grad–grad complexes on cuboid meshes (Hu et al., 2023), exact Hessian complexes on Worsey–Farin splits (Gong et al., 2023), and distributional Hessian and divdiv complexes on triangulations whose cohomology is isomorphic to the continuous versions and hence to de Rham cohomology with coefficients (Hu et al., 2023). Taken together, these works indicate that finite element conformal Hessian complexes are part of a larger FEEC/BGG synthesis in which exactness, cohomology, and operator-commuting discretization are treated as primary design principles rather than secondary consequences.