Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2205.06442 · Physical review. B./Physical review. B

Migdal-Eliashberg theory as a classical spin chain

Emil A. Yuzbashyan, B. L. Altshuler

Published Jul 18, 2022• 21 citations
arXivHTMLPDF

Abstract

We formulate the Migdal-Eliashberg theory of electron-phonon interactions in terms of classical spins by mapping the free energy to a Heisenberg spin chain in a Zeeman magnetic field. Spin components are energy-integrated normal and anomalous Green's functions and sites of the chain are fermionic Matsubara frequencies. The Zeeman field grows linearly with the spin coordinate and competes with ferromagnetic spin-spin interaction that falls off as the square of the inverse distance. The spin-chain representation makes a range of previously unknown properties plain to see. In particular, infinitely many new solutions of the Eliashberg equations both in the normal and superconducting states emerge at strong coupling. These saddle points of the free-energy functional correspond to spin flips. We argue that they are also fixed points of kinetic equations and play an essential role in far from equilibrium dynamics of strongly coupled superconductors. Up to an overall phase, the frequency-dependent gap function that minimizes the free energy must be non-negative. There are strong parallels between our Eliashberg spins and Anderson pseudospins, though the two sets of spins never coincide.

DOI
10.1103/physrevb.106.014512

Related papers

Monte Carlo solver and renormalization of Migdal-Eliashberg spin chain

Yang-Zhi Chou, Zhentao Wang, S. Das Sarma

Motivated by the recently developed classical spin model for the Migdal-Eliashberg theory, we develop numerical and analytical methods based on this spin-chain representation and apply these methods to the Bogoliubov-Tomachov-Morel-Anderson pairing potential, which incorporates the phonon-mediated attraction and Coulomb repulsion. We show that the Monte Carlo method with heat bath updates can efficiently obtain the gap functions even for the situations challenging for the iterative solvers, suggesting an unprecedented robust approach for solving the full nonlinear Migdal-Eliashberg theory. Moreover, we derive the renormalization of all the couplings by tracing out the high-frequency spins in the partition function. The derived analytical renormalization equations produce the well-known ${\ensuremath{\mu}}^{*}$ effect for the Bogoliubov-Tomachov-Morel-Anderson pairing potential and can be generalized to other superconductivity problems. We further point out that several interesting features (e.g., sign changing in the frequency-dependent gap function) can be intuitively understood using the classical spin-chain representation for Migdal-Eliasherg theory. Our results show the advantage of using the spin-chain representation for solving Migdal-Eliashberg theory and provide new ways for tackling general superconductivity problems.

arXiv:2311.12094•Physical review. B./Physical review. B•Feb 20, 2024•0 citations

Exact single and two particle non-equilibrium Green's functions for the time-dependent Holstein-Hubbard model in the atomic limit

R. D. Nesselrodt, James K. Freericks

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales -- for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections -- but strongly correlated electrons require a more complex description and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We derive nonequilibrium spectral moment sum rules for the first four moments of the single-particle Green's function. We then analytically confirm that this Green's function exactly satisfies the zeroth, first, and second moment sum rules, and provide numerical evidence that the third spectral moment sum rule also holds. Finally we calculate the two-particle Green's function in the atomic limit. These two functions are required inputs to create a self-consistent strong-coupling expansion for the lattice Green's functions of the Hubbard-Holstein model.

•Mar 26, 2021•0 citations

Migdal’s theorem and the pseudogap

Philippe H. Monthoux

We study a model of quasiparticles on a two-dimensional square lattice coupled to Gaussian-distributed dynamical molecular fields. We consider two types of such fields, a vector molecular field that couples to the quasiparticle spin density and a scalar field coupled to the quasiparticle number density. The model describes quasiparticles coupled to spin or charge fluctuations, and is solved by a Monte Carlo sampling of the molecular-field distributions. When the molecular-field correlations are sufficiently weak, the corrections to the self-consistent Eliashberg theory do not bring about qualitative changes in the quasiparticle spectrum. But for a range of model parameters near the magnetic boundary, we find that Migdal's theorem does not apply and the quasiparticle spectrum is qualitatively different from its Eliashberg approximation. In the range of model parameters studied, we find that the transverse spin-fluctuation modes play a key role. While a pseudogap opens when quasiparticles are coupled to antiferromagnetic fluctuations, such a pseudogap is not observed in the corresponding charge-fluctuation case for the range of parameters studied, where vertex corrections are found to effectively reduce the strength of the interaction. This suggests that one has to be closer to the border of long-range order to observe pseudogap effects in the charge-fluctuation case than for a spin-fluctuation-induced interaction under otherwise similar conditions. An important feature of the magnetic pseudogap found in the present calculations is that it is strongly anisotropic. It vanishes along the diagonal of the Brillouin zone and is large near the zone boundary. In the case of ferromagnetic fluctuations, we also find a range of model parameters with qualitative changes in the quasiparticle spectral function not captured by the one-loop approximation, that is, the quasiparticle peak splits into two. We find that one needs to be closer to the magnetic boundary to observe the pseudogap effects in the nearly ferromagnetic case relative to the nearly antiferromagnetic one, under otherwise similar conditions. We provide intuitive arguments to explain the physical origin of the breakdown of Migdal's theorem.

