In equilibrium systems amplitude and phase collective modes are decoupled, as they are mutually orthogonal excitations. The direct detection of these Higgs and Leggett collective modes by linear-response measurements is not possible, because they do not couple directly to the electromagnetic field. In this work, using numerical exact simulations we show for the case of two-gap superconductors, that optical pump–probe experiments excite both Higgs and Leggett modes out of equilibrium. We find that this non-adiabatic excitation process introduces a strong interaction between the collective modes, which is absent in equilibrium. Moreover, we propose a type of pump–probe experiment, which allows to probe and coherently control the Higgs and Leggett modes, and thus the order parameter directly. These findings go beyond two-band superconductors and apply to general collective modes in quantum materials.

We theoretically study the pump-probe response of nonequilibrium BCS superconductors coupled to optical phonons. For ultrashort pump pulses a nonadiabatic regime emerges, which is characterized by oscillations of the superconducting order parameter as well as by the generation of coherent phonons. Using the density-matrix formalism, we compute the pump-probe response in the nonadiabatic regime of the coupled Bogoliubov quasiparticle-phonon system and determine the signatures of the order parameter and the phonon oscillations in the pump-probe conductivity. We find that the nonadiabatic dynamics of the BCS superconductor reflects itself in oscillations of the pump-probe response as functions of delay time δt between pump and probe pulses. We argue that from the analysis of this oscillatory behavior both frequency and decay time of the algebraically decaying order-parameter oscillations can be inferred. Similarly, the coherent phonons are evidenced in the pump-probe conductivity by oscillations with the frequency of the phonons. Remarkably, we find that the oscillatory response in the pump-probe conductivity is resonantly enhanced when the frequency of the order-parameter oscillations is tuned to the phonon energy.

Unitary transformations are an essential tool for the theoretical understanding of many systems by mapping them to simpler effective models. A systematically controlled variant to perform such a mapping is a perturbative continuous unitary transformation (pCUT) among others. So far, this approach required an equidistant unperturbed spectrum. Here, we pursue two goals: First, we extend its applicability to non-equidistant spectra with the particular focus on an efficient derivation of the differential flow equations, which define the enhanced perturbative continuous unitary transformation (epCUT). Second, we show that the numerical integration of the flow equations yields a robust scheme to extract data from the epCUT. The method is illustrated by the perturbation of the harmonic oscillator with a quartic term and of the two-leg spin ladders in the strong-rungcoupling limit for uniform and alternating rung couplings. The latter case provides an example of perturbation around a non-equidistant spectrum.

Recent findings of new Higgs modes in unconventional superconductors require a classification and characterization of the modes allowed by nontrivial gap symmetry. Here we develop a theory for a tailored nonequilibrium quantum quench to excite all possible oscillation symmetries of a superconducting condensate. We show that both a finite momentum transfer and quench symmetry allow for an identification of the resulting Higgs oscillations. These serve as a fingerprint for the ground state gap symmetry. We provide a classification scheme of these oscillations and the quench symmetry based on group theory for the underlying lattice point group. For characterization, analytic calculations as well as full scale numeric simulations of the transient optical response resulting from an excitation by a realistic laser pulse are performed. Our classification of Higgs oscillations allows us to distinguish between different symmetries of the superconducting condensate.

The density matrix formalism and the equation of motion approach are two semi-analytical methods that can be used to compute the non-equilibrium dynamics of correlated systems. While for a bilinear Hamiltonian both formalisms yield the exact result, for any non-bilinear Hamiltonian a truncation is necessary. Due to the fact that the commonly used truncation schemes differ for these two methods, the accuracy of the obtained results depends significantly on the chosen approach. In this paper, both formalisms are applied to the quantum Rabi model. This allows us to compare the approximate results and the exact dynamics of the system and enables us to discuss the accuracy of the approximations as well as the advantages and the disadvantages of both methods. It is shown to which extent the results fulfill physical requirements for the observables and which properties of the methods lead to unphysical results. PACS. 02.30.Mv Approximations and expansions -02.60.Cb Numerical simulation; solution of equations -05.30.Jp Boson systems arXiv:1701.07282v3 [cond-mat.str-el]