Anisotropic second-harmonic generation in superconducting nanostructures

Circuits based on superconducting nanostructures are among the most promising platforms for quantum computing. Understanding how device geometry governs nonlinear electrodynamics is crucial for implementing superconducting quantum technologies. However, to date, research has largely been limited to superconducting nanostructures with collinearly aligned static and dynamic applied magnetic fields. Here, we analyze the dynamics of Meissner currents and Abrikosov vortices in a superconducting nanocube exposed to combined static and microwave magnetic fields, extending the analysis to a more general excitation geometry. We demonstrate that, in a noncollinear configuration,the magnetization component parallel to the static field develops a dominant second-harmonic response under the microwave driving. This effect is strongly enhanced when Meissner currents saturate at static fields just below the thresholds for successive vortex nucleation. By numerically solving the time-dependent Ginzburg-Landau equations, we show that the response originates from Meissner-current saturation combined with the nonlinear oscillations of normal-phase indentations, yielding an anisotropic second-harmonic signal that is directionally separated from, and not overshadowed by, the first-harmonic component of the dynamic magnetization. These findings are relevant for superconducting devices that require controllable high-frequency nonlinearity.

💡 Research Summary

This paper investigates the nonlinear electromagnetic response of a superconducting nanocube when subjected to combined static (DC) and microwave (AC) magnetic fields that are not collinear. Using the time‑dependent Ginzburg‑Landau (TDGL) formalism, the authors numerically solve the coupled equations for the order parameter ψ and the vector potential A in a three‑dimensional type‑II superconductor with dimensions a = h = 250 nm, Ginzburg‑Landau parameter κ = 3, London penetration depth λ = 60 nm, and coherence length ξ ≈ 20 nm. The static field B_DC is applied along the z‑axis, while the AC field B_AC(t) = B_AC cos(ω_B t) is applied either parallel (z‑direction) or perpendicular (y‑direction) to B_DC after the system has reached a steady state under B_DC.

The static magnetization ⟨M_z,DC⟩ versus B_DC exhibits sharp discontinuities at 240 mT, 400 mT, and 550 mT, each corresponding to the nucleation of four Abrikosov vortices. Just before each nucleation event, the Meissner screening current reaches a critical saturation value, producing a maximal diamagnetic response that then drops abruptly when vortices appear. This behavior sets the stage for a strongly nonlinear response when an AC field is added.

When the AC field is parallel to B_DC, the vortices undergo a breathing mode: they contract and expand symmetrically without changing shape. The resulting magnetization M_z oscillates primarily at the drive frequency ω_B, with only weak higher‑harmonic content. The symmetry of the current distribution remains largely intact, so inversion symmetry is not strongly broken and second‑harmonic generation (SHG) is modest.

In contrast, when the AC field is applied perpendicular to B_DC (along y), the vortex lines swing within the y‑z plane while normal‑phase indentations at the top and bottom surfaces oscillate in size. This motion breaks inversion symmetry and redistributes magnetic flux, leading to a pronounced second‑harmonic component in the z‑component of the magnetization (M_z ∝ cos 2ω_B t). Simultaneously, the y‑component M_y exhibits a strong first‑harmonic response, while the x‑component remains negligible. Importantly, the SHG signal in M_z is directionally separated from the dominant first‑harmonic response in M_y, allowing clean detection.

The authors identify two mechanisms responsible for the observed SHG: (1) saturation of Meissner currents near the critical field, which creates normal‑phase regions that oscillate nonlinearly under the AC drive, and (2) collective oscillations of Abrikosov vortices that modify the current pathways. The SHG amplitude peaks when B_DC is just below a vortex‑nucleation threshold (e.g., B_DC ≈ 370 mT), where the Meissner current is maximally saturated and a small AC perturbation produces large changes in the normal‑phase volume. Increasing B_AC enhances nonlinearity up to a point, after which excessive drive depins vortices and reduces the SHG strength.

The study demonstrates that the angle between DC and AC fields controls the efficiency of SHG, with a sin θ dependence (θ being the angle between the fields). This anisotropic nonlinear response is absent in collinear configurations and provides a new design knob for superconducting quantum devices. Because the second‑harmonic signal is not masked by the first harmonic and can be spatially or directionally separated, it can be exploited for high‑frequency signal processing, nonreciprocal transmission, and controlled nonlinearities in qubit control circuits.

Overall, the work extends the understanding of nonlinear superconducting electrodynamics beyond the traditional collinear paradigm, showing that geometric confinement at the nanoscale, combined with non‑collinear magnetic excitation, yields a robust, anisotropic SHG. The findings suggest practical routes to engineer controllable high‑frequency nonlinearity in superconducting nanostructures, with potential applications in quantum computing, microwave photonics, and superconducting metamaterials.