Collective excitations and divergent spin currents in non-centrosymmetric superconductors

Basic Information

• Paper ID: 2511.02032
• Title: Collective excitations and divergent spin currents in non-centrosymmetric superconductors
• Authors: Markus Lysne, Philipp Werner, Nikolaj Bittner (University of Fribourg, Switzerland)
• Classification: cond-mat.supr-con (Superconductivity), cond-mat.str-el (Strongly Correlated Electron Systems)
• Publication Date: November 5, 2025
• Paper Link: https://arxiv.org/abs/2511.02032

Abstract

This paper investigates collective modes in non-centrosymmetric superconductors (NCS) with Rashba spin-orbit coupling under laser irradiation. Utilizing the Anderson pseudospin resonance concept, the authors reveal that laser driving not only leads to resonant enhancement of the known third harmonic generation (THG) response in charge currents, but also induces resonant enhancement of second harmonic generation (SHG) in spin currents. The authors propose a theoretical framework without interband transitions to explain this phenomenon and verify it through numerical simulations including interband effects, thereby elucidating the role of collective modes in the long-time dynamics of superconductors.

Research Background and Motivation

Research Questions

This study addresses the detection and characterization of collective excitation modes in non-centrosymmetric superconductors (NCS), specifically:

1. How to observe Higgs and Leggett modes in NCS with spin-orbit coupling
2. How nonlinear spin current responses reflect these collective modes
3. The correlation mechanism between spin dynamics under optical driving and superconducting order parameter oscillations

Importance of the Problem

1. Fundamental Physics: Non-centrosymmetric superconductors allow mixing of singlet and triplet pairing, serving as an important platform for studying unconventional superconductivity
2. Experimental Challenge: Detailed characterization of superconducting gaps through traditional probes remains difficult, necessitating new nonequilibrium detection methods
3. Technological Potential: Understanding spin current dynamics is crucial for spintronics applications
4. Theoretical Gap: While collective modes in conventional superconductors have been extensively studied, related phenomena in NCS remain relatively unexplored

Limitations of Existing Approaches

1. Traditional third harmonic generation (THG) studies focus primarily on charge current responses
2. Insufficient investigation of nonlinear responses in spin currents
3. Incomplete theoretical descriptions of Higgs and Leggett modes in NCS
4. Lack of systematic analytical theoretical framework for handling mixed pairing scenarios

Research Motivation

The authors are driven by:

• Development of ultrafast terahertz spectroscopy enabling direct detection of collective modes
• Rich spin dynamics exhibited by spin-orbit coupled systems
• Need to establish a unified theoretical framework connecting optical driving, collective modes, and spin currents

Core Contributions

1. Theoretical Framework: Establishes an analytical theory based on the Anderson pseudospin formalism that describes collective modes in NCS without requiring interband transitions
2. Novel Physical Phenomena: First reveals divergent enhancement of second harmonic generation (SHG) in spin currents under Anderson pseudospin resonance conditions, analogous to third harmonic generation in charge currents
3. Linearized Equations: Derives closed linearized equations of motion for order parameter components Eq. (33) and susceptibility expressions Eq. (34) for predicting collective mode frequencies
4. Numerical Verification: Validates analytical theory through time-dependent mean-field simulations and investigates interband coupling effects and long-time dynamics
5. Mechanism Clarification: Systematically decomposes contributions from Higgs, phase, and charge density fluctuation modes to spin current response, finding charge density fluctuations as the dominant contribution
6. Universality Proof: Demonstrates that photoinduced spin currents are universal features of spin-orbit coupled systems, not limited to specific superconducting configurations

Detailed Methods

Task Definition

Study the response of non-centrosymmetric superconductors driven by an external periodic electric field $E(t) = -\partial_t A(t)$, where:

• Input: Laser pulse parameters (frequency Ω, field strength $E_0$, number of cycles M), system parameters (Rashba coupling α, pairing strength $V_0$, singlet/triplet ratios $e_s, e_{tr}$)
• Output: Spin current $J_{yx}(t)$, order parameter dynamics $\Delta_s(t), \Delta_{tr}(t)$, collective mode spectrum
• Constraints: Weak-field approximation (linearized theory) or self-consistent mean-field (numerical simulation)

Model Architecture

1. Hamiltonian Construction

Total Hamiltonian $\hat{H} = \hat{H}_0 + \hat{H}_I$, where:

