✓✓ Beats tuned baseline 2026

Reciprocal-Lattice Gauge-Covariant Bloch Network

Implementation & benchmark of arXiv:2608.14348 — Full-field and Bloch-periodic-factor discretizations: Accuracy and phantom modes

Usefulness7/10

Difficulty4/10

Novelty8/10

Architecture Accuracy Stability Generalization Embedding Mlp Regularization Pde Harmonic analysis Symmetry

Source paper: Full-field and Bloch-periodic-factor discretizations: Accuracy and phantom modes arXiv:2608.14348 ⓘ · analyzed Sep 1, 2026

AI-generated research hypothesis, automatically tested. Not peer-reviewed.

Idea description

Build a Bloch-conditioned neural model whose periodic-factor representation transforms covariantly when the supplied Bloch wavenumber is shifted by a reciprocal lattice vector. Either canonicalize q to the first Brillouin zone or augment training with mathematically paired examples whose outputs differ by the exact phase gauge. This prevents the network from learning inconsistent predictions for physically identical Bloch modes.

Formulas

$$E(x)=\widetilde{E}_q(x)e^{iqx},\qquad \widetilde{E}_q(x+a)=\widetilde{E}_q(x),$$

$$L(q)\widetilde{E}=\widetilde{E}''+2iq\widetilde{E}'+(k^2-q^2)\widetilde{E}=0,$$

$$\widetilde{E}_{q+m\kappa}(x)=\widetilde{E}_q(x)e^{-im\kappa x},\qquad \kappa=\frac{2\pi}{a},\quad m\in\mathbb{Z},$$

$$\widehat{F}(q+m\kappa)[x]=e^{-im\kappa x}\widehat{F}(q)[x]$$

Mathematical statement

For a one-dimensional cell of length a, the physical field satisfies E''(x)+k^2E(x)=0 with Bloch boundary conditions E(a)=zE(0), E'(a)=zE'(0), where z=exp(iqa). Writing E(x)=\widetilde{E}_q(x)exp(iqx) gives the periodic-factor operator L(q)\widetilde{E}=\widetilde{E}''+2iq\widetilde{E}'+(k^2-q^2)\widetilde{E}=0 and \widetilde{E}_q(x+a)=\widetilde{E}_q(x). Let \kappa=2\pi/a. Since exp(i(q+m\kappa)x)=exp(iqx)exp(im\kappa x), the same physical field has the gauge transformation \widetilde{E}_{q+m\kappa}(x)=\widetilde{E}_q(x)exp(-im\kappa x). The homogeneous-cell aliases are \widetilde{E}_m(x)=exp(im\kappa x), q_{m,\pm}=\pm k-m\kappa. The implementation exploits this exact reciprocal-lattice covariance rather than treating q as an ordinary scalar feature.

Implementation notes

Integrate this at the input/output interface of a Bloch-conditioned MLP or Fourier neural operator that maps a periodic material profile and Bloch wavenumber q to a complex periodic field, band frequency, or mode coefficient tensor. Represent complex fields as two real channels sampled at cell locations x_j. First implement the canonicalized baseline: compute m=round(q/\kappa), q_0=q-m\kappa, feed q_0 to the network, predict u_0(x_j), and return u(x_j)=u_0(x_j)exp(iq_0x_j)exp(im\kappa x_j). If the target is a periodic factor, transform it as exp(-im\kappa x_j)u_0(x_j) according to the gauge convention. For the stronger version, retain raw q and add a covariance loss. Use m in {-2,-1,1,2}, generate q_shift=q+m\kappa, and transform the reference prediction with y_shift=exp(-im\kappa x)y(q). Optimize task_loss(y(q),target(q))+lambda_cov mean_j|y(q_shift)-exp(-im\kappa x_j)y(q)|^2. The quantities \kappa and the phase factor are exact; prediction errors and lambda_cov are empirical. Start with a 1D periodic coefficient dataset and a small Fourier neural operator, comparing raw-q conditioning, Brillouin-zone clipping, and gauge covariance. Measure error at unseen reciprocal shifts and consistency of reconstructed physical fields. Success means near-zero post-gauge discrepancy and lower error on aliased queries without increasing model size.

