We investigate the dynamics of a phonon-mediated superconductor driven out of equilibrium. The electronic hopping amplitude is ramped down in time, resulting in an increased electronic density of states. The dynamics of the coupled electron-phonon model is investigated by solving Migdal-Eliashberg equations for the double-time Keldysh Green's functions. The increase of the density of states near the Fermi level leads to an enhancement of superconductivity when the system thermalizes to the new state at the same temperature. We provide a time- and momentum-resolved view on this thermalization process and show that it involves fast processes associated with single-particle scattering and much slower dynamics associated with the superconducting order parameter. The importance of electron-phonon coupling for the rapid enhancement and the efficient thermalization of superconductivity is demonstrated, and the results are compared to a BCS time-dependent mean-field approximation.
We investigate the dynamics of a phonon-mediated superconductor driven out of equilibrium. The electronic hopping amplitude is ramped down in time, resulting in an increased electronic density of states. The dynamics of the coupled electron-phonon model is investigated by solving Migdal-Eliashberg equations for the double-time Keldysh Green's functions. The increase of the density of states near the Fermi level leads to an enhancement of superconductivity when the system thermalizes to the new state at the same temperature. We provide a time- and momentum-resolved view on this thermalization process and show that it involves fast processes associated with single-particle scattering and much slower dynamics associated with the superconducting order parameter. The importance of electron-phonon coupling for the rapid enhancement and the efficient thermalization of superconductivity is demonstrated, and the results are compared to a BCS time-dependent mean-field approximation.
Motivated by the experiment by M. Budden {\em et al.} [Nature Physics {\bf 17}, 611 (2021)], who observed signatures of long-lived photo-induced superconductivity, we develop an accurate analytical/computational approach to non-equilibrium superconductivity following a quench. We consider the BCS-Holstein model, which includes both integrable local electron-electron interactions and integrability-breaking electron-phonon coupling. We develop Keldysh-Eliashberg theory on the Kadanoff-Baym contour, which enables us to describe non-equilibrium dynamics of the superconductor. We consider a quench in interactions, which results in a dynamic transition from the initial superconducting state to a normal thermal state in the end of the evolution. It is shown that the dynamics contain two stages: The early-time integrable behavior, involving coherent oscillations of the superconducting order parameter, crosses over to the late-time ergodic dynamics exhibiting a thermal decay into an equilibrium state. In the former regime, our computational approach both reproduces exact analytical results on the integrable dynamics of the order parameter and generalizes those to the case of an initial thermal state. The method also succeeds for the first time in describing both integrable-to-thermalizing crossover and the late-time thermal decay, which is shown to be consistent with the time-dependent Ginzburg-Landau theory (with the exponential decay time dependent on the density of quasiparticle excitations). We observe the electron distribution function approaching the Fermi-Dirac thermal distribution at final stages. The details of two-time non-equilibrium dynamics depend on the density of quasiparticles in the initial state and the integrability-breaking parameters, which under certain conditions may result in a long-lived transient superconductivity consistent with experiment.
Motivated by the observation of superconductivity in SrTiO3 and Bi, we analyze phonon-mediated superconductivity in three-dimensional systems at low carrier density, when the chemical potential $\ensuremath{\mu}$ is comparable to or even smaller than the characteristic phonon frequency ${\ensuremath{\omega}}_{L}$. We consider the attractive part of the Bardeen-Pines pairing interaction, in which the frequency-dependent electron-phonon interaction is dressed by the Coulomb potential. This dressing endows the pairing interaction with momentum dependence. We argue that the conventional Migdal-Eliashberg approximation becomes invalid when $\ensuremath{\mu}\ensuremath{\le}{\ensuremath{\omega}}_{L}$, chiefly because the dominant contribution to pairing comes from electronic states away from the Fermi surface. We obtain the pairing onset temperature, which is equal to Tc in the absence of phase fluctuations, as a function of $\ensuremath{\mu}/{\ensuremath{\omega}}_{L}$. We find both analytically and numerically that Tc increases as the ratio $\ensuremath{\mu}/{\ensuremath{\omega}}_{L}$ becomes smaller. In particular, in the dilute regime, $\ensuremath{\mu}\ensuremath{\rightarrow}0$, it holds that ${T}_{c}\ensuremath{\propto}{\ensuremath{\omega}}_{L}{(\frac{Ry}{{\ensuremath{\omega}}_{L}})}^{\ensuremath{\eta}}$, where Ry is the Rydberg constant and $\ensuremath{\eta}\ensuremath{\sim}0.2$.
We investigate the nonequilibrium dynamics of the weak-coupling Hubbard-Holstein model after a sudden switch-on of the electron-phonon interaction within nonequilibrium dynamical mean-field theory (DMFT). Using the self-consistent Migdal approximation for the electron-phonon coupling together with second-order perturbation theory for the electron-electron interaction, we show that the relaxation dynamics exhibits a crossover between electron-dominated and phonon-dominated regimes, extending to finite Hubbard interaction the scenario previously identified in the Holstein model. To investigate the microscopic buildup of the thermal state, we analyze the dynamics within the Step-by-Step DMFT framework. In the plane of real time and DMFT iteration number, thermalization is marked by a sharp propagating front. This front appears in electronic observables already for weak quenches within the simulated time window, whereas the phononic sector exhibits a visible front only at sufficiently strong coupling. Thus, at weak coupling the local dispersionless phonons show a delayed onset of front formation, while near and beyond the crossover the front develops on comparable timescales in both the electronic and phononic sectors. Whenever both fronts are resolved, they propagate with the same velocity, showing that thermalization spreads coherently through the coupled electron-phonon system.
We study the relaxation of the Holstein model after a sudden switch-on of the interaction by means of the nonequilibrium dynamical mean field theory, with the self-consistent Migdal approximation as an impurity solver. We show that there exists a qualitative change in the thermalization dynamics as the interaction is varied in the weak-coupling regime. On the weaker interaction side of this crossover, the phonon oscillations are damped more rapidly than the electron thermalization time scale, as determined from the relaxation of the electron momentum distribution function. On the stronger interaction side, the relaxation of the electrons becomes faster than the phonon damping. In this regime, despite long-lived phonon oscillations, a thermalized momentum distribution is realized temporarily. The origin of the ``thermalization crossover'' found here is traced back to different behaviors of the electron and phonon self-energies as a function of the electron-phonon coupling. In addition, the importance of the phonon dynamics is demonstrated by comparing the self-consistent Migdal results with those obtained with a simpler Hartree-Fock impurity solver that neglects the phonon self-energy. The latter scheme does not properly describe the evolution and thermalization of isolated electron-phonon systems.
Optimizing various properties of materials have been under intensive focus for many decades. In recent years, the exposure of ultrafast techniques enhance the possibility of seizing this goal. As by employing these techniques, one can expose the system under strong short fields, and investigate the transient response of various low-lying excitations. Such studies will strengthen our understanding of interplays among various degrees of freedom in correlated materials, and thereby facilitate the chance to manipulate different properties of these systems. Of our particular attention is the electron-lattice dynamics, where modulating a particular vibrational mode of the crystal lattice, will lead to the emergence of new quasi stable'' lattice structures as well as novel transient phenomena. Aside from these experimental prospects, the above mentioned goals were also the subject of an extensive theoretical studies. In this thesis, we investigated the in and out of equilibrium dynamics of electronic systems which are interacting with lattice degrees of freedom. We employed the two-time Keldysh Green's function within the framework of the dynamical mean field theory to conduct our researches on this subject. In our first attempt, we studied the dynamical response of a one-electron system under a linear electron-lattice interaction. We explored the real-time formation of a bound quasiparticle, known as a polaron, under various parameter regimes of the model. We presented that when the carrier is suddenly coupled to the lattice degrees of freedom, the long-time response of the system in the strong coupling limit exhibits both bound, polaronic, and unbound, delocalized, states. We characterized the nature of such a mixed state using the adiabatic picture of the model. We also showed that the formed polaron is dressed by excited phononic states, in contrast to its equilibrium counterpart, where the ground-state phononic states are bounding the electron. We also addressed the response of a many-body electronic system, weakly coupled to a bosonic bath, under various time dependent excitations, which brings the system to a bad-metallic phase. We employed a model which treated the Coulomb repulsive interaction locally. We presented that the relaxation of this excited state varies by the strength of the electron-electron interaction. We demonstrated that, for both insulating and metallic regimes of the system, a slow relaxation dynamics can be observed close to the phase transition from the metallic to insulating regimes of the system. Our analysis revealed that in the metallic side, this slow evolution is rooted in the spin-related physics of the problem. However, in the insulating regimes of the model, our results showed that presence of a gap in the spectral density of the insulating system prohibits the relaxation of the excited particles. We also studied the associated timescale for the adiabatic dynamics near the Mott transition. Our results exhibited that there is a nonadiabatic window'' in the metallic side of the phase transition, in which the adiabatic timescale is extremely long. We pointed out that the presence of such a window is intertwined with the slow relaxation timescale of the spin-related dynamics close to the Mott transition. We also captured an induced metallic phase near metal-to-insulator transition out of equilibrium. We initiated this study using the slave-rotor impurity solver, which is been coupled to the self-energies obtained from the first-order Luttinger-Ward functional for the Phonon degrees of freedom. We presented the emergence of the larger spectral density at small frequencies both in metallic and insulating phases close to the phase transition, which is surrounded by pseudo gaps. We explored the physics of such spectral densities as well as their dependencies on the phononic parameters. We postulated that this response of the system is due to the evolution of the system, in a non-thermal mechanism, towards a temperature which is less than the initial temperature of the system, before coupling to the bath. We also pointed out that the associated very slow relaxation dynamics is the result of the formation of psuedo-gaps around the quasiparticle peak.
We investigate an exact nonequilibrium solution of a two-site electron-phonon model, where an infrared-active phonon that is nonlinearly coupled to the electrons is driven by a laser field. The time-resolved electronic spectrum shows coherence-incoherence spectral weight transfer, a clear signature of light-enhanced electron-phonon coupling. The present study is motivated by recent evidence for enhanced electron-phonon coupling in pump-probe terahertz and angle-resolved photoemission spectroscopy in bilayer graphene when driven near resonance with an infrared-active phonon mode [E. Pomarico et al., Phys. Rev. B 95, 024304 (2017)], and by a theoretical study suggesting that transient electronic attraction arises from nonlinear electron-phonon coupling [D. M. Kennes et al., Nat. Phys. 13, 479 (2017)]. We show that a linear scaling of light-enhanced electron-phonon coupling with the pump field intensity emerges, in accordance with a time-nonlocal self-energy based on a mean-field decoupling using quasiclassical phonon coherent states. Finally, we demonstrate that this leads to enhanced double occupancies in accordance with an effective electron-electron attraction. Our results suggest that materials with strong phonon nonlinearities provide an ideal playground to achieve light-enhanced electron-phonon coupling and possibly light-induced superconductivity.