arXiv:cond-mat/0301228•Physical review. B, Condensed matter•Aug 7, 2003•18 citations

Exact solution of two simple non-equilibrium electron-phonon and electron-electron coupled systems: The atomic limit of the Holstein-Hubbard model and the generalized Hatsugai-Komoto model

R. D. Nesselrodt, J. K. Freericks

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales---for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We then examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.

arXiv:2103.14771•Physical review. B./Physical review. B•Oct 4, 2021•7 citations

Gor'kov-Hedin-Baym equations for quantum many-body systems with spin-dependent interactions

Christopher Lane

Driven by the need to understand and determine the presence of non-trivial superconductivity in real candidate materials, we present a generalized set of self-consistent Gor'kov-Hedin-Baym equations with spin dependent electron-electron and electron-phonon interactions. This extends Hedin's original equations to treat quantum many-body systems where electronic and lattice correlations along with relativistic effects coexist on the same footing with superconductivity. The leading order self-energies yields a generalization of the Migdal-Eliashberg theory and by iterating this set of equations generalized ladder vertex corrections naturally emerge.

arXiv:2506.07302•Physical review. B./Physical review. B•Mar 20, 2026•1 citation

Breakdown of Migdal-Eliashberg theory in the infinite-dimensional Holstein model

J. P. Hague, N. d’Ambrumenil

We investigate the applicability of Migdal--Eliashberg (ME) theory by revisiting Migdal's analysis within the dynamical mean-field theory framework. First, we compute spectral functions, the quasi-particle weight, the self energy, renormalised phonon frequency and resistivity curves of the half-filled Holstein model. We demonstrate how ME theory has a phase-transition-like instability at intermediate coupling, and how the Engelsberg--Schrieffer (ES) picture is complicated by low-energy excitations from higher order diagrams (demonstrating that ES theory is a very weak coupling approach). Through consideration of the lowest-order vertex correction, we analyse the applicability of ME theory close to this transition. We find a breakdown of the theory in the intermediate coupling adiabatic limit due to a divergence in the vertex function. The region of applicability is mapped out, and it is found that ME theory is only reliable in the weak coupling adiabatic limit, raising questions about the accuracy of recent analyses of cuprate superconductors which do not include vertex corrections.

