---
title: Discrete Magnetic Laplacian
url: https://www.emergentmind.com/topics/discrete-magnetic-laplacian
type: topic
---

# Discrete Magnetic Laplacian

Searching arXiv for recent and foundational work on discrete magnetic Laplacians and closely related spectral settings.
Search query: "discrete magnetic Laplacian periodic graphs triangulations Sierpinski gasket"
A discrete magnetic Laplacian is a phase-twisted Laplace operator on a discrete geometric object—most commonly a weighted graph, but also a triangulation, a finite quotient of a periodic graph, or a graph approximation to a fractal—in which the magnetic field is encoded by antisymmetric edge phases or, in higher-dimensional simplicial settings, by holonomy through faces. In the graph-theoretic setting, the operator acts by replacing ordinary differences across an edge with covariant differences of the form \(f(x)-e^{i\theta_{x,y}}f(y)\), or equivalently by a twisted derivative \(d_\alpha\) and the factorization \(\Delta_\alpha=d_\alpha^*d_\alpha\) [1710.01157, 1507.02638]. Across the literature considered here, the subject develops along four tightly connected lines: gauge covariance and holonomy on graphs, Floquet-theoretic reductions of periodic operators to finite magnetic matrices, higher-degree magnetic Hodge theory on triangulations, and approximation frameworks in which graph magnetic operators converge to fractal or continuum magnetic operators [1808.07762, 2105.10171, 1604.01340].

## 1. Graph-theoretic definitions and gauge structure

On weighted graphs, the magnetic data are encoded by antisymmetric phases on oriented edges. In the weighted framework of discrete cusps and funnels, a graph is \(\mathcal G=(\mathcal E,\mathcal V,m)\) or \((\mathcal E,\mathcal V,m,\theta)\), where \(m:\mathcal V\to(0,\infty)\) is a vertex weight, \(\mathcal E:\mathcal V\times\mathcal V\to[0,\infty)\) is symmetric, and \(\theta_{x,y}\in \mathbb R/2\pi\mathbb Z\) satisfies \(\theta_{x,y}=-\theta_{y,x}\) and vanishes off the edge set. The associated quadratic form is
\[
Q_{\mathcal G,\theta}(f)=\frac12\sum_{x,y\in\mathcal V}\mathcal E(x,y)\,\left|f(x)-e^{i\theta_{x,y}}f(y)\right|^2,
\]
and the magnetic Laplacian is the nonnegative self-adjoint operator associated with this form; on finitely supported functions,
\[
(\Delta_{\mathcal G,\theta}f)(x)=\frac1{m(x)}\sum_{y\in\mathcal V}\mathcal E(x,y)\bigl(f(x)-e^{i\theta_{x,y}}f(y)\bigr)
\]
[1507.02638, 2507.05766].

A closely related formulation uses a twisted derivative. For an oriented weighted graph \(\mathbf W=(\mathbf G,m)\) with vector potential \(\alpha:E\to\mathbb T\), the twisted derivative is
\[
(d_\alpha f)_e = e^{i\alpha_e/2}f(\partial_+e)-e^{-i\alpha_e/2}f(\partial_-e),
\]
and the discrete magnetic Laplacian is
\[
\Delta_\alpha=d_\alpha^*d_\alpha.
\]
On vertices it is written
\[
\left( \Delta_\alpha f\right) (v) = \rho(v)f(v)-\frac1{m(v)}\sum_{e\in E_v} e^{\,i\alpha_e(v)}f(v_e)\,m_e,
\]
with \(\rho(v)=m(E_v)/m(v)\) [1710.01157].

Gauge invariance is a structural principle in every setting discussed here. On weighted graphs, two magnetic potentials with the same holonomy are gauge equivalent, and the corresponding magnetic Laplacians are unitarily equivalent [1507.02638]. In the periodic-graph setting, if \(\alpha'=\alpha+d\varphi\), then \(\Delta_\alpha\) and \(\Delta_{\alpha'}\) are unitarily equivalent via multiplication by \(e^{i\varphi}\) [1710.01157]. In the continuous magnetic Schrödinger framework, the same principle appears as
\[
\mathbf A'=\mathbf A+\nabla \phi,\qquad H_{\mathbf A'}=e^{i\phi}H_{\mathbf A}e^{-i\phi},
\]
a fact explicitly identified as transferable to the discrete case [1405.7912].

