We show that the gauge-invariant kinetic equation of superconductivity provides an efficient approach to study the electromagnetic response of the gapless Nambu-Goldstone and gapful Higgs modes on an equal footing. We prove that the Fock energy in the kinetic equation is equivalent to the generalized Ward's identity. Hence, the gauge invariance directly leads to the charge conservation. Both linear and second-order optical responses are analytically investigated. The linear response of the Higgs mode vanishes in the long-wave limit, whereas the linear response of the Nambu-Goldstone mode interacts with the long-range Coulomb interaction, causing the original gapless energy spectrum effectively lifted up to the plasma frequency as a result of the Anderson-Higgs mechanism, consistent with previous work. The second-order response exhibits interesting physics. On the one hand, a finite second-order optical response of the Higgs mode is obtained in the long-wave limit. We reveal that this response, which has been experimentally observed, is attributed solely to the drive effect rather than the widely considered Anderson-pump effect. On the other hand, the second-order optical response of the Nambu-Goldstone mode, free from the influence of the long-range Coulomb interaction and hence the Anderson-Higgs mechanism, is predicted. We find that both Anderson-pump and drive effects play important role in this response. A tentative scheme to detect this second-order response is proposed.
We theoretically investigate collective excitations in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states of Pauli-limited superconducting films. When the long-range Coulomb interaction is absent, excitation spectra consist of two gapless and three gapped modes. The gapless modes are the Nambu-Goldstone modes associated with the spontaneous breaking of the
U(1) symmetry and the translational symmetry. The gapped modes include the Higgs mode and the twofold degenerate modes that cause the oscillation of the domain width and grayness of FFLO nodal planes. We find that the long-range Coulomb interaction only gaps out the gapless phase mode through the Anderson-Higgs mechanism, while the other modes remain unaffected. Furthermore, the field evolution of the dispersion of the gapless elastic mode, the Nambu-Goldstone mode associated with the translational symmetry breaking, is associated with that of the bandwidth of the mid-gap Andreev bound states. We demonstrate that the signature of the elastic mode can be detected by measuring the density-density response function.
Within a gauge-invariant microscopic kinetic theory, we study the electromagnetic response in the superconducting states. Both superfluid and normal-fluid dynamics are involved. We predict that the normal fluid is present only when the excited superconducting velocity
vs is larger than a threshold ${v}_{L}=|\mathrm{\ensuremath{\Delta}}|/{k}_{F}$. Interestingly, with the normal fluid, we find that there exists friction between the normal-fluid and superfluid currents. Due to this friction, part of the superfluid becomes viscous. Therefore a three-fluid model, normal fluid and nonviscous and viscous superfluids, is proposed. For the stationary magnetic response, at
vs<vL with only the nonviscous superfluid, the Meissner supercurrent is excited and the gap equation can be reduced to the Ginzburg-Landau equation. At ${v}_{s}\ensuremath{\ge}{v}_{L}$, with the normal fluid and nonviscous and viscous superfluids, in addition to the directly excited Meissner supercurrent in the superfluid, a normal-fluid current is also induced through the friction drag with the viscous superfluid current. Due to the normal-fluid and viscous-superfluid currents, the penetration depth is influenced by the scattering effect. In addition, a modified Ginzburg-Landau equation is proposed. We predict an exotic phase in which both the resistivity and superconducting gap are finite. As for the optical response, the excited
vs oscillates with time. When
vs<vL, only the nonviscous superfluid is present, whereas at ${v}_{s}\ensuremath{\ge}{v}_{L}$, normal fluid and nonviscous and viscous superfluids are present. We show that the excited normal-fluid current exhibits the Drude-model behavior, while the superfluid current consists of the Meissner supercurrent and Bogoliubov quasiparticle current. Due to the friction between the superfluid and normal-fluid currents, the optical conductivity is captured by the three-fluid model. Finally, we also study the optical excitation of the Higgs mode. By comparing the contributions from the drive and Anderson-pseudospin pump effects, we find that the drive effect is dominant at finite temperature whereas at zero temperature, both effects contribute.
This thesis reports on theoretical investigations of the ultrafast nonequilibrium dynamics in superconductors. In particular, two nonequilibrium scenarios are modeled by combining several analytical and numerical methods: transient superconducting phase transition and collective excitations of a superconducting condensate. In the first part of the thesis the possibility to induce singlet superconducting correlations far from equilibrium is explored by applying time-dependent Lanczos algorithm to the one-dimensional extended Hubbard model. Here, the system is driven into a nonequilibrium state by (i) interaction quench and (ii) quench by pulse. In both cases the buildup of the superconducting correlations is studied in view of several correlation functions. Moreover, a transient Meissner effect, which reflects fingerprints of the induced superconductivity, is observed in calculations of the important for experiments time-dependent optical conductivity. In addition, the differences between nonequilibrium dynamics of the extended Hubbard model initiated by interaction quench and by pulse quench are discussed in detail. The second part is devoted to the excitation of collective modes in a superconductor out of equilibrium. First of all, focus is given to the properties of collective modes close to equilibrium conditions, where a linearized response can be expected. Within the framework of the (Nambu-) matrix kinetic theory in the clean limit a comprehensive gauge-invariant description of the collective modes associated with phase fluctuations of the order parameter is analytically formulated for the generalized case of two- band superconducting systems. Particular emphasis is on the application to the recently discovered noncentrosymmetric superconductors. A detailed analysis of the collective modes includes specification of the gauge mode (common to all superconductors) and the Leggett mode (collective fluctuations of the interband phase difference), and investigation of their subtle interplay in view of validity of the charge conservation law and the participation in the Anderson-Higgs mechanism. Further, the nonequilibrium excitation of the superconducting gaps for the case of two-band superconductors with singlet pairing correlations is investigated in a view of optical pump-probe experiments. For this purpose a semi-numerical method based on the density- matrix theory is formulated. In this case additional amplitude (Higgs) modes were identified, which allow to directly detect the properties of the superconducting condensates as a function of time. Moreover, a new interplay between the phase (Leggett) and amplitude (Higgs) modes is studied in detail. In particular it was found that the nonadiabatic excitation process leads to strong coupling between both modes. Finally, based on calculations of the time-dependent optical conductivity the prediction is made that the coupled Higgs and Leggett modes can be clearly visible in the pump-probe absorption spectra.
In superconductors gauge U(1) symmetry is spontaneously broken. According to Goldstone theorem this breaking of a continuous symmetry establishes the existence of the Bogoliubov phase mode while the gauge-invariant response also includes the amplitude fluctuations of the order parameter. The latter, which are also termed ’Higgs’ modes in analogy with the standard model, appear at the energy of the spectral gap 2Δ when the superconducting ground state is evaluated within the weak coupling BCS theory and are therefore damped. Previously we have shown that in the time-dependent Gutzwiller approximation (TDGA) Higgs modes appear inside the gap with a finite binding energy relative to the quasiparticle continuum. Here we show that the binding energy of the Higgs mode becomes exponentially small in the weak coupling limit converging to the BCS solution. On the other hand, well-defined undamped amplitude modes exist in strongly coupled superconductors when the interaction energy becomes of the order of the bandwidth.
A diagrammatic method of obtaining exact gauge-invariant response functions in strongly correlated Fermi superfluids is implemented for several example condensed matter systems of current interest. These include: topological superfluids, high temperature superconductors, and superfluids with finite center-of-mass momentum pairing known as Fulde-Ferrell superfluids. Much of the literature on these systems has focused on single-particle properties or alternatively has invoked simple approximations to treat response functions. The goal is to show that, for this wide class of topical problems, one can compute exact response functions. This enables assessment of the validity of different physical scenarios and allows a very broad class of experiments to be addressed.,The method developed is based on deriving the full electromagnetic vertex, which satisfies the Ward-Takahashi identity, and determining the collective modes in a manner compatible with the self-consistent gap equation. In the condensed phase of a superfluid and a superconductor, where gauge invariance is spontaneously broken, it is crucial to determine the collective modes from the gap equation in a manner which restores gauge invariance. Our diagrammatic framework provides a very general and powerful method for obtaining these collective modes in a variety of strongly correlated Fermi superfluids. We show that a full electromagnetic vertex satisfying the Ward-Takahashi identity ensures the
f-sum rule is satisfied and thus charge is conserved. ,This diagrammatic method is implemented for both normal and superfluid phases. While there are no collective modes in the normal phase, the Ward-Takahashi identity plays a similarly important role. In particular, for the normal phase we study Rashba spin-orbit coupled Fermi gases with intrinsic pairing in the absence and presence of a magnetic field. Exact density and spin response functions are obtained, even in the absence of a spin conservation law, providing signatures of topological order. Another set of examples studied in the normal phase include a broad collection of models of the cuprate pseudogap. For these cases too, the exact response functions are obtained and signatures of each model are identified. A detailed discussion of the superfluid phase is provided by the study of Fulde-Ferrell superfluids. The generally omitted collective modes are computed and the exact gauge-invariant response is obtained. The first complete calculation of the superfluid density, incorporating the collective mode contribution and without any assumptions on the size of the gap, is provided. It is found that the amplitude mode causes the superfluid to become unstable at temperatures much lower than the mean-field transition temperature.,As an alternative to the diagrammatic approach, the functional path integral technique for Fermi superfluids is presented. While this approach is widely studied in the literature, it is of restricted utility, in comparison to the diagrammatic approach presented here, due to the inability to perform the functional integral except in very limited cases. A short coming in the current literature is identified, in that the standard method of deriving gauge-invariant electrodynamics leads to a violation of the compressibility sum rule. It is found that to satisfy this sum rule the amplitude mode of the order parameter must be incorporated. A reformulation of the path integral method for Fermi superfluids is derived, which has both gauge-invariant electrodynamics and thermodynamics compatible with the compressibility sum rule, to all orders of approximation.
The Higgs mode is a key component in the spontaneous breaking of a continuous symmetry along with the Nambu-Goldstone mode, and has been studied extensively for homogeneous systems. We consider it for inhomogeneous systems, using the superfluid of harmonically trapped ultracold Fermi atomic gas. The Fermionic field operators are expanded in a complete set of wave functions corresponding to inhomogeneous situation. Within the Hatree-Fock approximation, we derive integral equations from the Bogoliubov-de Gennes equations, which lead to the frequencies of the collective modes, including the Higgs and Nambu-Goldstone modes. The results show that the frequency of the Higgs mode equals twice the absolute value of the order parameter at the center of trap. This feature is robust against variations in the interaction strength, trap potential, and temperature. These results are consistent with previous theoretical and experimental studies.