Migdal Approximation Research Articles

Correlated Mott insulators in a strong electric field: The effects of phonon renormalization

We characterize the response of a Mott insulating system to a static electric field in terms of its conducting and spectral properties. Dissipation is included by a coupling to fermionic baths and to either optical or acoustic phonons. This paper extends and completes the analysis made in a previous work by the authors [Mazzocchi et al., Phys. Rev. B 106, 125123 (2022)]. In the present work, phonons are included diagrammatically within the Migdal approximation by also including self-consistency from the electronic feedback. The nonequilibrium steady state is addressed by means of the dynamical mean-field theory based on the nonequilibrium Green's function approach, while the so-called auxiliary master equation approach is employed as impurity solver. With optical phonons the self-consistency suppresses the steady-state current for field strengths comparable to the band gap with respect to the non-self-consistent case. This is due to the interaction of phonons with the hot electrons of the lattice which increases their temperature, thus providing a less effective relaxation channel for the current-induced Joule heat. In addition, in the case of optical phonons, the results are essentially independent of the temperature of the fermionic baths, as the latter is sensibly smaller than their characteristic frequency. On the other hand, with acoustic phonons the steady-state current is slightly suppressed by the self-consistent treatment only at field strengths close to half of the gap and especially at very small phonon frequency. Also, in this case, the results seem to slightly depend on the temperature of the fermionic baths.

Cumulant expansion in the Holstein model: Spectral functions and mobility

We examine the range of validity of the second-order cumulant expansion (CE) for the calculation of spectral functions, quasiparticle properties, and mobility of the Holstein polaron. We devise an efficient numerical implementation that allows us to make comparisons in a broad interval of temperature, electron-phonon coupling, and phonon frequency. For a benchmark, we use the dynamical mean-field theory which gives, as we have recently shown, rather accurate spectral functions in the whole parameter space, even in low dimensions. We find that in one dimension, the CE resolves well both the quasiparticle and the first satellite peak in a regime of intermediate coupling. At high temperatures, the charge mobility assumes a power law $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}2}$ in the limit of weak coupling and $\ensuremath{\mu}\ensuremath{\propto}{T}^{\ensuremath{-}3/2}$ for stronger coupling. We find that, for stronger coupling, the CE gives slightly better results than the self-consistent Migdal approximation (SCMA), while the one-shot Migdal approximation is appropriate only for a very weak electron-phonon interaction. We also analyze the atomic limit and the spectral sum rules. We derive an analytical expression for the moments in CE and find that they are exact up to the fourth order, as opposed to the SCMA where they are exact to the third order. Finally, we analyze the results in higher dimensions.

Impact ionization processes in a photodriven Mott insulator: Influence of phononic dissipation

We study a model for photovoltaic energy collection consisting of a Mott insulating layer in presence of acoustic phonons, coupled to two wide-band fermion leads at different chemical potentials and driven into a nonequilibrium steady state by a periodic electric field. We treat electron correlations with nonequilibrium dynamical mean-field theory (DMFT) using the so-called auxiliary master equation approach as impurity solver and include dissipation by acoustic phonons via the Migdal approximation. For a small hybridization to the leads, we obtain a peak in the photocurrent as a function of the driving frequency which can be associated with impact ionization processes. For larger hybridizations the shallow peak suggests a suppression of impact ionization with respect to direct photovoltaic excitations. Acoustic phonons slightly enhance the photocurrent for small driving frequencies and suppress it at frequencies around the main peak at all considered hybridization strengths.

Correlated Mott insulators in strong electric fields: Role of phonons in heat dissipation

We study the spectral and transport properties of a Mott insulator driven by a static electric field into a non-equilibrium steady state. For the dissipation, we consider two mechanisms: Wide-band fermion reservoirs and phonons included within the Migdal approximation. The electron correlations are treated via non-equilibrium dynamical mean field theory with an impurity solver suitable for strong correlations. We find that dissipation via phonons is limited to restricted ranges of field values around Wannier-Stark resonances. To cover the full range of field strengths, we allow for a small coupling with fermionic baths, which stabilizes the solution. When considering both dissipation mechanisms, we find that phonons enhance the current for field strengths close to half of the gap while lowering it at the gap resonance as compared to the purely electronic dissipation used by Murakami and Werner [arXiv:1804.08257]. Once phonons are the only dissipation mechanism, the current in the metallic phase is almost one order of magnitude smaller than the typical values obtained by coupling to a fermionic bath. In this case, the transport regime is characterized by an accumulation of charge in the upper Hubbard band described by two effective chemical potentials.

