arxiv: 2604.15116 · v1 · submitted 2026-04-16 · 🧮 math.NA · cs.NA· physics.comp-ph

Recognition: unknown

Asymptotic gauge-invariant Hybrid High-Order method for magnetic Schr\"odinger equations

Joubine Aghili

Authors on Pith no claims yet

Pith reviewed 2026-05-10 10:18 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.comp-ph

keywords Hybrid High-Order methodmagnetic Schrödinger equationgauge invariancepolyhedral meshescovariant gradientGårding inequalityAharonov-Bohm effectconvergence analysis

0 comments

The pith

A Hybrid High-Order method for magnetic Schrödinger equations constructs a discrete covariant gradient to achieve asymptotic gauge covariance on polyhedral meshes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical scheme for the Schrödinger equation coupled to a magnetic vector potential. Physical observables must remain unchanged under gauge transformations, and the work ensures this property holds at the discrete level without unphysical artifacts. It does so by defining a covariant discrete gradient operator that works on general polyhedral meshes. The resulting bilinear form is shown to be asymptotically gauge-invariant, to deliver optimal convergence rates, and to satisfy a discrete Gårding inequality that stabilizes the ground-state computation.

Core claim

By constructing a discrete covariant gradient operator on arbitrary polyhedral meshes, the Hybrid High-Order method produces a bilinear form that is asymptotically gauge-invariant. The scheme attains optimal convergence rates in appropriate norms and preserves a discrete Gårding inequality, which guarantees a stable discrete ground state. These properties are confirmed by experiments that recover the Fock-Darwin energy levels and reproduce the Aharonov-Bohm effect.

What carries the argument

The discrete covariant gradient operator, which incorporates the magnetic vector potential into the Hybrid High-Order reconstruction to enforce asymptotic gauge covariance in the discrete bilinear form.

If this is right

Optimal convergence rates hold for the approximation of solutions and energies.
The discrete Gårding inequality ensures stability of the computed ground state.
The scheme reproduces gauge-dependent physical phenomena such as the Aharonov-Bohm effect without artifacts.
The method applies directly to general polyhedral meshes without requiring additional regularity assumptions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The asymptotic gauge covariance may allow reliable long-time simulations of time-dependent magnetic Schrödinger problems on complex geometries.
The construction could be adapted to other gauge-invariant quantum models, such as those with electromagnetic coupling in higher dimensions.
Mesh-independent error bounds for gauge-related quantities would follow if the discrete operator satisfies uniform stability estimates beyond the paper's analysis.

Load-bearing premise

A discrete covariant gradient operator can be built on arbitrary polyhedral meshes that preserves the asymptotic gauge covariance property without extra mesh regularity or restrictions.

What would settle it

Numerical experiments on irregular polyhedral meshes that exhibit either loss of asymptotic gauge covariance or failure to recover optimal convergence rates would disprove the central claims.

read the original abstract

We introduce a Hybrid High-Order (HHO) method for the Schr\"odinger equation in the presence of a magnetic vector potential. In quantum mechanics, physical observables are invariant under continuous gauge transformations, which must be kept at the discrete level to avoid unphysical artifacts. To address this, we construct a discrete covariant gradient operator on arbitrary polyhedral meshes. We prove that the resulting discrete bilinear form guarantees gauge covariance asymptotically at the discrete level. The resulting scheme achieves optimal convergence rates and preserves a discrete Garding inequality, guaranteeing a stable ground state. The theoretical properties of the scheme are corroborated by numerical experiments, including the computation of the Fock-Darwin fundamental energy and replicating the Aharonov-Bohm effect.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a new discrete covariant gradient for HHO on polyhedral meshes that delivers asymptotic gauge invariance for the magnetic Schrödinger equation, plus optimal rates and stability.

read the letter

The paper introduces an asymptotically gauge-invariant Hybrid High-Order method for the magnetic Schrödinger equation. The central new element is a discrete covariant gradient operator defined on arbitrary polyhedral meshes that makes the discrete bilinear form gauge covariant in the limit. They do a solid job laying out the construction and then proving the key properties. The asymptotic gauge covariance is established, along with optimal convergence rates in the appropriate norms and a discrete Gårding inequality that supports stability of the computed ground state. The numerical experiments replicate known physical behaviors like the Fock-Darwin spectrum and the Aharonov-Bohm effect, which is good to see. This work fits into the line of research on structure-preserving discretizations for quantum problems. The use of HHO allows for high-order accuracy and flexibility with mesh types, which is a strength compared to more standard finite element approaches that might require special handling for the magnetic term. The main limitation is that the invariance is asymptotic rather than exact. This is probably unavoidable in the HHO setting without additional constraints, but it does mean that for very coarse meshes or certain gauges there could be small unphysical effects. The assumption that the operator can be built without extra regularity on the potential or mesh seems reasonable, and the paper discusses practical aspects adequately. The error analysis and numerical results hold up without inconsistencies. This paper is useful for researchers in numerical analysis who work on discretizations for Schrödinger-type equations with electromagnetic fields. Someone looking for a method that balances accuracy, stability, and physical fidelity on general meshes would find it relevant. It has enough substance in the construction and analysis to warrant sending it out for peer review rather than a desk reject.

