Defect-Mediated Skyrmion Nucleation

• Defect-mediated skyrmion nucleation is the process where material imperfections reduce energy barriers to form stable, topologically nontrivial spin textures.
• It exploits modifications in exchange stiffness, Dzyaloshinskii–Moriya interactions, and anisotropy to enable precise control in spintronic and data storage applications.
• Experimental protocols using TEM, STM, and nanofabrication demonstrate that engineered defects significantly lower nucleation energy, enhancing skyrmion stability and device performance.

Defect-mediated skyrmion nucleation refers to the process by which local structural or atomic imperfections facilitate the formation of topologically nontrivial spin textures—skyrmions—by reducing the nucleation energy barriers, locally distorting the exchange and Dzyaloshinskii–Moriya interaction (DMI) landscape, or providing pinning sites for topological transition pathways. This mechanism is relevant across a wide class of chiral magnets and synthetic multilayers and is central to the controlled generation, stability, and manipulation of skyrmions for applications in information storage and spintronic devices.

1. Theoretical Framework and Free-energy Landscapes

Skyrmion nucleation is governed by the interplay among exchange stiffness, DMI, uniaxial anisotropy, and Zeeman coupling (as well as demagnetizing fields in thin-film geometries). For bulk cubic chiral magnets such as β-Mn–CoZn–Mn, the micromagnetic continuum free-energy functional is given by

$$F[\mathbf{m}] \;=\; \int d^3r \left\{ A\,|\nabla \mathbf{m}|^2 + D\,\mathbf{m}\cdot(\nabla\times\mathbf{m}) - \mu_0 M_s\,\mathbf{H}\cdot\mathbf{m} \right\},$$

where $\mathbf{m}(\mathbf{r})$ is the unit vector of the local magnetization, $A$ the exchange stiffness ($\sim$10 pJ/m), $D$ the DMI constant ($\sim$0.5 mJ/m²), $M_s$ the saturation magnetization, and $\mathbf{H}$ the applied field. In multilayers with interfacial DMI and perpendicular magnetic anisotropy (PMA), an additional term $-K_{\mathrm{eff}} m_z^2$ is important (Kim et al., 2019, Yin et al., 2020).

Defects, which may take the form of atomic substitutions, interfacial roughness, local strain, nonmagnetic vacancies, or nanofabricated pinning structures, modify the local parameters $A$, $D$, $K_{\mathrm{eff}}$, or $M_s$, thereby lowering the nucleation barriers for nontrivial spin structures. The activation energies, critical radii for nucleation, and attempt frequencies can be modeled via minimum-energy path (MEP) methods such as the string or geodesic nudged elastic band (GNEB) approaches (Kim et al., 2019, Uzdin et al., 2017).

2. Defect Types and Microscopic Mechanisms

Structural Defect Classes:

• Disclination Dipoles (e.g., 5–7 defects): In skyrmion clusters evolving from the conical phase, Lorentz-TEM reveals local coordination defects where one skyrmion is surrounded by five nearest neighbors (pentagon), adjacent to one with seven (heptagon). These wedge defects act as energetically favorable sites for skyrmion "self-splitting," analogous to mitotic division, facilitating rapid defect-mediated nucleation at lattice imperfections (Kim et al., 2019).
• Nonmagnetic Vacancies: In atomistic models, missing magnetic atoms at or near the skyrmion boundary result in the removal of exchange and anisotropy energy costs precisely where the nucleation core forms, reducing the saddle-point energy by 5–7 meV per atom. A three-atom defect can decrease the barrier by up to 70% compared to defect-free systems (Uzdin et al., 2017).
• Substitutional Doping / Strain Centers: Substitutional defects (e.g., Cl for I in monolayer CrI₃) produce local lattice distortions, breaking inversion symmetry and inducing or locally enhancing DMI, which stabilizes spin-canted bubbles or, at sufficient strength, nucleates true skyrmions (Beck et al., 2021).
• Interface-induced Nanocrystallites: In ultrathin films (CoFeB/Ta), atomic-scale defects such as small nanocrystallites or vertical pinhole-like features, controllable via inverted growth protocols, create locally reduced PMA and inhomogeneous DMI, correlating directly with skyrmion nucleation centers (Yin et al., 2020).
• Lithographically Defined Pinning Sites: Engineered notches or high-anisotropy barriers at nanowire edges prevent the escape or expulsion of intermediate topological defects (e.g., vertical Bloch lines), enabling deterministic, single-skyrmion writing via localized current injection (Dutta et al., 2018).

