Symmetric Deep Neural Networks

• Symmetric deep neural networks are architectures that enforce permutation invariance, enabling effective high-dimensional function approximation.
• They utilize symmetric Korobov spaces and squared-ReLU subnets to achieve dimension-free approximation rates and mitigate the curse of dimensionality.
• The design improves computational efficiency and generalization in fields like physics and finance by integrating sparse grid symmetrization and Vandermonde-inverse aggregation.

Symmetric deep neural networks are architectures designed to exploit permutation symmetry inherent in function classes encountered in scientific and mathematical modeling, particularly for high-dimensional tasks. These models enforce invariance under permutations of input coordinates, leading to substantial computational advantages and rigorous improvements in both approximation and generalization for functions possessing such symmetry. The paradigm offers dimension-free rates, avoiding the curse of dimensionality previously endemic to neural approximations of symmetric functions, as established by the dimension-free approximation and learning guarantees for symmetric Korobov spaces (Lu et al., 16 Nov 2025).

1. Symmetric Korobov Spaces and Function Classes

Symmetric Korobov spaces are a central construct for analyzing permutation-symmetric functions in multiple dimensions. Let $r \geq 1$ and $d \geq 1$; the periodic Korobov space $H^r_{\mathrm{Kor}}(d)$ is defined on $[0,1]^d$ as the set of periodic functions $f$ admittting a Fourier expansion $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$, equipped with norm

$$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$

Equivalently, the zero-boundary "hat-basis" formulation establishes

$$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$

with semi-norm $$|f|_{2,2} := \bigl\| \partial_{x_1}^2 \cdots \partial_{x_d}^2 f \bigr\|_{L^2(\Omega)}$$. Functions $f$ are called symmetric if $d \geq 1$0 for all $d \geq 1$1; accordingly, the symmetric subspace $d \geq 1$2 is defined by restriction to symmetric functions, and in Fourier coordinates requires $d \geq 1$3 for all coordinate permutations.

2. Dimension-Free Approximation by Deep Symmetric Networks

The main theorem establishes that any $d \geq 1$4 can be approximated, for integer $d \geq 1$5, by a symmetric squared-ReLU network $d \geq 1$6 of the form

$d \geq 1$7

where each $d \geq 1$8 is itself a squared-ReLU network of width $d \geq 1$9 and depth $H^r_{\mathrm{Kor}}(d)$0. The total number $H^r_{\mathrm{Kor}}(d)$1 of summands satisfies

$H^r_{\mathrm{Kor}}(d)$2

The energy-norm error satisfies

$H^r_{\mathrm{Kor}}(d)$3

where $H^r_{\mathrm{Kor}}(d)$4 depends polynomially on $H^r_{\mathrm{Kor}}(d)$5 but not exponentially. The approximation rate $H^r_{\mathrm{Kor}}(d)$6 is thus dimension-free; to drive the energy-norm below $H^r_{\mathrm{Kor}}(d)$7 requires $H^r_{\mathrm{Kor}}(d)$8, achievable with network depth $H^r_{\mathrm{Kor}}(d)$9, width $[0,1]^d$0, and weights bounded by $[0,1]^d$1.

3. Permutation-Invariant Network Architecture

Symmetry is imposed by grouping tensor-product sparse grid basis functions $[0,1]^d$2 into symmetrized blocks

$[0,1]^d$3

While direct summation over $[0,1]^d$4 permutations is intractable, Lemma 4.1 represents $[0,1]^d$5 as a linear combination of only $[0,1]^d$6 exponentials of inner-product features $[0,1]^d$7, $[0,1]^d$8, with $[0,1]^d$9. Recovery of the symmetrized output is performed through a Vandermonde-inverse linear layer.

Each $f$0 is approximated by feeding the scalar hat function $f$1 into a product-of-exponentials, using shallow squared-ReLU subnets and an $f$2-deep binary-tree of ReLU-based bilinear blocks for the $f$3-fold product, requiring $f$4 neurons. A final linear layer of width $f$5 combines these channels, with global weight-sharing across combinatorial block types to ensure permutation invariance.

4. Mathematical Framework Underpinning Dimension-Free Rates

Key ingredients for dimension-free results include:

• The energy-based sparse grid index set $f$6, which replaces total-degree sets to reduce the dominant $f$7 term in the error estimate. This produces cardinality $f$8.
• Exploiting permutation symmetry by aligning with ordered multi-indices ($f$9) and symmetrizing bases, resulting in $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$0 distinct symmetric blocks (exponential in $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$1 only).
• Realizing each symmetrized block $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$2 via a squared-ReLU subnet of width $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$3, depth $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$4, and $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$5 parameters, attaining $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$6-accuracy $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$7.
• By truncating to $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$8 blocks and approximating each to error $$f(x) = \sum_{k \in \mathbb{Z}^d} \hat{f}_k e^{2\pi i k \cdot x}$$9, total $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$0 error is $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$1, yielding an algebraic, truly dimension-free rate.

The relevance lies in reducing the exponential cost normally expected in $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$2 for generic approximators, making the approach scalable to high-dimensional symmetric problems.

5. Sample Complexity and Generalization Guarantees

For supervised learning of symmetric Korobov functions, let the target $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$3 satisfy $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$4. Observed i.i.d. samples $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$5 are distributed so that $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$6 and $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$7 almost surely, and the hypothesis class $$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$8 is the set of symmetric networks described above.

The empirical risk minimizer

$$\|f\|_{H^r_{\mathrm{Kor}}}^2 = \sum_{k \in \mathbb{Z}^d} |\hat{f}_k|^2 \prod_{j=1}^d (1 + |k_j|^2)^r.$$9

admits the bound

$$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$0

where $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$1 is polynomial in $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$2, $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$3, $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$4. By choosing $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$5 and $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$6, one achieves $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$7. High-probability bounds are also established: for any $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$8, with probability at least $$X^{2,2}(\Omega) := \{ f \in L^2(\Omega) : f|_{\partial\Omega} = 0, D^\alpha f \in L^2(\Omega) \;\forall\, |\alpha|_\infty \leq 2 \},$$9,

$$|f|_{2,2} := \bigl\| \partial_{x_1}^2 \cdots \partial_{x_d}^2 f \bigr\|_{L^2(\Omega)}$$0

A plausible implication is that learning symmetric function classes with deep networks can achieve sample and approximation efficiency competitive with classical statistical rates, with dimension-independent leading factors.

6. Implications and Significance for High-Dimensional Learning

The dimension-free results obtained for symmetric deep neural networks represent a substantial advance over previous approximation and generalization bounds, as both the convergence rates and constant prefactors scale at most polynomially with ambient dimension, as opposed to classical exponential dependencies. This suggests a scalable pathway for approximating physically or mathematically symmetric models, such as those in computational physics, finance, and chemistry.

The architectural insights—enforcing permutation invariance via sparse grid symmetrization and Vandermonde-based aggregation—may generalize to other domains requiring strict invariance under variable permutation, such as set-based models or particle-interaction networks. Broadly, the approach expands the class of feasible problems for neural approximation and learning in high-dimensional symmetric settings, and demonstrates that by carefully matching neural architecture to underlying function symmetry, one can eliminate a principal bottleneck traditionally faced by generic deep learning models (Lu et al., 16 Nov 2025).