A Gauge-Invariant Finite Element Method for the Time-Dependent Ginzburg--Landau Equations
Abstract
It is well known that the Ginzburg--Landau equations are gauge invariant. The most popular gauge choice in practical simulations is the temporal gauge, also known as the zero electric potential gauge. However, the commonly used Nédélec approximations may produce spurious giant vortices on domains with re-entrant corners. To address this issue while retaining the Nédélec approximation, we propose a finite element scheme that preserves the gauge invariance of the time-dependent Ginzburg--Landau equations at the discrete level. Under the temporal gauge, the method employs lowest-order Nédélec edge elements for the magnetic potential and continuous piecewise-linear elements for the order parameter on triangular and tetrahedral meshes, and yields a linear, decoupled formulation that requires only the solution of linear systems at each time step. We prove that the scheme is gauge invariant, preserves a discrete maximum bound principle for the order parameter, |ψhn|L∞(Ω)≤ 1, and satisfies a discrete energy law under mild step-size and mesh conditions. Numerical experiments confirm these theoretical properties. In particular, on domains with re-entrant corners, the scheme computes the correct physical solution, whereas conventional temporal-gauge discretizations with Nédélec elements produce an unphysical solution.
Related
- Bounds-Constrained Finite Element Approximation of Time-Dependent Partial Differential Equations
- Grad-div stabilized finite element Picard-Newton iterative methods for the Navier-Stokes equations
- A dual-field structure-preserving mixed finite element discretization for incompressible Hall MHD equations
- Robust preconditioners for finite element approximations of control-state constrained optimal control problems governed by elliptic equations
- GridapGeosciences.jl: A Julia finite element package for partial differential equations on general manifolds
BibTeX
@misc{gao2026gaugeinvariant,
title = {{A Gauge-Invariant Finite Element Method for the Time-Dependent Ginzburg--Landau Equations}},
author = {Huadong Gao and Yuhui Huang and Wen Xie},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.24470},
eprint = {2610.24470},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.24470},
}
RIS
TY - UNPB
TI - A Gauge-Invariant Finite Element Method for the Time-Dependent Ginzburg--Landau Equations
AU - Huadong Gao
AU - Yuhui Huang
AU - Wen Xie
PY - 2026
DA - 2026/10/17
PB - arXiv
DO - 10.48550/arxiv.2610.24470
UR - https://arxiv.org/abs/2610.24470
AB - It is well known that the Ginzburg--Landau equations are gauge invariant. The most popular gauge choice in practical simulations is the temporal gauge, also known as the zero electric potential gauge. However, the commonly used Nédélec approximations may produce spurious giant vortices on domains with re-entrant corners. To address this issue while retaining the Nédélec approximation, we propose a finite element scheme that preserves the gauge invariance of the time-dependent Ginzburg--Landau equations at the discrete level. Under the temporal gauge, the method employs lowest-order Nédélec edge elements for the magnetic potential and continuous piecewise-linear elements for the order parameter on triangular and tetrahedral meshes, and yields a linear, decoupled formulation that requires only the solution of linear systems at each time step. We prove that the scheme is gauge invariant, preserves a discrete maximum bound principle for the order parameter, |ψ_h^n|_L^∞(Ω)≤ 1, and satisfies a discrete energy law under mild step-size and mesh conditions. Numerical experiments confirm these theoretical properties. In particular, on domains with re-entrant corners, the scheme computes the correct physical solution, whereas conventional temporal-gauge discretizations with Nédélec elements produce an unphysical solution.
ER -