Excitation and ionization of a particle in a one-dimensional potential well of zero radius by an extremely short light pulse

The Migdal sudden perturbation approximation is used to solve the problem of excitation and ionization particles in a one-dimensional potential of zero radius with an extremely short pulse. There is only one energy level in such a one-dimensional the delta-shaped potential well. It is shown that for pulse durations shorter than the characteristic period of oscillations of the wave function of the particle in the bound state, the population of the level (and the probability of ionization) is determined by the ratio of the electric area of the pulse to the characteristic "scale" of the area inversely proportional to the size of localization of the particle in a bound state. Keywords: unipolar pulses, electrical pulse area, atomic scale of electric area.

Возбуждение и ионизация частицы в одномерной потенциальной яме нулевого радиуса предельно коротким световым импульсом

The Migdal sudden perturbation approximation is used to solve the problem of excitation and ionization particles in a one-dimensional potential of zero radius with an extremely short pulse. There is has only one energy level in such a one-dimensional the delta-shaped potential well. It is shown that for pulse durations shorter than the characteristic period of oscillations of the wave function of the particle in the bound state, the population of the level (and the probability of ionization) is determined by the ratio of the electric the area of ​​the pulse to the characteristic “scale” of the area inversely proportional to the area of ​​localization of the particle in a bound state.

Unconventional superconductivity mediated solely by isotropic electron-phonon interaction

Unconventional superconductivity is commonly linked to electronic pairing mechanisms, since it is believed that the conventional electron-phonon interaction (EPI) cannot cause sign-changing superconducting gap symmetries. Here, we show that this common understanding needs to be revised when one considers a more elaborate theory of electron-phonon superconductivity beyond standard approximations. We self-consistently solve the full-bandwidth, anisotropic Eliashberg equations including vertex corrections beyond Migdal's approximation assuming the usual isotropic EPI for cuprate, Fe-based, and heavy-fermion superconductors with nested Fermi surfaces. In the case of the high-${T}_{c}$ cuprates we find a $d$-wave order parameter, as well as a nematic state upon increased doping. For Fe-based superconductors, we obtain ${s}_{\ifmmode\pm\else\textpm\fi{}}$ gap symmetry, while for heavy-fermion ${\mathrm{CeCoIn}}_{5}$ we find unconventional $d$-wave pairing. These results provide a proof of concept that EPI cannot be excluded as a mediator of unconventional and of high-${T}_{c}$ superconductivity.

Exploring multichannel superconductivity in ThFeAsN

We investigate theoretically the superconducting state of the undoped Fe-based superconductor ThFeAsN. Using input from $ab~initio$ calculations, we solve the Fermi-surface based, multichannel Eliashberg equations for Cooper-pair formation mediated by spin and charge fluctuations, and by the electron-phonon interaction (EPI). Our results reveal that spin fluctuations alone, when coupling only hole-like with electron-like energy bands, can account for a critical temperature $T_c$ up to $\sim7.5\,\mathrm{K}$ with an $s_{\pm}$-wave superconducting gap symmetry, which is a comparatively low $T_c$ with respect to the experimental value $T_c^{\mathrm{exp}}=30\,\mathrm{K}$. Other combinations of interaction kernels (spin, charge, electron-phonon) lead to a suppression of $T_c$ due to phase frustration of the superconducting gap. We qualitatively argue that the missing ingredient to explain the gap magnitude and $T_c$ in this material is the first-order correction to the EPI vertex. In the noninteracting state this correction adopts a form supporting the $s_{\pm}$ gap symmetry, in contrast to EPI within Migdal's approximation, i.e., EPI without vertex correction, and therefore it enhances tendencies arising from spin fluctuations.

Functional-integral approach to Gaussian fluctuations in Eliashberg theory

