Abstract
Topologically nontrivial spin textures such as vortices, skyrmions, and monopoles are promising candidates as information carriers for future quantum information science. Their controlled manipulation including creation and annihilation remains an important challenge toward practical applications and further exploration of their emergent phenomena. Here, we report controlled evolution of the helical and skyrmion phases in thin films of multiferroic Te-doped Cu2OSeO3 as a function of material thickness, dopant, temperature, and magnetic field using in situ Lorentz phase microscopy. We report two previously unknown phenomena in chiral spin textures in multiferroic Cu2OSeO3: anisotropic scaling and channeling with a fixed-Q state. The skyrmion channeling effectively suppresses the recently reported second skyrmion phase formation at low temperature. Our study provides a viable way toward controlled manipulation of skyrmion lattices, envisaging chirality-controlled skyrmion flow circuits and enabling precise measurement of emergent electromagnetic induction and topological Hall effects in skyrmion lattices.
INTRODUCTION
Skyrmions, i.e., smooth and topologically stable spin textures with a nontrivial topology (1, 2), potentially offer notable advancement for the regime of next-generation electronics (3, 4), spintronics (5), and condensed matter physics (6–10). Magnetic skyrmions have been reported in several chiral magnetic single crystals (11–13), polycrystalline alloys (14), multilayered thin films (15–18), and multiferroics (19, 20) and are envisaged as potential candidates for high-density racetrack memories (21), cache memory (22), ultralow energy spintronic switches (23), and skyrmion reshufflers for probabilistic computing (24). For the realization of these potential applications, the energy efficient creation, annihilation, and controlled manipulation of skyrmion lattices (SkLs) (5, 25, 26) in designed architectures are essential. Manipulation of SkLs has been demonstrated using spin-polarized currents (27), electric field (28), uniaxial strains (29), boundary conditions (30, 31), and magnetoelectric coupling (32–34). However, deterministic control of SkLs still remains elusive, thus necessitating systematic studies for tailored boundary conditions.
Insulating Cu2OSeO3 (CSO) has gained notable interest due to its magnetoelectric coupling and is considered as an ideal platform to achieve electric field control of SkLs. Recently, an intriguing discovery of a new skyrmion phase (35) and of low-temperature skyrmions (36, 37) has been made in CSO crystals, implying richer phase diagrams. The existence of two distinct skyrmion phases in a multiferroic material could be a generic character of underlying physics and may assist in skyrmion manipulation (36). In chiral cubic CSO crystals, the structure and prevalence of the spin textures are determined by the competing magnetic Hamiltonian terms: exchange correlation (J), relativistic Dzyaloshinskii-Moriya (DM) interactions, and crystal anisotropy energy. The magnetic ground state in CSO crystals at low temperature is represented by helical spin arrangements with the modulation period (λ) proportional to the ratio of J and DM interactions (15). Within a range of an applied external magnetic field, a hexagonally packed SkL along the direction of the applied magnetic field is predictably formed.
Chemical doping has been used to enhance DM interactions in B20 compounds due to the resulting variation in DM vector or electronic band filling (38). Similarly, it has been shown that the energy bandgap values change with the Te-dopant concentration in CSO crystals that affects the magnetic properties of the system (39). For one-dimensional (1D) or 2D confined structures, helical spin textures are often observed to be distorted, and geometric anisotropies develop because of important roles played by surface/interface and demagnetization (or stray field) effects. In general, the skyrmion phase can be stabilized over wider ranges of temperature (T) and magnetic field (H) in confined structures, compared to the narrow A-phase in bulk crystals. Thus, a systematic study on chemical doping and sample geometry is necessary to understand the skyrmion phase dynamics in multiferroic CSO crystals toward their deterministic manipulation.
Here, we report the study of systematic control of helical-to-skyrmion phase transitions in thin films of multiferroic Te-doped CSO single-crystals as a function of material thickness, doping, temperature, and magnetic field. By in situ Lorentz phase microscopy, nanoscale magnetic spin evolutions in tailored thin films were directly imaged to obtain Te-doping and thickness-dependent phase maps obtained over a range of T and H. With increasing magnetic field, we find an anisotropic nonlinear change in the spiral helix periodicity accompanied by reduced Lorentz imaging contrast that are associated with out-of-plane tilting of the helix based on our theoretical model. In addition, we report gradual skyrmion injection from thinner to thicker sections with an edge-locked wave vector, termed “skyrmion channeling.” Our systematic study demonstrates that thickness and Te doping are viable parameters for controlling skyrmion H-T phase spaces in multiferroic CSO. The thickness-dependent skyrmion channeling envisages chirality-controlled skyrmion injections in skyrmion flow circuits for designing future skyrmionic devices.
RESULTS AND DISCUSSION
Lorentz phase images (Fig. 1, B and C) taken at 25 K and under a small residual magnetic field (~11.7 mT) show two helical spin modulations (h[110] and h[001]) with wave vectors Q oriented along [110] and [001] directions, respectively. The h[110] is observed in the area near the original surface, while the h[001] modulation is found in the area far from the Pt/CSO surface in both samples. Using the transport-intensity-equation–based phase-retrieval method (40), we reconstructed the magnetic induction map of helical phases for both samples, as shown in Fig. 1 (D and E). The modulation periods (λ) at 11.7 mT measured from the magnetic induction map were 58 and 61 nm for the undoped and doped samples, respectively, implying that Te-doping effects result in a slight change in the ratio between J and DM interactions. Shibata et al. (41) extensively studied the doping effect in Mn1−xFexGe and found a marked impact on skyrmion size depending on crystal composition and point to competing interactions such as magnetic anisotropy and exchange interactions for the compensation of a weakened DM interaction. In our study, the DM interaction becomes weakened by the Te substitution with respect to the exchange interaction J due to the increase in the lattice constant by substitutional Te doping (ionic radius of Te4+ is larger than that of Se4+ by ~10%). However, the exact impact of doping on the nature of DMI and other relevant energies in CSO requires further theoretical investigation. At the applied field of 82.6 mT, both samples showed SkLs with lattice constants of 67 and 72 nm, respectively (Fig. 1, F and G). Eventually, the SKLs destabilize and reach the poled ferromagnetic state at large magnetic fields.
(A to C) Sample schematic and Lorentz images of Te-doped and undoped CSO transmission electron microscopy samples at 25 K under residual magnetic field (H ~ 11 mT) along the [
The magnetic-field dependence of the spin textures at 25 K is shown in Fig. 2 (a series of Lorentz images obtained with the increase of applied magnetic field perpendicular to the sample). Corresponding diffractograms from each thickness section show that the evolution of helical-skyrmion phases is clearly thickness dependent. The h[110] phase shows a notable expansion in its modulation period and a clockwise rotation up to ~15o with increasing magnetic field (for example, compare h[110] between 11.7 and 62 mT in the t = 100 nm section).
Helical-to-skyrmion phase transition in Te-doped CSO observed by Lorentz microscopy at 25 K under various external magnetic fields (H) along the [110] direction (A to E). Lorentz phase images from three different thicknesses (69, 100, and 141 nm, respectively) and their corresponding diffractograms (right) showing field and thickness dependence of the spin textures. Anisotropic scaling of two orthogonal helical phases, edge-induced skyrmion nucleation, and skyrmion channeling are observed upon magnetic field application.
In contrast, the modulation period of the h[001] phase does not change notably, while its area gradually decreases with increasing magnetic field. In neutron scattering experiments, the helical states at low magnetic fields are observed along the cubic axes in CSO due to the cubic magnetic anisotropy. We note that the h[110] phase deviates from the cubic axis and only shows a scaling behavior with applied magnetic field. We discuss this anisotropic scaling behavior of helical phases later in more detail.
At the applied field of 41.5 mT, skyrmions begin to nucleate at the edge of the thinnest section (t = 69 nm), as shown in Fig. 2B. Similar nucleation at the sample edge has been reported in metallic FeGe crystals and attributed to preferential nucleation at distorted edge states in the helix (42). Two sets of SkLs (SkL1 and SkL2) are initially formed independently (see the t = 69 nm section with 62 mT and its Fourier transform in Fig. 2C), of which wave vectors are tilted 30o away from each other. With increasing magnetic field, SkL1 becomes dominant over the other, eventually the sample is encompassed by a single set of SkL1. The evolution of magnetic spin textures in the Te-doped and undoped CSO samples at 25 K shows sequential “skyrmion filling” from thinnest sections to thicker ones with increasing magnetic field from 11.7 to 184.3 mT (see the Supplementary Materials for data of helical-skyrmion phase dynamics for the undoped sample). With further increase in magnetic field, the SkL destabilizes and rotates locally that indicates the wave vector coherence has been lost between neighboring sections as demonstrated by the presence of Moiré fringes in the t = 100 nm section at 163.8 mT. Subsequently, SkLs expand with increased rotation up to 30o clockwise before eventual transitioning to the uniform ferromagnetic phase.
The evolution of magnetic spin textures in both undoped and doped samples at 25 K shows sequential skyrmion channeling, in which the SkL nucleates from the helical phase at the edge of the thinnest sample section and propagates to the thicker section, through the helical phase, with increasing magnetic field (Fig. 1, B and C). We term skyrmion channeling as the sequential filling of different thickness sections with skyrmions through transformation of the h[001] helices. We demonstrate that the channeling process for both the doped and undoped samples at lower temperatures, where the helical-skyrmion phase dynamics are less affected by thermal fluctuations, proceeds rather slowly. Figure 3 shows Lorentz images from the two thinnest sections of the doped and undoped samples near the onset of channeling at 15 K. Two sets of SkLs (SkL1 and SkL2) are observed in both samples upon skyrmion nucleation. SkL1 is bound to the top and bottom of the sample edges. Channeling occurs predominantly along the bottom edge into the thicker section [the top edge has Au/Pt from focused ion beam (FIB) preparation], as indicated with blue arrows in Fig. 3 (B and D). As channeling progresses, the SkL1 expands at the expense of SkL2 in both doped and undoped samples. This channeling process at 15 K occurs over a time of about 20 s, as shown in Fig. 3. Effects of geometric confinement and the presence of discrete boundaries have been studied theoretically and experimentally in other skyrmion systems (43–45). Our findings of skyrmion channeling and anisotropic scaling provide important insight toward skyrmion manipulation. Note that the channeling SkL (SkL1) in both the doped and undoped samples have a wave vector component normal to the bottom sample edge. It should be noted that the bottom edge of the doped sample is along the cubic axes ([100]), but that of the undoped sample is 45o tilted from the cubic axes. Similar edge-locked skyrmion wave vectors are reported in FeGe nanostripes (43). We suggest that optimally packed skyrmions at the sample edge minimize fringing fields leading to formation of skyrmion chains along the edges with the wave vector along the normal to the edge. The Moiré fringes found in Fig. 2E occur when the edge-locked skyrmions are detached from the edges and thus are free to rotate. We note that edge geometries, specifically the presence of Au/Pt on one surface, may influence the evolution of skyrmion phases and is worth further investigation. Effects of different edge termination and local strains play important roles on phase transitions and spin texture evolution upon external magnetic field. These edge-locked SkL wave vectors, independent of doping, thickness, temperature, and underlying crystallographic directions, provide a viable way to realize skyrmion flow circuits in thin-film geometry for practical applications.
Skyrmions nucleated at edges in the thinnest section gradually channel across thicker sections. Both doped (A to C) and undoped (D to F) samples show skyrmion channeling. Red dashed lines separate SkL1 and SkL2.
The anisotropic scaling behaviors (dependency of the helical modulation period on external magnetic field) in the doped and undoped samples are plotted in Fig. 4. The h[110] phases (domain II) show a marked expansion of up to 35%, particularly in the t = 100 nm section. In contrast, the h[001] phase (domain I) expansion is notably less prominent and below 5%. Similar scaling behavior is observed in SkLs (initial lattice constant of ~70 nm under 11 mT to 130 nm under 200 mT, measured with diffractograms shown in Fig. 2). In addition to the wavelength expansion, the modulation amplitude of the spin spiral in the domain II drops compared to the domain I by a factor of ~5 to 10 with increasing magnetic field H. The observed effects of the helix modulation period expansion strongly depend on the thickness of the CSO film. These experimental observations can be naturally explained if we assume that, at finite magnetic field H, spin spiral in the domain II tilts out of plane toward the direction of the magnetic field, as shown Fig. 5A. The in-plane component of wave vector of the spiral tilted at the angle θ decreases as Q|| = Q cos(θ); hence, the period of the spiral measured in Lorentz microscopy λ = 2π/(Q cos(θ)) increases with the tilting angle θ. A similar phenomenon of tilting of the spin spiral in CSO crystals toward the direction of the magnetic field was observed in bulk measurements with small-angle neutron scattering (35, 36) and has been attributed to a competition between either cubic and exchange anisotropies (35) or cubic anisotropy and dipolar interactions (36).
(A) Anisotropic contribution to the energy of the spin spiral in domains I and II at H/(ρQ2) = 0.5, β/(ρQ2) = 0.1, keff/(ρQ2) = 0.5. (B and C) Comparison of experimental data on the scaling of the measured spin spiral wave vector Q|| versus perpendicular magnetic field H in (A) undoped CSO, (C) Te-doped CSO. Squares show experimental data, and blue, red, and cyan solid lines are the fit with the theoretical formula, eq. S7. The open (filled) markers correspond to spirals in domain I, the h[001] phase (domain II, the h[110] phases). The black solid line is the theoretical prediction for the domain I, λ = λ0. (D) Extracted value of the phenomenological parameter β versus sample thickness t. The solid lines show β (t) ∝ 1/t.
Left: Out-of-plane tilting of the spin h[110] spiral in the external magnetic field (domain II) H||z. Right: Orientation of the crystal axes in the CSO sample.
From a theoretical point of view, we can understand this behavior as follows. The spin spiral in CSO crystals is generated by the DM interaction, the energy density reads
In very thin films, the energy of the out-of-plane spin spiral with
At sufficiently small magnetic fields, the minimum energy in the domain I corresponds to θ = 0; see Fig. 4A. Therefore, the wave vector of the spin spiral is a constant as a function of the perpendicular magnetic field H. Conversely, for the domain II, the energy is minimal at a nonzero tilting angle θ(H) (see Fig. 4A), and the tilting angle evolves with the magnetic field. The corresponding theoretical curves for
The out-of-plane tilting of the spiral results in the decrease of the relative ratio of helix amplitude A in domains II and I measured with Lorentz microscopy. Since Lorentz microscopy probes the averaged in-plane spin component over the thickness of the sample, we obtain the amplitude suppression
The tilting angle
Last, on the basis of a series of temperature dependent measurements, the H-T phase diagrams for the both doped and undoped samples are shown in Fig. 6. Figure 6 (A to D) represents measurements on the Te-doped sample, whereas Fig. 6 (E to H) was obtained on the undoped sample. The colored regions in both cases denote the appearance of the skyrmion phase. Thickness dependences are clearly seen for both samples; the stable H-T space of SkLs increased with decreasing film thickness. In addition, the skyrmion phase diagrams are quantitatively different for the Te-doped CSO compared to the undoped CSO. For example, the critical magnetic fields (dashed lines on phase diagrams in Fig. 6) to nucleate/annihilate skyrmions in the temperature range of 20 to 40 K are larger for the Te-doped CSO compared to those of undoped CSO. We note that the SkL stability is more strongly dependent on sample thickness compared to that for the undoped sample. In addition, when the film thickness is above 100 nm, a slight tendency of skyrmion phase splitting is observed (see Fig. 6, C and D for t = 141 and 174 nm) for the doped sample, which is consistent with the recently reported two independent skyrmion phases (35) and the skyrmion phase splitting reported for Zn-doped CSO crystals (47). This indicates that the Te doping promotes the separation between the thermal fluctuation–induced skyrmion phase at high temperature and the second skyrmion phase at low temperature. Thus, the Te doping in CSO crystals provides an advantageous platform to investigate the second skyrmion phase in thin films. Previously, chemical doping was reported to stabilize the metastable skyrmion phases at low temperature. We note that there are differences in our obtained phase diagrams compared to chemical doping effects in bulk crystals studied by neutron scattering studies (11, 48, 49), which highlight the importance of boundary conditions and the thin-film geometry on the evolution of skyrmion phases.
Magnetic field–temperature phase maps of the both doped (A to D) and undoped (E to H) CSO crystals with increasing film thickness, respectively. The doped sample shows the SkL stability is more strongly dependent on film thickness compared to the undoped sample. Both doped and undoped samples show that SkLs are more stable with decreasing film thickness.
In summary, the dynamics of the helical and skyrmion phases in thin films of multiferroic Te-doped and undoped CSO with discrete thickness sections are investigated systematically using in situ Lorentz microscopy. We find that, with increasing applied field, the helical-to-skyrmion phase transition are characterized by anisotropic scaling of spiral helical-phase periodicities, from edge-enhanced skyrmion nucleation to gradual skyrmion channeling to thicker sections. Our theoretical study shows that the anisotropic scaling behavior of helices is the result of out-of-plane tilting of the helix. Thickness-dependent gradual skyrmion channeling and enhanced dependence on boundary conditions by Te doping reported in this study shed new light on the SkL manipulation for skyrmion-based spintronics.
MATERIALS AND METHODS
Single-crystalline undoped and Te-doped CSO crystals were synthesized by the chemical vapor transport method (19) yielding millimeter-sized crystals. Detailed information on synthesis and crystal characterization can be found in the Supplementary Materials. To investigate the effects of chemical doping and sample thickness on the skyrmion phase dynamics, transmission electron microscopy (TEM) samples with various thicknesses were prepared by FIB. Schematics of TEM samples and Lorentz images of undoped and Te-doped samples are shown in Fig. 1 (A to C). Each undoped and Te-doped CSO samples have four sections of different thickness, which are measured by the low-loss function of electron energy-loss spectroscopy. The thin films prepared in [
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/13/eaax2138/DC1
X-ray diffraction
Energy-dispersive x-ray spectroscopy
Second skyrmion phase at lower T and higher magnetic field
Micromagnetic simulation
Theoretical analysis
Table S1. Crystal data and structure refinement for Cu2OSe(Te)O3.
Table S2. Fractional atomic coordinates and equivalent isotropic displacement parameters (Å2 × 103).
Table 3. Anisotropic displacement parameters (Å2 × 103).
Fig. S1. Energy-dispersive x-ray spectroscopy maps of Cu, Se, O, and Te.
Fig. S2. Skyrmion phase diagrams obtained from Te-doped and undoped TEM samples with relatively uniform thickness (~150 nm).
Fig. S3. Micromagnetic simulations.
Movie S1. Video of Te-doped sample under increasing magnetic field at 15 K.
Movie S2. Video of Te-doped sample under increasing magnetic field at 25 K.
Movie S3. Video of Te-doped sample under increasing magnetic field at 40 K.
Movie S4. Video of undoped sample under increasing magnetic field at 15 K.
Movie S5. Video of undoped sample under increasing magnetic field at 25 K.
Movie S6. Video of undoped sample under increasing magnetic field at 40 K.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).