V. Plastovets,A. S. Mel'nikov,A. I. Buzdin
Abstract:We study the coherent dynamic interaction of a time-dependent spin-splitting field with the homogeneous superconducting order parameter $\Delta(t)$ mediated by spin-orbit coupling using the time-dependent Bogoliubov-de Gennes theory. In the first part of the work we show that linear response of the superconductor is strongly affected by the Zeeman field and spin-flip processes, giving rise to multiple resonant frequencies of the superconducting Higgs modes. These modes can be excited either by a quench, or by an additional non-stationary component of the spin-splitting field, which couples linearly to the Higgs modes. In the second part, we analyze the nonadiabatic dynamics of quasiparticle states arising from the intersection of spectral branches from different spin subbands, which can be provoked by a linearly growing Zeeman field. We provide insights into the dependence of the order parameter $\Delta(t)$ on this field and interference effects caused by tunneling of states at the avoided crossing points. We also show that since the nonadiabatic tunneling is related to spin-flip processes, the quasiparticle gas experiences a dynamic magnetization that contributes to its spin susceptibility.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the collisionless dynamics of the superconducting energy gap caused by the spin - splitting field in superconductors. Specifically, the research focuses on the following aspects:
1. **Linear response analysis**: The paper first explores the time evolution of a small perturbation $\delta\Delta(t)$ of the superconducting energy gap in the presence of a static spin - splitting field $h_0$. Through linear response theory, the authors find that the superconducting Higgs mode can be excited by the spin - splitting field, and these modes have multiple resonance frequencies, including the standard $2\Delta_0$ and two additional frequencies $2(\Delta_0\pm h_0)$.
2. **Non - adiabatic dynamics**: The paper also analyzes the non - adiabatic dynamics of quasiparticle states under the action of a time - dependent spin - splitting field $h(t)$. In particular, under a linearly increasing spin - splitting field $h(t)=\gamma t$, when the field strength exceeds the superconducting energy gap $\Delta_0$, the quasiparticle spectral branches of different spin sub - bands will cross. The non - adiabatic transition at this crossing point is equivalent to a dynamic spin - flip process, which has a significant impact on the superconducting order parameter $\Delta(t)$.
3. **The role of spin - orbit coupling**: Spin - orbit coupling (SOC) plays a crucial role in the interaction between different spin sub - bands. The paper discusses in detail how spin - orbit coupling affects the dynamics of the superconducting energy gap, especially in the presence of a spin - splitting field.
4. **Dynamic magnetization**: Since non - adiabatic tunneling is related to the spin - flip process, the quasiparticle gas undergoes a dynamic magnetization, which contributes to its spin susceptibility.
### Main conclusions
- **Higgs mode with multiple resonance frequencies**: In the presence of a static spin - splitting field $h_0$, the superconducting Higgs mode can be excited and has multiple resonance frequencies, including $2\Delta_0$ and $2(\Delta_0\pm h_0)$.
- **Non - adiabatic transition**: Under a linearly increasing spin - splitting field $h(t)=\gamma t$, when the field strength exceeds the superconducting energy gap $\Delta_0$, the quasiparticle spectral branches of different spin sub - bands will cross, leading to non - adiabatic transitions and further affecting the superconducting order parameter $\Delta(t)$.
- **Dynamic magnetization**: The non - adiabatic tunneling process leads to the dynamic magnetization of the quasiparticle gas, which contributes to the spin susceptibility of the system.
### Formula summary
- **Self - consistent equation of the superconducting order parameter**:
\[
\Delta(t)=-\lambda\sum_{i.c.}\langle\psi_k^\dagger(t)\tau_\Delta\psi_k(t)\rangle
\]
- **TDBdG equations**:
\[
i\frac{\partial}{\partial t}\psi_k(t)=\left[H_0 + V(t)\right]\psi_k(t)
\]
where,
\[
H_0=\begin{pmatrix}
\xi_k\sigma_0 - h_0\sigma_z&i\sigma_y\Delta_0\\
- i\sigma_y\Delta_0&-\xi_k\sigma_0 + h_0\sigma_z
\end{pmatrix}
\]
\[
V(t)=\begin{pmatrix}
-\delta h(t)\sigma_z&i\sigma_y\delta\Delta(t)\\
- i\sigma_y\delta\Delta(t)&\delta h(t)\sigma_z
\end{pmatrix}
\]
- **Linear response function**:
\[
\chi_{\Del