The Eliashberg theory of superconductivity is based on a dynamical electron-phonon interaction as opposed to a static interaction present in BCS theory. The standard derivation of Eliashberg theory is based on an equation of motion approach, which incorporates certain approximations such as Migdal's approximation for the pairing vertex. In this paper we provide a functional-integral-based derivation of Eliashberg theory and we also consider its Gaussian-fluctuation extension. The functional approach enables a self-consistent method of computing the mean-field equations, which arise as saddle-point conditions, and here we observe that the conventional Eliashberg self energy and pairing function both appear as Hubbard-Stratonovich transformations. An important consequence of this fact is that it provides a systematic derivation of the Cooper and density-channel interactions in the Gaussian fluctuation response. We also investigate the strong-coupling fluctuation diamagnetic susceptibility near the critical temperature.

Influence of phonon renormalization in Eliashberg theory for superconductivity in two- and three-dimensional systems

Eliashberg's foundational theory of superconductivity is based on the application of Migdal's approximation, which states that vertex corrections to first order electron-phonon scattering are negligible if the ratio between phonon and electron energy scales is small. The resulting theory incorporates the first Feynman diagrams for electron and phonon self-energies. However, the latter is most commonly neglected in numerical analyses. Here we provide an extensive study of full-bandwidth Eliashberg theory in two and three dimensions, where we include the full back reaction of electrons onto the phonon spectrum. We unravel the complex interplay between nesting properties, Fermi surface density of states, renormalized electron-phonon coupling, phonon softening, and superconductivity. We propose furthermore a scaling law for the maximally possible critical temperature $T_c^{\textrm{max}}\propto\lambda (\Omega ) \sqrt{\Omega_0^2-\Omega^2}$ in 2D and 3D systems, which embodies both the renormalized electron-phonon coupling strength $\lambda(\Omega)$ and softened phonon spectrum $\Omega$. Also, we analyze for which electronic structure properties a maximal $T_c$ enhancement can be achieved.

Relative importance of nonlinear electron-phonon coupling and vertex corrections in the Holstein model

Determining the range of validity of Migdal’s approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal’s approximation holds. Our results demonstrate that questions regarding the validity of Migdal’s approximation go hand in hand with questions of the validity of a linear e-ph interaction.

Combining Eliashberg Theory with Density Functional Theory for the Accurate Prediction of Superconducting Transition Temperatures and Gap Functions.

We propose a practical alternative to Eliashberg equations for the abinitio calculation of superconducting transition temperatures and gap functions. Within the recent density functional theory for superconductors, we develop an exchange-correlation functional that retains the accuracy of Migdal's approximation to the many-body electron-phonon self-energy, while having a simple analytic form. Our functional is based on a parametrization of the Eliashberg self-energy for a superconductor with a single Einstein frequency, and enables density functional calculations of experimental excitation gaps. By merging electronic structure methods and Eliashberg theory, the present approach sets a new standard in quality and computational feasibility for the prediction of superconducting properties.

Full-bandwidth Eliashberg theory of superconductivity beyond Migdal's approximation

We solve the anisotropic, full-bandwidth and non-adiabatic Eliashberg equations for phonon-mediated superconductivity by fully including the first vertex correction in the electronic self-energy. The non-adiabatic equations are solved numerically here without further approximations, for a one-band model system. We compare the results to those that we obtain by adiabatic full-bandwidth, as well as Fermi-surface restricted Eliashberg-theory calculations. We find that non-adiabatic contributions to the superconducting gap can be positive, negative or negligible, depending on the dimensionality of the considered system, the degree of non-adiabaticity, and the coupling strength. We further examine non-adiabatic effects on the transition temperature and the electron-phonon coupling constant. Our treatment emphasizes the importance of overcoming previously employed approximations in estimating the impact of vertex corrections on superconductivity and opens a pathway to systematically study vertex correction effects in systems such as high-$T_c$, flat band and low-carrier density superconductors.

Eliashberg theory of phonon-mediated superconductivity — When it is valid and how it breaks down

