---
title: Weak Localization of Bosons
url: https://www.emergentmind.com/topics/weak-localization-of-bosons
type: topic
---

# Weak Localization of Bosons

Searching arXiv for the cited works to ground the article in current records.
Weak localization of bosons denotes interference-driven localization phenomena in disordered bosonic systems, with manifestations that depend strongly on the platform and observable. In quasi-two-dimensional superconducting films, the localized carriers can be Cooper pairs of charge \(2e\), and the characteristic low-temperature correction appears as a logarithm in the **resistance**, \(R(T)=R_0+A_b\ln(T/T_0)\), rather than the familiar fermionic correction in the conductance [1002.1720]. In disordered cold-atom settings, the same interference physics is formulated through crossed diagrams, coherent backscattering, and disorder-averaged transport equations, while in quasiperiodic lattice systems it appears through localization diagnostics such as the inverse participation ratio, superfluid fraction, participation ratio, and decay of one-body correlations [1307.4940] [1506.07692] [1610.05398] [2502.14440]. Across these realizations, weak localization is distinguished from purely classical trapping by its dependence on phase coherence, reversed-path interference, and disorder-controlled crossover or transition scales.

## 1. Conceptual basis and distinction from fermionic weak localization

In the conventional two-dimensional fermionic case, weak localization gives a correction to the sheet conductance of the form
\[
G(T)=G_0-\frac{e^2}{2\pi^2\hbar}\ln(T/T_0),
\]
or more generally
\[
G(T)=G_0+A\ln T,
\]
with \(A<0\) set by electron-electron interaction and spin-orbit coupling in the Abrahams-Aronov framework [1002.1720]. Das and Doniach predicted that a dilute “gas” of Cooper pairs subject to disorder and a perpendicular field should instead show a weak-localization correction appearing directly in the resistance,
\[
R(T)=R_0+A_b\ln(T/T_0),
\]
with \(A_b>0\) arising from interference of bosonic paths; in the strong-field bosonic regime one also expects a field dependence of \(A_b(B)\) and a localization length \(\xi_b\) set by pair size and disorder [1002.1720].

In microscopic transport theory for interacting bosons in a weak random potential, the disorder-averaged dynamics can be reduced to a nonlinear integral equation of Boltzmann type for the single-particle diffusive flux. The relevant building blocks are the propagation kernel \(P_E\), an elastic collision kernel \(g_{E_1,E_2}\), and an inelastic collision kernel \(f_{E_1,E_2,E_3}\), while weak localization enters through crossed, or Cooperon, corrections rather than the ladder sector that describes incoherent diffusion [1307.4940]. In a complementary semiclassical treatment of nonlinear bosonic matter waves in billiard geometries, ladder pairings produce no coherent backscattering, whereas crossed pairings of a trajectory and its time reverse generate the weak-localization peak [1112.5603].

These formulations establish a central distinction: for bosons, weak localization is not tied to a single universal observable. In superconducting films it appears as \(R\sim\ln T\); in disordered matter-wave transport it appears in the coherent backscattering cone; and in lattice realizations it is tracked by fragmentation and localization diagnostics. This suggests that “weak localization of bosons” is best understood as a family of interference phenomena rather than a single transport law.

## 2. Quasi-two-dimensional superconducting films and the \(R\sim\ln T\) signature

A direct experimental realization was reported for ultrathin amorphous Bi films grown by in-situ quench-condensation at \(7\,\mathrm K\) under UHV on a polished \((100)\) SrTiO\(_3\) crystal with ex-situ deposited sputtered Pt four-probe contacts of \(100\,\text{\AA}\). A \(10\,\text{\AA}\) amorphous Sb wetting underlayer was deposited first, followed by amorphous Bi to a total thickness of \(11.2\,\text{\AA}\). From the Hall-effect value \(n_{2D}\approx1.5\times10^{16}\,\mathrm{cm}^{-2}\), one obtains \(k_F\ell\approx6\) and \(\ell\approx2.4\,\text{\AA}\), with normal-state sheet resistance just above the superconducting transition \(R_\square\approx3384\,\Omega\). Measurements were performed in DC four-probe geometry in a dilution refrigerator down to \(T\approx100\,\mathrm{mK}\), with perpendicular field up to \(10\,\mathrm T\) and extensive lead filtering to minimize electromagnetic noise [1002.1720].