Mechanistically, defects act by concentrating the magnetization twist, locally reducing effective energy barriers, and providing sites for nucleation, division, or stabilization of topological spin textures.

3. Energetics and Quantitative Scaling Relations

Defect-mediated nucleation is characterized by a significant reduction of the required energy to transition from a homogeneous or weakly modulated magnetic background to a skyrmion or skyrmion cluster state.

• Classical nucleation picture: The critical nucleus work is given by $W_c = \gamma\,2\pi R_c - \Delta f\,\pi R_c^2$ with $\gamma$ the line tension and $\Delta f$ the free-energy density difference. The critical radius $R_c = \gamma/\Delta f$.
• Activation energy ($\Delta E$): In β-Mn–CoZn–Mn, the nucleation barrier for an isolated skyrmion is $\Delta E_{\text{iso}}\approx 1.2\times10^{-19}$ J; at a 5–7 defect, this is reduced to $\Delta E_{\text{defect}}\approx 3.0\times10^{-20}$ J ($\sim$25% of the isolated value). In atomistic tracks, a single missing atom reduces the barrier by 15 meV, while a 3-atom cluster yields a $\sim$47 meV lower barrier (Kim et al., 2019, Uzdin et al., 2017).
• Phenomenological scaling: For nonmagnetic defects,

$$\Delta E_{\text{nuc}}(D,N_d,r_d) = \Delta E_0(D) - \gamma(N_d)\exp\left[-\frac{r_d^2}{2\xi^2}\right],$$

where $D$ is the track width, $N_d$ the number of defect atoms, $r_d$ the nucleus-defect distance, and $\xi$ the healing length ($\sim$skyrmion radius) (Uzdin et al., 2017).

• Density scaling with defect concentration: In CoFeB/Ta films,

$$K_{\mathrm{eff}} = K_0\,\left(1 + \epsilon\,t_d/t_f\right)$$

links effective anisotropy to defect fraction $\epsilon$. The skyrmion density exhibits a linear dependence $$\rho_{\text{sk}}(H_r) \approx C\epsilon\,[H_{\text{sat}} - |H_r|]$$ (Yin et al., 2020).

• Attempt frequencies and nucleation rates: Arrhenius-type rates with $\Gamma_0 \sim 10^8$–$10^{10}\,\text{s}^{-1}$, nucleation times $t_{\text{inc}}\sim 1.8$ s followed by rapid growth $\tau_{\text{grow}}\sim 0.1$ s for conical-phase skyrmion crystals (Kim et al., 2019).

Quantitative performance metrics—such as deterministic nucleation time ($\sim$0.2–1 ns), energy cost (a few fJ per skyrmion), and reliability under material inhomogeneity (>90% up to 30% disorder)—define practical limits for device engineering (Dutta et al., 2018).

4. Experimental Protocols and Observational Signatures

Defect-mediated skyrmion nucleation has been verified through diverse experimental approaches:

• In-situ Lorentz TEM: Reveals direct formation pathways from precursor “embryos” to full skyrmions at defect sites, as well as mitosis-like division events within skyrmion clusters (Kim et al., 2019).
• First-order reversal curve (FORC) measurements: In perpendicularly magnetized ultrathin films, FORC protocols select field intervals that expose or trap skyrmions at hard-center defect sites, permitting room-temperature control of skyrmion populations up to $2\times 10^2$ mm$^{-2}$ (Yin et al., 2020).
• Spin-polarized scanning tunneling microscopy (STM): Preferential nucleation of skyrmions at surface defects in Pd/Fe/Ir(111), with activation energies quantitatively matched to theoretical predictions ($\sim20$–$60$ meV) (Uzdin et al., 2017).
• Atomistic engineering of defects: Single-atom substitution in 2D magnets such as CrI₃ can nucleate spin-canted bubbles; implantation of 3d/4d elements in PdFe/Ir(111) bilayers enables deterministic seeding, pinning, or deletion of skyrmions with programmable magnetization (Beck et al., 2021, Fernandes et al., 2023).
• Nanofabrication: Patterned notches or high-PMA edge regions in nanowire geometries enables fully deterministic, single-skyrmion writing via localized current pulses (Dutta et al., 2018).

