Polaron/QCD Correspondence

• Polaron/QCD correspondence is an analogy that maps strong-coupling dressing in condensed matter systems onto non-perturbative QCD phenomena such as confinement and meson structures.
• It employs doped antiferromagnets and Fröhlich-polaron models to simulate analogues of quark confinement, mesons, tetraquark-like states, and dressing via collective bosonic modes.
• The framework is supported by analytical methods (e.g., Born–Oppenheimer, Lee–Low–Pines) and numerical techniques (DMRG, ARPES), suggesting experimental tests in ultracold atom setups.

Polaron/QCD correspondence denotes a family of analogical constructions that map strong-coupling dressing phenomena in condensed-matter systems onto non-perturbative structures familiar from quantum chromodynamics. In current arXiv usage, the term covers at least two distinct but related programs. One treats a doped antiferromagnetic spin chain in a staggered field as a controlled analog simulator of confinement, mesons, and tetraquark-like bound states, with magnons providing a polaronic dressing bath for confined partons (Čubela et al., 2022). The other translates established Fröhlich-polaron results into low-energy hadron physics by identifying the polaron with a constituent quark, phonons with pions, and the polaron ground-state energy with the nucleon mass scale (Afonin et al., 14 Jul 2025).

1. Conceptual scope and correspondence rules

In the doped-antiferromagnet formulation, the correspondence is organized around emergent partons and strings. A mobile hole fractionalizes into a spinon and a chargon, the hole motion leaves a geometric string of misaligned spins, and a staggered field turns that string into a linear confining potential. The resulting spinon–chargon bound state is identified as a meson, magnons play the role of gluonic collective modes, and a meson bound to a magnon is interpreted as a tetra-parton state analogous to a tetraquark-like bound state (Čubela et al., 2022).

In the nucleon-focused formulation, the correspondence is stated as an explicit dictionary: photons ↔ gluons, phonons ↔ pions, electrons ↔ quarks, polaron ↔ constituent quark, electron mass in the conduction band $m$ ↔ current quark mass $m_q$, effective electron mass $m^*$ ↔ constituent quark mass $m_q^*$, and the absolute value of the electron ground-state energy $|E_0|$ in the lattice ↔ the ground-state energy of a dressed quark in the QCD vacuum, identified with the nucleon mass $M_N$. The dimensionless electron–phonon coupling $\alpha$ is mapped to an effective quark–pion coupling $\alpha_{\text{eff}}$, while the phonon energy $\hbar\omega_0$ is mapped to $m_\pi c^2$, with $m_q$0 used thereafter so that $m_q$1 (Afonin et al., 14 Jul 2025).

The common principle is strong dressing of a fermionic degree of freedom by a bosonic medium. In one case the bosonic modes are magnons in a doped antiferromagnet; in the other they are pions viewed as collective excitations of the QCD vacuum. This suggests a shared organizing language of self-energy, effective mass, binding, and confinement-like structure, while leaving the microscopic content of the two correspondences sharply distinct.

2. Microscopic realization in a doped antiferromagnetic chain

The condensed-matter realization studies a single mobile dopant in a one-dimensional antiferromagnetic spin chain subject to a staggered Zeeman field. Two microscopic Hamiltonians are used as starting points, the $m_q$2–$m_q$3 model and the Hubbard model,

$m_q$4

$m_q$5

At large $m_q$6, a Schrieffer–Wolff transformation maps the Hubbard model to the $m_q$7–$m_q$8 model with $m_q$9. The setup uses a chain of length $m^*$0, one hole so that $m^*$1, and representative numerical parameters $m^*$2 with $m^*$3 in the range $m^*$4 (Čubela et al., 2022).

The strong-coupling description employs a parton decomposition $m^*$5, together with the single-occupancy constraint $m^*$6. The hole displaces spins along its path, producing a geometric string $m^*$7 of length $m^*$8 connecting a neutral spinon and a spinless chargon. In the Ising limit, retaining $m^*$9 and $m_q^*$0 and initially setting $m_q^*$1, the staggered field generates a linear confining potential,

$m_q^*$2