The phase identification was operational. Superconducting behavior was defined by \(dR/dT>0\) at the lowest \(T\) with \(R\to0\); insulating behavior by \(dR/dT<0\) at the lowest \(T\) with \(R\to\infty\); and a mixed or nonmonotonic regime by a resistance maximum at \(T_{\max}(B)\) followed by a low-temperature downturn, interpreted as superconducting islands in an insulating matrix [1002.1720].

At zero field and for \(B\le1\,\mathrm T\), \(T>2\,\mathrm K\), the sheet conductance obeys
\[
G(T)=G_0+A_f\ln T,
\]
with \(A_f\approx5.23\times10^{-6}\,\Omega^{-1}\,\square^{-1}\), in quantitative agreement with weak-localization plus electron-electron interaction theory for strong spin-orbit scattering [1002.1720]. In the intermediate temperature window \(2\,\mathrm K<T<4\,\mathrm K\) and fields \(0.6\,\mathrm T<B<5\,\mathrm T\), the sheet resistance collapses under Fisher’s quantum-critical form,
\[
R(B,T)/R_c = F(|B-B_c|/T^{1/\nu z}),
\]
with \(R_c=3384\,\Omega\), \(B_c=1.58\,\mathrm T\), and \(\nu z\approx0.4\). This scaling breaks down for \(T<2\,\mathrm K\), where \(R(T)\) becomes strongly nonmonotonic [1002.1720].

Deep on the insulating side, for \(B\ge2.5\,\mathrm T\) and \(T<1\,\mathrm K\), the key result is that the **resistance itself**, not the conductance, is a clear logarithm of temperature over \(0.1\)–\(1\,\mathrm K\) and up to \(10\,\mathrm T\):
\[
R(T)=R_0+A_b\ln T.
\]
At \(B=10\,\mathrm T\), the best fit gives \(A_b\approx1.13\times10^3\,\Omega\,\square^{-1}\) per \(\ln\) decade with \(\chi^2\approx1.3\times10^3\), whereas fitting \(G\sim\ln T\) gives \(\chi^2\approx8.2\times10^3\). Converting to bulk resistivity using film thickness \(\approx11.2\,\text{\AA}\) yields \(A_b^{3D}\approx-6.9\times10^{-5}\,\Omega\,\mathrm{cm}\) per \(\ln T\), close to the value \(-7.2\times10^{-5}\,\Omega\,\mathrm{cm}\) reported for underdoped cuprates. By contrast, forcing \(G\sim\ln T\) in this regime yields \(A_f'\approx1.63\times10^{-5}\,\Omega^{-1}\,\square^{-1}\), about three times the zero-field value and therefore inconsistent with electronic quantum corrections [1002.1720].

The interpretation is specific: \(R\sim\ln T\), rather than \(G\sim\ln T\), with a coefficient of order \(10^{-5}\,\Omega\,\mathrm{cm}\) over a broad field and temperature range, is taken as the hallmark of bosonic weak localization. In this picture, the localized carriers are Cooper pairs, and their interference produces a positive back-reflection correction to conductance, which appears experimentally as increasing \(R(T)\) on cooling. The absence of a large magnetoresistance peak was noted to argue against simple boson-pair hopping or Josephson-junction-array physics and to suggest weakly localized bosons rather than a true Bose glass [1002.1720].

## 3. Microscopic transport theory and coherent backscattering of interacting bosons

A microscopic scattering theory for interacting bosons in a three-dimensional weakly disordered potential begins from
\[
H=H_0+V+U,
\]
with \(H_0=\int d^3k\,E_k\,a_k^\dagger a_k\), \(E_k=k^2\) in units \(\hbar^2/2m\equiv1\), a Gaussian white-noise disorder potential satisfying \(\langle V(r)V(r')\rangle=(4\pi/\ell_{\mathrm{dis}})\delta(r-r')\), and a two-body interaction \(U(r_1-r_2)\) [1307.4940]. The disorder-averaged single-particle Green function is diagonal in momentum,
\[
\langle k'|G_V(E)|k\rangle=(2\pi)^3\delta(k-k')G_E(k),\qquad
G_E(k)=\frac{1}{E-k^2+i/(2\ell_{\mathrm{dis}})},
\]
where \(\ell_{\mathrm{dis}}\) is the disorder mean free path [1307.4940]. At low energy, the two-body \(T\)-matrix satisfies
\[
\langle k'|T_U(E)|k\rangle \simeq 8\pi a_s\left[1-i\sqrt{E/2}\,a_s+\dots\right],
\]
with \(a_s\) the \(s\)-wave scattering length [1307.4940].