Referee Report

0 major / 3 minor

Summary. The manuscript introduces a Hybrid High-Order (HHO) method for the magnetic Schrödinger equation. It constructs a discrete covariant gradient operator on arbitrary polyhedral meshes to enforce asymptotic gauge covariance at the discrete level. The authors prove that the resulting bilinear form achieves this property, delivers optimal convergence rates, and satisfies a discrete Gårding inequality that guarantees stability of the ground state. These claims are supported by numerical experiments on the Fock-Darwin fundamental energy and the Aharonov-Bohm effect.

Significance. If the central proofs hold, the work provides a useful contribution to structure-preserving discretizations for quantum problems with magnetic potentials. Gauge invariance is a fundamental physical requirement, and extending HHO methods to preserve it asymptotically on general meshes addresses a practical need. The combination of the discrete Gårding inequality, optimal rates, and targeted numerical tests on physically relevant examples strengthens the case for the method's utility in computational quantum mechanics.

minor comments (3)

[Abstract] Abstract: the phrase 'guarantees gauge covariance asymptotically at the discrete level' is central but left imprecise; the introduction or §2 should state explicitly whether this refers to the mesh-size limit h→0, a specific scaling with the gauge parameter, or another regime.
[Method section] The construction of the discrete covariant gradient (presumably in §3) is stated to work on arbitrary polyhedral meshes without extra restrictions; a brief remark confirming that no hidden quasi-uniformity or shape-regularity assumption enters the gauge-covariance proof would strengthen the claim.
[Numerical experiments] Numerical section: the experiments on the Fock-Darwin energy and Aharonov-Bohm effect are described only at a high level; adding a short table or plot caption that lists the observed convergence orders for different polynomial degrees would make the 'optimal rates' claim immediately verifiable.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript on the asymptotically gauge-invariant Hybrid High-Order method for magnetic Schrödinger equations and for recommending minor revision. The referee's summary accurately reflects the construction of the discrete covariant gradient on polyhedral meshes, the proof of asymptotic gauge covariance, optimal convergence rates, the discrete Gårding inequality, and the numerical tests on the Fock-Darwin energy and Aharonov-Bohm effect. As no specific major comments were provided, we will incorporate minor revisions to enhance clarity, notation, and presentation where appropriate.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper constructs a discrete covariant gradient operator within the Hybrid High-Order framework for the magnetic Schrödinger equation on polyhedral meshes, then proves asymptotic gauge covariance of the bilinear form, optimal convergence rates, and preservation of a discrete Gårding inequality. These steps derive directly from the operator definition and standard finite-element analysis techniques without reducing any claimed result to a fitted parameter, self-referential definition, or load-bearing self-citation chain. The derivation remains self-contained against the continuous problem and external benchmarks for discretization stability.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence and properties of a newly constructed discrete covariant gradient operator on polyhedral meshes; no free parameters or invented physical entities are mentioned.

axioms (1)

domain assumption A discrete covariant gradient operator exists on arbitrary polyhedral meshes that enables asymptotic gauge covariance for the bilinear form.
Invoked to define the discrete scheme and prove its properties.

pith-pipeline@v0.9.0 · 5416 in / 1115 out tokens · 29669 ms · 2026-05-10T10:18:10.802327+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 2 canonical work pages

[1]

Aharonov and D

[1]Y. Aharonov and D. Bohm,Significance of electromagnetic potentials in the quantum theory, Physical Review, 115 (1959), p. 485–491, https://doi.org/10.1103/physrev.115.485, http: //dx.doi.org/10.1103/PhysRev.115.485. [2]F. Alouges and V. Bonnaillie-No ¨el,Numerical computations of fundamental eigenstates for the schr¨ odinger operator under constant mag...

work page doi:10.1103/physrev.115.485 1959
[2]

[14]I. W. Sudiarta and D. J. W. Geldart,Solving the schr¨ odinger equation for a charged particle in a magnetic field using the finite difference time domain method, Journal of Physics A: Mathematical and Theoretical, 40 (2007), p

2007
[3]

Performance and Analysis of the Alchemical Transfer Method for Binding-Free-Energy Predictions of Diverse Ligands

[15]N. Yenugu, A. K. Tiwari, and S. Sen,A grid-based gauge-invariant non-perturbative solution of the schr¨ odinger equation for diatomic molecules in strong magnetic fields, Journal of Chemical Theory and Computation, 21 (2025), p. 9753–9771, https://doi.org/10.1021/acs. jctc.5c00972, http://dx.doi.org/10.1021/acs.jctc.5c00972. [16]N. Yi, Y. Huang, and H...

work page doi:10.1021/acs 2025