The gauge-invariant content is holonomy or flux through cycles. On periodic graphs, for a cycle \(c\) in the cycle space \(\mathcal C\) of the fundamental graph,
\[
\Phi_\alpha(c)=\left(\sum_{\mathbf e\in c}\alpha(\mathbf e)\right)\bmod 2\pi\in(-\pi,\pi]
\]
[1808.07762]. On weighted graphs with cusp geometry, if
\[
\gamma=(x_0,x_1)+\cdots +(x_{N-1},x_N),\qquad x_N=x_0,
\]
then
\[
\operatorname{Hol}_\theta(\gamma)=\theta_{x_0,x_1}+\cdots+\theta_{x_{N-1},x_0}\in \mathbb R/2\pi\mathbb Z,
\]
and the magnetic potential is determined up to gauge by this holonomy map [1507.02638]. This identifies the discrete magnetic Laplacian as an operator depending on cycle fluxes rather than on a particular edge-wise representative.

## 2. Discrete geometric settings: graphs, triangulations, and fractals

The most elementary setting is the graph magnetic Laplacian on vertices. For periodic discrete graphs \(\mathcal G=(\mathcal V,\mathcal E)\), the combinatorial magnetic Laplacian is
\[
(\Delta_\alpha f)(v)=\sum_{e=(v,u)\in\mathcal A}\Bigl(f(v)-e^{i\alpha(e)}f(u)\Bigr),
\]
or equivalently
\[
(\Delta_\alpha f)(v)=\varkappa_v f(v)-\sum_{e=(v,u)\in\mathcal A}e^{i\alpha(e)}f(u),
\]
where \(\varkappa_v\) is the degree of \(v\) [1808.07762]. In weighted settings the vertex weight \(m(x)\) appears explicitly, as above [1507.02638, 2507.05766]. These formulations are not identical as operators, but they share the same phase-twisted neighbor-interaction mechanism.

A higher-dimensional simplicial version is developed for weighted triangulations. A magnetic triangulation is
\[
T_\alpha=(K_\alpha,F),
\]
where \(K_\alpha=(V,E,\alpha)\) is a weighted magnetic graph and \(F\) is a symmetric set of triangular faces. The cochain spaces are weighted \(\ell^2\)-spaces on vertices, skew-symmetric edge functions, and skew-symmetric face functions, and the total Hilbert space is
\[
l^2(T_\alpha)=l^2(V)\oplus l^2(E)\oplus l^2(F).
\]
The magnetic derivative on \(0\)-forms is
\[
d_{\alpha}^{0}g(x,y):=e^{\frac{i\alpha(y,x)}{2}}g(y)-e^{\frac{i\alpha(x,y)}{2}}g(x),
\]
and the magnetic exterior derivative on \(1\)-forms is
\[
d_{\alpha}^{1}(\varphi)(x,y,z):= e^{\frac{i}{6}(\alpha(x,z)+\alpha(y,z))} \varphi(x,y)+ e^{\frac{i}{6}(\alpha(y,x)+\alpha(z,x))} \varphi(y,z)+ e^{\frac{i}{6}(\alpha(z,y)+\alpha(x,y))} \varphi(z,x).
\]
From these operators one obtains the magnetic Gauss–Bonnet operator
\[
T_\alpha= \begin{pmatrix} 0&\delta^0_{\alpha}&0 \\ d^{0}_{\alpha}&0&\delta^{1}_{\alpha} \\ 0&d^{1}_{\alpha}&0 \end{pmatrix},
\qquad \Delta_\alpha:=T_\alpha^2
\]
[2105.10171].

