Relaxation of nonlinear oscillations in BCS superconductivity

The diagonal case of the $sl(2)$ Richardson-Gaudin quantum pairing model \cite{Richardson1,Richardson2,Richardson3,Richardson4,Richardson5,Richardson6,G audin76} is known to be solvable as an Abel-Jacobi inversion problem \cite{SOV,Kuznetzov,Kuz1,Kuz2,Kuz3,Kuz4,Kuz5,YAKE04}. This is an isospectral (stationary) solution to a more general integrable hierarchy, in which the full time evolution can be written as isomonodromic deformations. Physically, the more general solution is appropriate when the single-particle electronic spectrum is subject to external perturbations. The asymptotic behavior of the nonlinear oscillations in the case of elliptic solutions is derived.