Non-interacting part: $\hat{H}_0 = \sum_{k,s_1s_2} [(\epsilon^0_k - \mu_c)\tau_0 + \mathbf{d}_k \cdot \boldsymbol{\tau}]_{s_1,s_2} \hat{c}^\dagger_{ks_1}\hat{c}_{ks_2}$

where $\mathbf{d}_k = \alpha(-\sin k_y, \sin k_x, 0)$ is the Rashba spin-orbit coupling term. Through unitary transformation $U_k$ it is diagonalized to: $\hat{H}_0 = \sum_{k\mu} \epsilon_{k\mu} \hat{a}^\dagger_{k\mu}\hat{a}_{k\mu}$ with band dispersion $\epsilon_{k\mu} = \epsilon^0_k + \mu|\mathbf{d}_k| - \mu_c$, where $\mu = \pm$ labels the two bands.

Interaction part Eq. (4): $V_{k,k';s_1s_2s'_2s'_1} = \frac{V_0}{4}\{e_s\phi_{s_1,s_2}\phi^\dagger_{s'_2,s'_1} + e_{tr}(\mathbf{g}_k \cdot \boldsymbol{\phi})_{s_1,s_2}(\mathbf{g}_{k'} \cdot \boldsymbol{\phi})^\dagger_{s'_2,s'_1} + \text{mixed terms}\}$

where $e_s^2 + e_{tr}^2 + e_m^2 = 1$ characterizes singlet, triplet, and mixed pairing strengths.

2. Anderson Pseudospin Representation

Introduce pseudospin operators Eq. (13): $\hat{\sigma}^i_{k\mu} = \frac{1}{2}(\hat{a}^\dagger_{k,\mu}, \hat{a}_{-k,\mu}) \cdot \tau^i \cdot \begin{pmatrix}\hat{a}_{k,\mu} \\ \hat{a}^\dagger_{-k,\mu}\end{pmatrix}$

The equilibrium Hamiltonian is written as: $\hat{H}^{eq}_{APS} = \sum_{k\mu} \boldsymbol{\sigma}_{k\mu} \cdot \mathbf{b}_{k\mu}$

where the pseudomagnetic field $\mathbf{b}_{k\mu} = (\Delta'_{k\mu}, -\Delta''_{k\mu}, \frac{1}{2}(\epsilon_{k\mu} + \epsilon_{-k\mu}))$.

3. Equations of Motion

Pseudospins satisfy the Bloch equation Eq. (16): $\partial_t \boldsymbol{\sigma}_{k\mu}(t) = 2\mathbf{b}_{k\mu}(t) \times \boldsymbol{\sigma}_{k\mu}(t)$

where the electric field is introduced through Peierls substitution as time-dependent terms: $b^z_{k\mu}(t) = \frac{1}{2}(\tilde{\epsilon}_{k,qA(t),\mu} + \tilde{\epsilon}_{-k,qA(t),\mu})$

Technical Innovations

1. Linearized Theory Framework

For weak fields, introduce deviations $\delta\Delta_{k\mu}(\omega), \delta\sigma_{k\mu}(\omega), \delta b^z_{k\mu}(\omega)$, deriving the response matrix Eq. (20): $\delta\sigma_{k\mu}(\omega) = \mathbf{M}_{k\mu}(\omega) \cdot \delta\mathbf{b}_{k\mu}(\omega)$

The key innovation is establishing closed equations Eq. (33): $\delta\boldsymbol{\Delta}(\omega) = (1 - \chi(\omega))^{-1}\left[\sum_{k\mu}\Gamma_{k\mu}\mathbf{M}_{k\mu}(\omega)\delta b^z_{k,\mu}(\omega)\mathbf{e}_z + \text{c.c.}\right]$

where the susceptibility matrix Eq. (34): $\chi(\omega) = \sum_{k\mu}\Gamma_{k\mu}\mathbf{M}_{k\mu}(\omega)\mathbf{F}_{k\mu} + \text{complex conjugate terms}$

Collective mode frequencies are determined by $\det(1-\chi(\omega)) = 0$.

2. Nonlinear Response of Spin Currents

Derive the second-order response of spin currents Eq. (45): $J^{A^2}_{yx}(\omega) = q\sum_{k\mu}\delta\sigma^z_{k\mu}(\omega)\left[\partial_{k_x}\epsilon^0_k\sin(\varphi_k)\tau^z_{\mu\mu} + \tau^0_{\mu\mu}\partial_{k_x}d_{k,y}\right]$

Key insight: $\delta\sigma^z_{k\mu}$ diverges as $\omega \to 2E_{k\mu}$, leading to resonant enhancement.

Further decomposition Eq. (46): $\delta\sigma^z_{k\mu}(\omega) = \frac{\tanh(\beta E_{k\mu}/2)}{E_{k\mu}(\omega^2 - 4E^2_{k,\mu})}\left[-2\epsilon_{k,\mu}\Delta'_{k\mu}\delta\Delta'_{k\mu} + i\omega\Delta'_{k\mu}\delta\Delta''_{k\mu} + 2(\Delta'_{k\mu})^2\delta b^z_{k\mu}\right]$

