🔬 physics

Inverse-Designed Annular Pinning of Giant-Vortex Transitions in Mesoscopic Superconducting Cylinders

This study demonstrates that inverse-designed annular pinning landscapes can effectively lower the magnetic field threshold for the transition from a single-vortex to a giant-vortex state in mesoscopic superconducting cylinders by approximately 21%, offering a symmetry-preserving method to tune discrete vortex-state transitions.

Original authors: Dwi Sabda Budi Prasetya, Pekik Nurwantoro², Muhammad Farchani Rosyid²

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Dwi Sabda Budi Prasetya, Pekik Nurwantoro², Muhammad Farchani Rosyid²

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world where electricity flows without any resistance at all, a phenomenon known as superconductivity. It's like a magic highway where cars (electrons) zoom forever without ever hitting a pothole or losing speed. But this magic only happens when things get incredibly cold. In these super-highways, there are tiny whirlpools of magnetic force called "vortices." Think of them as little tornadoes spinning inside the material. Usually, these tornadoes are messy and chaotic, but in very small, perfectly shaped cylinders of superconductor, they behave like disciplined soldiers. They don't just appear randomly; they line up in specific, numbered groups. One tornado, two tornadoes, three tornadoes. Scientists call these groups "vortex states."

The big question for researchers is: How do we control when these soldiers change their formation? If you have a tiny cylinder, you can't just wait for the magnetic field to get strong enough to force a second tornado to join the first one. You want to be the architect, deciding exactly when that switch happens. This paper explores a clever trick to force that change to happen earlier, using a concept called "inverse design." Instead of guessing what shape to make, they started with the goal (make the switch happen sooner) and worked backward to find the perfect shape to get there. It's like trying to figure out exactly how to bend a slide so that a toy car reaches the bottom in exactly three seconds, rather than just building a slide and seeing how fast it goes.

The Story of the Magic Ring

In this study, the scientists looked at a tiny, superconducting cylinder. In a normal, perfect cylinder, the system likes to stay in a "single-vortex" state (one tornado) until the magnetic field gets quite strong. Only then does it switch to a "giant-vortex" state (two tornadoes sharing the same spot). The researchers wanted to know: Can we build a tiny, artificial defect inside the cylinder to make that switch happen much sooner?

They decided to try a very specific shape: a ring. Imagine a donut-shaped weakness in the superconductor, placed not in the very center, but in a circle around the middle. They used a powerful computer simulation to test thousands of different rings, changing their strength, their distance from the center, and how thick they were. They were looking for the "Goldilocks" ring—not too weak, not too strong, and not in the wrong spot.

The Big Discovery

The team found a winning design. By creating a specific ring-shaped weakness with a strength of 1.20, a radius of 1.90 ξT, and a width of 0.36 ξT, they managed to change the rules of the game.

In a normal, perfect cylinder, the switch from one vortex to two happens at a magnetic field strength of 0.485411 Hc2. However, with their specially designed ring, that switch happened much earlier, at just 0.383876 Hc2. That is a reduction of 20.92%.

To prove this wasn't just a fluke, they tested a specific magnetic field strength of 0.434644 Hc2. At this exact level, the normal cylinder was still stubbornly holding onto its single vortex. But the cylinder with the magic ring had already happily switched to the two-vortex state. It was like the ring-pinned cylinder was a sprinter who had already crossed the finish line while the normal cylinder was still tying its shoes.