After disorder averaging and tracing over \(N-1\) atoms, the transport problem is expressed in terms of a propagation kernel
\[
P_E(r,r')=\frac{e^{-|r-r'|/\ell_{\mathrm{dis}}}}{4\pi \ell_{\mathrm{dis}}|r-r'|^2},
\]
together with elastic and inelastic collision kernels \(g\) and \(f\). The resulting nonlinear transport equation for the disorder-averaged spectral density \(I_E(r)\) contains an incoming source term, the propagation term \(\int dr'P_E(r,r')I_E(r')\), and quadratic collision terms in \(I_E\) built from \(g\) and \(f\) [1307.4940]. The formalism reduces to the Gross-Pitaevskii mean-field description when only elastic collisions are retained, but even weak interactions generate inelastic energy redistribution and drive thermalization of the single-particle current [1307.4940].

Weak localization is encoded in the crossed sector. The crossed propagator is
\[
P_E^{(C)}(r,r')=(4\pi/\ell_{\mathrm{dis}})\int d^3k_1d^3k_2\,e^{-i(k_1-k_2)\cdot(r-r')}
G_E(k_1)G_{\tilde E}^*(k_2),
\]
with \(E+\tilde E=E_i+E_d\), and the crossed density \(C_E(r)=C_E^{(1)}(r)+C_E^{(2)}(r)\) obeys coupled integral equations whose solution determines the interference part of the backscattered flux [1307.4940]. In slab far field, the coherent backscattering contribution takes the form
\[
\Delta I(\theta)=\int \frac{d^3Q}{(2\pi)^3}\,C(Q)e^{iQ\cdot(\hat n_d-\hat n_i)L},
\]
with static Cooperon
\[
C(Q)=\frac{\ell_{\mathrm{dis}}^2}{4\pi \ell_{\mathrm{tr}}(Q^2+1/\ell_{\mathrm{tr}}^2)},
\qquad \ell_{\mathrm{tr}}\simeq \ell_{\mathrm{dis}}.
\]
This is the explicit bosonic counterpart of weak-localization interference in the language of coherent backscattering [1307.4940].

The inelastic channel modifies the coherent backscattering cone in a highly structured way. Inelastic scattering broadens and changes the amplitude of the cone in a narrow spectral window around \(E_d\approx E_i\), with width \(\Delta E/E_i\sim1/(\ell_{\mathrm{dis}}\sqrt{E_i})\) and enhancement greater than \(2\) at the peak energy. At \(\theta=0\), purely elastic nonlinearity yields a crossover from constructive to destructive coherent backscattering as the dephasing parameter \(\beta\) grows, whereas inelastic collisions slow the suppression of the cone and leave a residual enhancement due to many-wave interference in the inelastic channel [1307.4940]. This makes clear that interactions do not merely wash out bosonic interference; they can also reshape it spectrally.

## 4. Nonlinear bosonic matter waves and dephasing of weak localization

For dilute Bose-Einstein condensates propagating through irregularly shaped billiards connected to waveguides, the relevant mean-field equation is the inhomogeneous Gross-Pitaevskii equation
\[
i\hbar\partial_t\Psi(r,t)=\left[-\frac{\hbar^2}{2m}\left(\nabla-\frac{iA(r)}{\hbar}\right)^2+V(r)+g\frac{\hbar^2}{2m}|\Psi(r,t)|^2\right]\Psi(r,t)+S(r,t),
\]
where \(A(r)\) is an artificial vector potential breaking time-reversal symmetry, \(V(r)\) is the hard-wall billiard potential, and \(g\) is the effective two-dimensional interaction strength [1112.5603]. For a monochromatic source \(S(r,t)=S_0\chi_i(y)\delta(x-x_L)e^{-i\mu t/\hbar}\) and stationary ansatz \(\Psi(r,t)=\psi(r)e^{-i\mu t/\hbar}\), one obtains the stationary nonlinear equation
\[
(\mu-H)\psi(r)-g\frac{\hbar^2}{2m}|\psi(r)|^2\psi(r)=S_0\chi_i(y)\delta(x-x_L),
\]
with \(H=[-i\hbar\nabla-A(r)]^2/(2m)+V(r)\) [1112.5603].

