Bond Bipolaron Model

• Bond bipolaron model is an electron–phonon framework where lattice distortions modulate hopping amplitudes rather than local charge, leading to pair formation.
• It employs off-diagonal (bond-centered) coupling, as seen in bond-SSH and Peierls models, to produce compact bipolarons with low effective mass and enhanced mobility.
• The approach integrates sign-problem-free Monte Carlo simulations and semiclassical analyses to elucidate binding energies, effective masses, and superconducting transition temperatures.

to=functions.list_tools 经彩票json 凤凰大参考 เงินฟรี 彩神争霸有 to=functions.list_tools 人人中彩票json to=functions.get_tool_schema 天天中彩票开奖json to=functions.arxiv_search 天天中彩票开奖json ￣第四色നം 聚利squery":"bond bipolaron SSH superconductivity dispersion Monte Carlo","max_results":10} Bond bipolaron model usually denotes a class of electron–phonon Hamiltonians in which lattice distortions modulate electronic hopping on bonds rather than couple primarily to on-site charge density. In the bond Su–Schrieffer–Heeger (SSH) or bond Peierls realization, oscillators live on nearest-neighbor links, or equivalently the interaction depends on differences of neighboring displacements, so pairing is generated by phonon-assisted modulation of kinetic energy. Across square, cubic, triangular, and one-dimensional lattices, this mechanism has been reported to produce compact bipolarons with comparatively low effective mass, a property central to dilute bipolaronic superconductivity scenarios (Zhang et al., 2021, Zhang et al., 2024, Zhang, 10 Jul 2025). The label is not fully unique, however: related literature extends it to intersite or bond-centered pairs in longer-range or materials-specific models, whereas continuum translation-invariant Fröhlich bipolarons belong to a distinct theoretical framework (0904.3678, Lakhno, 2013).

1. Conceptual scope and defining features

The defining feature of a bond bipolaron model is off-diagonal electron–phonon coupling: the phonon field modifies a hopping matrix element, not an on-site potential. In the square-lattice bond-SSH model, for example, the electronic term is nearest-neighbor hopping plus on-site Hubbard repulsion, while the interaction couples the bond operator $(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})$ to a bond phonon coordinate $(b_b^\dagger+b_b)$ (Zhang et al., 2024). In one-dimensional Peierls formulations, the same structure appears as a coupling to displacement differences $(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$ (Nocera et al., 2020). In both cases, the phonon field reshapes the kinetic network itself.

This distinguishes bond bipolarons from Holstein-type bipolarons, where phonons couple to local density and strong coupling usually produces heavy, site-centered objects. The bond-SSH literature instead emphasizes a kinetic-energy/lattice-dynamical pairing mechanism that can keep the pair both compact and mobile (Kim et al., 2023). In the screened Hubbard–Fröhlich chain, the same broad bond-centered interpretation appears in the intersite singlet $S1$ regime, where the relevant bound state is bond-like rather than on-site (0904.3678).

Formulation	Defining interaction	Characteristic result in cited work
Bond SSH / bond Peierls lattice model	Phonons modulate nearest-neighbor hopping on bonds	Small, relatively light bipolarons and high dilute-limit $T_c$ estimates (Zhang et al., 2024)
1D Peierls liquid	Bond coupling plus optical phonons, optionally dispersive	Stable dilute bipolaron liquid up to quarter filling (Nocera et al., 2020)
Hard-core SSH model	Bond coupling generates repulsion plus pair hopping	Repulsively bound excited bipolaron above the continuum (Sous et al., 2017)
Translation-invariant Fröhlich bipolaron	Continuum Fröhlich coupling, delocalized ground state	Explicitly not a lattice bond model (Lakhno, 2013)

The term therefore has a narrow usage and a broader usage. In the narrow sense it refers to bond-centered SSH/Peierls lattice Hamiltonians. In the broader sense it may include intersite or bond-disproportionated paired states in organic semiconductors, bismuthates, or screened Fröhlich models, provided the pairing is not reducible to a purely local negative-$U$ picture (Perroni et al., 2012, Jiang et al., 2020).

2. Hamiltonians and parameterizations

A standard bond-SSH Hamiltonian on a square or cubic lattice is written as

$H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$

with

$$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$

$$H_{\rm ph}=\omega_0\sum_b\left(b_b^\dagger b_b+\frac12\right)-t_{\rm ph}\sum_{\langle b,b'\rangle}(b_b^\dagger b_{b'}+{\rm h.c.}),$$