arXiv:cond-mat/0106355•Jun 18, 2001•0 citations
We formulate the Migdal-Eliashberg theory of electron-phonon interactions in terms of classical spins by mapping the free energy to a Heisenberg spin chain in a Zeeman magnetic field. Spin components are energy-integrated normal and anomalous Green's functions and sites of the chain are fermionic Matsubara frequencies. The Zeeman field grows linearly with the spin coordinate and competes with ferromagnetic spin-spin interaction that falls off as the square of the inverse distance. The spin-chain representation makes a range of previously unknown properties plain to see. In particular, infinitely many new solutions of the Eliashberg equations both in the normal and superconducting states emerge at strong coupling. These saddle points of the free-energy functional correspond to spin flips. We argue that they are also fixed points of kinetic equations and play an essential role in far from equilibrium dynamics of strongly coupled superconductors. Up to an overall phase, the frequency-dependent gap function that minimizes the free energy must be non-negative. There are strong parallels between our Eliashberg spins and Anderson pseudospins, though the two sets of spins never coincide.
Motivated by the recently developed classical spin model for the Migdal-Eliashberg theory, we develop numerical and analytical methods based on this spin-chain representation and apply these methods to the Bogoliubov-Tomachov-Morel-Anderson pairing potential, which incorporates the phonon-mediated attraction and Coulomb repulsion. We show that the Monte Carlo method with heat bath updates can efficiently obtain the gap functions even for the situations challenging for the iterative solvers, suggesting an unprecedented robust approach for solving the full nonlinear Migdal-Eliashberg theory. Moreover, we derive the renormalization of all the couplings by tracing out the high-frequency spins in the partition function. The derived analytical renormalization equations produce the well-known ${\ensuremath{\mu}}^{*}$ effect for the Bogoliubov-Tomachov-Morel-Anderson pairing potential and can be generalized to other superconductivity problems. We further point out that several interesting features (e.g., sign changing in the frequency-dependent gap function) can be intuitively understood using the classical spin-chain representation for Migdal-Eliasherg theory. Our results show the advantage of using the spin-chain representation for solving Migdal-Eliashberg theory and provide new ways for tackling general superconductivity problems.
One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales -- for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections -- but strongly correlated electrons require a more complex description and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We derive nonequilibrium spectral moment sum rules for the first four moments of the single-particle Green's function. We then analytically confirm that this Green's function exactly satisfies the zeroth, first, and second moment sum rules, and provide numerical evidence that the third spectral moment sum rule also holds. Finally we calculate the two-particle Green's function in the atomic limit. These two functions are required inputs to create a self-consistent strong-coupling expansion for the lattice Green's functions of the Hubbard-Holstein model.
We study a model of quasiparticles on a two-dimensional square lattice coupled to Gaussian-distributed dynamical molecular fields. We consider two types of such fields, a vector molecular field that couples to the quasiparticle spin density and a scalar field coupled to the quasiparticle number density. The model describes quasiparticles coupled to spin or charge fluctuations, and is solved by a Monte Carlo sampling of the molecular-field distributions. When the molecular-field correlations are sufficiently weak, the corrections to the self-consistent Eliashberg theory do not bring about qualitative changes in the quasiparticle spectrum. But for a range of model parameters near the magnetic boundary, we find that Migdal's theorem does not apply and the quasiparticle spectrum is qualitatively different from its Eliashberg approximation. In the range of model parameters studied, we find that the transverse spin-fluctuation modes play a key role. While a pseudogap opens when quasiparticles are coupled to antiferromagnetic fluctuations, such a pseudogap is not observed in the corresponding charge-fluctuation case for the range of parameters studied, where vertex corrections are found to effectively reduce the strength of the interaction. This suggests that one has to be closer to the border of long-range order to observe pseudogap effects in the charge-fluctuation case than for a spin-fluctuation-induced interaction under otherwise similar conditions. An important feature of the magnetic pseudogap found in the present calculations is that it is strongly anisotropic. It vanishes along the diagonal of the Brillouin zone and is large near the zone boundary. In the case of ferromagnetic fluctuations, we also find a range of model parameters with qualitative changes in the quasiparticle spectral function not captured by the one-loop approximation, that is, the quasiparticle peak splits into two. We find that one needs to be closer to the magnetic boundary to observe the pseudogap effects in the nearly ferromagnetic case relative to the nearly antiferromagnetic one, under otherwise similar conditions. We provide intuitive arguments to explain the physical origin of the breakdown of Migdal's theorem.
One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales---for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. We then examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.
Driven by the need to understand and determine the presence of non-trivial superconductivity in real candidate materials, we present a generalized set of self-consistent Gor'kov-Hedin-Baym equations with spin dependent electron-electron and electron-phonon interactions. This extends Hedin's original equations to treat quantum many-body systems where electronic and lattice correlations along with relativistic effects coexist on the same footing with superconductivity. The leading order self-energies yields a generalization of the Migdal-Eliashberg theory and by iterating this set of equations generalized ladder vertex corrections naturally emerge.
We investigate the applicability of Migdal--Eliashberg (ME) theory by revisiting Migdal's analysis within the dynamical mean-field theory framework. First, we compute spectral functions, the quasi-particle weight, the self energy, renormalised phonon frequency and resistivity curves of the half-filled Holstein model. We demonstrate how ME theory has a phase-transition-like instability at intermediate coupling, and how the Engelsberg--Schrieffer (ES) picture is complicated by low-energy excitations from higher order diagrams (demonstrating that ES theory is a very weak coupling approach). Through consideration of the lowest-order vertex correction, we analyse the applicability of ME theory close to this transition. We find a breakdown of the theory in the intermediate coupling adiabatic limit due to a divergence in the vertex function. The region of applicability is mapped out, and it is found that ME theory is only reliable in the weak coupling adiabatic limit, raising questions about the accuracy of recent analyses of cuprate superconductors which do not include vertex corrections.