In this work, we demonstrate all-optical, stochastic magnetisation modulation of a topological spin structure at room temperature. Specifically, we show that the magnetisation winding of Néel spin twists in dipolar skyrmions within kagome-type Fe3Sn2 thin films can be changed by femtosecond optical light pulses. The ultrafast switching phenomena were monitored by in situ Lorentz TEM imaging at magnetic remanence and at room temperature. A statistical evaluation of the electron micrographs uncovered a stochastic mechanism of spin texture variation within the optically excited dipolar skyrmions, while retaining the Bloch-type internal circular domain wall chirality. In addition, utilising an off-axis magnetic field as an additional experimental degree-of-freedom, the optically induced interconversion of type II bubbles and type I skyrmions was observed, involving a substantial spin reorganisation in the Bloch- and Néel structures. Theoretical calculations, based on micromagnetic modelling and simulations of the electron phase shift, provide insights into the magnetisation dynamics of the dipolar skyrmions induced by the ultrafast light pulse. The kagome-type Fe3Sn2 material thus constitutes an intriguing platform for the development of all-optical magnetisation switching of topological particles.
For investigating the optical response of FeSn by in situ Lorentz TEM (Fig. 1a, b), a single-crystalline, electron-transparent lamella was prepared (for details see Methods), which exhibits a stripe-like domain pattern in the magnetic ground state with alternating out-of-plane magnetic field alignment. By applying an external magnetic field of 0.53 T perpendicular to the kagome plane (magnetic easy axis), dipolar skyrmions were generated and subsequently imaged at magnetic remanence. A typical Lorentz TEM micrograph recorded with 1.2 mm image defocus is presented in Fig. 1c, showing regions characterised by stripe-like magnetic domain states alongside numerous circular-shaped particles displaying different contrast patterns. The majority of the particles have been identified as dipolar skyrmions of type I, which can be described as a cylindrical domain with a magnetisation antiparallel to the surrounding out-of-plane spin orientations, and continuous closure domain walls. The mid section of the dipolar skyrmions shows a Bloch-type texture with clockwise (cw) or counterclockwise (ccw) field rotation, while close to the surfaces the field arrangement aligns in a twisted Néel structure. Since type I bubbles can be described with the same integer topological charge Q as chiral Bloch skyrmions, they are termed dipolar skyrmions. As marked by a rectangle in Fig. 1c, some particles show contrast variations in their circular shape, indicating domain walls with a spin alignment oriented towards the in-plane direction. These particles are identified as type II magnetic bubbles, which have a vanishing topological charge.
We observe different helicities of type I dipolar skyrmions in FeSn, characterised by bright or dark Lorentz contrast. The dipolar skyrmions have a hybrid magnetic state containing both Bloch- and Néel-type magnetic textures with a smooth transition, moving from the mid-section of the film to the top and bottom surfaces. Whereas the helicity of the Bloch domain wall determines the Lorentz image contrast in the outer parts of the dipolar skyrmion, the presence of twisted Néel caps and their respective helicities is the contributing factor to the spot-like contrast observed in the middle of dipolar skyrmions and type II bubbles. Figure 1d shows the spin texture of a type I dipolar skyrmion, obtained by a micromagnetic simulation (for details, see Methods and Section 2.4). The hybrid Bloch-Néel structure generates a characteristic Fresnel image contrast sensitive to in-plane components of the magnetic field. Specifically, a bright or dark central dot appears when the rotation directions of the bottom and top Néel caps are both ccw or cw, respectively. In cases where the upper and lower surfaces exhibit opposite spin rotation directions, the central image intensity is neutral. All three cases are also visible in the Lorentz TEM micrograph displayed in Fig. 1c. To further disentangle the individual contributions of the different parts of the dipolar skyrmion on the total Lorentz contrast, we calculated the electron phase shift maps for thin film slices for the case of equal Néel cap rotation at both surfaces, as shown in Fig. 1e. As visible in the slices, the image contrast in the centre of the skyrmion originates from both surfaces, whereas the outer Bloch-type contrast is dominated from mid-section slices. For the co-rotating Néel caps, the contributions from both surfaces add up in the total phase shift (Fig. 1e, right panel). In contrast, for counter-rotating Néel caps (Fig. 1f) the surface contributions cancel out, leaving a neutral image contrast in the centre of the dipolar skyrmion (see Supplementary Fig. 1 for further combinations of Bloch wall and Néel cap helicities).
In parallel, an in situ heating experiment inside the TEM was carried out to study the temperature dependence of dipolar skyrmions in FeSn near the Curie temperature (see Supplementary Fig. 5). The recorded Lorentz TEM images showed that type II bubbles transformed to type I structure, then they tend to be replaced by stripe domains the temperature increased. The magnetic contrast has gradually decreased and disappeared near 643 K. After reaching T = 670 K and slightly higher temperature, the specimen was slowly cooled. Interestingly, no formation of bubbles or dipolar skyrmions occurred upon cooling. The emergence of stripe domains was only observed, thereby indicating the detrimental effect of temperature on the stability of dipolar skyrmions.
As an experimental platform to study the optical control of internal degrees-of-freedom in dipolar skyrmions, we utilize the ultrafast transmission electron microscope (UTEM) at the Regensburg Center for Ultrafast Nanoscopy (RUN). Upon illumination of the FeSn lamella by a single, fs optical pulse, as depicted in Fig. 1a (800 nm central wavelength, 169 fs pulse duration, p-polarized 100 nJ pulse energy, corresponding to a fluence of 14 mJ/cm), Lorentz TEM images in the equilibrium state and in magnetic remanence were recorded. A series of such micrographs is presented in Fig. 2. After each light pulse, we observe slight variations in the position and shape of individual particles, and, more strikingly, significant contrast fluctuations are observed at the centres of dipolar skyrmions, where the image contrast can invert, diminish, or reappear, showing either bright or dark contrast. However, the outer ring of the circular contrast features remains unaffected by optical excitation. The observed phenomenon is attributed to light-induced switching of the helicity of the Néel caps of individual dipolar skyrmions, while the core Bloch domain wall maintains its configuration. An assembled video of the FeSn specimen subjected to 100 consecutive incident optical pulses, each with a fluence of 14 mJ/cm, is provided as Supplementary Movie 1.
In our experiments, the incident optical pulse fluence was varied between 10.6 mJ/cm and 17.7 mJ/cm, which does not seem to have a significant effect on the core contrast variations (see Supplementary Figs. 2 and 3). The temperature increase of the specimen in this range remains below 100 degree (see Supplementary Methods 4), which is much lower than the Curie temperature of FeSn. In the limit of low pulse fluences (below 10.6 mJ/cm), a substantial decline in switching probability was found. Conversely, for pulse fluences exceeding 35.4 mJ/cm, structural damage and non-reversible magnetic contrast changes were observed.
In order to extract correlations within the stochastic light-induced switching process, subsequent Lorentz TEM images were analysed, each recorded after excitation with a single optical pulse for a total of up to 100 incident pulses (for details see Supplementary Method 1). Histograms of the distribution of the image intensity within a circular region of the dipolar skyrmions show three distinct peaks (see Supplementary Fig. 2) allowing for a classification of the Néel caps into one of three states: the bright contrast state, corresponding to ccw rotating Néel caps at the top and bottom of the FeSn lamella (ccw state); the mixed state, associated with opposite helicities of the two Néel caps (m state); or the dark contrast state, representing a cw rotation of the Néel caps (cw state). To clearly distinguish between these states, upper and lower image contrast threshold values were assigned to each individual dipolar skyrmion based on the minima of the image intensity distribution histograms (see Supplementary Methods 1 and Fig. 2). Figure 3a and b presents the results of two selected analyses for two dipolar skyrmion with ccw and cw rotation of the Bloch wall, respectively. By averaging over all incident pulses, the probability of an individual dipolar skyrmion adopting the ccw, m, or cw state after excitation by a single light pulse is determined. These probabilities are shown in Fig. 3c and d for a total of 16 dipolar skyrmions. Notably, the Néel caps of all investigated dipolar skyrmions adopt the mixed state in approximately 50 % of the events, with relatively small variation (m state standard deviation (std): 0.07) compared to the inter-skyrmion variation observed for the other states (ccw state std: 0.14, cw state std: 0.13). Furthermore, it has been observed that dipolar skyrmions with a ccw Bloch wall exhibit a slightly higher probability for the Néel caps to adopt the ccw state, while the reverse is true for skyrmions with a cw Bloch wall. However, this asymmetry is relatively minor compared to the fluctuations observed between individual skyrmions. Measurements conducted at varying incident optical fluences reveal an overall similar switching behaviour, with a slightly enhanced asymmetry in the probabilities for ccw and cw Bloch walls at incident optical fluences of 17.7 mJ/cm (see the Supplementary Fig. 3). This small probability asymmetry is tentatively ascribed to a magnetic coupling between the rotation of the Néel caps and the internal Bloch wall. However, the substantial inter-skyrmion variations imply that the micromagnetic environment of the dipolar skyrmion exerts a more significant influence on the probability of the Néel caps adopting a particular magnetic state. The switching probabilities from an initial state (i) to a final state (f) after excitation with a single light pulse with a fluence of 14 mJ/cm is presented as a probability matrix in Fig. 3f, which is averaged over dipolar skyrmions with ccw and cw Bloch walls. It is hypothesised that optical excitation will lead to an increased spin temperature, levelling out the corrugation of the spin-free-energy landscape and allowing for an interconversion of individual Néel cap states, as illustrated in Fig. 3g. As the system cools by heat transfer to the environment, the magnetic configuration relaxes into one of the local minima of the free magnetic energy landscape, a process potentially influenced by the local magnetic environment. If the Néel cap configuration is adopted randomly, the final Néel cap state should be independent of the initial state and thus all rows of the probability matrix should be equal. Given that the mixed state can be formed by two Néel cap configurations, the rows of the model probability matrix are given by [1/4 1/2 1/4]. In contrast, the experimentally determined probability matrix demonstrates a dependence on the initial state, suggesting that the switching process follows Markovian behaviour. Only the switching from the mixed state proceeds with the expected probabilities of a random process. Starting from either the ccw or cw twisted Néel state, the dipolar skyrmion tends to remain in the original configuration, demonstrating that the simultaneous reversal of both Néel cap helicities is less likely. At the same time, the probability of switching to the mixed state remains approximately 50 %, indicating that at the experimentally adopted optical excitation level, the helicity of one Néel cap is randomly selected. Notably, applying a magnetic field perpendicular to the sample surface using the electron microscope's objective lens does not alter the overall switching behaviour but significantly enhances the previously mentioned asymmetry (see Supplementary Methods 2 and Fig. 3 for data acquired under an external field of 133 mT). Furthermore, it appears that extending the pulse length to 1 ps does not modify the switching characteristics, as shown in Supplementary Methods 3 and Fig. 4. Potential approaches for controlling the Néel cap switching behaviour by fine-tuning the external magnetic field are currently under investigation.
Following the implementation of a protocol that involved the application of an external magnetic field of 0.4 T to the slightly tilted FeSn lamella, the predominant type of produced magnetic particles was found to be of type II structure. In particular, Lorentz TEM images recorded at magnetic remanence reveal the presence of Bloch wall features with alternating contrast (Fig. 4a). In the analysed image region, 65 particles were present, 4 of which had a type I and 61 had type II structure. After 17 optical pulses with a fluence of 10.6 mJ/cm, a transition to a type I was observed in all type II bubbles. The number of particles of the two types observed after each light pulse is plotted in Fig. 4b. Following the conversion of all bubbles from type II to type I, the newly formed dipolar skyrmions exhibit 43 ccw and 22 cw Bloch magnetisation rotations. Figure 4c illustrates the step-by-step changes of five type II bubbles, which exhibit uniform Bloch helicities. In particle 1, the segment between its two internal domain walls shortened after shots #2 and #4, forming a type I structure with ccw field rotation after shot #5. After shot #9, three out of five particles were transformed into type I structures with the same Bloch helicity. Interestingly, particle number 4, formed after shot #12, displays a cw Bloch rotation despite having the same initial internal state and similar magnetic environment as the others. This observation suggests that the light-triggered type II to type I transition is, at least in part, a random process and is likely connected to the mobility of the Bloch domain defects and their annihilation dynamics in the partially melted spin configuration shortly after the optical pulse. In the absence of an applied field, only the interconversion from type II to type I structures was observed, but not in the opposite direction. This finding suggests that the type II structure possesses a higher free energy compared to the type I bubbles, and that the activation free energy barrier in between these two spin configurations can be surpassed in the hot-spin state. In particular, at an applied field of approximately 133 mT, rare evidence for the light-induced generation of a type II bubble from a type I Bloch structure was observed at an incident optical fluence of 14 mJ/cm (see Supplementary Movie 2).
To understand the fundamental processes involved in the formation of the dipolar skyrmion spin texture and the switching induced by fs light pulses, micromagnetic simulations were performed based on the Landau-Lifshitz-Gilbert (LLG) equation (see Methods) to reveal details of the underlying spin texture variations. Taking values for the saturation magnetisation and the uniaxial magnetocrystalline anisotropy from magnetometry experiments, we performed detailed micromagnetic simulations. For a direct comparison between experiment and theory, in particular with respect to long-range dipolar interactions, the simulation geometry was chosen to be comparable to the sample used in experiments: a thin film with dimensions of 2000 nm × 2000 nm × 200 nm, as shown in Fig. 5. For the top and bottom surfaces of the film we considered soft ferromagnetic properties, in order to approximate the non-ideal structure resulting from the ion milling fabrication and possible oxidation of the TEM specimen. For the soft layers, saturation magnetisation and exchange stiffness were assumed to be the same as in the ideal case, but with a vanishing crystalline anisotropy (K = K = 0). Additionally, thin layers without modified surfaces were tested as well in order to investigate the magnetic state in the ideal case.
The influence of the fs optical pulse was simulated as a temperature spike in the electron system, via a random temperature-dependent thermal field, by raising T to a peak value within 20 ps and subsequently gradually lowering it to room temperature within 100 ps. Various values of T were tested to determine the effective peak temperature that best reproduces the experimental observations, and this was found to be T = 12000 K. In principle, the relevant temperature in the experiments corresponds to the temperature of the spin system, which in metals is expected to increase after fs optical excitation on a few-100-fs time-scale, resulting from electron-spin coupling. Due to the simultaneously occurring electron-phonon coupling on a timescale of about 1 ps, we expect a peak spin temperature on the order of several hundred K. This non-equilibrium heating of the spin system is anticipated to have a small but finite influence on the magnetisation switching dynamics. However, experimental studies employing optical pulses with durations of about 1 ps suggest that the transient lattice temperature spike and the subsequent re-freezing of the spin configuration, play a decisive role in determining the final magnetic state (see Supplementary Methods 3 and Fig. 4). However, in the micromagnetic simulation, macrospins are considered, neglecting spin fluctuations within a simulation cell, which become important at high spin temperatures. Therefore, larger effective temperatures are required in the simulations to induce comparable spin dynamics as observed in the experiments.
Figure 5a shows the top view of the simulated spin field in the thin-film mid-section depicting the Bloch-type domain wall surrounding each magnetic particle. Based on the three-dimensional spin field, we simulated the expected electron phase shift, as shown in Fig. 5b1. All relevant magnetic features observed in the experiments, such as type I dipolar skyrmions, elongated bubble shapes, cw and ccw Bloch helicities, and the presence of Néel caps, are well reproduced in the simulation. The correspondingly simulated magnetic structure after applying a temperature peak and subsequent relaxation is presented in Fig. 5c1. Comparing the magnetic structure before and after excitation (Fig. 5b1,c1), the experimentally observed switching of Néel caps can be clearly reproduced in the simulations. For example, the magnetic texture for the dipolar skyrmion highlighted in Fig. 5b1,c1 changes at the top surface from cw to a ccw rotating Néel cap (Fig. 5b2,c2). Together with the non-changing cw rotating Néel cap at the bottom surface (not shown), the contrast within the centre of the dipolar skyrmion switches from dark to neutral.
In the simulations, we have observed four possible scenarios for different T (not shown here). In the low-temperature limit, the pulse introduces a small degree of randomization but the micromagnetic state remains overall unchanged during and after the pulse. In the temperature range 11 000 K <T < 18 000 K, the temperature pulse provides enough thermal energy to randomly reverse the twisting of the Néel caps, while the body of the Bloch dipolar skyrmions remains unchanged. In the range 18 000 K <T < 24 000 K, the thermal pulse stochastically switches the configuration of both the Néel caps and the Bloch dipolar skyrmions. In the high-temperature limit, the Zeeman energy of the random field overcomes the ferromagnetic exchange, leading to a completely demagnetized state that does not return to its original configuration with dipolar skyrmions after the pulse is removed. The excellent agreement between simulation and experimental observation provides substantial support for our interpretation of the light-induced processes and captures the stochastic nature and diversity of the involved switching phenomena.