Empirical scaling laws and field-mapping protocols systematically correlate defect distributions to skyrmion densities and nucleation thresholds.

5. Universal Patterns, Scaling Laws, and Materials Design Principles

Recent ab-initio studies reveal universal patterns in skyrmion magnetizations, spin, and orbital moments as functions of the defect species and positions:

• Emergent field descriptors: Skyrmion energies and magnetization patterns are governed by the three-spin scalar chirality $$\chi = (\mathbf{S}_i \times \mathbf{S}_j)\cdot\mathbf{S}_k$$ and the two-spin vector chirality $\kappa = \mathbf{S}_i \times \mathbf{S}_j$. Defect-induced variations in local $\chi$ or $\kappa$ modulate nucleation thresholds and stability (Fernandes et al., 2023).
• Defect species dependence: Implantation of late 3d elements (Fe, Co) near the skyrmion core enhances spin and chiral orbital magnetization, while early elements (Ti, Cr, Mo) suppress or delete skyrmions. The dependence is “W-shaped” (3d series) or “V-shaped” (4d) in defect atomic number.
• Design guidelines: Control of skyrmion stability and read-out contrast can be systematically engineered by choosing defect type, location (within 0.5 nm of the core for maximal effect), and density (≤1% monolayer to prevent collective pinning) (Fernandes et al., 2023).
• Ab-initio scaling: For PdFe/Ir(111) with single-atom defects,

| Defect Z | $S_\text{imp}$ ($\mu_B$) | $J_{\mathrm{Fe-imp}}/J_{\mathrm{Fe-Fe}}$ | $\Delta M^\text{spin}_\text{sk}$ ($\mu_B$) | $\Delta M^\text{orb}_{1\sigma}$ ($\mu_B$) | $\Delta M^\text{orb}_{3-\text{spin}}$ (%) | |:--------:|:-----------------------:|:---------------------------------------:|:------------------------------------------:|:------------------------------------------:|:------------------------------------------:| | Ti | –0.8 | 0.2 | –25 | –0.8 | –5 | | Cr | 2.1 | 0.5 | –15 | –0.3 | –2 | | Fe | 3.2 | 1.0 | +20 | +0.9 | +6 | | Co | 2.8 | 1.1 | +22 | +0.7 | +5 | | Mo | –0.2 | –0.1 | +15 | –0.2 | –3 |

(Fernandes et al., 2023)

• Implementation pseudocode: Select target magnetization or nucleation threshold, choose defect type/position, implant, and apply just-above-threshold field or current for deterministic nucleation.

6. Comparison with Classical Nucleation and Broader Context

Skyrmion nucleation at defects exhibits both analogies and departures from classical crystallization theory:

• Analogies: Defect-assisted lower barriers yield classical-looking nucleation bursts, with defined incubation times and embryo formation. Cluster growth is governed by effective potentials exhibiting both attraction and repulsion, promoting densely packed hexagonal lattices analogous to atomic crystal growth.
• Critical distinction: Skyrmion growth is non-conservative—unlike atoms, skyrmions can “divide” (mitosis-like self-splitting at 5–7 defects) and need not rely on particle diffusion. The "defect-mediated mitosis" process provides a fundamentally new, non-conservative channel for internal generation of topological charge without any atomic analogue (Kim et al., 2019).

Defect-mediated nucleation is therefore a unifying mechanism enabling both rapid, low-barrier skyrmion creation at deterministic locations and tunable control of skyrmion populations and properties in real materials, bridging experimental, theoretical, and technological domains.