Strong light-matter and electron-phonon coupling generate ground states dressed by virtual bosonic excitations, making bare-state truncations and perturbative treatments unreliable in the ultrastrong-coupling regime. We introduce a nonperturbative variational ground-state framework based on a state-dependent polaron transformation, combined with a product-state ansatz and a second-order perturbative correction for residual matter-boson entanglement. We show that the optimized transformed frame becomes asymptotically decoupled at infinite coupling, because the leading linear coupling is canceled while off-diagonal matter transitions are suppressed by displaced-oscillator overlaps. The approach is asymptotically correct in both weak- and strong-coupling limits and remains accurate in the intermediate regime, where fixed polaron transformations are least reliable. Dicke-model benchmarks reproduce ground-state energies, fidelities, and the superradiant transition, with second-order energy errors below 0.2%. Holstein-model benchmarks yield errors below 0.5% and clarify how translational symmetry affects wave-function quality. This dressed-basis framework enables nonperturbative modeling of strongly coupled light-matter and electron-phonon systems.
Strong light-matter and electron-phonon coupling generate ground states dressed by virtual bosonic excitations, making bare-state truncations and perturbative treatments unreliable in the ultrastrong-coupling regime. We introduce a nonperturbative variational ground-state framework based on a state-dependent polaron transformation, combined with a product-state ansatz and a second-order perturbative correction for residual matter-boson entanglement. We show that the optimized transformed frame becomes asymptotically decoupled at infinite coupling, because the leading linear coupling is canceled while off-diagonal matter transitions are suppressed by displaced-oscillator overlaps. The approach is asymptotically correct in both weak- and strong-coupling limits and remains accurate in the intermediate regime, where fixed polaron transformations are least reliable. Dicke-model benchmarks reproduce ground-state energies, fidelities, and the superradiant transition, with second-order energy errors below 0.2%. Holstein-model benchmarks yield errors below 0.5% and clarify how translational symmetry affects wave-function quality. This dressed-basis framework enables nonperturbative modeling of strongly coupled light-matter and electron-phonon systems.
Strong light-matter and electron-phonon coupling generate ground states dressed by virtual bosonic excitations, making bare-state truncations and perturbative treatments unreliable in the ultrastrong-coupling regime. We introduce a nonperturbative variational ground-state framework based on a state-dependent polaron transformation, combined with a product-state ansatz and a second-order perturbative correction for residual matter-boson entanglement. We show that the optimized transformed frame becomes asymptotically decoupled at infinite coupling, because the leading linear coupling is canceled while off-diagonal matter transitions are suppressed by displaced-oscillator overlaps. The approach is asymptotically correct in both weak- and strong-coupling limits and remains accurate in the intermediate regime, where fixed polaron transformations are least reliable. Dicke-model benchmarks reproduce ground-state energies, fidelities, and the superradiant transition, with second-order energy errors below 0.2%. Holstein-model benchmarks yield errors below 0.5% and clarify how translational symmetry affects wave function quality. This dressed-basis framework enables nonperturbative modeling of strongly coupled light-matter and electron-phonon systems.
The Holstein model often serves as an archetype for electron-phonon interactions and polaron formation in solids. However, precise descriptions of the Holstein polaron are difficult when the phonon frequency is small and the electron-phonon coupling is strong, due to the presence of many phonons in the ground-state. We present a semi-analytical approximation that consists of a variational ansatz with clouds of phonons surrounding the electron in the form of coherent states. This becomes particularly simple and exact in the Lang-Firsov limit. We determine the domain of validity away from this limit, and further explore the improvement achieved with a removal of the requirement that the phonon clouds form coherent states. Both approximations work extremely well for the ground-state energy at strong coupling, and both work surprisingly well also at weak-coupling. For the effective mass, however, both approximations have difficulty, particularly in the crossover region. Moreover, the coherent-state ansatz (CSA) retains an unphysical jump in the crossover region. Nonetheless, the CSA provides a simple and intuitive picture of the polaron ground-state wavefunction, and in addition predicts accurate values for the ground-state energy.
Recent development of theoretical method for dissipative quantum transport has achieved notable progresses in the weak or strong electron-phonon coupling regime. However, a generalized theory for dissipative quantum transport in a wide parameter space had not been established. In this work, a variational polaron theory for dissipative quantum transport in a wide range of electron-phonon coupling is developed. The optimal polaron transformation is determined by the optimization of the Feynman-Bogoliubov upper bound of free energy. The free energy minimization ends up with an optimal mean-field Hamiltonian and a minimal interaction Hamiltonian. Hence, second-order perturbation can be applied to the transformed system, resulting in an accurate and efficient method for the treatment of dissipative quantum transport with different electron-phonon coupling strength. Numerical benchmark calculation on a single site model coupled to one phonon mode is presented.
Strong coupling between mobile charge carriers or excitons and lattice vibrations gives rise to polarons: quasiparticles formed by charge carriers dressed by a cloud of phonon excitations. Polaron formation determines transport properties of systems ranging from transition-metal oxides and organic semiconductors to light-harvesting complexes. In excitonic systems that admit a hard-core boson description, mapping the excitations to pseudospins provides a natural framework for studying polaron physics, including the nonperturbative regime in which hopping amplitudes, phonon energies, and coupling strengths are comparable. Here, we propose a trapped-ion blueprint for realizing spin-phonon Hamiltonians with programmable spatial coupling between excitation density and phonon displacements, encompassing local Holstein and long-range Fröhlich-type interactions. To access strongly dressed quasiparticles whose bare spectral weight is suppressed, we introduce displaced Green's functions implemented through spin-dependent phonon displacements. Their spectral response reveals the polaron band spectrum, including binding energies and phonon-excitation sidebands, while spatially resolved measurements probe the structure of the phonon dressing cloud, providing a route to polaron tomography. We further complement this spectroscopy with direct and ancilla-assisted measurements of spin-phonon correlations that track the coupled dynamics of the excitation and its environment. These approaches connect the emergent quasiparticle spectrum to the microscopic buildup of its dressing cloud, enabling controlled studies of how the spatial range of coupling shapes polaron formation and mobility.
Polaron ground states and finite-temperature dynamics remain challenging to simulate because existing methods struggle to combine nonperturbative accuracy, systematic improvability, and scalability from models to materials-specific Hamiltonians. We introduce a translationally invariant variational coupled-cluster (CC) theory for polarons, termed delocalized CC (dCC), with closed-form energies at cost as low as O(N3) and no phonon-number cutoff. dCC accurately describes the ground states of the one- and two-dimensional Holstein and Su--Schrieffer--Heeger (optical and bond) models and the Fr{ö}hlich model, in close agreement with density matrix renormalization group (DMRG) and diagrammatic Monte Carlo benchmarks. A projected tangent-space response formalism built on the same ansatz yields electron-addition spectral functions and optical conductivities at zero and finite temperature. The resulting spectra agree well with DMRG, Lanczos, and neural-network quantum-state benchmarks while retaining a physically interpretable excitation hierarchy, and extend to two-dimensional lattices at finite temperature beyond the practical reach of these methods. The same framework applies directly to \textit{ab initio} electron--phonon matrix elements, yielding LiF electron- and hole-polaron binding energies that match state-of-the-art many-body calculations. These results establish dCC as a unified variational framework for polaron ground states and dynamics, from model systems to real materials.
We develop and compare several analytical approximations for the ground state and zero-temperature polaron dispersion in finite-width, nonparabolic conduction bands. The main focus of the work is an extension of the Feynman variational method to a tight-binding lattice, where the effective-mass approximation is no longer applicable. The resulting variational formulation is not restricted to a specific phonon dispersion or electron-phonon interaction and provides a uniform description across weak-, intermediate-, and strong-coupling regimes. In addition, we revisit and generalize other analytical approaches traditionally formulated for continuum polarons, including canonical transformations and self-consistent Wigner-Brillouin-type approximations. For lattice polarons, these methods exhibit qualitative features absent in the continuum case, such as a nontrivial connection between weak- and strong-coupling limits. We show that an improved Wigner-Brillouin scheme yields a momentum-dependent polaron self-energy free of resonances and in good agreement with numerically exact results. All methods are applied to the Holstein model on a tight-binding lattice and are benchmarked against numerically exact calculations, including diagrammatic Monte Carlo (both our calculations and preceding works) and exact diagonalization results. Furthermore, the analytical approaches are extended to polarons with Rashba-type spin-orbit coupling, providing a stringent test of their applicability in systems with nontrivial band structure. Our results demonstrate that the modified Feynman variational method yields ground-state energies and dispersions with accuracy comparable to, and in many cases exceeding, that of other established analytical approaches. The developed framework offers a versatile and reliable analytical description of lattice polarons beyond the continuum approximation.