We analyze the validity of Eliashberg theory of phonon-mediated superconductivity in 2D systems in light of recent extensive Monte-Carlo studies of the Holstein model. Conventional wisdom says that Eliashberg theory is applicable as long as vertex corrections remain small. For small ratio of the phonon energy Ω0 and the Fermi energy EF, this condition is supposed to hold even when the dimensionless electron–phonon coupling λ is larger than one, i.e., in the strong coupling regime. A comparison between various quantities computed in the Migdal approximation and those computed by Quantum Monte Carlo prove that this belief is wrong, and we identify analytically some of the ways in which this breakdown occurs for various “normal state” properties at λ=λcr, where λcr=O(1). The breakdown occurs at temperatures high enough that neither superconducting nor charge–density wave correlations extend over any significant range of distances, so it cannot be associated with the onset of an instability toward any of the relevant ordered ground-states — rather it is associated with the local physics of classical bipolaron formation. Still, we show that certain properties, including the superconducting Tc and the superconducting gap structure below Tc, can be accurately inferred from the strong-coupling limit of Eliashberg theory at λ≤λcr.

Temperature-filling phase diagram of the two-dimensional Holstein model in the thermodynamic limit by self-consistent Migdal approximation

We study the temperature-filling phase diagram of the single-band Holstein model in two dimensions using the self-consistent Migdal approximation, where both the electron and phonon self-energies are treated on an equal footing. By employing an efficient numerical algorithm utilizing fast Fourier transforms to evaluate momentum and Matsubara frequency summations, we determine the charge-density-wave (CDW) and superconducting transition temperatures in the thermodynamic limit using lattice sizes that are sufficient to eliminate significant finite size effects present at lower temperatures. We obtain the temperature-filling phase diagrams for a range of coupling strengths and phonon frequencies for the model defined on a square lattice with and without next-nearest neighbor hopping. We find the appearance of a superconducting dome with a critical temperature that decreases before reaching the $\mathbf{q}_{\text{max}} = (\pi,\pi)$ CDW phase boundary. For very low phonon frequencies, we also find an incommensurate CDW phase with the ordering vector $\mathbf{q}_{\text{max}} \approx (\pi,\pi)$ appearing between the commensurate CDW and superconducting phases. Our numerical implementation can be easily extended to treat momentum-dependent electron-phonon coupling, as well as dispersive phonon branches, and has been made available to the public.

Fully gapped s -wave superconductivity enhanced by magnetic criticality in heavy-fermion systems

In heavy fermion systems, higher-rank multipole operators are active thanks to the strong spin-orbit interaction (SOI), and the role of diverse multipole fluctuations on the pairing mechanism attracts a lot of attention. Here, we study a mechanism of superconductivity in heavy fermion systems, by focusing on the impact of vertex corrections (VCs) for the pairing interaction going beyond the Migdal approximation. In heavy fermion systems, strong interference between multipole fluctuations cause significant VCs, that represent many-body effects beyond mean-field-type approximations. Especially, the coupling constants between electrons and charged-bosons, including the electron-phonon coupling constant, are strongly magnified by the VCs. For this reason, moderate even-rank (=electric) multipole fluctuations give large attractive interaction, and therefore s-wave superconductivity can emerge in heavy-fermion systems. In particular, phonon-mediated superconductivity is expected to be realized near the magnetic criticality, thanks to the VCs due to magnetic multipole fluctuations. The present mechanism may be responsible for the fully gapped s-wave superconducting state realized in CeCu2Si2.

Nonempirical Calculation of Superconducting Transition Temperatures in Light-Element Superconductors.

Recent progress in the fully nonempirical calculation of the superconducting transition temperature (Tc ) is reviewed. Especially, this study focuses on three representative light-element high-Tc superconductors, i.e., elemental Li, sulfur hydrides, and alkali-doped fullerides. Here, it is discussed how crucial it is to develop the beyond Migdal-Eliashberg (ME) methods. For Li, a scheme of superconducting density functional theory for the plasmon mechanism is formulated and it is found that Tc is dramatically enhanced by considering the frequency dependence of the screened Coulomb interaction. For sulfur hydrides, it is essential to go beyond not only the static approximation for the screened Coulomb interaction, but also the constant density-of-states approximation for electrons, the harmonic approximation for phonons, and the Migdal approximation for the electron-phonon vertex, all of which have been employed in the standard ME calculation. It is also shown that the feedback effect in the self-consistent calculation of the self-energy and the zero point motion considerably affect the calculation of Tc . For alkali-doped fullerides, the interplay between electron-phonon coupling and electron correlations becomes more nontrivial. It has been demonstrated that the combination of density functional theory and dynamical mean field theory with the ab initio downfolding scheme for electron-phonon coupled systems works successfully. This study not only reproduces the experimental phase diagram but also obtains a unified view of the high-Tc superconductivity and the Mott-Hubbard transition in the fullerides. The results for these high-Tc superconductors will provide a firm ground for future materials design of new superconductors.