The triangulation setting differs from the graph setting in a decisive respect: the magnetic field becomes face holonomy. For a face,
\[
\widehat{\alpha}_{xyz}=Hol_\alpha(x,y,z,x)=\alpha_{xy}+\alpha_{yz}+\alpha_{zx}.
\]
Then \(d_\alpha^1d_\alpha^0\) need not vanish, and this failure is controlled explicitly by \(\sin(\widehat\alpha/6)\). The magnetic Hodge Laplacian therefore acquires off-diagonal terms \(d_\alpha^1d_\alpha^0\) and \(\delta_\alpha^0\delta_\alpha^1\), absent in the non-magnetic case [2105.10171]. This is the simplicial analogue of curvature entering the de Rham complex.

A different extension appears on the Sierpiński gasket. There the magnetic operator is defined on a Hilbert space \(H\) of \(1\)-forms with derivation \(\partial:F\to H\), using a real-valued \(1\)-form \(a\) as magnetic potential:
\[
(\partial+ia):F\to H,\qquad E^a(f)=\|(\partial+ia)f\|_H^2.
\]
The associated Dirichlet magnetic Laplacian \(M^a\) is defined through the closed form \((E^a,F_0)\). At graph level, the finite approximants carry magnetic energies
\[
E_m^{a_m}(f)=\sum_{x,y:x\sim_m y}\Bigl|f(x)-f(y)e^{i a_m(e_{xy})}\Bigr|^2
\]
and magnetic Laplacians
\[
M_m^{a_m}f(x) = -\sum_{y:y\sim_m x}\Bigl(f(x)-f(y)e^{i a_m(e_{xy})}\Bigr),
\]
with renormalized convergence to the fractal operator [1604.01340].

## 3. Periodic graphs and Floquet reduction

On periodic graphs, the discrete magnetic Laplacian is inseparable from Floquet theory. For a connected, locally finite \(\Gamma\)-periodic graph \(\mathcal G\subset\mathbb R^d\) with finite quotient \(G_*=\mathcal G/\Gamma\), the periodic magnetic operator admits a direct-integral decomposition over the quasimomentum torus \(\mathbb T^d\) [1808.07762]. In the standard fiber representation,
\[
UH_\alpha U^{-1}=\int_{\mathbb T^d}^{\oplus} H_{\tau,\alpha}(\vartheta)\,\frac{d\vartheta}{(2\pi)^d},
\]
with fiber magnetic Laplacians
\[
(\Delta_{\tau,\alpha}(\vartheta)f)(v)
= \varkappa_v f(v)-\sum_{\mathbf e=(v,u)\in A_*} e^{\,i(\alpha(\mathbf e)+(\tau(\mathbf e),\vartheta))}f(u).
\]

A refined version replaces the edge-index form \(\tau\) and the magnetic form \(\alpha\) by minimal forms \(m\in\mathscr F(\kappa)\) and \(\phi\in\mathscr F(\alpha)\), giving
\[
\mathscr U_{m,\phi}H_\alpha \mathscr U_{m,\phi}^{-1}
= \int_{\mathbb T^d}^{\oplus} H_{m,\phi}(\vartheta)\,\frac{d\vartheta}{(2\pi)^d},
\]
with
\[
(\Delta_{m,\phi}(\vartheta)f)(v)
= \varkappa_v f(v)-\sum_{\mathbf e=(v,u)\in A_*} e^{\,i(\phi(\mathbf e)+(m(\mathbf e),\vartheta))}f(u).
\]
Minimality here means support minimality among all \(1\)-forms with the same cycle fluxes [1808.07762]. The resulting integers
\[
\mathcal I=\frac12\#\operatorname{supp}m,\qquad
\mathcal I_\alpha=\frac12\#\operatorname{supp}\phi,
\]
are invariants of the periodic graph and of the magnetic Laplacian, respectively. They measure how many coefficients in the fiber matrices genuinely depend on quasimomentum and on magnetic potential [1808.07762].