corresponding respectively to Higgs, phase, and charge density fluctuation contributions.

3. Distinction from Traditional Methods

• vs. Diagrammatic techniques 23: This work provides more intuitive physical pictures and analytical expressions
• vs. Pure numerical methods: Linearized theory is computationally efficient and rapidly predicts resonance frequencies
• Innovation: First application of Anderson pseudospin resonance concept to spin current SHG, establishing parallel theoretical framework with THG

Experimental Setup

Model Parameters

Band structure Eq. (50):

• Hopping integral: $\epsilon^0_k = -2t_0(\cos k_x + \cos k_y)$, $t_0 = 1$ (energy unit)
• Rashba coupling: $\alpha = 1.2$ (ensures clear Fermi surface separation)
• Chemical potential: $\mu = -1.5$
• Debye frequency: $\omega_D = 0.3$

Pairing parameters:

• Interaction strength: $V_0 = -16$
• Three pairing ratio configurations:
  1. Nearly pure singlet: $e_s = 0.99, e_{tr} = 0.141$
  2. Dominant triplet: $e_s = 0.05, e_{tr} = 0.999$
  3. Mixed state: $e_s = 0.7, e_{tr} = 0.714$

Laser parameters:

• Frequency: $\Omega = 0.1$ (matched with gap for resonance conditions)
• Field strength: $E_0 = 0.02$ (weak field ensures linear approximation validity)
• Number of cycles: $M = 9$ (sufficient to excite collective modes while avoiding excessive heating)

Numerical Methods

1. k-space grid: 300-800 k-points per dimension (300 for analytical calculations, 800 for numerical simulations)
2. Broadening parameter: $\eta = 0.01$ (artificial broadening for frequency domain calculations)
3. Temperature range: $\beta = 5$ to $\beta = 100$ (spanning the superconducting transition)
4. Time evolution: Commutator-free expansion algorithm 37 for solving time-dependent density matrix equations Eq. (49)

Comparison Methods

1. Analytical vs. Numerical: Linearized theory vs. full time-dependent mean-field
2. Intraband vs. Interband: Intraband pairing approximation only vs. including interband coupling
3. Self-consistent vs. Frozen: Self-consistent order parameter updates vs. fixed equilibrium order parameter

Evaluation Metrics

1. Collective mode frequencies: Extracted from Fourier transform peaks of order parameter oscillations
2. Spin current SHG intensity: $|J_{yx}(2\Omega)|$
3. Order parameter amplitude: $|\Delta_s(2|\Delta_-|)|, |\Delta_{tr}(2|\Delta_-|)|$
4. Leggett mode intensity: Phase difference $L(t) = \frac{1}{2\pi}[\theta_+(t) - \theta_-(t)]$

Experimental Results

Main Results

1. Identification of Collective Modes (Figures 3-4)

Time-domain behavior:

• After excitation by single-cycle broadband pulse ($\Omega = 1.8, E_0 = 0.2$), order parameter oscillations exhibit $1/\sqrt{t}$ decay
• $|\Delta_s(t)|$ oscillation amplitude ~0.4%, $|\Delta_{tr}(t)|$ ~0.03% (due to small triplet component)
• Leggett mode $L(t)$ remains observable even at $e_s = 0.99$

Frequency-domain characteristics (Figure 4):

• Main peaks at $2|\Delta_+| \approx 0.396$ and $2|\Delta_-| \approx 0.372$
• Analytical theory (gold curve) agrees well with numerical results, though not fully capturing the double-peak structure
• Triplet contribution peak is weak, requiring 13.1× magnification to be visible

