Extended time-dependent Ginzburg-Landau theory Konstantin V. Grigorishin∗ Bogolyubov Institute for Theoretical Physics of the National Academy of Sciences of Ukraine, 14-b Metrolohichna str. Kiev, 03143, Ukraine. We formulate the gauge invariant Lorentz covariant Ginzburg-Landau theory which describes nonstationary regimes: relaxation of a superconducting system accompanied by eigen oscillations of internal degrees of freedom (Higgs mode and Goldstone mode), and also forced oscillations under the action of an external gauge field. The theory describes Lorentz covariant electrodynamics of superconductors where Anderson-Higgs mechanism occurs, at the same time the dynamics of conduction electrons remains non-relativistic. It is demonstrated that Goldstone oscillations cannot be accompanied by oscillations of charge density and they generate the transverse field only. In addition, we consider Goldstone modes and features of Anderson-Higgs mechanism in two-band superconductors. We study dissipative processes, which are caused by movement of the normal component of electron liquid and violate the Lorentz covariance, on the examples of the damped oscillations of the order parameter and the skin-effect for electromagnetic waves. An experimental consequence of the extended time-dependent Ginzburg-Landau theory regarding the penetration of the electromagnetic field into a superconductor is proposed. PACS numbers: 74.20.De, 74.25.-q, 74.25.N-, 74.25.Nf Keywords: Ginzburg-Landau theory, Lorentz covariance, gauge invariance, Anderson-Higgs mechanism, Lon- don penetration depth, wave skin-effect I. INTRODUCTION The system undergoing the second-order phase transition (on example of superconductivity) is described with a Landau functional Z Z 2 b F = d3rF = d3r ~ (rΨ) rΨ+ + a jΨj2 + jΨj4 ; (1) 4m 2 iθ(r) where F is density of free energy, a = α(T − Tc), Ψ = jΨ(r)j e is a two-component order parameter (the wave 2 function of condensate of Cooper pairs) so that ns = 2 jΨj is density of superconducting (SC) electrons, a term 2 ~ 2 4m jrΨj can be understood as density of kinetic energy of Cooper pairs of mass 2m each. Configuration of the field Ψ(r) which minimizes the functional (1) is obtained from equation: δF @F @F 2 = − r = 0 ) ~ ∆Ψ − aΨ − b jΨj2 Ψ = 0: (2) δΨ+ @Ψ+ @(rΨ+) 4m 2 b 4 This configuration is shown in Fig.1a: they say that the field is a string laying in a valley of the potential a jΨj + 2 jΨj . At the same time, out the equilibrium the field Ψ depends on time. For small deviations from the equilibrium it is natural to assume [1, 2], that the time derivative @Ψ=@t is proportional to the variational derivative of the free energy functional δF/δΨ+ which is equal to zero at the equilibrium. Thus, one can write the time-dependent Ginzburg-Landau equation [1{10] (TDGL equation) in a form: 2 @Ψ δF @ ~ = − ) τ = ξ2∆ + − j j2 ; (3) 4mD @t δΨ+ @t arXiv:1906.02097v8 [cond-mat.supr-con] 27 May 2021 p where = Ψ=Ψ0 is the dimensionless order parameter, Ψ0 = −a=b is an equilibrium value of the spatial homo- geneous order parameter, D is a diffusion coefficient of electrons in the normal state. The temperature-dependent q 2 2 ~ ~ coherence length ξ(T ) = 4mja(T )j and the temperature-dependent relaxation time τ(T ) = 4mDja(T )j can be found from the microscopic theory for the case of a gapless SC alloy containing a high concentration of paramagnetic im- purities [3, 4] and for dirty superconductors in the Ginzburg-Landau (GL) regime jTc − T j << Tc [2, 5]. However, in ∗Electronic address: [email protected] 2 Figure 1: Variation of Landau free energy Fs − Fn (for example, SC state s versus normal state n) with the order parameter Ψ. (a) - equilibrium state of the field, (b) - monotonous relaxation to the equilibrium, (c) - relaxation to the equilibrium with damped oscillation, (d) - forced oscillation under the action of an external field. the general case, we have to treat these parameters phenomenologically. For T > Tc the TDGL equation has a form @Ψ 2 π τ0 = ξ ∆Ψ − Ψ, where τ0 = ~ is the relaxation time for a homogeneous mode (here and further kB = 1) [6]. @t 8(T −Tc) In this case the equilibrium value is hΨi = 0 but Ψ2 6= 0. The corresponding relaxation processes are illustrated in Fig.1b. Thus, Eq.(3) describes relaxation of the order parameter at small deviations from the equilibrium. At the same time, in the system undergoing the second-order phase transition the collective excitations can exist which are resonant oscillations. When a continuous symmetry is spontaneously broken, there emerge two types of collective modes in general: massive Higgs mode, which is oscillation of modulus jΨj of the order parameter, and Goldstone mode, which is oscillation of the phase θ. As it has been demonstrated in [11{13] in collisionless approximation the Higgs mode in SC system can exist as perturbation of amplitude of the gap ∆ which takes the form of oscillations having a frequency ∼ j∆j. At the same time, the oscillations of order parameter damp in some time τ. Thus, the system is out from the equilibrium then the relaxation process depends on relation between period of the eigen oscillations ∼ 1=! and the damping time: if 1=! τ then aperiodic relaxation occurs as shown in Fig.1b and it is described with TDGL equation (3). In general case the parameters 1=! and τ can be in arbitrary relation, hence the relaxation process can have more complicated form, for example, if 1=! τ then an oscillatory process with small damping occurs as shown in Fig.1c. Moreover, an external field can swing the system, that is the undamped oscillations occur, while heat is released (for example, electromagnetic wave, falling on a superconductor, induces Foucault currents both normal jn and superconducting js). Such situation is shown in Fig.1d. In presence of electric ' and magnetic A potentials the replacements @ @ i2e i2e ! + '; r ! r − A (4) @t @t ~ c~ must be done in free energy functional (1) and in Eqs.(2,3) for a gauge invariance. Moreover, the total current is a sum 1 @A ie~ + + of normal current and supercurrent: j = σ −∇' − c @t +js, where σ is conductivity, js = − 2m (Ψ rΨ − ΨrΨ )− 3 2e2 2 mc jΨj A. Thus, the TDGL equations determine dynamics of both order parameter and electromagnetic field Aµ ≡ ('; A). On the other hand, equations for electromagnetic field should be Lorentz covariant both in vacuum and within any medium, despite dynamics of particles (medium) is non-relativistic, for example, Maxwell equations in dielectrics. However the light speed in the medium is less than the speed in vacuum c and is determined with dynamic properties of the system. At the same time, the dissipative mechanisms give terms which violate Lorentz covariance since the dissipation distinguishes a time direction, i.e., violates the time symmetry t $ −t which is symmetry of the Lorentz boost. For example, Maxwell equation in medium with conductivity σ has a form: 1 @D 4π curlH = + σE: (5) c @t c 4π Here, the first two terms are an equation from Lorentz covariant Maxwell equations, the last term c σE is a dissipative part. Depending on material and processes occurring in it, some terms in the equations can be neglected. So, in metals the dissipative term dominates, for example, the strength of electrostatic field must be E = 0 inside metal: the nonzero strength would lead to a current, meanwhile the propagation of the current is associated with energy dissipation according to Joule-Lenz law Q = j2/σ = σE2 and, therefore, it cannot be supported in an equilibrium state by itself [14]. For not very large frequencies of electromagnetic waves (less than plasma frequency) a condition σ 1 @A ! "(!) is satisfied for good metals, then the conduction current −σ c @t gives main contribution in electromagnetic " @E response that gives the skin-effect, at the same time the displacement current 4π @t can be neglected [14, 15]. On the contrary, in dielectrics the dissipative processes can be neglected in the first approximation (since conductivity is zero); hence, Lorentz covariant electrodynamics remains only. In bad conductors (semiconductors, weak electrolytes, weakly ionized plasma) both conduction and displacement currents can be important. Analogously to dielectrics and metals electrodynamics of superconductors must be composed of both Lorentz covariant part and dissipative part (due to friction of normal electrons, damping of collective excitations and breaking of Cooper pairs). Which term is dominant is determined with physical conditions and processes. Unlike the normal metal phase, in the SC phase the conductivity τph 2 essentially depends on temperature due to density of normal electrons nn: σ(T ) = m e nn(T ) ) σ(T ! 0) ! 0 (here τph is the mean free time of electron caused by electron-phonon interaction), that is at low temperatures the contribution of the dissipative part decreases. It should be noted that superconductivity is a thermodynamic effect and is not electrodynamic one (superconductor is not ideal conductor). Hence equations for electromagnetic field in superconductor are result of variation of some functional of action (or free energy functional for equilibrium δS case): = 0. The action S [Ψ(r; t);Aµ(r; t)] must be Lorentz invariant in order to ensure Lorentz covariant δAµ electrodynamics of superfluid component. The dissipative terms are introduced in the equations of motion by means of Rayleigh dissipation function. At the same time, the Lorentz invariance of the action should have consequences for dynamics of the scalar field Ψ: the collective pseudo-relativistic excitations (Higgs mode and Goldstone mode) occur. In [2, 5] the general equations for the dynamic behavior of dirty superconductors in GL regime jTc − T j << Tc are derived from microscopic theory.