A complementary formulation appears in the spectral-gap theory of periodic graphs. There the Floquet fibers of a periodic Laplacian on the infinite graph are identified with discrete magnetic Laplacians on the finite quotient graph. If \(\widetilde{\mathbf G}\to \mathbf G\) is a periodic covering and \(\mathcal A_{\mathbf H}\) is the set of quotient vector potentials with the lifting property
\[
e^{i\alpha_{[e]}}=\chi(\operatorname{ind}_{\mathbf H}(e)),
\]
then
\[
\sigma(\Delta^{\widetilde{\mathbf W}})=\bigcup_{\alpha\in\mathcal A_{\mathbf H}}\sigma(\Delta_\alpha^{\mathbf W}).
\]
In this sense, the vector potential on the quotient is a Floquet parameter [1710.01157]. The finite-dimensional magnetic Laplacian on the quotient is therefore the decisive object in band localization and gap detection.

The periodic viewpoint also yields quantitative spectral information. For the fiber eigenvalues
\[
\lambda_{\alpha,1}(\vartheta)\le \cdots \le \lambda_{\alpha,\nu}(\vartheta),
\]
the spectrum is the union of bands
\[
\sigma(H_\alpha)=\bigcup_{n=1}^{\nu}\sigma_n(H_\alpha),\qquad
\sigma_n(H_\alpha)= [\lambda_{\alpha,n}^-,\lambda_{\alpha,n}^+].
\]
The total spectral measure satisfies
\[
|\sigma(H_\alpha)|\le 4\mathcal I,
\]
and each band obeys
\[
\sigma_n(H_\alpha)\subset [\lambda_n^0,\lambda_n^0+2\varkappa_+^m]
\]
[1808.07762]. These estimates are sharper than bounds using only the Betti number \(\beta=\#E_*-\#V_*+1\).

## 4. Spectral phenomena: gaps, localization, cusps, and funnels

One major theme is the relation between magnetic flux and spectral gaps. In the periodic-gap framework, edge and vertex virtualization produce comparison operators \(\Delta^{\mathbf W^-}\) and \(\Delta^{\mathbf W^+}\) satisfying
\[
\Delta^{\mathbf W^-}\preccurlyeq \Delta_\alpha^{\mathbf W}\preccurlyeq \Delta^{\mathbf W^+},
\]
provided the magnetic potential is supported on a selected edge set and a suitable neighboring vertex set is virtualized [1710.01157]. Consequently,
\[
\sigma(\Delta_\alpha^{\mathbf W}) \subset \bigcup_{k=1}^{|V|}[\lambda_k(\Delta^{\mathbf W^-}),\lambda_k(\Delta^{\mathbf W^+})].
\]
This yields a geometric criterion for nonempty magnetic spectral-gap sets. If \(v_0\) is a center vertex with cycle edges \(A(v_0)\), then the quantity
\[
\delta = \rho(v_0) - \sum_{e\in A(v_0)}\frac{m_e}{m((v_0)_e)} - \frac{m(A(v_0))}{m(v_0)}
\]
satisfies: the Lebesgue measure of the magnetic spectral-gap set \(MS^{\mathbf W}\) is at least \(\delta\); in particular, \(\delta>0\) implies \(MS^{\mathbf W}\neq\emptyset\) [1710.01157].

On discrete cusps, the magnetic field can change the global spectral type. A discrete cusp is a twisted product
\[
\mathcal G=\mathcal G_1\times_{\mathcal I}\mathcal G_2
\]
with \(m_1(x)\to 0\) at infinity, \(\mathcal G_2\) finite, and \(\Delta_{\mathcal G_1,\theta_1}\) bounded. In this setting,
\[
\Delta_{\mathcal G,\theta}
= \Delta_{\mathcal G_1,\theta_1}\otimes \frac{1_{\mathcal I}(\cdot)}{m_2(\cdot)}
+\frac1{m_1(\cdot)}\otimes \Delta_{\mathcal G_2,\theta_2}.
\]
The key theorem states that
\[
\Delta_{\mathcal G,\theta}\text{ has compact resolvent}
\quad\Longleftrightarrow\quad \operatorname{Hol}_{\theta_2}\neq 0.
\]
When the transversal holonomy is nontrivial,
\[
\mathcal D(\Delta_{\mathcal G,\theta}^{1/2})
=\mathcal D(\deg_{\mathcal G}^{1/2}(\cdot)),
\qquad
\lim_{n\to\infty}
\frac{\lambda_n(\Delta_{\mathcal G,\theta})}
{\lambda_n(m_1^{-1}(\cdot)\otimes \Delta_{\mathcal G_2,\theta_2})}
=1.
\]
When \(\operatorname{Hol}_{\theta_2}=0\), the form domain can differ from the degree form domain [1507.02638]. The magnetic field therefore affects not only the spectrum but also the energy space.

Discrete funnels exhibit a different asymptotic regime. The model graph is a twisted product \(G=G_1\times_3 G_2\) with
\[
\mathcal V_1=\mathbb N,\qquad m_1(n)=e^n,\qquad \mathcal E_1(n,n+1)=e^{(2n+1)/2}.
\]
The funnel Laplacian decomposes as
\[
\Delta_G = \Delta_{G_1}\otimes \frac1{m_2(\cdot)} + \frac1{m_1(\cdot)}\otimes \Delta_{G_2},
\]
so the transverse term is asymptotically negligible because \(1/m_1(n)=e^{-n}\to0\) [2507.05766]. The essential spectrum is determined by the radial channel:
\[
\alpha=e^{1/2}+e^{-1/2}-2,\qquad
\beta=e^{1/2}+e^{-1/2}+2,
\]
and
\[
\sigma_{\mathrm{ess}}(\Delta_G)=[\alpha,\beta]\cdot \frac1{m_2(V_2)}.
\]
Under long-range perturbations of weights, phases, and an electric potential, the perturbed operator satisfies a Mourre estimate and a Limiting Absorption Principle away from thresholds and embedded eigenvalues; singular continuous spectrum is empty [2507.05766]. In this model, the radial magnetic potential is gauge removable because the radial graph is a half-line.

The fractal case shows another spectral pattern. For the magnetic Laplacian \(M^a\) on the Sierpiński gasket, the operator has compact resolvent and pure point spectrum accumulating at \(+\infty\). Its eigenvalue counting asymptotics are unchanged by the magnetic field:
\[
\lim_{x\to\infty} \rho^a(x)x^{-\log 3/\log 5} -\chi(\log x)=0,
\]
with the same periodic function \(\chi\) as for the ordinary Laplacian [1604.01340]. At the same time, low-generation eigenvalues can shift with the flux, as illustrated by
\[
M_1^{\beta b}f_k
=\Bigl(4-2\cos\Bigl(\frac{2\pi k}{3}+\frac{2\beta}{\sqrt{30}}\Bigr)\Bigr)f_k
\]
on the first graph approximation [1604.01340].

## 5. Self-adjointness, domains, and higher-degree structure

Essential self-adjointness is a central analytical issue for discrete magnetic operators on noncompact discrete geometries. On weighted triangulations, the decisive completeness hypothesis is \(\chi\)-completeness. It requires cutoff functions \(\chi_n\in C^c(V)\), \(0\le \chi_n\le 1\), equal to \(1\) on an exhaustion \(B_n\subset V\), together with the uniform estimates
\[
\frac{1}{c(x)}\sum_{e\in E,\ e^+=x} r(e)|d^0\chi_n(e)|^2\le C \qquad (C_1),
\]
and
\[
\frac{1}{r(x,y)}\sum_{t\in F_{xy}} s(x,y,t)\left|d^0\chi_n(t,x)+d^0\chi_n(t,y)\right|^2\le C \qquad (C_2).
\]
Under this hypothesis, the magnetic Gauss–Bonnet operator \(T_\alpha\) is essentially self-adjoint on compactly supported cochains, and therefore so is the magnetic Laplacian \(\Delta_\alpha=T_\alpha^2\) [2105.10171].

The same paper introduces a bounded-curvature condition on the magnetic potential:
\[
\forall x\in V,\qquad
\frac{1}{c(x)}\sum_{e\in F_x}s(x,e)\sin^2\left(\frac{\widehat{\alpha}(x,e)}{6}\right)\le C.
\]
Under \(\chi\)-completeness and bounded curvature, the closure of \(T_\alpha\) has the domain decomposition
\[
Dom (T_{\alpha,\min})=Dom (d^{0}_{\alpha,\min})\oplus\left(Dom (\delta^{0}_{\alpha,\min})\cap Dom (d^{1}_{\alpha,\min})\right) \oplus Dom (\delta^{1}_{\alpha,\min}),
\]
which separates the vertex, edge, and face components [2105.10171]. The appearance of the face holonomy in the curvature bound is intrinsically two-dimensional.

In the cusp setting, domain theory is more delicate because the magnetic field can eliminate or preserve a transversal zero mode. The paper on discrete cusps proves that
\[
\mathcal D(\deg_{\mathcal G}^{1/2}(\cdot))\subset \mathcal D(\Delta_{\mathcal G,\theta}^{1/2}),
\]
but equality may fail when the transversal holonomy vanishes [1507.02638]. This establishes that magnetic and non-magnetic form domains can differ. The result is unusual enough that the paper explicitly identifies it as one of its most striking conclusions.

In the funnel setting, the unperturbed form domain is given explicitly as
\[
\mathcal G=D(\Delta_G^{1/2})
= \ell^2(\mathcal V_1,m_1)\,\widehat\otimes\, D\!\left(\frac1{\sqrt{m_2(\cdot)}}\right),
\]
and the perturbed analysis is carried out in a form-commutator framework \(C^{1,1}(A;\mathcal G,\mathcal G^*)\), rather than through a purely operator-domain computation [2507.05766]. This suggests that in weighted magnetic graph models with long-range perturbations, form methods can be more robust than direct operator-domain arguments.

## 6. Conceptual boundaries and relation to continuous magnetic Laplacians

The expression “discrete magnetic Laplacian” does not refer to a single theory. In the sources considered here it can denote at least three distinct but related constructions. First, it denotes graph magnetic Laplacians with phase-twisted differences on vertices [1710.01157, 1808.07762, 1507.02638, 2507.05766]. Second, it denotes higher-degree magnetic Hodge Laplacians on simplicial complexes, where face holonomy becomes curvature and the magnetic Gauss–Bonnet operator is the fundamental object [2105.10171]. Third, it can denote graph approximants to a non-Euclidean limit operator, as on the Sierpiński gasket, where finite graph magnetic Laplacians converge after renormalization to a fractal magnetic Laplacian [1604.01340].

By contrast, some closely related literature concerns continuous magnetic Laplacians rather than discrete ones. The monograph “Little Magnetic Book” studies the continuous magnetic Schrödinger operator
\[
H_{\mathbf A}=(-i\nabla-\mathbf A)^2,\qquad H_h(\mathbf A)=(-ih\nabla-\mathbf A)^2,
\]
with an emphasis on discrete spectrum, essential spectrum, localization, effective Hamiltonians, Born–Oppenheimer approximation, and norm resolvent convergence [1405.7912]. This work is not about graph magnetic Laplacians, but it supplies a continuous spectral framework that informs discrete models.

Likewise, the study of curved waveguides in a strong uniform magnetic field analyzes the continuous operator
\[
\mathscr P_h = (-ih\nabla-\mathbf A)^2-h
\]
on a curved strip with Dirichlet boundary conditions, proving existence of discrete spectrum below the essential spectrum under a curvature condition [2205.12733]. Here “discrete” refers to the discrete spectral component, not to a discrete underlying space. The distinction is terminologically important: a discrete magnetic Laplacian in graph theory is a magnetic operator on a discrete combinatorial or weighted geometry, whereas a continuous magnetic Laplacian may possess discrete spectrum without being discrete in this combinatorial sense [2205.12733].

A plausible implication is that the common language of gauge invariance, model operators at infinity, holonomy, and localization allows continuous and discrete magnetic theories to inform one another without collapsing their distinction. In the discrete setting, this has already produced a broad landscape: phase-twisted graph Laplacians, Floquet-reduced finite magnetic matrices on periodic graphs, magnetic Hodge theory on triangulations, and renormalized graph approximations to fractal magnetic operators [1808.07762, 1710.01157, 2105.10171, 1604.01340].

Source: https://www.emergentmind.com/topics/discrete-magnetic-laplacian