2. Resonant Enhancement of Spin Currents (Figure 5)

Temperature dependence Figure 5(b):

• Sharp peak appears at $\beta \approx 11$ ($T \approx T_c$)
• Peak intensity enhanced by approximately one order of magnitude compared to normal state
• Comparison of three methods:
  • $J_{yx}$ (analytical): peak ~0.0010
  • $J^{SCMF,intra}_{yx}$ (numerical without interband): peak ~0.0007, closer to analytical
  • $J^{SCMF}_{yx}$ (numerical with interband): peak ~0.0014, approximately 2× without interband

Contribution decomposition Figure 5(c):

• Charge density fluctuation contribution (blue) dominates, ~90% of total signal
• Higgs contribution (red) ~8%
• Phase mode contribution (green) ~2%
• Far from resonance, density term contribution increases further

3. Comparison of Different Pairing Ratios

Nearly pure singlet ($e_s = 0.99$, Figure 6):

• Single resonance peak at $\beta \approx 10$
• Order parameter component amplitudes: $|\Delta_s| \sim 0.003$, $|\Delta_{tr}| \sim 0.0001$
• Spin current SHG peak value ~0.00015

Dominant triplet ($e_s = 0.05$, Figure 7):

• Resonance peak shifts to $\beta \approx 8$
• Triplet component dominates: $|\Delta_{tr}| \sim 0.004$ >> $|\Delta_s| \sim 0.0001$
• Spin current response characteristics similar but different resonance conditions

Mixed state ($e_s = 0.7$, Figure 8):

• Double resonance structure: $\beta \approx 10$ and $\beta \approx 15$
• Leggett mode relative intensity maximum (~0.0003)
• Sharper spin current peak, suggesting Leggett mode enhancement effect

4. Phase Evolution (Figure 9)

Fixed $\beta = 100$, scanning laser frequency:

• Spin current phase jumps $\pi/2$ at $\Omega \approx |\Delta_-|$
• Self-consistent calculation shows weaker enhancement (peak ~0.00015)
• Frozen order parameter calculation shows stronger peak (~0.00025)
• Confirms dominant role of charge density fluctuations

Ablation Studies

While not explicitly labeled as "ablation studies," the paper systematically investigates various factors:

1. Interband coupling effects:
   • Including interband coupling enhances spin currents by ~2×
   • Qualitative behavior remains consistent
2. Self-consistency impact:
   • Self-consistent order parameter updates reduce effective gap
   • Explains why analytical theory overestimates peak values
3. Pairing type dependence:
   • Pure singlet/triplet lacks Leggett mode
   • Mixed state shows most prominent Leggett mode
4. Broadening parameter:
   • $\eta = 0.01$ provides reasonable spectral resolution
   • Smaller η yields sharper peaks but higher computational cost

Experimental Findings

1. Universality: Spin current SHG resonant enhancement is a universal feature of spin-orbit coupled superconductors
2. Mechanism correspondence: Spin current SHG and charge current THG follow the same Anderson pseudospin resonance mechanism
3. Dominant contribution: Charge density fluctuations rather than Higgs mode are the primary cause of resonant enhancement (consistent with reference 23)
4. Leggett mode role: Leggett mode may enhance spin current response in mixed pairing
5. Decay law: Order parameter oscillations follow $1/\sqrt{t}$ decay, same as traditional s-wave superconductors 17
6. Phase characteristics: Spin current SHG exhibits $\pi/2$ phase jump similar to THG

Related Work

Theory of Non-centrosymmetric Superconductors

1. Gor'kov & Rashba 1: First prediction of singlet-triplet mixing due to spin degeneracy lifting
2. Frigeri et al. 3-6: Developed phenomenological theory of NCS, introduced Dzyaloshinskii-Moriya interaction
3. Samokhin & Mineev 7: Studied gap structure
4. Vorontsov et al. 8: Predicted surface bound states and spin currents

Nonequilibrium Superconductivity

5. Tsuji & Aoki 17: Pioneering work establishing Anderson pseudospin resonance theory explaining THG
6. Matsunaga et al. 14,40: Experimental direct observation of Higgs mode
7. Murotani et al. 18: Extended to Leggett modes in multiband superconductors
8. Cea et al. 36: Analyzed Cooper pair vs. Higgs mode contributions in THG

Spin Current Physics