Verification

Beats tuned baseline

Stage 1 · Toy mechanism gate: Passed ✓

Stage 2 · Mechanism transferred to benchmark: Confirmed ✓

Stage 2 · Practical benchmark result: Significant win ✓

Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works

Stage 1 — Mechanism check agent confidence 9/10

Built a reciprocal-lattice gauge-aware Bloch MLP wrapper and numerical verification suite. For shifts m=-3..3, phase reconstruction error was 3.43e-16 and operator covariance residual was 6.60e-15, matching machine-precision predictions; the naive untransformed alias mismatch was exactly 1.4142 as predicted. In the toy MLP, mean nonzero-shift RMSE fell from 0.8810 for raw-q conditioning to 0.02344 with canonicalization plus the exact gauge phase, demonstrating the promised mechanism.

Agent confidence: 9/10

Baseline: Raw-q MLP: nonzero reciprocal-shift mean RMSE 0.88098; m=0 RMSE 0.02692.
Idea: Canonicalized q with exact gauge transformation: nonzero reciprocal-shift mean RMSE 0.02344; identical 0.02344 RMSE across m=-3..3.
Math check: Maximum phase identity error 3.43e-16; maximum operator covariance relative error 6.60e-15; predicted naive alias RMS sqrt(2)=1.41421, observed 1.41421.

Limitations: This is a small synthetic 1D experiment with a hand-designed periodic target and an MLP, not a full material-profile dataset or Fourier neural operator. It tests canonicalization and exact gauge wrapping, but not the stronger covariance-loss training variant, physical band-structure accuracy, FLOPs, or large-scale generalization.

How to run: python3 experiment.py

Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.14348) Sep 2, 2026

Beats the tuned baseline by 94.8% (p=0.0081, wins 8 / 8 paired seeds; metric: lower is better).

Benchmark: Bloch_gauge_pde
Model: mlp_tiny
Paired seeds: 8
Baseline mean: 0.3918
Idea mean: 0.0205
Effect (Δ): -0.3713 (−94.8%; negative = idea better)
Wins: 8 / 8 paired seeds
p-value: 0.0081 (permutation test, 20 000 shuffles)
Smallest detectable effect: ±3.2%
Mechanism: Confirmed ✓
Practical effect: Helps
Baseline tuning: swept over 3 configs

Limitations:

The custom task is a fast synthetic 1D Bloch-factor surrogate rather than a discretized heterogeneous Maxwell/Helmholtz solver or full Fourier neural operator. Covariance-loss augmentation, physical band-structure quantities, larger models, speed/FLOP tradeoffs, and real material-profile data were not tested.

How to run: python3 stage2_bench.py

Verdict computed by deterministic test code from paired-seed statistics — not by the language model.

Stage 2 — Benchmark Sep 1, 2026 latest Worked ✓

Agent confidence: 9/10

Built and ran a structurally matched custom 1D Bloch/PDE benchmark because no built-in track contained reciprocal-lattice gauge structure. With identical mlp_tiny systems and an equal three-point learning-rate sweep, canonicalized q plus the exact gauge phase achieved significantly lower test MSE than raw-q conditioning. The trained-model mechanism signature also confirmed the predicted reciprocal-shift covariance quantitatively.

Baseline: Raw-q mlp_tiny: tuned lr=0.006, 25 epochs; full 8-seed test MSE mean 0.3917703, std 0.0135264.
Idea: Canonicalized q plus exact gauge phase: lr=0.006, 25 epochs; full 8-seed test MSE mean 0.0205089, std 0.0010668; paired delta -0.3712614, 8/8 wins, permutation p=0.0081.

Benchmark result