$$H_{\rm e-ph}=g\sum_{b=\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})(b_b^\dagger+b_b).$$

Here $(b_b^\dagger+b_b)$0 is bare hopping, $(b_b^\dagger+b_b)$1 the on-site Hubbard repulsion, $(b_b^\dagger+b_b)$2 the bond electron–phonon coupling, and $(b_b^\dagger+b_b)$3 the hopping of optical phonons between parallel nearest-neighbor bonds (Zhang et al., 2024). The dispersive phonon band is characterized by

$(b_b^\dagger+b_b)$4

and a common dimensionless coupling is

$(b_b^\dagger+b_b)$5

with $(b_b^\dagger+b_b)$6 for $(b_b^\dagger+b_b)$7 and $(b_b^\dagger+b_b)$8 for $(b_b^\dagger+b_b)$9 (Zhang et al., 2024).

A widely used Einstein-phonon version omits $(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$0 and keeps only

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$1

On the triangular lattice this yields the same formal bond-centered structure, with $(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$2 and hence

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$3

The triangular geometry enters not through a new interaction term but through the coordination number $(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$4 and the denser bond network (Zhang, 10 Jul 2025).

In one dimension, the Peierls or SSH form is often written instead as

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$5

or, equivalently, in terms of displacement differences

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$6

A common 1D convention is

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$7

These notational differences reflect geometry and phonon placement, not a change in the core principle: the interaction remains off-diagonal in the electronic basis (Nocera et al., 2020, Marijanović et al., 18 Feb 2025).

With phonon dispersion, the phonon propagator becomes spatially nonlocal. In the dispersionless case,

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$8

whereas for dispersive optical phonons,

$(b_i^\dagger+b_i-b_{i+1}^\dagger-b_{i+1})$9

with $S1$0 the modified Bessel function. This nonlocality is the formal locus at which phonon bandwidth enters exact Monte Carlo sampling (Zhang et al., 2024).

3. Binding mechanism, quasiparticle structure, and diagnostics

The central physical mechanism is not a local lattice trap but bond-assisted coherent motion. In the bond-SSH picture, two electrons lower their energy by dressing the hopping network, so the effective attraction is generated primarily by phonon-assisted modulation of hopping rather than by density coupling (Zhang et al., 2024). The semiclassical square-lattice analysis makes this point explicitly: in the strong-coupling limit the bipolaron does not sit in a deep isolated minimum but in a continuous degenerate manifold of bond distortions satisfying

$S1$1

which allows the pair to slide with only a small tunneling barrier (Kim et al., 2023).

The standard diagnostics are the bipolaron binding energy, effective mass, and size. In much of the lattice literature,

$S1$2

$S1$3

$S1$4

A positive $S1$5 indicates a bound pair, a small $S1$6 favors mobility, and a small $S1$7 raises the density at which non-overlapping bosonic pairs can still be realized (Zhang et al., 2024). Some authors instead define the binding energy with the opposite sign, for example

$S1$8

so the sign convention must be checked paper by paper (Marijanović et al., 18 Feb 2025).

A recurring result is that bond coupling alters the usual mass–binding relation. In dispersionless and dispersive bond-SSH models, stronger coupling can make the pair more compact without automatically producing the extreme mass growth familiar from Holstein physics (Zhang et al., 2024). In the sign-free Monte Carlo study, the superconducting route is described as lying between the “Scylla of large size of moderately light bipolarons and Charybdis of large mass of compact bipolarons,” which states the optimization problem directly (Zhang et al., 2021).

Not all bond-bipolaron branches are attraction-driven. In the hard-core 1D SSH model, virtual phonons generate an effective nearest-neighbor repulsion

$S1$9

together with a next-nearest-neighbor hopping

$T_c$0

The result is a repulsively bound nearest-neighbor pair stabilized by hard-core statistics, pair hopping, and lattice-band kinematics; it lies above the two-polaron continuum but below higher phonon sectors, and can have a negative effective mass near $T_c$1 (Sous et al., 2017).

The regime dependence is substantial. A variational study of the 1D SSH model found that bipolarons remain light in the anti-adiabatic regime but become exponentially heavy in the adiabatic regime (Marijanović et al., 18 Feb 2025). By contrast, the 2D semiclassical bond-SSH analysis in the adiabatic limit found only weak mass enhancement and a much larger upper bound on $T_c$2 than in the Holstein model (Kim et al., 2023). This suggests strong sensitivity to dimensionality, phonon implementation, and approximation scheme.

4. Computational and analytical frameworks

The modern bond-bipolaron literature is methodologically notable because several of its core results were obtained with sign-problem-free or numerically exact techniques. A foundational development was the sign-free Monte Carlo formulation combining a lattice path integral for the particle sector with either real-space diagrammatics or a Fock-path-integral treatment of the phonons. In the chosen basis, the elementary hopping, phonon-creation, and phonon-annihilation matrix elements are non-negative, permitting efficient sampling of the projected propagator and direct extraction of energy, effective mass, and pair size (Zhang et al., 2021).

This framework was extended to dispersive phonons in a diagrammatic Monte Carlo treatment that keeps electron worldlines on the lattice and handles the phonon sector by real-space diagrammatics. Because the dispersive propagator is sampled exactly, nonlocal optical-phonon propagation does not require uncontrolled approximations (Zhang et al., 2024). A related exactly sign-problem-free QMC approach was used to incorporate long-range Coulomb repulsion in the two-dimensional square lattice bond-SSH model (Zhang, 2024).

On the triangular lattice, diagrammatic Monte Carlo with a lattice path-integral formulation for both electron and phonon degrees of freedom was used to extract the long-imaginary-time propagator

$T_c$3

and the diffusive estimator

$T_c$4

which gives the effective mass from the center-of-mass broadening (Zhang, 10 Jul 2025).

Analytical control has come mainly from semiclassical and effective-model reductions. In the adiabatic regime, instanton theory yields the effective tunneling amplitude

$T_c$5

with

$T_c$6

For the square-lattice bond-SSH model, the fitted strong-coupling asymptotics

$T_c$7

imply a much smaller instanton action than in the Holstein case (Kim et al., 2023). In the anti-adiabatic hard-core SSH limit, the bound-state spectrum can instead be obtained from the pole of the two-particle propagator via an exact BBGKY equation-of-motion treatment (Sous et al., 2017).

At finite density in one dimension, DMRG established that Peierls bipolaron liquids remain stable over a broad range of coupling and filling, with phase separation appearing only in a narrow regime where the linearized model becomes unphysical unless phonon dispersion is carefully tuned (Nocera et al., 2020).

5. Bipolaronic superconductivity and transition-temperature estimates

The superconducting interpretation of the bond bipolaron model is usually dilute and bosonic: tightly bound pairs are treated as composite bosons undergoing either a Berezinskii–Kosterlitz–Thouless transition in two dimensions or Bose–Einstein condensation in three dimensions. In two dimensions,

$T_c$8

whereas in three dimensions

$T_c$9

Imposing a non-overlap condition for the pairs gives the practical upper-bound-type estimate

$U$0

so the competition is explicitly between compactness and mobility (Zhang et al., 2024).

In the square-lattice bond-SSH model with dispersive optical phonons, the central result is that phonon dispersion makes bipolarons more extended but lighter, so $U$1 retains its dome shape while the near-optimal region broadens. In two dimensions the broader dome extends toward larger $U$2; in three dimensions the peak can shift and, in deep adiabatic regimes with sufficiently large phonon bandwidth, $U$3 may continue increasing where the dispersionless model would already be suppressed by very large mass (Zhang et al., 2024). This led to the conclusion that phonon dispersion is not a limiting factor for high-temperature bipolaronic superconductivity in the bond SSH model.

Long-range Coulomb repulsion suppresses but does not necessarily eliminate the mechanism. In the 2D square-lattice study with

$U$4

the practical estimate becomes

$U$5

for long-range Coulomb repulsion. For $U$6 and $U$7, the optimal $U$8 is reduced by about $U$9 for $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$0 and by about $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$1 for $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$2 relative to $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$3, while only very strong repulsion such as $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$4 drives the peak below $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$5 and effectively removes it (Zhang, 2024).

Geometry also matters. On the triangular lattice, bond-centered coupling with $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$6 was found to sustain relatively high $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$7 over a wide range of $H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$8, with peak values around

$H=H_{\rm e}+H_{\rm ph}+H_{\rm e-ph},$9

for representative parameters and still

$$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$0

at $$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$1. A notable result is that moderate on-site repulsion,

$$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$2

enhances $$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$3 by preventing overly localized and therefore overly heavy on-site pairs (Zhang, 10 Jul 2025). This accords with the earlier sign-free Monte Carlo observation that on-site repulsion can help $$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$4-wave superconductivity in this setting, contrary to conventional intuition (Zhang et al., 2021).

6. Extensions, nonstandard variants, and conceptual boundaries

The bond-bipolaron concept extends beyond the dilute two-electron square-lattice problem. In one dimension, DMRG studies found “Peierls bipolaron liquids” that are stable against phase separation in the dilute limit and remain stable up to quarter filling, except in a narrow region where the linearized bond-coupling model becomes unphysical. Phonon dispersion can further stabilize the liquid and push the phase-separation boundary to larger coupling (Nocera et al., 2020). This places bond bipolarons within a many-body liquid framework rather than only a two-particle bound-state framework.

Materials-oriented models broaden the term further. In organic semiconductors at dielectric interfaces, interface phonons provide the dominant binding tendency while a bulk SSH-like coupling modifies the spatial structure and transport; the conclusions were described as supporting a bond-bipolaron interpretation based on extended distortions rather than a purely local trap (Perroni et al., 2012). In hole-doped bismuthates, a semiclassical three-orbital model with off-diagonal coupling to Bi–O bond distortions led to a picture in which the bond-disproportionated insulator is a frozen bipolaron crystal that melts into a dynamic bipolaron/polaron liquid under heating or doping (Jiang et al., 2020).

The concept also includes qualitatively nonstandard bound states. The 1D hard-core SSH model supports a repulsively bound excited bipolaron above the two-polaron continuum, protected by energy and momentum conservation (Sous et al., 2017). This shows that a bond bipolaron need not be a low-energy attractive bound state.

Two conceptual boundaries are especially important. First, “bond bipolaron” is not synonymous with all light intersite bipolarons. The screened Hubbard–Fröhlich chain supports robust intersite singlets under strong Coulomb repulsion and is often interpreted as bond-like, but its microscopic interaction is long-range density coupling rather than canonical bond-SSH coupling (0904.3678). Second, the translation-invariant Fröhlich bipolaron of Tulub–Lakhno theory is explicitly a continuum, delocalized, translation-invariant object and not a lattice bond model. Its variational estimate

$$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$5

and its stability thresholds in $$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$6 belong to a different framework (Lakhno, 2013).

A persistent theme across the literature is that “light and compact” is a regime statement, not a theorem. Some bond-SSH formulations yield weak mass enhancement and unusually large superconducting scales (Kim et al., 2023), whereas others find exponential mass growth in the adiabatic regime (Marijanović et al., 18 Feb 2025). What is robust is narrower: bond-centered coupling systematically changes the tradeoff between binding, size, and mobility relative to Holstein physics, and it opens parameter regions in which mobile bipolarons, stable bipolaron liquids, or broadened high-$$H_{\rm e}=-t\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{j,\sigma}c_{i,\sigma}+{\rm h.c.})+U\sum_i n_{i,\uparrow}n_{i,\downarrow},$$7 domes become possible (Zhang et al., 2024, Nocera et al., 2020).