9. Bercioux & Lucignano 31: Review of quantum transport in Rashba systems
10. Hamamoto et al. 32: Nonlinear spin current generation
11. Silaev et al. 25: Higgs mode-driven spin and charge currents

Relationship of This Work to Related Work

• Relative to 17,18: Extended to NCS and spin current responses
• Relative to 23: Provides more intuitive analytical framework, verifies density fluctuation dominance conclusion
• Relative to 25: Demonstrates photoinduced spin currents are universal rather than configuration-specific
• Advantages: Unified treatment of singlet-triplet mixing, establishes parallel SHG-THG theory, systematic numerical verification

Conclusions and Discussion

Main Conclusions

1. Theoretical Contribution: Establishes complete theoretical framework for NCS based on Anderson pseudospin formalism, including analytical linearized equations and numerical time-dependent mean-field methods
2. Novel Physical Phenomena: First reveals divergent enhancement of spin current second harmonic generation under Anderson pseudospin resonance conditions ($\Omega \approx |\Delta|$)
3. Mechanism Clarification: Spin current SHG follows the same resonance mechanism as charge current THG, but charge density fluctuations are the dominant contribution (~90%), with Higgs mode contribution relatively small (~8%)
4. Collective Mode Characteristics:
   • Higgs mode frequency: $\omega = 2|\Delta_\pm|$
   • Leggett mode most prominent in mixed pairing
   • Oscillations exhibit $1/\sqrt{t}$ decay
5. Experimental Significance: Spin current SHG can serve as a new tool for detecting collective modes in low-$T_c$ NCS, as it generates less heating

Limitations

1. Theoretical Assumptions:
   • Only considers intraband pairing (neglecting interband pairing may underestimate certain effects)
   • Weak-field approximation limits description of strong driving regimes
   • Neglects impurity scattering (may enhance Higgs contribution 23)
2. Model Simplifications:
   • Uses simple 2D Rashba model; actual materials are more complex
   • Does not include phonon retardation effects 41
   • Neglects dipole coupling 29,38
3. Numerical Limitations:
   • k-space discretization may affect accuracy
   • Artificial broadening parameter $\eta$ affects peak shape
   • Computational resources limit parameter space exploration
4. Experimental Feasibility:
   • Spin current detection techniques (e.g., magnetic moment precession in adjacent ferromagnets 31,35) remain challenging
   • Requires extremely low temperatures and strong terahertz fields

Future Directions

1. Dirty limit effects: Study how impurity scattering enhances Higgs contribution (analogous to reference 23 THG study in NCS)
2. Multiorbital extension: Include orbital degrees of freedom and interorbital dipole coupling
3. Competing order parameters: Study coexistence with charge density waves, spin density waves, etc.
4. Strong-field nonlinearity: Go beyond linearized approximation, explore higher-order harmonics
5. Real materials: Apply to specific NCS materials like CePt₃Si, UIr
6. Dynamical phase transitions: Study possible photoinduced phase transitions under strong driving
7. Topological effects: Combine with topological superconductor edge states

In-Depth Evaluation

Strengths

1. Theoretical Innovation:

• First systematic establishment of theoretical framework for spin current nonlinear response in NCS
• Cleverly extends Anderson pseudospin resonance concept from THG to SHG
• Derives compact closed equations Eq. (33) with good physical transparency

2. Method Completeness:

• Analytical theory and numerical simulations mutually verify, enhancing credibility
• Systematically compares different approximations (intraband/interband, self-consistent/frozen)
• Detailed mathematical derivations (Appendices A-D) ensure reproducibility

3. Deep Physical Insights:

• Reveals dominant role of charge density fluctuations, correcting potential misconceptions
• Clarifies Leggett mode enhancement in mixed pairing
• Establishes parallel theoretical structure between SHG and THG

4. Clear Result Presentation:

• Well-designed figures (Figure 1 conceptual diagram, Figures 2-9 data plots)
• Quantitative comparison of different methods (three curves in Figure 5)
• Multi-angle verification (time domain Figure 3, frequency domain Figure 4, temperature scan Figures 5-8, frequency scan Figure 9)

5. High Writing Quality:

• Clear logic: introduction → methods → results → discussion
• Complete technical details without redundancy
• Appropriate citation of 42 references, fully reflecting academic context

Weaknesses

1. Theoretical Limitations:

