We present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

• Received 16 June 2020
• Accepted 7 August 2020

DOI:https://doi.org/10.1103/PhysRevB.102.054517