Abstract
Realizing room-temperature magnetic skyrmions in two-dimensional van der Waals ferromagnets offers unparalleled prospects for future spintronic applications. However, due to the intrinsic spin fluctuations that suppress atomic long-range magnetic order and the inherent inversion crystal symmetry that excludes the presence of the Dzyaloshinskii-Moriya interaction, achieving room-temperature skyrmions in 2D magnets remains a formidable challenge. In this study, we target room-temperature 2D magnet Fe3GaTe2 and unveil that the introduction of iron-deficient into this compound enables spatial inversion symmetry breaking, thus inducing a significant Dzyaloshinskii-Moriya interaction that brings about room-temperature Néel-type skyrmions with unprecedentedly small size. To further enhance the practical applications of this finding, we employ a homemade in-situ optical Lorentz transmission electron microscopy to demonstrate ultrafast writing of skyrmions in Fe3-xGaTe2 using a single femtosecond laser pulse. Our results manifest the Fe3-xGaTe2 as a promising building block for realizing skyrmion-based magneto-optical functionalities.
Similar content being viewed by others
Introduction
Magnetic skyrmions, which are topological swirling spin configurations stabilized by DzyaloshinskiiâMoriya interaction (DMI)1,2, have garnered significant interest over the past decade because of their nanometric scale and exotic magnetoelectronic properties3, such as the topological Hall effect (THE)4, skyrmion Hall effect5, and ultra-low current density for motion6,7, making them ideal information carriers for future high-density and fast-speed data storage8, quantum and neuromorphic computation9,10. Since the existence of skyrmion crystal was first verified in helimagnet MnSi11, various skyrmion-host three-dimensional (3D) bulk materials have been discovered with abundant magnetic and electronic features12,13. Compared to the 3D bulk ferromagnets, two-dimensional (2D) van der Waals (vdW) ferromagnets have inherent superiorities for practical applications in spintronic devices due to their unique atomic layered structure14,15, such as long-range magnetic order down to atomic thickness16,17, wide flexibility for stacking artificial heterostructures18, high sensitivity to external field perturbations19, and possible compatibility with modern integrated circuit process20. Therefore, exploring magnetic skyrmions with small size, especially at room temperature (RT), in 2D vdW ferromagnetic materials with facile tunability have become a focal point of magnetic and topological order of matters, as well as spintronic applications21.
To realize RT magnetic skyrmions within the isolated 2D vdW ferromagnets, DMI induced by inversion crystal symmetry breaking and RT ferromagnetism are crucial prerequisites22. Nevertheless, for most of the intrinsic 2D vdW ferromagnets discovered hitherto, such as Fe3GeTe223, Fe5GeTe224, CrGeTe325, and CrI326, etc., on one hand, their natural centrosymmetric crystal structures exclude the DMI;27 on the other hand, due to the strong MerminâWagner fluctuations that suppress the intrinsic magnetic order in the 2D limit28, their Curie temperatures (Tc) are typically below RT. Therefore, on the way pursing RT magnetic skyrmions in 2D vdW ferromagnets, how to simultaneously improve Tc and introduce DMI is an internationally recognized difficulty. Although several methods, such as elemental doping and substitution etc.29,30, have shown the possibility to increase Tc while breaking the centrosymmetric structure for creating DMI in intrinsic 2D vdW ferromagnets31,32, the elemental controllability and raising the Tc above RT are still challenging33,34. One better choice is to directly introduce DMI into intrinsic 2D vdW ferromagnets with Tc above RT. The iron-based ternary telluride Fe3GaTe235 and chromium-based binary telluride CrTe236 are the only two 2D vdW intrinsic magnets that exhibit RT ferromagnetism discovered hitherto, but merely the Fe3GaTe2 exhibits both above Tc and large RT perpendicular magnetic anisotropy (PMA) up to the order of magnitudes of 105âJ/m3â35,. Recently, a Fe3GaTe2-based magnetic tunnel junction (MTJ) achieved a large tunnel magnetoresistance (TMR) of 85% at RT37, highlighting the great potential of Fe3GaTe2 for developing spintronic devices. However, the intrinsically inversion-symmetric crystal structure of Fe3GaTe2 possesses an inherent obstacle. Hitherto, the universal strategy that enables directly breaking the inversion symmetry and achieving substantial DMI for stabilizing RT skyrmions in the Fe3GaTe2 have not been realized.
In this study, we discover that the iron deficiency in Fe3GaTe2 can lead to a pronounced displacement of the Fe atoms within the crystal structure. Based on systematic structural analysis and first-principles calculations, we find that this atomic displacement causes a transformation from the original centrosymmetric crystal structure to a non-centrosymmetric structure, resulting in a significant DMI. Combined Lorentz transmission electron microscopy (LTEM), magnetic force microscopy (MFM) and magneto-transport measurements demonstrate that the non-stoichiometric Fe3-xGaTe2 could accommodate Néel-type skyrmions together with a prominent topological Hall effect over a broad temperature range of 330âK to 100âK. Moreover, the size of the skyrmions decreases as the sample thickness becomes thinner, and field-free sub-100ânm skyrmions can be obtained at RT when the thickness falls below a threshold ranging from 40 to 60ânm, which are the smallest skyrmions achieved hitherto in 2D vdW magnets. More intriguingly, with the use of a homemade in-situ optical LTEM, we realize an ultrafast writing of RT skyrmions in the non-stoichiometric Fe3-xGaTe2 thin flakes by a single femtosecond (fs) laser pulse, which offers a possible avenue for the realization of ultrafast and energy-efficient skyrmion-based logic and memory devices.
Results and discussions
To control the Fe content, we systematically grew a series of Fe3-xGaTe2 single crystals by varying the Fe content in the raw material composition, utilizing a Te-flux method (see Methods section and Supplementary Note 1). To determine the chemical composition of the as-grown crystals, energy dispersive X-ray spectroscopy (EDX) analyses were conducted on the surfaces of Fe3-xGaTe2 nanoflakes (Fig. 1a) that were exfoliated and placed onto the Si3N4 membrane (see Methods). The ratio of raw materials and the corresponding final crystal composition are listed in Table S1, Supplementary Fig. S1 and Fig. 1b. We found that the Fe deficiencies always exist in these crystals, while the minimum and maximum Fe contents correspond to Fe2.84±0.05GaTe2 and Fe2.96±0.02GaTe2, respectively. This result implies the feasibility of inducing Fe deficiency in the samples. To highlight the existence of Fe vacancies, the subsequent studies were focused on the minimum Fe content sample Fe2.84±0.05GaTe2. The Raman measurements revealed a clear red-shift of the A1 peak in the non-stoichiometric Fe2.84±0.05GaTe2 to 126âcmâ1, in contrast to 130âcmâ1 observed for the stoichiometric Fe3GaTe2 reported previously (Fig. 1c)35, which can be attributed to the localized effect of Fe-vacancy defects38. We conducted the magnetic characterization of the Fe2.84±0.05GaTe2 crystal with the external magnetic field (B) applied along the out-of-plane direction. Figure 1d displays the temperature-dependent magnetization (M-T) curve measured under a small magnetic field of 30 mT using a field-cooled protocol. By calculating the first derivative (dM/dT) of the M-T curve, we determined the Curie temperature (Tc) of the Fe2.84±0.05GaTe2 crystal to be ~ 350âK, which is slightly lower than that of the stoichiometric Fe3GaTe2. Additionally, an obvious magnetization kink was observed in the M-T curve at around 290âK, as indicated by the dashed box in Fig. 1d. This particular kink is typically regarded as a signature of rotation of the magnetic easy axis, and is often observed in skyrmion-hosting magnetic systems23,39. Moreover, the field-dependent magnetization curves for the Fe2.84±0.05GaTe2 sample with minimum Fe content reveal an out-of-plane easy magnetization direction at room-temperature (Fig. S2). These curves exhibit magnetic anisotropy almost identical to that of Fe2.96±0.02GaTe2 sample with high Fe content.
The crystal structure of the stoichiometric Fe3GaTe2 was confirmed to have the centrosymmetric space group of P63/mmc, and can be visualized as a series of Te-Fe3Ga-Te monolayers stacked along the c-axis (Supplementary Fig. S3). Within each monolayer, the central Feii-Ga slice flanked by two adjacent Fei slices is sandwiched between two outer Te slices, indicating a protected \(c\,\to \,-c\) mirror symmetry. Regarding the non-stoichiometric Fe2.84±0.05GaTe2 in this work, we analyzed the crystal structure on a microscopic scale (high-resolution scanning transmission electron microscopy, HR-STEM) and a macroscopic scale (single-crystal X-ray diffraction, XRD). Figure 1e presents a typical high-angle annular dark field (HAADF) image of the Fe2.84±0.05GaTe2 nanoflake along the [0001] zone axis. The image demonstrates that the Fei atom columns (dark brown circles) are hexagonally surrounded by six Te-Feii-Te-Ga atomic columns (marked with light brown and pale green balls), matching well with the standard Fe3GaTe2 crystal structure. Moreover, its lattice parameters (aâ=âbâ=â4.08âà ) are slightly enlarged compared to those of the stoichiometric Fe3GaTe2 (aâ=âbâ=â3.99âà )35. These observations suggest that the Fe-deficiency does not cause noticeable lattice distortion in the ab plane. However, the HAADF image along the [\(11\bar{2}0\)] zone axis (Fig. 1f) reveals that Feii atoms deviate clearly from the center position of the Te slices along the c-axis, which is also supported by the annular bright-field (ABF-STEM) image in Fig. S4. By referencing the center of the two Te atoms in a magnified ABF-STEM image (Fig. S5), an averaged Feii deviation is calculated as â0.16â±â0.06âà over an area of 2 à 17 unit cells (see Supplementary Note 2). Notably, this deviation can be also identified in the bottom layer of the Fe2.84±0.05GaTe2 unit cell, and both the offset and deviation direction are the same as those observed in the top layer. To gain additional insights into the structural features of the Fe2.84±0.05GaTe2, we performed single-crystal XRD on a single crystal with dimensions of approximately 3 à 3 mm2. A total of 534 symmetry-independent reflections, corresponding to Miller indices of â5 ⤠hââ¤â5, â5 ⤠kââ¤â5, â20 ⤠lââ¤â21, were collected for the structural determination and refinement. It is well known that for the centrosymmetric space group P63/mmc of stoichiometric Fe3GaTe2, the reflection patterns, such as (\({hh}\overline{2h}{l}\)) and (000âl), are permitted only for even values of l (lâ=â2n, where n are integers), whereas they are forbidden for the odd values (lâ=â2nâ+â1)33,40. In the case of Fe2.84±0.05GaTe2, however, a series of weak (\(11\bar{2}l\)) and (000âl) reflection patterns, such as (\(11\bar{2}3\)) and (\(11\bar{2}5\)) (see Fig. S6), were detected for lâ=â2nâ+â1, suggesting a substantial deviation from the original crystal structure of Fe3GaTe2. By fitting the reflections while allowing for the relaxation of z-positions and occupation ratio of Fe atoms, we obtained a symmetry-lowering structural model with the non-centrosymmetric space group P3m1 (Supplementary Note 2 and Fig. S7). The refined crystal structure revealed that the Feii atoms are deficient and located at the Wyckoff site 1c and 1b at (x, y, z) = (2/3, 1/3, 0.7482) and (1/3, 2/3, 0.2486) with a slight deviation of Feii atoms from the Ga-Ga plane, as shown in Fig. 1g. Based on the experimentally established crystal structure (Fig. S8a), we simulated the HR-STEM images along the [0001] and [\(11\bar{2}0\)] zone axis (Fig. S8b, c), which agree well with our experimental observations41. Moreover, we also simulated the selected area electron diffraction (SAED) along the \([10\bar{1}0]\) and \([11\bar{2}0]\) zone axes. The presence of odd l values in the (000âl) diffractions, evident in both simulated and experimental SAED results (Fig. 1h and Fig. S9), further confirms that the crystal structure of the Fe2.84±0.05GaTe2 belong to the non-centrosymmetric P3m1 space group, rather than the centrosymmetric P63/mmc structure of Fe3GaTe2. These simulations agree well with our experimental observations, providing solid evidential support for our structural modelâs reliability.
In comparison to the stoichiometric Fe3GaTe2 with a centrosymmetric structure, the presence of Fe deficiency in Fe2.84±0.05GaTe2 should exert a pivotal influence on the Feii deviation for the asymmetric structure. Our refined single-crystal XRD indicates that Fe deficiency is predominantly concentrated at the Feii positions with an occupancy ratio of 0.8467. Additionally, the upper-layer Fei-a sites have an occupancy ratio of 0.9688, while the under-layer Fei-b sites are nearly fully occupied (Fig. 1g). As observed in the line profile of Fei-a and Fei-b atoms in the ABF-STEM image (Fig. S5c), it is apparent that the image intensity of Fei-a above Feii is weaker than that of Fei-b below Feii. Since the ABF imaging intensity is generally proportional to the number of projected atoms42, the contrast difference between Fei-a and Fei-b indicates asymmetric site occupations, suggesting a small quantity of Fe vacancies in the Fei-a site, which is consistent with the results of single-crystal XRD. To assess the influence of Fei-a and Feii vacancies on Feii deviation, we further conducted first-principles calculations involving structure relaxation under three scenarios: no vacancy, Fei-a vacancy, and Feii vacancy (Supplementary Note 3). The electron density of Fe3-xGa atoms is depicted in Fig. S10 to facilitate a comparison of the alterations in Feii chemical bonding: (a) The perfect Fe3GaTe2, with no Fe vacancy, showcases a hexagonally bonded Feii-Ga plane. In this arrangement, the centrally positioned Fei-a and Fei-b dimers do not form direct bonds with Feii and Ga atoms. Thus, the overall chemical bonding is mirror-symmetric along the Feii-Ga plane with no Feii deviation. (b) The presence of Feii vacancy induces deformation of the Ga electron density within the ab plane. Nevertheless, no bonding is established between the Feii-Ga plane and the Fei atoms, which remain a mirror-symmetric electron density with no Feii deviation. (c) In case of Fei-a vacancy, there is additional electron-density overlapping between the lower Fei-b atom and its three nearest Feii atoms, while no overlapping between Fei-b and its three nearest Ga atoms. As a result, chemical bonding between the Feii and Fei-b atoms induces a substantial Feii deviation, with a calculated δc (Feii â Ga) of about â0.0554âà , which compares favorably to the XRD result of â0.0871âà . Furthermore, the calculated formation energy value for Fei vacancy (2.96âeV/Fe) is higher than Feii vacancy (2.86âeV/Fe), indicating that the formation of Feii vacancies is more favorable. The above Fe vacancy model and calculated Feii deviation align well with the analysis from single-crystal XRD and ABF-STEM image. Therefore, we conclude that the asymmetric vacancy of Fei-a induces a displacement of Feii atoms towards the âc direction, which results in the symmetry breaking of the Fe2.84±0.05GaTe2 crystal structure.
Due to the asymmetric crystal structure, a pronounced DMI is naturally expected in the non-stoichiometric Fe2.84±0.05GaTe2. In quantitative terms, the DMI vector can be expressed as
where D is the DMI constant, \({\hat{{{{{{\bf{u}}}}}}}}_{{{\mbox{ij}}}}\) represents the unit vector from Fei atom to Feii atom, and \(\hat{{{{{{\bf{z}}}}}}}\) represents the unit vector from magnetic Feii atom to heavy Te atom. For Fe2.84±0.05GaTe2, since its structural symmetry is broken by the Feii atom deviation, the DMI is proposed to originate from the interactions between the neighboring Fei-Feii pair and the adjacent Te atoms. As schematically shown in Fig. 1i, each Fe2.84±0.05GaTe2 unit cell has two distinct DMI sources: (a) the interaction between the Fei-Feii atom pair and the upper Te atom (corresponding to D1 vector, represented by the red arrow) and (b) the interaction between the Fei-Feii pair and the lower Te atom (corresponding to D2 vector, represented by the blue arrow). Equation (1) indicates that the two DMI vectors D1 and D2 are perpendicular to the Fei-Feii-Te triangle of atoms. This is consistent with the Fert-Levy DMI observed at heavy metal/ferromagnet interfaces, which can lead to the formation of Néel-type skyrmions2. Moreover, we found that the directions of D1 and D2 are opposite due to the opposite directions of \(\hat{{{{{{\bf{z}}}}}}}\) for D1 and D2. For the centrosymmetric Fe3GaTe2, because the Feii atom is located at the center of two adjacent Te atoms, the upper D1 and lower D2 vectors would always cancel out with each other, resulting in an effective net DMI vector (Deff) of zero. In the case of Fe2.84±0.05GaTe2, however, the deviation of Feii atoms shows â0.0871âà atom displacement towards âc directions, which breaks the inversion symmetry, and thus makes the nonequal D1 and D2 yield a nonzero Deff within each monolayer. As for the total DMI (Dtotal) of the unit cell, since the Deff vectors in top and bottom layers have the same direction, the magnitude of Dtotal is the sum of the two vectors. Moreover, Fig. S11b illustrates the effective net DMI vectors Deff viewed from [\(0001\)] zone axis, which are perpendicular to the Fei-Feii-Te atom cross sections and exhibit threefold rotational symmetry within the ab plane. Based on the model depicted in the aforementioned illustration, we quantitatively investigated the relationship between the Feii deviation value δc(Feii â Ga) and the DMI constant D (see Supplementary Note 4 and Fig. S12 for details). For a centrosymmetric structure (δcâ=â0), the absence of Feii deviation yields Dâ=â0âmJ/m2. Conversely, a non-centrosymmetric structure (δcâ=ââ0.0871âà ), determined by single-crystal XRD, corresponds to Dâ=â0.91âmJ/m2. This significant DMI constant meets the crucial requirement for the formation of skyrmions43,44,45,46, suggesting that the non-stoichiometric Fe2.84±0.05GaTe2 has the potential to exhibit topological magnetism.
For a magnetic system hosting topological spin configurations, the total Hall resistivity (\({\rho }_{{{\mbox{xy}}}}\)) typically comprises three components4,47:
where \({\rho }_{{xy}}^{N}{=R}_{0}H\) is the normal Hall resistivity, \({R}_{0}\) is the normal Hall coefficient; \({\rho }_{{xy}}^{A}={S}_{A}{\rho }_{{xx}}^{2}M(H)\) is the anomalous Hall resistivity,\(\,{S}_{A}\) is the scaling coefficient, \({\rho }_{{xx}}\) is the longitudinal resistivity, M(H) is the magnetic field-dependent magnetization; and \({\rho }_{{xy}}^{T}\) is the topological Hall resistivity. Of these components, \({\rho }_{{xy}}^{T}\), which is driven by the local spin chirality, is widely regarded as a key transport signature of topological spin configurations. To investigate the potential existence of \({\rho }_{{xy}}^{T}\) in Fe2.84±0.05GaTe2, we fabricated a Hall device with a sample thickness of 250ânm (Fig. 2a) and conducted measurements of \({\rho }_{{xy}}\) with the external magnetic field applied along the normal direction of the device over the temperature range of 350 â 10âK, as shown in Fig. 2b. The results indicate that \({\rho }_{{xy}}\) varies in a nonlinear manner with sweeping the magnetic field, suggesting the presence of a ferromagnetic-order induced anomalous Hall effect. We subsequently fitted \({\rho }_{{xy}}\) using the M(H) and \({\rho }_{{xx}}(H)\) curves in the high-field region based on Eq. (2). Figure 2c displays both the fitting and experimental \({\rho }_{{xy}}(H)\) curves at 300âK. It is clear that there is a discrepancy at the low-field region, indicating the presence of a pronounced THE component in \({\rho }_{{xy}}\). The extracted magnetic field-dependent \({\rho }_{{xy}}^{T}\) curves over the temperature range of 350 â 10âK are summarized in Fig. 2d. It is important to note that the maximum value of \({\rho }_{{xy}}^{T}\) (~1.13âμΩ cm) at 300âK is one order of magnitude higher than that of the material systems hosting RT skyrmions, such as Co-doped 2D magnets (Fe0.5Co0.5)5GeTe2 (\(\sim\)0.8âμΩcm)31, Kagome ferromagnets Fe3Sn2 (~ â0.4âμΩâcm)47, and Ir/Fe/Co/Pt multilayers (~0.03âμΩâcm)48, and is comparable to that of noncoplanar ferromagnet Cr5Te6 (~1.6âμΩâcm) at 90 K49. Furthermore, we investigated thickness- and magnetic field-dependent Hall resistivity \(\rho_{xy}\) and topological Hall resistivity \(\rho^{T}_{xy}\) under room temperature (Fig. S13). As the sample thickness increases, the topological Hall signals gradually strengthen and shift towards higher magnetic field. More importantly, the THE signals persist over a broad temperature range and various thickness, suggesting the existence of topological spin configurations in Fe2.84±0.05GaTe2.
To directly visualize the possible topological spin textures associated with the topological Hall effect, we conducted cryo-LTEM experiments on an exfoliated [0001]-oriented Fe2.84±0.05GaTe2 nanoflake with a thickness of approximately 100ânm, as schematically illustrated in Fig. 3a. For the cryo-LTEM measurements, we heated the sample above Tc (380âK) and then cooled it to the desired temperature with a liquid nitrogen TEM holder under a 30 mT out-of-plane magnetic field. Following this field-cooling procedure, no magnetic contrast in the Lorentz phase images was observed when the electron beam was injected along the normal direction of the nanoflake at 300âK (Fig. 3b and Supplementary Fig. S14). However, if we titled the nanoflake an angle (θ) of ±20° away from the horizontal plane, bubble-like domains with half-dark/bright contrast were detected. It is well known that the magnetic contrast discernible in a Lorentz phase image acquired via Fresnel imaging mode is predominantly caused by the deflection of the electron beam due to the in-plane magnetic field as it passes through a magnetic sample. Since in the Bloch-type domain walls the magnetization, or spin, direction is out-of-plane, the Lorentz deflection is thus in-plane and perpendicular to the walls, yielding clear domain walls contrast. In contrast, in the Néel-type domain walls, the magnetization direction in the walls is in-plane and the Lorentz deflection is along the wall length, resulting in no contrast. However, once the sample is tilted, this is no longer the case as an out-of-plane magnetization component is generated50,51,52. The angle-dependent contrast modulation seen in our experiments unambiguously suggests that the bubble-like domains we imaged are Néel-type. We further analyzed their in-plane spin structures using the transport-of-intensity equation (TIE)53, as displayed in the right panel of Fig. 3b. We found that the in-plane magnetic induction was composed of a pair of conjoined clockwise and counterclockwise spin swirls, which agrees well with the calculated magnetic induction map for the Néel-type skyrmions31,52. On the basis of the deduced double in-plane swirls and Néel-type spin arrangements, we conclude that the observed bubble-like domains are indeed Néel-type skyrmions.
Due to the high Tc of Fe2.84±0.05GaTe2, its magnetic phase is expected to remain stable well above room temperature. The upper panel of Fig. 3c presents the Lorentz phase images taken at a zero magnetic field over the temperature range of 300 to 340âK after the field-cooling operation (see Methods). The results indicate that the skyrmion phase remain stable up to 330âK, which represents a record-high value compared with that of the skyrmion-hosting magnetic vdW materials reported to date. However, when the temperature exceeds 330âK, the magnetic skyrmions become elongated and gradually transformed into stripe domains due to the presence of thermal fluctuations, which are proposed to be large enough to overcome the topologically protected energy barrier between skyrmions and stripe domains above 330âK. Furthermore, the magnetic field-dependent domain evolution process after the field-cooling operation was studied within a temperature range of 340â100âK (lower panel of Fig. 3c and Supplementary Fig. S15). To illustrate the correlation between the magnetic states and both the magnetic field and temperature, we construct a magnetic phase diagram (Fig. 3d), in which the red region denotes that the sample hosts high-density skyrmions at the corresponding T-B plane, while the blue region indicates the absence of skyrmions. As displayed in Fig. 3d, high-density, field-free skyrmions can be stabilized in a wide temperature range of 100â330âK, demonstrating that Fe2.84±0.05GaTe2 is a promising material platform for use in spintronic devices.
We further use magnetic force microscopy (MFM), which is sensitive to the out-of-plane magnetic field in the sample, to study the effect of nanoflake thickness on the stabilization of the skyrmion phase at room temperature after field-cooling manipulation. All the nanoflakes used for the MFM measurements were freshly exfoliated from the same batch to ensure consistency (Supplementary Fig. S16). Additionally, we used a low-moment MFM tip (<10âmT) to minimize the influence of the tipâs magnetic field on the domain structures during scanning. Figure 4a illustrates a series of MFM images that record the magnetic field-dependent domain evolution processes at six typical thicknesses. Notably, the vertically arranged MFM images represent the same sample area under different out-of-plane fields, while the horizontal direction represents the variation in thickness for different samples. For the 250ânm thick sample, the skyrmions are arranged densely into honeycomb lattices with a large skyrmion size (dsk) of approximately 180ânm at zero magnetic field. As the external magnetic field increases, there is little change in dsk, accompanied, however, by a decrease in density with the sudden annihilation of skyrmions at 210âmT. When the sample reaches saturation magnetization, the MFM images with uniform contrast represent the fully ferromagnetic (FM) state. Varying the sample thickness could significantly affect the skyrmion size due to the change in the strength of dipole-dipole interaction33,54. Upon decreasing the sample thickness, the skyrmion size at a zero magnetic field reduces correspondingly, reaching a minimum skyrmion size of 87ânm when the sample thickness is decreased below 100ânm (Fig. 4b). This value is much smaller than that of the skyrmions in other skyrmion-hosting vdW magnets reported to date, such as Cr1+xTe2 (~400ânm)34, Fe3GeTe2 (~250ânm)55, Fe5GeTe2 (~200ânm)24, and (Fe0.5Co0.5)5GeTe2 (~150ânm)31, as displayed in Fig. 4c (see Supplementary Fig. S17 for detailed size distribution of field-free skyrmions in Fe2.84±0.05GaTe2 at varied temperatures, and Fig. S18 for the corresponding Micromagnetic simulations). Figure 4d provides a comprehensive overview of the RT skyrmion phase diagram, depicting the relationship between thickness, magnetic field, and the occurrence of skyrmion states. The field-free skyrmions are stable in a broad thickness range and the highest density appears to be between 46 and 60ânm. As the thickness of Fe2.84±0.05GaTe2 nanoflake decreases to 25ânm, the stripe domains remain while the field-cooling process no longer generates skyrmions. Compared with the previously reported skyrmion-hosting 2D material (Fe0.5Co0.5)5GeTe231, Fe2.84±0.05GaTe2 exhibits smaller magnetic parameters such as DMI constant D, saturation magnetization Ms, threshold of sample thickness t and etc. (see Fig. S19 and Supplementary Note 5 for the determination of magnetic parameters), which contribute to the reduction in skyrmion size (see Table S3, Fig. S20 and Supplementary Note 6 for the corresponding micromagnetic simulations). Moreover, our micromagnetic simulations demonstrate the influence of in-plane magnetic field B on skyrmion shape, revealing a progressive transformation from a circular to an elliptical configuration (see Fig. S21 and Supplementary Note 7).
Achieving ultrafast and energy-efficient writing of skyrmions is a critical requirement in the pursuit of their practical applications in high-speed and low-power spintronic devices. Conventional field-cooling operations are unsuitable as they involve heating the entire sample above the Tc and slow cooling to the desired temperature55,56,57. In contrast, fs laser pulses can demagnetize a localized region on the micrometer scale within fs timescale followed by a picosecond (ps) thermal quenching58,59,60,61,62,63,64. This unique approach enables the creation of metastable magnetic states at extremely short timescales, offering a promising pathway for the ultrafast writing of skyrmions. Several studies have reported successful writing of skyrmions using short laser pulses in various materials, including magnetic multilayer films60,61,62, FeGe63, and Co9Zn9Mn264. To further investigate the possibility of ultrafast writing of RT skyrmions in the 2D van der Waals Fe2.84±0.05GaTe2 nanoflakes, we performed Lorentz phase microscopy measurements under fs laser excitation using our homemade in-situ optical LTEM setup that enables single-shot fs laser pulse excitation on the sample (520ânm wavelength, 50 μm focal spot size, and 300âfs pulse duration), as illustrated in Fig. 5a. We recorded the fs laser writing processes by taking snapshots of Lorentz phase images before and after the laser pulse excitation (fluence of ~ 11âmJ/cm2), while the electron beam was blanked during laser pulse excitation to avoid unexpected sample damage. Upon laser pulse excitation of the initial domains, ultrafast heating with electron temperature (Telec) above Tc would temporarily melt the long-range magnetic order (Fig. 5b) on the fs time scale59,60,63. Subsequently, the paramagnetic state is rapidly quenched on the ps timescale due to the electron-phonon coupling, where the phonon temperature (Tphon) follows the variation of Telec (see Supplementary Note 8 and Fig. S22 for detailed calculations of temperature evolutions based on a two-temperature model). Such quenching process leads to a rapid decrease of the temperature to below Tc and initiates the remagnetization process of the spin system. Since only skyrmion is a stable solution for the Fert-Levy DMI observed in Fe2.84±0.05GaTe2, it is expected that skyrmion states could be created during the thermal relaxation in the quenching process.
To demonstrate the feasibility of using a fs laser pulse to write RT skyrmions in the 2D van der Waals Fe2.84±0.05GaTe2 material, we conducted a step-by-step single-shot fs laser pulse excitation measurement on different initial field-polarized magnetic states, including stripe and single domains. These initial field-polarized domain structures were achieved by varying the out-of-plane magnetic field from a saturation state (±2âT), as indicated by the normalized Hall hysteresis curve in Fig. 5c. Moreover, the corresponding Lorentz phase images captured the magnetic states before and after the laser pulse excitation at six representative magnetic fields (see Supplementary Figs. S23, S24 for more details). These images provide compelling evidence that the presence of an external assisting magnetic field is vital for the successful fs laser writing of skyrmions in the material60,61,62,63. At Bâ=â0âmT, only stripe domains were observed after the laser excitation (fluence of ~11âmJ/cm2). As B is swept towards positive saturation, the width of the stripe domains gradually decreases. After the fs laser excitation, the magnetic states accessible by the laser pulse initially exhibit a mixed state of stripe domains and skyrmions at 11âmT and then transition to fully-formed skyrmions with the highest density at 46âmT. As B is increased to 80âmT, the initial stripe domains start to fade out, and the laser-accessible skyrmions completely vanish, transforming into a single-domain state. Compared to the single domain state stabilized by a pure magnetic field (Bâ=â104âmT), the assisting magnetic field decreases significantly, suggesting that for the laser-writing operation the laser functions as an efficient field that can effectively shift the magnetic states, possessing a higher energy state in the magnetization process to a lower energy state. When the magnetic field is increased to 104âmT, the initial stripe domains transform into a single domain state, and the laser pulse with the used fluence can no longer induce a magnetic phase transition. However, as the magnetic field is decreased from the positive saturation field, the single domain state persists at a much lower magnetic field of 46âmT due to the hysteresis effect. Interestingly, we find that there exists a magnetic field window (69âmT to 46âmT) where laser-accessible skyrmions can be created from the initial single domain state, even though it is in a higher energy state than that of the skyrmions in the magnetization process. As the magnetic field continues to decrease, the behavior under laser excitation in the low-field region (35âmT to 0âmT) is similar to that at the beginning of the positive field sweep (0âmT to 35âmT). These results indicate that the laser-writing operation is independent of the initial domain structures but rather dependent on the strength of the assisting magnetic field. Furthermore, conventional zero-field cooling can only result in the formation of interconnected, relatively long stripe domains, but does not spontaneously lead to the creation of skyrmions57. To demonstrate the differences with in-site fs laser quenching approach, we conducted fluence-dependent laser pulse excitation without magnetic field (see Fig. S25 and Supplementary Note 9). After a single laser pulse with fluence of 1.3 mJ/cm2, stripe domains show slightly domain wall movement. Upon increasing to 9.4 mJ/cm2, a hybrid state with both stripes and skyrmions are formed, while at 11 mJ/cm2 only stripe domains observed. This indicates that a higher fluence of fs laser can completely demagnetize the sample during the laser-writing process, regardless of the initial magnetic state.
To further understand the underlying magnetization dynamic process of the ultrafast laser writing of skyrmions in Fe2.84±0.05GaTe2 under external magnetic field assistance, we performed finite-element micromagnetic simulations on the subsequent magnetic structure evolution of an initial stripe domain after a fs laser pulse excitation at a certain external magnetic field based on the Landau-Lifshitz-Gilbert equation with Langevin dynamics65,66,67. We considered the following scenarios for the simulations: (i) the fs laser pulse interacts with the magnetic structures through photothermal effect; and (ii) the fs pulse heats the sample above the Tc and melts the electronic spin structures, but without changing the atomic lattice. Specifically, the laser quenching-induced magnetization dynamics was achieved by relaxing the magnetic system from the laser-induced paramagnetic state under an out-of-plane magnetic field of 46 mT (see Supplementary Note 5 and Fig. S19 for more details about the simulations). As shown in Fig. 5d (see details in Supplementary Movie 1), following the excitation by the femtosecond (fs) laser pulse, the initial melted spin state (snapshot at t0) rapidly evolves into numerous nanoscale spin clusters. These clusters contain topological defects, including skyrmionic and anti-skyrmionic nucleation centers (snapshot at t1). This transformation occurs due to the ultrafast cooling, achieved at a quenching rate of up to 1012âK/s58,68. Because only skyrmion is a stable solution for the Fert-Levy DMI observed in Fe2.84±0.05GaTe2, in the further cooling process the anti-skyrmionic nucleation centers merge and annihilate with the nearby skyrmionic nucleation centers and appear less frequently until completely disappear, namely, the topological fluctuations60,69, resulting in formation of a pure skyrmionic state (snapshot at t2). The skyrmion nucleation process is complete at t2, but the skyrmion size has not reached the energy equilibrium state. In order to minimize the magnetostatic energy, the domain walls move and gradually evolve into a uniform skyrmion lattice (snapshot at t3). The micromagnetic simulations reproduce well the magnetic phase transitions observed in our experiments, confirming the reliability of our established relationship between the laser-writing process, magnetic field strength, and magnetic domain evolution. Such evidenced controllable ultrafast laser writing of skyrmions in 2D van der Waals magnetic materials provides opportunities for both fundamental researches and device applications towards magneto-optical control of spin topologies.
We report the discovery of a field-free sub-100ânm Néel-type skyrmion state in non-stoichiometric Fe2.84±0.05GaTe2 over a broad temperature range from 330âK to 100 K. Using HR-STEM and single-crystal XRD, we determine that the deviation of Feii atoms from the center position of the Te slices due to the asymmetric Fei-a vacancies induces a transformation of the crystal structure from centrosymmetric to non-centrosymmetric, enabling the formation of skyrmions through the in-plane isotropic DMI. LTEM (along with MFM) shows the size of the skyrmions decreases with sample thickness, and a field-free sub-100ânm skyrmion state was achieved at RT within a specific sample thickness range of 40ânm to 60ânm. Furthermore, we demonstrate that a single fs-laser pulse can rapidly generate field-free sub-100ânm skyrmions from both stripe domains and single domains. Our study demonstrates not only the non-stoichiometric Fe2.84±0.05GaTe2 to be a promising material platform for exploring magnetic skyrmions, but also the fs-laser can be a powerful tool to manipulate and control topological chiral spin textures to realize skyrmion-based high-speed logic and memory applications.
Methods
Single crystal growth and structure characterization
Single crystals of Fe3-xGaTe2 were grown by the self-flux method. The mixtures of Fe (99.99%), Ga (99.99%), and Te (99.99%) elements were mounted in an alumina crucible and sealed inside a quartz tube under high vacuum (~10â4âPa). The mixtures were firstly heated at 1150â°C for 24âh, then followed by slow-cooling down to 850°C for 3 weeks. Finally, excessive molten flux was centrifuged to separate the single crystals. Its chemical composition was determined by energy dispersive x-ray spectroscopy mapping (EDX, Bruker Nano GmbH Berlin). A comprehensive overview of the raw material composition and the final crystal composition are listed in Supplementary Note 1 and Table S1. The single-crystal X-ray diffraction was carried out with a four-circle diffractometer (Bruker D8 venture). The refined crystal structure was solved and refined by using the Bruker SHELXTL Software Package. The HR-STEM images were acquired by high-resolution transmission electron microscope (HRTEM, JEOL ARM200F). Magnetometry measurements were carried out with the Quantum Design PPMS.
LTEM and MFM measurements
The fresh Fe2.84±0.05GaTe2 nanoflakes utilized for Hall devices, LTEM and MFM experiments were prepared through an all-dry mechanical-transfer method within an argon-filled glovebox. In an argon-protected environment, the nanoflakes were first produced on PDMS stamp by micromechanical cleavage, and then transferred onto Si3N4 membrane or SiO2/Si substrate, with or without pre-patterned Au electrodes. To maintain the integrity of the samples, these freshly prepared Hall devices, LTEM, and MFM specimens were promptly placed into a plastic box within the argon-filled glovebox and securely sealed with parafilm before removal. During subsequent measurements outside the glovebox, the exposure of the samples to air was minimized, and all the sample transfers were conducted within a confined timeframe of no more than 10âminutes. LTEM measurements were performed on Thermo Fisher Talos F200i at an acceleration voltage of 200âkV. The objective lens was used to apply a magnetic field perpendicular to the sample plane by controlling the excitation current. The specimen was in-suit warmed and cooled by a liquid nitrogen double-tilt sample holder. In order to perform field cooling with accurate temperature control, the sample was initially cooled to +20âK above the desired temperature at a rate of â5 K/min, and then slowly cooled down to the desired temperature at a rate of â1âK/min. The in-plane magnetization distribution map was reconstructed from the under- and over-focused images using the transport-of-intensity equation (TIE) approach. MFM measurements were performed by scanning probe microscopy (MFP-3D, Asylum Research), which equipped a low-moment magnetic tip (PPP-LM-MFMR, Nanosensors) and VFM3 component (Asylum Research). To protect samples from air degradation, Hall and LTEM measurements were conducted under vacuum conditions, while MFM measurements were carried out in an environment continuously flushed with argon gas to ensure effective protection.
In-situ optical LTEM experiments
To directly visualize the fs laser writing of RT skyrmions in the 2D van der Waals Fe2.84±0.05GaTe2 nanoflakes, we performed the in-situ optical LTEM experiments in our homemade 4D-electron microscopy (Thermo Fisher Talos F200i), which enables single-shot fs laser pulse excitation on the sample under the Lorentz phase imaging mode. The Lorentz phase images were acquired under the Fresnel mode, in which the external perpendicular magnetic field was applied by the objective lens with controlled lens current. The fs laser system was triggered externally with a digital delay generator which outputted single-shot fs laser pulses with 520ânm wavelength, 300âfs pulse duration, and 11âmJ/cm2 fluence, where the laser spot size was adjusted to be 50 μm to ensure homogeneous illumination on the sample.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The single crystal X-ray diffraction, Lorentz transmission electron microscopy, and characterization data generated in this study are provided in the Supplementary Information file. The data that support the findings of this study are available from the corresponding authors upon reasonable request.
References
Fert, A., Reyren, N. & Cros, V. Magnetic skyrmions: advances in physics and potential applications. Nat. Rev. Mater. 2, 17031 (2017).
Fert, A. & Levy, P. M. Role of anisotropic exchange interactions in determining the properties of spin-glasses. Phys. Rev. Lett. 44, 1538â1541 (1980).
Kang, W., Huang, Y., Zhang, X., Zhou, Y. & Zhao, W. Skyrmion-electronics: an overview and outlook. Proc. IEEE 104, 2040â2061 (2016).
Neubauer, A. et al. Topological Hall effect in the A phase of MnSi. Phys. Rev. Lett. 102, 186602 (2009).
Litzius, K. et al. Skyrmion Hall effect revealed by direct time-resolved X-ray microscopy. Nat. Phys. 13, 170â175 (2017).
Jiang, W. et al. Blowing magnetic skyrmion bubbles. Science 349, 283â286 (2015).
Iwasaki, J., Mochizuki, M. & Nagaosa, N. Current-induced skyrmion dynamics in constricted geometries. Nat. Nanotechnol. 8, 742â747 (2013).
Wang, Y. et al. Electric-field-driven non-volatile multi-state switching of individual skyrmions in a multiferroic heterostructure. Nat. Commun. 11, 3577 (2020).
Song, K. M. et al. Skyrmion-based artificial synapses for neuromorphic computing. Nat. Electron. 3, 148â155 (2020).
Xia, J., Zhang, X., Liu, X., Zhou, Y. & Ezawa, M. Universal quantum computation based on nanoscale skyrmion helicity qubits in frustrated magnets. Phys. Rev. Lett. 130, 106701 (2023).
Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. Science 323, 915â919 (2009).
Tokura, Y. & Kanazawa, N. Magnetic skyrmion materials. Chem. Rev. 121, 2857â2897 (2021).
Kanazawa, N., Seki, S. & Tokura, Y. Noncentrosymmetric magnets hosting magnetic skyrmions. Adv. Mater. 29, 1603227 (2017).
Li, H., Ruan, S. & Zeng, Y.-J. Intrinsic Van Der Waals magnetic materials from bulk to the 2D limit: new frontiers of spintronics. Adv. Mater. 31, 1900065 (2019).
Gibertini, M., Koperski, M., Morpurgo, A. F. & Novoselov, K. S. Magnetic 2D materials and heterostructures. Nat. Nanotechnol. 14, 408â419 (2019).
Gong, C. et al. Discovery of intrinsic ferromagnetism in two-dimensional van der Waals crystals. Nature 546, 265â269 (2017).
Huang, B. et al. Layer-dependent ferromagnetism in a van der Waals crystal down to the monolayer limit. Nature 546, 270â273 (2017).
Gong, C. & Zhang, X. Two-dimensional magnetic crystals and emergent heterostructure devices. Science 363, eaav4450 (2019).
Burch, K. S., Mandrus, D. & Park, J.-G. Magnetism in two-dimensional van der Waals materials. Nature 563, 47â52 (2018).
Zhang, L., Zhou, J., Li, H., Shen, L. & Feng, Y. P. Recent progress and challenges in magnetic tunnel junctions with 2D materials for spintronic applications. Appl. Phys. Rev. 8, https://doi.org/10.1063/5.0032538 (2021).
Yao, X., Hu, D. & Dong, S. Modulation of skyrmionic magnetic textures in two dimensional vdW materials and their heterostructures. iScience https://doi.org/10.1016/j.isci.2023.106311 (2023).
Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions. Nat. Nanotechnol. 8, 899â911 (2013).
Ding, B. et al. Observation of Magnetic Skyrmion Bubbles in a van der Waals Ferromagnet Fe3GeTe2. Nano Lett. 20, 868â873 (2020).
Schmitt, M. et al. Skyrmionic spin structures in layered Fe5GeTe2 up to room temperature. Commun. Phys. 5, 254 (2022).
Han, M.-G. et al. Topological magnetic-spin textures in two-dimensional van der Waals Cr2Ge2Te6. Nano Lett. 19, 7859â7865 (2019).
Ghosh, S., StojiÄ, N. & Binggeli, N. Comment on âMagnetic skyrmions in atomic thin CrI3 monolayer. Appl. Phys. Lett. 114, 232402 (2019).
Laref, S., Kim, K.-W. & Manchon, A. Elusive Dzyaloshinskii-Moriya interaction in monolayer Fe3GeTe2. Phys. Rev. B 102, 060402 (2020).
Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models. Phys. Rev. Lett. 17, 1133â1136 (1966).
Chen, X. et al. Pervasive beyond room-temperature ferromagnetism in a doped van der Waals magnet. Phys. Rev. Lett. 128, 217203 (2022).
Zhang, H. et al. A room temperature polar magnetic metal. Phys. Rev. Mater. 6, 044403 (2022).
Zhang, H. et al. Room-temperature skyrmion lattice in a layered magnet (Fe0.5Co0.5)5GeTe2. Sci. Adv. 8, eabm7103 (2022).
Meisenheimer, P. et al. Ordering of room-temperature magnetic skyrmions in a polar van der Waals magnet. Nat. Commun. 14, 3744 (2023).
Chakraborty, A. et al. Magnetic skyrmions in a thickness tunable 2D ferromagnet from a defect driven DzyaloshinskiiâMoriya interaction. Adv. Mater. 34, 2108637 (2022).
Saha, R. et al. Observation of Néel-type skyrmions in acentric self-intercalated Cr1+δTe2. Nat. Commun. 13, 3965 (2022).
Zhang, G. et al. Above-room-temperature strong intrinsic ferromagnetism in 2D van der Waals Fe3GaTe2 with large perpendicular magnetic anisotropy. Nat. Commun. 13, 5067 (2022).
Freitas, D. C. et al. Ferromagnetism in layered metastable 1T-CrTe2. J. Phys. Condens. Matter 27, 176002 (2015).
Zhu, W. et al. Large room-temperature magnetoresistance in van der Waals ferromagnet/semiconductor junctions. Chin. Phys. Lett. 39, 128501 (2022).
Cong, X., Liu, X.-L., Lin, M.-L. & Tan, P.-H. Application of Raman spectroscopy to probe fundamental properties of two-dimensional materials. NPJ 2D Mater. Appl. 4, 13 (2020).
Hou, Z. et al. Observation of various and spontaneous magnetic skyrmionic bubbles at room temperature in a frustrated Kagome magnet with uniaxial magnetic anisotropy. Adv. Mater. 29, 1701144 (2017).
Flack, H. & Shmueli, U. International Tables for Crystallography A, 114â128 (2016).
Barthel, J. Probe: a software for high-resolution STEM image simulation. Ultramicroscopy 193, 1â11 (2018).
Kourkoutis, L. F., Song, J. H., Hwang, H. Y. & Muller, D. A. Microscopic origins for stabilizing room-temperature ferromagnetism in ultrathin manganite layers. Proc. Natl Acad. Sci. 107, 11682â11685 (2010).
Yang, H., Thiaville, A., Rohart, S., Fert, A. & Chshiev, M. Anatomy of Dzyaloshinskii-Moriya Interaction at Co/Pt Interfaces. Phys. Rev. Lett. 115, 267210 (2015).
Yang, J. H. et al. Strong Dzyaloshinskii-Moriya interaction and origin of ferroelectricity in Cu2OSeO3. Phys. Rev. Lett. 109, 107203 (2012).
Xiang, H. J., Kan, E. J., Wei, S.-H., Whangbo, M. H. & Gong, X. G. Predicting the spin-lattice order of frustrated systems from first principles. Phys. Rev. B 84, 224429 (2011).
Cui, Q., Liang, J., Shao, Z., Cui, P. & Yang, H. Strain-tunable ferromagnetism and chiral spin textures in two-dimensional Janus chromium dichalcogenides. Phys. Rev. B 102, 094425 (2020).
Li, H. et al. Large topological Hall effect in a geometrically frustrated kagome magnet Fe3Sn2. Appl. Phys. Lett. 114 https://doi.org/10.1063/1.5088173 (2019).
Soumyanarayanan, A. et al. Tunable room-temperature magnetic skyrmions in Ir/Fe/Co/Pt multilayers. Nat. Mater. 16, 898â904 (2017).
Chen, Y. et al. Observation of colossal topological Hall effect in noncoplanar ferromagnet Cr5Te6 thin films. Adv. Funct. Mater. n/a, 2302984 (2023).
Garlow, J. A. et al. Quantification of mixed Bloch-N\âeel topological spin textures stabilized by the Dzyaloshinskii-Moriya interaction in Co/Pd multilayers. Phys. Rev. Lett. 122, 237201 (2019).
Pollard, S. D. et al. Observation of stable Néel skyrmions in cobalt/palladium multilayers with Lorentz transmission electron microscopy. Nat. Commun. 8, 14761 (2017).
Peng, L. et al. Tunable Néelâbloch magnetic twists in Fe3GeTe2 with van der Waals structure. Adv. Funct. Mater. 31, 2103583 (2021).
Yao, Y. et al. Chirality flips of skyrmion bubbles. Nat. Commun. 13, 5991 (2022).
Wang, H. et al. Characteristics and temperature-field-thickness evolutions of magnetic domain structures in van der Waals magnet Fe3GeTe2 nanolayers. Appl. Phys. Lett. 116, 192403 (2020).
McCray, A. R. C. et al. Thermal hysteresis and ordering behavior of magnetic skyrmion lattices. Nano Lett. 22, 7804â7810 (2022).
Ding, B. et al. Tuning the density of zero-field skyrmions and imaging the spin configuration in a two-dimensional Fe3GeTe2 magnet. NPG Asia Mater. 14, 74 (2022).
Birch, M. T. et al. History-dependent domain and skyrmion formation in 2D van der Waals magnet Fe3GeTe2. Nat. Commun. 13, 3035 (2022).
Fu, X. et al. Optical manipulation of magnetic vortices visualized in situ by Lorentz electron microscopy. Sci. Adv. 4, eaat3077 (2018).
Khela, M. et al. Laser-induced topological spin switching in a 2D van der Waals magnet. Nat. Commun. 14, 1378 (2023).
Büttner, F. et al. Observation of fluctuation-mediated picosecond nucleation of a topological phase. Nat. Mater. 20, 30â37 (2021).
Je, S.-G. et al. Creation of magnetic skyrmion bubble lattices by ultrafast laser in ultrathin films. Nano Lett. 18, 7362â7371 (2018).
Ogasawara, T., Iwata, N., Murakami, Y., Okamoto, H. & Tokura, Y. Submicron-scale spatial feature of ultrafast photoinduced magnetization reversal in TbFeCo thin film. Appl. Phys. Lett. 94, 162507 (2009).
Berruto, G. et al. Laser-induced skyrmion writing and erasing in an ultrafast cryo-lorentz transmission electron microscope. Phys. Rev. Lett. 120, 117201 (2018).
Shimojima, T. et al. Nano-to-micro spatiotemporal imaging of magnetic skyrmionâs life cycle. Sci. Adv. 7, eabg1322 (2021).
Ostler, T. A. et al. Ultrafast heating as a sufficient stimulus for magnetization reversal in a ferrimagnet. Nat. Commun. 3, 666 (2012).
Kuzâmin, M. D., Skokov, K. P., Diop, L. V. B., Radulov, I. A. & Gutfleisch, O. Exchange stiffness of ferromagnets. Eur. Phys. J. 135, 301 (2020).
Lichtenberg, T. et al. Anisotropic laser-pulse-induced magnetization dynamics in van der Waals magnet Fe3GeTe2. 2d Mater. 10, 015008 (2023).
Eggebrecht, T. et al. Light-induced metastable magnetic texture uncovered by in situ Lorentz microscopy. Phys. Rev. Lett. 118, 097203 (2017).
Zheng, F. et al. Skyrmionâantiskyrmion pair creation and annihilation in a cubic chiral magnet. Nat. Phys. 18, 863â868 (2022).
Acknowledgements
This work was supported by the National Key Research and Development Program of China at grant No. 2020YFA0309300, Science and Technology Projects in Guangzhou (Grant No. 202201000008), the National Natural Science Foundation of China (NSFC) at grant No. 12304146, 11974191, 12127803, 52322108, 52271178, U22A20117, and 12241403, China Postdoctoral Science Foundation (2023M741828), Guangdong Basic and Applied Basic Research Foundation (Grant nos. 2021B1515120047 and 2023B1515020112), the Natural Science Foundation of Tianjin at grant No. 20JCJQJC00210, the 111 Project at grant No. B23045, and the âFundamental Research Funds for the Central Universitiesâ, Nankai University (Grant Nos. 63213040, C029211101, C02922101, ZB22000104, and DK2300010207). This work was supported by the Synergetic Extreme Condition User Facility (SECUF). The work at Brookhaven National Laboratory was supported by the Materials Science and Engineering Divisions, U.S. DOE-BES under Contract No. DE-SC0012704.
Author information
Authors and Affiliations
Contributions
X.W.F., Z.P.H. and Z.F.L. conceived the research project. Z.F.L., H.Z. and J.T.G. synthesized the crystals and characterized the structure. Z.F.L. and H.Z. performed the Lorentz phase microscopy, magnetic force microscopy, and magneto-transport measurements. G.Q.L. performed the first-principles calculations. Z.F.L. performed the in-situ optical Lorentz phase microscopy. Z.F.L. and Q.P.W. did the micromagnetic simulation. The manuscript was drafted by Z.F.L., Z.P.H., Y.M.Z. and X.W.F. with contributions from H.Z., G.Q.L., J.T.G., Q.P.W., Y.D., Y.H., X.G.H., C.L., M.H.Q., X. S., R.C.Y., X.S.G., Z.M.L. and J.M.L. All the authors contributed to the discussion and revision of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Peter Meisenheimer, Yan Xin, and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, Z., Zhang, H., Li, G. et al. Room-temperature sub-100ânm Néel-type skyrmions in non-stoichiometric van der Waals ferromagnet Fe3-xGaTe2 with ultrafast laser writability. Nat Commun 15, 1017 (2024). https://doi.org/10.1038/s41467-024-45310-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-45310-2
This article is cited by
-
Distinct skyrmion phases at room temperature in two-dimensional ferromagnet Fe3GaTe2
Nature Communications (2024)
-
High-temperature Néel skyrmions in Fe3GaTe2 stabilized by Fe intercalation into the van der Waals gap
npj Spintronics (2024)