Breakdown of the Migdal approximation at Lifshitz transitions with giant zero-point motion in the H3S superconductor.

While 203 K high temperature superconductivity in H3S has been interpreted by BCS theory in the dirty limit here we focus on the effects of hydrogen zero-point-motion and the multiband electronic structure relevant for multigap superconductivity near Lifshitz transitions. We describe how the topology of the Fermi surfaces evolves with pressure giving different Lifshitz-transitions. A neck-disrupting Lifshitz-transition (type 2) occurs where the van Hove singularity, vHs, crosses the chemical potential at 210 GPa and new small 2D Fermi surface portions appear with slow Fermi velocity where the Migdal-approximation becomes questionable. We show that the neglected hydrogen zero-point motion ZPM, plays a key role at Lifshitz transitions. It induces an energy shift of about 600 meV of the vHs. The other Lifshitz-transition (of type 1) for the appearing of a new Fermi surface occurs at 130 GPa where new Fermi surfaces appear at the Γ point of the Brillouin zone here the Migdal-approximation breaks down and the zero-point-motion induces large fluctuations. The maximum Tc = 203 K occurs at 160 GPa where EF/ω0 = 1 in the small Fermi surface pocket at Γ. A Feshbach-like resonance between a possible BEC-BCS condensate at Γ and the BCS condensate in different k-space spots is proposed.

Effect of Van Hove singularities on high- Tc superconductivity in H3S

One of interesting open questions for the high transition temperature (Tc) superconductivity in sulfur hydrides is why high pressure phases of H3S have extremely high Tc's. Recently, it has been pointed out that the presence of the van Hove singularities (vHs) around the Fermi level is crucial. However, while there have been quantitative estimates of Tc based on the Migdal-Eliashberg theory, the energy dependence of the density of states (DOS) has been neglected to simplify the Eliashberg equation. In this study, we go beyond the constant DOS approximation and explicitly consider the electronic structure over 40eV around the Fermi level. In contrast with the previous conventional calculations, this approach with a sufficiently large number of Matsubara frequencies enables us to calculate Tc without introducing the empirical pseudo Coulomb potential. We show that while H3S has much higher Tc than H2S for which the vHs is absent, the constant DOS approximation employed so far seriously overestimates (underestimates) Tc by ~ 60K (~ 10K) for H3S (H2S). We then discuss the impact of the strong electron-phonon coupling on the electronic structure with and without the vHs and how it affects the superconductivity. Especially, we focus on (1) the feedback effect in the self-consistent calculation of the self-energy, (2) the effect of the energy shift due to the zero-point motion, and (3) the effect of the changes in the phonon frequencies due to strong anharmonicity. We show that the effect of (1)-(3) on Tc is about 10-30K for both H3S and H2S. Eventually, Tc is estimated to be 181K for H3S at 250GPa and 34K for H2S at 140GPa, which explains the pressure dependence of Tc observed in the experiment. In addition, we evaluate the lowest order vertex correction beyond the Migdal-Eliashberg theory and discuss the validity of the Migdal approximation for sulfur hydrides.

Multiple amplitude modes in strongly coupled phonon-mediated superconductors

We study collective amplitude modes of the superconducting order parameter in strongly-coupled electron-phonon systems described by the Holstein model using the nonequilibrium dynamical mean-field theory with the self-consistent Migdal approximation as an impurity solver. The frequency of the Higgs amplitude mode is found to coincide with the superconducting gap even in the strongly-coupled (beyond BCS) regime. Besides the Higgs mode, we unravel another collective mode involving the dynamics of both the phonons and the superconducting order parameter. The frequency of this mode, higher than twice the renormalized phonon frequency in the superconducting phase, is shown to reflect a strong electron-mediated phonon-phonon interaction. Both types of collective excitations contribute to time-resolved photoemission spectra after a strong laser pump as vertex corrections to produce resonance peaks, which allows one to distinguish them from quasiparticle excitations.