The contact term encodes a short-range attraction at zero separation. In the anisotropic case with couplings $m_q^*$3, the contact term becomes $m_q^*$4 for spin $m_q^*$5, with $m_q^*$6 used throughout (Čubela et al., 2022).

The same staggered field breaks SU(2), induces Néel order along $m_q^*$7, and gaps the spin-wave sector. Linear spin-wave theory yields

$m_q^*$8

For the isotropic chain,

$m_q^*$9

A crucial caveat is that, in the weak-field limit $|E_0|$0, the linear spin-wave form becomes unreliable near low energies; the spectrum approaches the spinon continuum with lower edge $|E_0|$1 (Čubela et al., 2022).

3. Mesons, tetra-partons, and polaronic solution methods

At strong coupling $|E_0|$2, a Born–Oppenheimer separation between a fast chargon and slow spinon/magnon degrees of freedom produces a meson with rovibrational quantum numbers $|E_0|$3, where $|E_0|$4 labels radial excitations and $|E_0|$5 is a parity quantum number. The meson dispersion in a fixed band is

$|E_0|$6

with a Franck–Condon renormalized hopping determined by overlaps of chargon wavefunctions at neighboring spinon positions. In the strong-coupling continuum limit, the vibrational energies take the Airy form

$|E_0|$7

so the level spacing scales as $|E_0|$8. The odd sector $|E_0|$9 has a node at the origin and lies above the corresponding even sector, so parity in one dimension plays the role that orbital angular momentum plays in two dimensions (Čubela et al., 2022).

A magnon can bind to a meson and thereby generate a composite containing one chargon and three confined spinons in total. The effective meson–magnon problem is written as

$M_N$0

where $M_N$1 is a meson operator in a fixed $M_N$2 band and $M_N$3 is a Bogoliubov magnon. The dominant couplings scale as $M_N$4 times overlap integrals of the string wavefunction. At weak coupling, the energy shift is

$M_N$5

while in the strong-field limit $M_N$6, a Schrieffer–Wolff calculation in the $M_N$7 sector gives a meson–magnon binding energy

$M_N$8

Numerically, stable binding is observed for $M_N$9 (Čubela et al., 2022).

Two complementary analytical tools are used. In the Lee–Low–Pines frame, the meson center of mass is removed, the Hamiltonian becomes block-diagonal in total momentum $\alpha$0, and quadratic expansion in magnons leads to a multi-mode Gaussian Hamiltonian diagonalized by a generalized Bogoliubov transform. In the Chevy-type approach, a two-magnon variational sector captures interaction-induced avoided quasiparticle decay, while a one-magnon projection produces a finite matrix whose lowest eigenvalue lies below the meson–magnon continuum and corresponds to the molecular tetra-parton state. In the weak-field, large-$\alpha$1 regime, the naive bare meson band would cross the meson–magnon continuum near $\alpha$2, but magnon zero-point fluctuations reduce the renormalized tunneling and the self-energy repels the bands, leaving the quasiparticle stable for all $\alpha$3 (Čubela et al., 2022).

The spectrum is identified numerically by DMRG and time-dependent MPS. The single-hole ARPES-like spectrum $\alpha$4 shows a sharp quasiparticle band with minimum at $\alpha$5 for $\alpha$6 and $\alpha$7, matched by the semi-analytical meson dispersion. Rotational ARPES at $\alpha$8 resolves even and odd vibrational resonances and demonstrates a well-defined internal meson structure. Spin-flip ARPES in the $\alpha$9 sector exhibits a weakly dispersing band below the meson–magnon continuum, and the binding energy extracted from ground-state DMRG becomes negative for $\alpha_{\text{eff}}$0. For $\alpha_{\text{eff}}$1 and $\alpha_{\text{eff}}$2, the ARPES band remains sharp and isolated over the entire Brillouin zone, consistent with interaction-stabilized avoided decay (Čubela et al., 2022).

4. QCD interpretation of the antiferromagnetic construction