• Intraband pairing approximation may miss important physics (though numerical includes interband, analytical theory does not)
• Does not consider impurities and retardation effects, may underestimate Higgs contribution
• Weak-field assumption limits description of experimentally relevant strong-field cases

2. Insufficient Numerical Verification:

• Only considers one Rashba model, lacks verification with other NCS models (e.g., Dresselhaus)
• Limited parameter space exploration (only three $e_s, e_{tr}$ values)
• No comparison with experimental data (possibly due to lack of relevant experiments)

3. Physical Interpretation Could Be Deeper:

• Microscopic mechanism for charge density fluctuation dominance not fully elucidated
• Leggett mode enhancement of spin current mechanism only speculative (Figure 8)
• Topological significance of phase jump not discussed

4. Practical Considerations:

• Challenges in spin current detection techniques not fully discussed
• Experimental feasibility of parameter choices ($\alpha=1.2, V_0=-16$) unclear
• Lacks predictions for specific materials

5. Technical Details:

• Choice of broadening parameter $\eta=0.01$ not well justified
• k-space grid convergence tests not shown
• Time step and total evolution time selection not specified

Impact

Contribution to the Field:

1. Opens new direction: Spin currents as new avenue for detecting collective modes
2. Theoretical tools: Provides methodology directly applicable to other NCS
3. Concept generalization: Successful extension of Anderson pseudospin resonance from THG to SHG

Practical Value:

1. Experimental guidance: Provides theoretical basis for designing experiments to detect collective modes in low-$T_c$ NCS
2. Spintronics: Understanding photoinduced spin currents meaningful for spintronic device design
3. Material screening: Theory can predict which materials exhibit strong spin current response

Reproducibility:

• High: Detailed formula derivations, explicit parameters, cited open-source algorithms
• Moderate: Numerical implementation details (e.g., time step) not fully disclosed
• Recommendation: Public code release would significantly enhance impact

Expected Citations:

• Nonequilibrium superconductor theory research will widely cite this work
• Spin current and spintronics fields may pay attention
• Experimental groups designing NCS detection experiments will reference this

Applicable Scenarios

1. Theoretical Research:

• Study collective modes in other unconventional superconductors (d-wave, p-wave)
• Extension to multiband, multi-orbital systems
• Combination with topological superconductor theory

2. Experimental Design:

• Guide terahertz pump-probe experiment design
• Optimize laser parameters to maximize signal
• Select appropriate NCS materials (e.g., CePt₃Si)

3. Material Screening:

• Predict which materials exhibit strong spin current SHG
• Identify mixed pairing systems with prominent Leggett modes
• Assess spin-orbit coupling strength effects

4. Device Applications:

• Design spintronic devices based on photoinduced spin currents
• Develop ultrafast spin switches
• Explore superconducting spintronics

Inapplicable Scenarios:

• Systems dominated by strong correlation effects (this work uses mean-field theory)
• Strong-field nonequilibrium states (beyond weak-field approximation)
• Room temperature or high-temperature applications (requires superconducting state)

References

Key Theoretical Foundations:

• 1 Gor'kov & Rashba (2001): Pioneering work on mixed pairing in NCS
• 17 Tsuji & Aoki (2015): Anderson pseudospin resonance theory
• 18 Murotani et al. (2017): Leggett modes in multiband superconductors

Experimental Background:

• 14 Matsunaga et al. (2013): First direct observation of Higgs mode
• 15 Giorgianni et al. (2019): Optical control of Leggett mode

Methodology:

• 37 Alvermann & Fehske (2011): Commutator-free expansion time evolution algorithm
• 29 Schüler et al. (2021): Gauge invariance in light-matter interaction in tight-binding models

Comparative Studies:

• 23 Klein & Manske (2024): Recent THG study in NCS
• 36 Cea et al. (2016): Analysis of Cooper pair vs. Higgs contribution in THG

Overall Assessment: This is a high-quality theoretical physics paper making substantive contributions to nonequilibrium dynamics of non-centrosymmetric superconductors. The theoretical framework is rigorous, numerical verification comprehensive, and physical picture clear. While there is room for improvement in experimental verification and certain technical details, the pioneering application of Anderson pseudospin resonance to spin current SHG opens new directions for the field. Expected to significantly impact nonequilibrium superconductor theory and spintronics experiments. Recommended for publication in high-level journals such as Physical Review B or Physical Review Letters.