The linear retarded Green function admits the Van Vleck form
\[
G_0(r,r';E)=\sum_\gamma A_\gamma
\exp\!\left[iS_\gamma(r,r';E)/\hbar-i\pi\mu_\gamma/2\right],
\]
where \(\gamma\) labels classical trajectories and \(S_\gamma\) is the action [1112.5603]. A weak vector potential contributes a dephasing phase \(\phi_\gamma\), while weak nonlinearity leads, after self-consistent iteration restricted to ladder pairings, to a dressed Green function with an additional phase
\[
\chi_\gamma=g(\hbar^2/m)\int_0^T |\psi^{(0)}(q_\gamma(t))|_d^2\,dt.
\]
In the diagrammatic interpretation, ladder pairings shift the density uniformly and do not produce weak localization, whereas crossed pairings of \(\gamma\) and its time reverse \(\tilde\gamma\) generate the coherent-backscattering correction [1112.5603].

The resulting weak-localization correction for retroreflection in the incident channel \(i\) is
\[
r_{ii}(B,g)=\tau_D/\tau_H+c_{ii}^{(g)},
\]
with
\[
c_{ii}^{(g)}=
\frac{\tau_D/\tau_H}{1+\tau_D/\tau_B+(g\,j^{\mathrm{in}}\tau_D)^2}.
\]
Equivalently, introducing the nonlinearity-induced dephasing rate \(\gamma_\phi(g)\equiv g\,j^{\mathrm{in}}\), one may write
\[
c_{ii}(B,g)=
\frac{\tau_D/\tau_H}{1+(\tau_D/\tau_B)+[\gamma_\phi(g)\tau_D]^2}.
\]
For total reflection and transmission,
\[
R=N_c\tau_D/\tau_H+c_{ii}^{(g)},\qquad
T=\tilde N_c\tau_D/\tau_H-c_{ii}^{(g)},
\]
or, in Lorentzian form,
\[
R=\frac{N_c}{N_c+\tilde N_c}+
\frac{\tilde N_c}{(N_c+\tilde N_c)^2}\,
\frac{1}{1+(B/B_0)^2+(g/g_0)^2},
\]
\[
T=\frac{\tilde N_c}{N_c+\tilde N_c}-
\frac{\tilde N_c}{(N_c+\tilde N_c)^2}\,
\frac{1}{1+(B/B_0)^2+(g/g_0)^2},
\]
where \(B_0\) sets the magnetic width and \(g_0=[j^{\mathrm{in}}\tau_D(N_c+\tilde N_c)]^{-1}\) sets the nonlinear dephasing scale [1112.5603].

Numerically, for hard-wall chaotic billiards with \(N_c=\tilde N_c=5\) at \(\mu=\mu_0\), finite-difference discretization at roughly \(30\) points per de Broglie wavelength, and nonlinearity confined to the cavity, the computed \(r_{ii}(B)\) and \(T(B)\) for \(g=0\ldots0.25\) show Lorentzian weak-localization peaks whose height scales as \(1/(1+(g/g_0)^2)\) [1112.5603]. A further nuance is that an apparent weak antilocalization regime can arise from short self-retracing paths, which form a non-universal reflection background and reduce the zero-field peak. In this formulation, the effect is attributed not to a universal sign reversal of bosonic interference, but to omitted short-path contributions in the Cooperon sum [1112.5603].

## 5. Quasiperiodic lattices, fragmentation, and localization diagnostics

In one-dimensional quasiperiodic lattices, the weakly interacting Bose gas is often modeled by the Bose-Hubbard Hamiltonian with Aubry-André onsite modulation,
\[
H=
-\sum_{\langle i,j\rangle}(a_i^\dagger a_j+\mathrm{h.c.})
+\lambda\sum_i\cos(2\pi\beta i)\hat n_i
+\frac{U}{2}\sum_i \hat n_i(\hat n_i-1),
\]
with \(\beta=(\sqrt5-1)/2\) and, in the weakly interacting high-filling regime, a condensate order parameter \(\psi_i=\langle a_i\rangle\) governed by the discrete nonlinear Schrödinger equation [1506.07692]. Two standard diagnostics are the inverse participation ratio,
\[
I=\sum_{i=1}^{N_s}|\psi_i|^4,
\]
and the superfluid fraction
\[
f_s=N_s^2\frac{E_{\rm cl}(\Theta)-E_{\rm cl}(0)}{\Theta^2}.
\]
For \(U=0\), the model is self-dual at \(\lambda_c=2\): all states are extended for \(\lambda<2\) and exponentially localized for \(\lambda>2\). With weak repulsion, the IPR becomes appreciable at \(\lambda\approx2\) but rises only gradually because the condensate fragments over several quasi-degenerate minima, a phenomenon described as multisite localization; at \(U N_b=0.5\), \(f_s\) drops to zero only for \(\lambda\gtrsim2.2\) [1506.07692]. In the associated classical dynamical-map formulation, the largest Lyapunov exponent equals the inverse localization length, and in the non-interacting Aubry-André model \(\xi^{-1}=\ln(\lambda/2)\) [1506.07692]. Bogoliubov-de Gennes analysis further shows that the quasiparticle spectrum develops gaps and a devil’s-staircase structure, the non-condensate density itself localizes in a multisite pattern, and the condensate fraction drops sharply around \(\lambda\simeq2\) [1506.07692].

A related few-body treatment of a finite one-dimensional quasiperiodic lattice uses
\[
H=\sum_j\left[-J(b_{j+1}^\dagger b_j+b_j^\dagger b_{j+1})
+V_j n_j+\frac{U}{2}n_j(n_j-1)\right],
\qquad
V_j=\Delta\cos(2\pi\alpha j+\phi),
\]
with \(V_p=10E_r\), \(J\sim0.03E_r\), \(\alpha\approx1.19721\), and weak interaction \(U\ll J\) [2502.14440]. Localization is characterized through the reduced one-body density matrix \(\rho_{ij}^{(1)}=\langle b_i^\dagger b_j\rangle\), its natural orbitals, the fragmentation parameter \(f=1-n_1/N\), the order parameter \(\Delta_{\rm coh}=\sum_\alpha (n_\alpha/N)^2\), the participation ratio \(PR=(\sum_i|\psi(i)|^4)^{-1}\), and the inverse localization length extracted from the asymptotic decay of \(\rho_{ij}^{(1)}\) [2502.14440]. In the \(U\to0\) limit, the sharp Aubry-André transition occurs at \(\Delta/J=2\), with
\[
\xi(\Delta)=\left[\ln\!\left(\frac{\Delta}{2J}\right)\right]^{-1},
\qquad
\xi(\Delta)\propto (\Delta/J-2)^{-1}
\quad \text{for } \Delta/J-2\ll1^+.
\]
Weak repulsion slightly increases the localization length through \(\delta\xi\sim +c(U/J)\xi^2\), \(c>0\), and for \(N=5\) the leading occupation \(n_1/N\) remains above \(95\%\) up to \(\Delta/J\approx1.8\), while for \(\Delta/J\gtrsim2.2\) the localization length agrees with \([\ln(\Delta/2J)]^{-1}\) within a few percent [2502.14440].

In two-dimensional optical lattices with uncorrelated on-site disorder \(\epsilon_i\in[-\Delta,+\Delta]\), mean-field computational experiments quantify localization through the inverse participation ratio,
\[
\mathrm{IPR}(\Delta)=
\frac{\sum_{i=1}^{\Omega}|\psi_i|^4}
{\left(\sum_{i=1}^{\Omega}|\psi_i|^2\right)^2},
\]
the density variance \(\mathrm{Var}[\rho]=\langle \rho_i^2\rangle-\langle \rho_i\rangle^2\), and the peak-fraction \(p_f(\Delta)\) defined by a density threshold \(\eta\) [1610.05398]. For square and triangular lattices with \(V_0\approx12E_R\) and \(U\approx0.01E_R\), weak disorder \(\Delta\lesssim0.3E_R\) leaves the condensate nearly uniform, while above \(\Delta_c\approx0.3E_R\) a smooth fragmentation begins: IPR grows exponentially in \(\Delta\), \(p_f\) decreases, and \(\mathrm{Var}[\rho]\) rises. The interval \(\Delta\sim0.3\)–\(0.6E_R\) is described as “weak-localisation,” while for \(\Delta\gtrsim0.6E_R\) the condensate has effectively collapsed onto a few isolated islands [1610.05398]. Honeycomb lattices behave differently because the effective tunnelling is larger at the same lattice depth: for \(V_0=12E_R\), no significant fragmentation is found up to \(\Delta\approx1E_R\), whereas increasing \(V_0\to24E_R\) produces a crossover at \(\Delta_c\approx0.4E_R\) [1610.05398]. The disorder-averaged momentum-space energy spectrum simultaneously evolves from sharp Bragg peaks at \(\Delta=0\) to a dense, almost flat spectrum at strong disorder [1610.05398].

| System | Control parameter | Localization signature |
|---|---|---|
| 1D Aubry-André condensate | \(\lambda_c=2\) | IPR rise, \(f_s\to0\), multisite localization [1506.07692] |
| 1D finite few-boson quasiperiodic lattice | \(\Delta/J=2\) | \(PR\) reduction, exponential decay of \(\rho_{ij}^{(1)}\), \(\xi=[\ln(\Delta/2J)]^{-1}\) [2502.14440] |
| 2D disordered optical lattices | \(\Delta_c\approx0.3E_R\) for square/triangular | IPR growth, \(\mathrm{Var}[\rho]\) increase, \(p_f\) decrease [1610.05398] |

These lattice studies do not use the same transport observable as the quasi-two-dimensional film work. Instead, they characterize weak localization through condensate fragmentation, loss of superfluid stiffness, finite-size crossover scales, and the spatial structure of the dominant orbital. A plausible implication is that, in bosonic systems, the operational meaning of weak localization is platform-specific even when the underlying interference physics is common.

## 6. Superconductor-insulator transitions, percolation, and interpretive issues

The quasi-two-dimensional film results place weak localization of bosons directly in the broader problem of magnetic-field-tuned superconductor-insulator transitions. In the intermediate regime, finite-size scaling with \(R_c=3384\,\Omega\), \(B_c=1.58\,\mathrm T\), and \(\nu z\approx0.4\) holds only for \(2\,\mathrm K<T<4\,\mathrm K\); below \(2\,\mathrm K\), the strongly nonmonotonic \(R(T)\) indicates a mixed regime rather than a direct continuous quantum critical point [1002.1720]. The interpretation given is that, at finite temperature, the film separates into superconducting regions embedded in an insulating matrix, forming a percolating network. In the \(T\to0\) limit, the system ends in one of two pure phases—superconducting or bosonic-localized insulator—depending on field, which is consistent with a percolative quantum phase transition rather than a direct \((2+1)\)-dimensional XY critical point [1002.1720].

The small effective exponent \(\nu z\approx0.4\) was explicitly compared with Monte Carlo studies of unfrustrated Josephson-junction arrays at rational frustration \(f=1/4\), where a first-order-like jump to a phase mixture is predicted for some parameters [1002.1720]. This comparison does not equate the systems, but it places the observed indirect transition route within a broader class of bosonic models where coexistence and percolation intervene between idealized critical regimes.

Several recurrent misconceptions are addressed by the cited literature. First, weak localization of bosons is not equivalent to any insulating bosonic phase: in the film work, the absence of a large magnetoresistance peak was used to distinguish weakly localized bosons from simple boson-pair hopping or a true Bose glass [1002.1720]. Second, interaction effects do not simply suppress interference monotonically: in three-dimensional disorder, inelastic collisions reduce coherence globally yet produce a residual enhancement of coherent backscattering in a narrow spectral band around the initial energy [1307.4940]. Third, apparent weak antilocalization in nonlinear bosonic matter waves need not signal a universal bosonic reversal of the interference correction; in the semiclassical billiard problem it is attributed to non-universal short self-retracing paths that modify the zero-field reflection background [1112.5603].

Taken together, these results define weak localization of bosons as a coherence-sensitive interference phenomenon whose observables include \(R\sim\ln T\), coherent backscattering lineshapes, condensate fragmentation, and decay of one-body correlations. The most direct condensed-matter signature remains the logarithmic temperature dependence of the **resistance** in the high-field insulating regime of ultrathin films, while cold-atom realizations provide microscopic control over collisions, spectral redistribution, geometry, and quasiperiodicity [1002.1720] [1307.4940] [1112.5603] [1506.07692] [1610.05398] [2502.14440].

Source: https://www.emergentmind.com/topics/weak-localization-of-bosons