The QCD reading of the doped spin chain is explicit. The spinon and chargon are mapped to quark and antiquark constituents; the geometric string with tension $\alpha_{\text{eff}}$3 is mapped to a QCD flux tube with string tension $\alpha_{\text{eff}}$4; magnons are mapped to gluonic collective modes; the spinon–chargon bound state is mapped to a hadronic meson; and the meson–magnon molecular state is mapped to a tetraquark-like bound state (Čubela et al., 2022).

Confinement is the central structural parallel. The linear potential $\alpha_{\text{eff}}$5 in the doped antiferromagnet mirrors the QCD confining potential, and the resulting rovibrational spectrum is governed by the same string-tension logic,

$\alpha_{\text{eff}}$6

In one dimension, parity splitting acts as an analog of a centrifugal barrier. In related two-dimensional studies cited in the same work, true rotational excitations generate Regge-like trajectories $\alpha_{\text{eff}}$7, whereas in the one-dimensional chain even/odd splitting provides only the parity remnant of that structure (Čubela et al., 2022).

The analogy extends from confinement to dressing. The meson’s magnon cloud is presented as a condensed-matter analog of gluonic or hadronic dressing clouds, and the interaction-induced avoided decay of a quasiparticle band into a continuum is compared to stabilization of hadronic resonances by strong coupling to multi-meson channels. The stated limits are equally important: magnons are not gauge bosons, the microscopic symmetry is SU(2) spin rather than SU(3) color, and the underlying system is a one-dimensional lattice rather than a three-dimensional continuum. The correspondence is therefore qualitative and structural rather than literal, even though the emergent $\alpha_{\text{eff}}$8 lattice gauge structure provides a controlled analog of confinement and dressing (Čubela et al., 2022).

The same work emphasizes experimental accessibility. One-dimensional Hubbard or $\alpha_{\text{eff}}$9–$\hbar\omega_0$0 chains can be realized with ultracold fermions such as Li or K in quantum gas microscopes, staggered fields can be implemented by local addressing or mixed-dimensional setups, and Rydberg dressing can engineer strong Ising couplings. The quoted parameter window is $\hbar\omega_0$1 and $\hbar\omega_0$2, with single-hole preparation by atom removal. ARPES-like spectroscopy, parity-resolved rotational ARPES, spin-flip ARPES, and direct imaging of spin–charge correlators and geometric strings are proposed as diagnostics of meson structure, magnon dressing, and tetra-parton binding (Čubela et al., 2022).

5. Nucleon physics from Fröhlich-polaron theory

The direct QCD application begins from the Fröhlich Hamiltonian,

$\hbar\omega_0$3

with dispersionless optical phonons, a linear electron–phonon coupling, and the usual dimensionless coupling $\hbar\omega_0$4. The paper emphasizes that early pion–nucleon Hamiltonians share the same formal feature of linear coupling to bosonic creation and annihilation operators, and that pions, like phonons, are spin-0 collective modes. The hadronic interaction is written in derivative form as

$\hbar\omega_0$5

The phenomenological pion–nucleon coupling is quoted as $\hbar\omega_0$6 with $\hbar\omega_0$7, and an illustrative spin–flavor SU(4) relation is $\hbar\omega_0$8 (Afonin et al., 14 Jul 2025).

The correspondence imports weak- and strong-coupling polaron results for the ground-state energy and effective mass. The quoted asymptotic forms are

$\hbar\omega_0$9

$m_\pi c^2$0

$m_\pi c^2$1

$m_\pi c^2$2

together with the improved strong-coupling value

$m_\pi c^2$3

These relations are then reinterpreted by the identifications $m_\pi c^2$4, $m_\pi c^2$5, $m_\pi c^2$6, $m_\pi c^2$7, and $m_\pi c^2$8 (Afonin et al., 14 Jul 2025).

The QCD inputs are $m_\pi c^2$9, $m_q$00, $m_q$01, and the Gell-Mann–Oakes–Renner relation

$m_q$02

Using $m_q$03, $m_q$04, and $m_q$05, the condensate is given as $m_q$06. The constituent-quark input is

$m_q$07

From the strong-coupling scaling $m_q$08 and $m_q$09, the paper derives

$m_q$10

and, after inserting the improved strong-coupling coefficient,