Why the Ring Works (and Why Other Shapes Don't)

You might think, "Why not just put a hole in the very center?" The scientists tested this idea, along with other shapes like a big, blurry patch of weakness or a random spot off to the side. They found that a hole in the center didn't help much because the first vortex already lives there; it doesn't need extra help. A random spot off to the side broke the symmetry and didn't work as well.

The magic of the ring is that it sits in the "recovery zone." Think of the single vortex as a small campfire in the center. The giant vortex (two tornadoes) needs a bit more space to breathe and spread out. The ring is placed exactly where the giant vortex needs to expand, making it easier for the second tornado to join the party without fighting against the material. It's a targeted move, not a general weakening of the whole system.

How Sure Are They?

The results come from detailed computer simulations, not a physical experiment in a lab yet. The authors ran their numbers through a rigorous "boundary-value problem" solver, which is like a super-precise calculator for physics equations. They also ran a secondary check using a 2D simulation to make sure their ring idea wasn't just an illusion caused by forcing the problem to look perfectly round. Even in this looser 2D test, the ring design still showed the smallest "energy penalty" for the two-vortex state compared to all other shapes they tested.

While they haven't built a physical cylinder with this ring yet, the math is very solid. They have shown that if you could manufacture this specific ring pattern (perhaps by zapping the material with ions or changing its thickness in a circle), you could program the superconductor to switch states at a lower magnetic field. It's a blueprint for engineering the future of superconducting devices, proving that with the right shape, you can tell these tiny magnetic tornadoes exactly when to change their dance formation.

Technical Summary: Inverse-Designed Annular Pinning of Giant-Vortex Transitions in Mesoscopic Superconducting Cylinders

Problem Statement
The study addresses the equilibrium vorticity selection in mesoscopic type-II superconducting cylinders, where the radius is comparable to the coherence length ( $\xi_T$ ). In such confined geometries, the system does not exhibit a continuous vortex density but rather discrete transitions between giant-vortex states characterized by integer vorticity $L$ (e.g., $L=0, 1, 2, \dots$ ). The specific problem is to determine whether an artificial, radially symmetric pinning landscape can be designed to quantitatively lower the critical magnetic field ( $H_{1\to2}$ ) at which the system transitions from a single-vortex state ( $L=1$ ) to a giant-vortex state ( $L=2$ ). Unlike prior works that focus on vortex motion, critical currents, or matching effects in antidot arrays, this work treats the transition field itself as an inverse-design objective.

Methodology
The research employs a two-stage computational workflow grounded in the axisymmetric Ginzburg–Landau (GL) theory:

Inverse Design via Variational Screening:
- The authors introduce an artificial annular pinning landscape defined by a Gaussian suppression of the linear GL coefficient: $V_{pin}(r) = V_0 \exp[-(r-r_p)^2 / (2\sigma^2)]$ .
- The design variables are the pinning strength ( $V_0$ ), ring radius ( $r_p$ ), and width ( $\sigma$ ).
- A fast variational model using a trial function $f_L(r) = A r^L \exp(-qr^2/2)$ is used to rapidly scan the parameter space and identify candidates that minimize the $H_{1\to2}$ transition field.
Full Boundary-Value Problem (BVP) Validation:
- The top candidate from the variational screen is re-evaluated by solving the full radial GL differential equation with Neumann boundary conditions.
- Equilibrium states are determined by calculating the Gibbs free energy ( $G_L$ ) for different vorticities ( $L=0, 1, 2, 3$ ) and identifying branch crossings where $G_L = G_{L+1}$ .
2D Pseudo-Time-Dependent Ginzburg–Landau (TDGL) Benchmark:
- To verify that the results are not artifacts of enforced axial symmetry, a fixed-field 2D relaxation benchmark is performed.
- This benchmark compares the optimized annular profile against four controls: a homogeneous cylinder, a central defect, a uniform suppression, and an off-axis defect.
- The metric used is the energy penalty $\Delta E = E(L=2) - E(L=1)$ at a fixed field where the homogeneous system is still in the $L=1$ state. Phase-winding diagnostics are used to confirm the topological nature of the states.

Key Results

Optimized Parameters: The inverse design identified an optimal annular profile with $V_0 = 1.20$ , $r_p = 1.90 \xi_T$ , and $\sigma = 0.36 \xi_T$ .
Transition Field Reduction: The optimized profile reduces the $L=1 \to L=2$ transition field from $H_{1\to2} = 0.485411 H_{c2}$ (homogeneous) to $H_{1\to2} = 0.383876 H_{c2}$ . This represents a 20.92% reduction.
Stability Window Preservation: The design shifts the transition to lower fields without destroying the stability of the $L=2$ state; the stability window ( $H_{2\to3} - H_{1\to2}$ ) remains comparable (0.166215 $H_{c2}$ vs. 0.162549 $H_{c2}$ ).
Same-Field State Reversal: At an applied field of $H_e/H_{c2} = 0.434644$ $H_{e} / H_{c 2} = 0.434644$ (between the two transition thresholds):
- The homogeneous cylinder remains in the $L=1$ ground state.
- The annularly pinned cylinder is stable in the $L=2$ ground state.
Benchmark Superiority: In the 2D relaxation benchmark at the same field, the annular ring produced the smallest $L=2$ energy penalty ( $\Delta E = 0.057359$ ) among all tested landscapes, outperforming the off-axis defect ( $\Delta E = 0.064084$ ) and significantly outperforming the homogeneous and central-defect controls.

Significance and Claims
The paper claims to demonstrate a symmetry-preserving route for tuning discrete vortex-state transitions in confined superconductors. The significance lies in the shift from predictive modeling (calculating states for a given sample) to inverse design (engineering a sample to select a specific state at a target field).

Mechanism: The annular suppression is positioned in the radial recovery zone of the order parameter. This location is energetically more effective for the broader $L=2$ core than for the $L=1$ core, thereby selectively lowering the condensation energy penalty for the giant-vortex state.
Distinction from Prior Work: The authors distinguish their work from studies on vortex clustering, critical current enhancement, or antidot matching effects. The novelty is the explicit optimization of a radial GL-coefficient landscape to program a specific thermodynamic branch crossing ( $H_{1\to2}$ ).
Limitations: The authors explicitly state that the study is a computational proof-of-concept within the static GL framework. The radial BVP enforces axisymmetry, and the 2D benchmark is a fixed-field relaxation test rather than a full dynamic simulation with self-consistent magnetic screening. Consequently, the results do not exhaustively address separated multivortex states, thermal fluctuations, or current-driven entry paths.

In conclusion, the work establishes that an inverse-designed annular pinning landscape can quantitatively program the onset of giant-vortex states, providing a design rule for where to weaken superconductivity to achieve specific topological transitions in mesoscopic cylinders.

Technical Summary: Inverse-Designed Annular Pinning of Giant-Vortex Transitions in Mesoscopic Superconducting Cylinders

More like this