$m_q$11

With the quoted inputs, this gives $m_q$12 (Afonin et al., 14 Jul 2025).

The same framework yields an effective coupling estimate,

$m_q$13

and, from the adapted scaling relation,

$m_q$14

For infrared freezing $m_q$15, this gives $m_q$16, consistent with the preceding estimate, and implies $m_q$17. The weak-coupling translation of the polaron energy produces

$m_q$18

so that $m_q$19–$m_q$20 at the nucleon scale gives $m_q$21–$m_q$22. An Appendix estimate instead gives

$m_q$23

which evaluates to approximately $m_q$24. The paper notes that this second value lies below standard phenomenology and suggests that missing strange-quark loop effects could account for part of the discrepancy. It also argues that the constituent-quark share of the nucleon mass is approximately one-third, with the remaining approximately two-thirds associated with gluons and trace anomaly (Afonin et al., 14 Jul 2025).

6. Limits, caveats, and contested points

The antiferromagnetic correspondence is explicitly constrained by dimensionality and symmetry. The model is a one-dimensional lattice system with SU(2) spin degrees of freedom, not a three-dimensional SU(3) gauge theory, and magnons are not gauge bosons. The analogy is therefore limited to an emergent description of partons, strings, bound states, and dressing clouds. In addition, the low-energy spin-wave description breaks down as $m_q$25, where the spectrum approaches a spinon continuum, so the clean separation into a meson band plus gapped magnon bath is field-dependent (Čubela et al., 2022).

The nucleon correspondence has a different set of limitations. It is intrinsically nonrelativistic, treats pions as noninteracting bosons in the spirit of the Fröhlich model, and does not include spin and relativistic effects in a full QCD sense. The identification $m_q$26 and the treatment of one “physical” constituent quark with the other two valence quarks folded into the environment are model assumptions. The mapping $m_q$27 becomes delicate in the chiral limit $m_q$28, and renormalization-scale issues enter through $m_q$29, $m_q$30, and $m_q$31 (Afonin et al., 14 Jul 2025).

A further point of sensitivity concerns numerical inputs and asymptotic coefficients. In the nucleon model, the quantitative success of the mass formula depends strongly on the improved strong-coupling coefficient $m_q$32; using the older value $m_q$33 would lower the predicted nucleon mass. The choice $m_q$34 is described as both an input and an output of the framework. Likewise, the sigma-term prediction depends strongly on whether one adopts the weak-coupling estimate $m_q$35–$m_q$36 or the Appendix ansatz giving approximately $m_q$37 (Afonin et al., 14 Jul 2025).

7. Extensions and open directions

In the doped-antiferromagnet program, the stated extension is to two-dimensional doped antiferromagnets, where Néel order is spontaneous and magnons become gapless Goldstones. The meson–magnon coupling is then stronger near $m_q$38, dressing is enhanced, the continuum is reshaped, and rotational excitations acquire true orbital angular momentum $m_q$39 and Regge-like trajectories. Tetra-parton binding can persist but must compete with the gapless magnon sector. The effective polaron framework is proposed as the analytic bridge that makes these higher-dimensional questions tractable (Čubela et al., 2022).

In the nucleon program, the proposed tests are scaling relations rather than analog quantum simulation. These include $m_q$40, the approximate relation $m_q$41 when $m_q$42, and $m_q$43 at low scales. The authors also pose open questions about the nucleon wave function in the polaron framework, the inclusion of relativity and spin, the existence of strong-interaction analogs of other polaron types, and a more rigorous derivation of the momentum dependence connecting pion–nucleon Hamiltonians to the Fröhlich kernel (Afonin et al., 14 Jul 2025).

Taken together, these lines of work use polaron physics as a technical language for non-perturbative QCD phenomena, but they do so at different levels. One provides a microscopically controlled condensed-matter analog of confinement and spectroscopy; the other proposes a phenomenological mass-and-dressing correspondence for nucleons. Their shared content is the treatment of strong dressing, emergent effective masses, bound-state formation, and boson-mediated self-energy as the central structures through which QCD-like behavior can be organized.