All-optical spin switching probability in [Tb/Co] multilayers
Long spin coherence length and bulk-like spin–orbit torque in ferrimagnetic multilayers
Jiawei Yu, Do Bang, … Hyunsoo Yang
Role of non-thermal electrons in ultrafast spin dynamics of ferromagnetic multilayer
Je-Ho Shim, Akbar Ali Syed, … Dong-Hyun Kim
Single-shot all-optical switching of magnetization in Tb/Co multilayer-based electrodes
L. Avilés-Félix, A. Olivier, … I. L. Prejbeanu
All-optical switching of magnetization in atomically thin CrI3
Peiyao Zhang, Ting-Fung Chung, … Xiang Zhang
Photoinduced spin dynamics in a uniaxial intermetallic heterostructure $$\hbox {TbCo}_2/\hbox {FeCo}$$ TbCo 2 / FeCo
Sergei Ovcharenko, Mikhail Gaponov, … Vladimir Preobrazhensky
Efficient spin–orbit torque in magnetic trilayers using all three polarizations of a spin current
Jeongchun Ryu, Ryan Thompson, … Byong-Guk Park
Single pulse all-optical toggle switching of magnetization without gadolinium in the ferrimagnet Mn2RuxGa
C. Banerjee, N. Teichert, … J. Besbas
Deterministic all-optical magnetization writing facilitated by non-local transfer of spin angular momentum
Youri L. W. van Hees, Paul van de Meugheuvel, … Reinoud Lavrijsen
All optical control of magnetization in quantum confined ultrathin magnetic metals
Saeedeh Mokarian Zanjani, Muhammad Tahir Naseem, … Mehmet Cengiz Onbaşlı
L. Avilés-Félix1 nAff4,
L. Farcis1,
Z. Jin1,
L. Álvaro-Gómez1,
G. Li2,
K. T. Yamada2,
A. Kirilyuk3,
A. V. Kimel2,
Th. Rasing2,
B. Dieny1,
R. C. Sousa1,
I. L. Prejbeanu1 &
L. D. Buda-Prejbeanu1
Scientific Reports volume 11, Article number: 6576 (2021) Cite this article
Magnetic properties and materials
Spintronics
Since the first experimental observation of all-optical switching phenomena, intensive research has been focused on finding suitable magnetic systems that can be integrated as storage elements within spintronic devices and whose magnetization can be controlled through ultra-short single laser pulses. We report here atomistic spin simulations of all-optical switching in multilayered structures alternating n monolayers of Tb and m monolayers of Co. By using a two temperature model, we numerically calculate the thermal variation of the magnetization of each sublattice as well as the magnetization dynamics of [\(\text {Tb}_n\)/\(\text {Co}_m\)] multilayers upon incidence of a single laser pulse. In particular, the condition to observe thermally-induced magnetization switching is investigated upon varying systematically both the composition of the sample (n,m) and the laser fluence. The samples with one monolayer of Tb as [\(\text {Tb}_1\)/\(\text {Co}_2\)] and [\(\text {Tb}_1\)/\(\text {Co}_3\)] are showing thermally induced magnetization switching above a fluence threshold. The reversal mechanism is mediated by the residual magnetization of the Tb lattice while the Co is fully demagnetized in agreement with the models developed for ferrimagnetic alloys. The switching is however not fully deterministic but the error rate can be tuned by the damping parameter. Increasing the number of monolayers the switching becomes completely stochastic. The intermixing at the Tb/Co interfaces appears to be a promising way to reduce the stochasticity. These results predict for the first time the possibility of TIMS in [Tb/Co] multilayers and suggest the occurrence of sub-picosecond magnetization reversal using single laser pulses.
All-optical switching (AOS) phenomena in magnetic materials have been object of thorough research over the last decade due to its fundamental interest and relevance as a possible new write scheme in spintronic devices1,2. This phenomenon, which was experimentally discovered for the first time in a rare earth (RE)-transition metal (TM) alloy1, has then been observed in different multilayers and alloys3,4,5. Two possible mechanisms for the magnetization reversal via ultrashort laser pulses were reported: a helicity-dependent one and another that is thermally-driven. In the helicity-dependent all-optical switching (HD-AOS), the asymmetry of the absorption coming from the magnetic circular dichroism is at the origin of the magnetization switching6. Additional mechanisms such as inverse Faraday effect have also been proposed to explain the features of HD-AOS1,7. On the other hand, the thermally induced magnetic switching (TIMS) generally observed in ferrimagnetic materials, is a process mediated by angular momentum transfer between magnetization sublattices during the remagnetization process, following the incidence of a single laser pulse. From these two mechanisms, HD-AOS is the most commonly encountered and it has been observed in a wide variety of materials. For instance, systematic experimental studies on the switching ability of RE-based and RE-free multilayered systems performed by El-Hadri et al.8, demonstrated that engineered synthetic ferrimagnets such as [Ho/Co], [Tb/Co] and [Co/Pt] based multilayers only show HD-AOS. Among all systems switched by HD-AOS, only a few multilayered systems exhibit the possibility of single-shot switching: Pt/Co/Pt/Co/GdFeCo9, Pt/Co/Gd10, FeCoB/Ta/[Tb/Co]\(_\text {N}\)11,12 and Pt/Co/Pt13. From a theoretical point of view, only two reports discuss the possibility of TIMS in \(\text {Tb}_{x} \text {Co}_{1-x}\) alloys14,15. The report from Moreno et al.14 based on atomistic spin dynamics modelling concludes that TIMS is possible in \(\text {Tb}_{x} \text {Co}_{1-x}\)14 alloys with compensation temperature higher than room temperature. However previous experiments only reported helicity-dependent AOS with multiple pulses or pure thermal demagnetization16,17,18. In this framework, it is important to mention that a recent observation of deterministic TIMS in [Tb/Co]-based multilayers evidences a limited fluence range for reversal in [Tb/Co] multilayers systems11.
Our objective is to model synthetic ferrimagnets composed of alternating ferromagnetic layers of Co and Tb with antiferromagnetic exchange coupling between them. The general structure of the multilayered sample is labeled as [\(\text {Tb}_n\)/\(\text {Co}_m\)], where the index n and m correspond to the number of monolayers of each element. In order to explore the probability of TIMS in these multilayers, an atomistic spin solver was used coupled with a two-temperature model (2TM) accounting for the laser heating. Our solver has been developed and benchmarked with respect to previous publications addressing the static magnetic properties and the magnetization dynamics of RE-TM-based systems, such as \(\text {Tb}_{x} \text {Co}_{1-x}\) alloys14 or GdFeCo alloys19. An array of localized magnetic moments on both Tb and Co sites is considered. The associated Heisenberg Hamiltonian reads:
$$\begin{aligned} \hat{{\mathcal {H}}}=-\frac{1}{2} \sum _{i\ne j} J_{ij} \mathbf {S}_{\textit{i}} \cdot \mathbf {S}_{\textit{j}}- \sum _{i} d_{i,z} S_{i,z}^2, \end{aligned}$$
with \(\mathbf {S}_{\textit{i}}\) the classical spin vectors with unit length associated to the magnetic moment direction on site i. The first term represents the Heisenberg exchange energy where the exchange interaction is either between spins of the Tb layer, spins of the Co layer or across the interface between Tb and Co spins. This term therefore contains three interactions \(J_{\text {Tb-Tb}}\), \(J_{\text {Co-Co}}\) and \(J_{\text {Tb-Co}}\), restricted in our model to the nearest neighbors. The second term represents the uniaxial anisotropy with the \(d_{i,z}\) constant per atom. The dynamics of the spins is governed by the stochastic Landau-Lifshitz-Gilbert equation of motion:
$$\begin{aligned} (1+\lambda _i^2)\mu _i\frac{d\mathbf {S}_{\textit{i}}}{dt}=-\gamma _i\mathbf {S}_{\textit{i}} \times [\mathbf {H_{\textit{i}}}(t)+\lambda _i \mathbf {S}_{\textit{i}} \times \mathbf {H_{\textit{i}}}(t)], \end{aligned}$$
where \(\gamma _i\), \(\mu _i\) and \(\lambda _i\) are the gyromagnetic ratio, the value of the magnetic moment and the damping constant of the spin i respectively. The effective field at each site i is derived from the Hamiltonian \(\mathbf {H_{\textit{i}}}(t)=-\frac{\delta \hat{{\mathcal {H}}}}{\delta \mathbf {S}_{\textit{i}}}+\zeta _i\) and includes a stochastic thermal contribution \(\zeta _i\), accounting for the coupling to the heat bath with the following white-noise properties:
$$\begin{aligned} \langle \zeta _i \rangle = 0, \langle \zeta _{i \eta } (0) \zeta _{j \theta }(t)\rangle = 2\delta _{ij}\delta _{\eta \theta }\delta (t) \frac{\lambda _i k_B T_e \mu _i}{\gamma _i}, \end{aligned}$$
where i, j correspond to the lattice sites and \(\eta \), \(\theta \) to the Cartesian components respectively. \(T_e\) is the electron temperature, which is related with the laser power during the excitation of a femtosecond long single laser pulse according to a standard 2TM expressed by two coupled equations:
$$\begin{aligned} C_e \frac{dT_e(t)}{dt}= & {} -G[T_e(t)-T_{ph}(t)]+P(t)-\frac{C_e}{t_0} \left [T_e(t)-T_0(t) \right] \end{aligned}$$
$$\begin{aligned} C_{ph} \frac{dT_{ph}(t)}{dt}= & {} G[T_e(t)-T_{ph}(t)]-\frac{C_{ph}}{t_0} \left[T_{ph}(t)-T_0(t)\right], \end{aligned}$$
where \(T_{ph}(t)\), \(C_e\) = \(\gamma T_e\) and \(C_{ph}\) correspond to the temperature of the phonon bath, the heat capacity of electrons and the heat capacity of the phonons respectively. G is an electron-phonon coupling parameter and P(t) is the external power injected by the laser pulses defined as \(P(t) = (I_0 F) e^{-(t/\tau _p)^2}\), where \(I_0\), F and \(\tau _p\) correspond respectively to the laser energy absorbed by the system, the fluence and the laser pulse duration. Equations (4) and (5) also include a Newton's cooling term which drives the system back to room temperature after the incidence of the laser pulse20,21,22. According to Newton's law the rate of heat loss is proportional to the temperature difference between the body and its surroundings and involves a characteristic time \(t_0\). It is important to mention that the heat produced by the laser pulse will slowly diffuse away from the electron system towards the heat sink fixed at room temperature after a few ps. In experiments, this is the required time for heat to dissipate to the substrate23.
We have considered a multilayered system consisting of \(30 \times 30 \times 30\) fcc cells with periodic boundary conditions. We assumed that there are no size differences between RE and TM atoms. Only exchange with the 12 nearest neighbors is taken into account. A schematic of the atom distribution is shown in Fig. S1 in the Supplementary data file. The local uniaxial anisotropy, the effective magnetic moment of Co and Tb sub-lattices and the interatomic exchange interactions values were taken from Moreno et al., which reports atomistic calculation of TIMS in \(\text {Tb}_{x} \text {Co}_{1-x}\) alloys. The values of various parameters are listed in Table 1. The magnetic moment as well as the intra-lattice exchange correspond to the bulk value of pure Co respectively Tb. Normally, it is expected that these parameters are sample dependent since the precise composition (number of multilayers, defects), the nature of the substrate, material deposition technique as well as annealing condition may influence the precise arrangement of magnetic moments in the lattice. We made the choice to keep these parameters constant upon varying the composition of the multilayers for comparison with a large panel of compositions. Calibration with respect to the real samples combined with ab-initio computation are desirable for a complete quantitative analysis but here we prefer to consider a perfect model system. The time integration of the above-mentioned equations was performed according to a Heun scheme using a time step below 0.1 fs. The parameters of Newton's law are only impacting the system behavior at longer timescales typically above 1 ns.
Table 1 List of parameters used for the numerically calculated static and dynamic properties of \(\text {Tb}_{n} \text {Co}_m\) multilayers.
Static magnetic properties of [\(\text {Tb}/\text {Co}\)] multilayers
The first aim of the study was to obtain the equilibrium properties of the multilayer as a function of their composition. Thus, the equilibrium magnetization was calculated as the spatial and time average of the magnetic moments over the volume of the sample. Figure 1 summarizes the results for three series of samples [\(\text {Tb}_4\)/\(\text {Co}_m\)] (Fig. 1a,d,g), [\(\text {Tb}_3\)/\(\text {Co}_m\)] (Fig. 1b,e,h) and [\(\text {Tb}_4\)/\(\text {Co}_8\)] with various degree of intermixing of the multilayered sample (Fig. 1c,f,i). The normalized \(M_s(T)\) curves for each individual Tb and Co sub-lattice (Fig. 1a,b,c), show typical decrease versus temperature due to increasing spin fluctuations.
Thermal variation of magnetization of [\(\text {Tb}_n\)/\(\text {Co}_m\)] multilayers. The subscripts indicate the number of atomic layers. (a–c) Normalized thermal variation of the magnetization of Tb and Co sublattices. (d–f) Reduced magnetic moment of sample [\(\text {Tb}_n\)/\(\text {Co}_m\)]. Compensation temperature as a function of the number of Co monolayers (g), (h) and intermixing (i). The insets from (i) show the composition of the sample with perfect interfaces and medium intermixing respectively. The grey lines indicate the room temperature.
In the case of ideal multilayered samples (Fig. 1a,b), the magnetization of the Co sublattice closely follows the typical behavior of a ferromagnet since the intra-layer exchange interaction is strong. The Curie temperature of the coupled multilayers is around or above 1000 K approaching the Curie temperature of the Co sub-lattice. The Tb sub-lattice having a weaker intra-layer exchange interaction exhibits magnetic ordering above its own bulk Curie temperature because of the interaction with the Co sublattice, the decrease of its magnetization follows an almost linear dependence with the temperature. The overall magnetization curves (Fig. 1d,e) indicate that the multilayers behave as synthetic ferrimagnets with a corresponding magnetic compensation point (\(\text {T}_{{comp}}\)). An exception is found for the sample [\(\text {Tb}_3\)/\(\text {Co}_3\)] (m = n = 3) having also much smaller Curie temperature around 705 K. Our simulations also confirm that the compensation temperature decreases almost linearly with the Co total magnetic moment as can be seen in Fig. 1g,h. Two samples [\(\text {Tb}_4\)/\(\text {Co}_8\)] and [\(\text {Tb}_3\)/\(\text {Co}_9\)] have a compensation point around room temperature and thus are suitable for TIMS experiments.
The impact of the interfacial intermixing has been investigated in the case of the [\(\text {Tb}_4\)/\(\text {Co}_8\)] sample. For that purpose, a soft, medium and strong intermixing were compared keeping constant the number of Co and Tb atoms. A moderate intermixing shifts the compensation point upwards, but also progressively reduces the Curie temperature. In the limit case of an homogeneous alloy, the Curie temperature is drastically reduced since the Co magnetic moments have less Co neighbors.
Dynamic magnetic properties of [\(\text {Tb}/\text {Co}\)] multilayers
The switching properties of the [\(\text {Tb}_n\)/\(\text {Co}_m\)] multilayers after the incidence of a single laser pulse were analyzed using the two temperature model (Eqs. 4 and 5). The total duration of the laser pulse is set at 50 fs and its fluence is varied up to 120 mJ/\(\text {cm}^2\). The magnetization of each sub-lattice was recorded at each time step and the value at 20 ps after illumination is included in the diagrams presented in Fig. 2 for two series of samples [\(\text {Tb}_3\)/\(\text {Co}_m\)] and [\(\text {Tb}_{1} \text {Co}_m\)] respectively.
Normalized magnetization of Tb and Co sublattices as a function of the laser fluence. \(M_z/M_s\) ratios are extracted from the magnetization curves of the [\(\text {Tb}_n\)/\(\text {Co}_m\)] multilayers 20 ps after the single laser pulse has interacted with the system. Each peak indicates switching, where Tb goes to a negative magnetization and Co goes to a positive magnetization. The noise observed in these graphs reflects the impact of random thermal fluctuations on the switching dynamics. A deterministic switching window is marked with light grey color for [\(\text {Tb}_1\)/\(\text {Co}_3\)] multilayers.
Figure 2a–c, corresponding to [\(\text {Tb}_{3}/ \text {Co}_m\)], shows random peaks of \(M_z/M_s\) after 20 ps, indicating that the magnetization reversal of each sub-lattice occurs randomly for n = 3. However, the multilayer [\(\text {Tb}_1\)/\(\text {Co}_3\)] in Fig. 2d presents a clear deterministic switching window for fluence values between 35 and 50 mJ/\(\text {cm}^2\) in which magnetization reversal occurs. For this composition, the exchange interaction between Tb and Co dominates the switching process. TIMS have been observed for [\(\text {Tb}_1\)/\(\text {Co}_3\)] and for [\(\text {Tb}_1\)/\(\text {Co}_2\)] with our model. However, Fig. 2d–f evidence that the fluence switching window rapidly disappears for \(m \ge \) 5. This is a strong indication that the magnetization switching in the [\(\text {Tb}_n\)/\(\text {Co}_m\)] requires a non-negligible number of interfaces with antiferromagnetic exchange in order to fully reverse the magnetization after each single laser pulse. As a matter of fact, the high number of antiferromagnetic interfaces present in [\(\text {Tb}_1\)/\(\text {Co}_3\)] system supports the theory according to which the switching originates at one interface and is followed by a propagation of the switching front towards the inner part of the layers comprising multilayer until a complete switching is achieved. The fluence threshold to trigger the switching increases with the number of Co monolayers, which is consistent with the increase of the Curie temperature. On the other hand, it is also important to notice that [\(\text {Tb}_1\)/\(\text {Co}_2\)] and [\(\text {Tb}_1\)/\(\text {Co}_3\)], whose Curie temperature and compensation temperature (around 750 K) are much higher than room temperature are the only systems that showed a clear fluence window where TIMS can be observed with our model. These results are in agreement with the observation of Moreno et al.14, in \(\text {Tb}_{x} \text {Co}_{1-x}\) alloys, in which systems with a compensation temperature higher than room temperature undergo thermally-induced magnetization switching in a wider fluence window14.
There is also a remarkable difference between the magnetization dynamics behavior of the sublattices in the [\(\text {Tb}_1\)/\(\text {Co}_3\)] and the [\(\text {Tb}_3\)/\(\text {Co}_9\)] systems as a function of the laser fluence shown in Fig. 2c,d. Although, these two systems have the same proportion of Tb and Co thicknesses, the switching window was fully suppressed, and additionally the response of the magnetization upon photo-excitation requires at least 60 mJ/\(\text {cm}^2\) to induce partial switching. This observation highlights and confirms the key role of the interface with antiferromagnetic exchange coupling in the switching mechanism of multilayer structures24.
Once the fluence switching window was obtained, further magnetization dynamic simulations were carried out in order to check whether thermally-induced magnetization switching can be fully deterministic. Although we have obtained a clear fluence switching window for [\(\text {Tb}_1\)/\(\text {Co}_3\)] and for [\(\text {Tb}_1\)/\(\text {Co}_2\)], additional simulated events showed that for some values between 35 and 50 mJ/\(\text {cm}^2\), thermally-induced magnetization switching is not fully deterministic. Figure 3a shows the normalized magnetization of Tb and Co as a function of the laser fluence for [\(\text {Tb}_1\)/\(\text {Co}_2\)], between 30 and 50 mJ/\(\text {cm}^2\) obtained from two different experiments. It can be seen that for fluences 42 mJ/\(\text {cm}^2\) and 44 mJ/\(\text {cm}^2\) the switching is not deterministic. Figure 3b plots the temporal dependence of the magnetization of the two different events for a fluence F = 40 mJ/\(\text {cm}^2\) and F = 42 mJ/\(\text {cm}^2\). The two events at each fluence value are indicated by the solid and dashed lines. On the magnetization curves in Fig. 3b, the single laser pulses were sent at time t = 7 ps. During the first picosecond after the heating, the single laser pulse generates a strong non-equilibrium state in the Tb and Co sublattices originating the decoupling of the sublattices and a subsequent demagnetization at different rates. It is observed in the inset of Fig. 3b that the magnetization of the Co sublattice crosses zero before the one of the Tb sublattice, due to the demagnetization of Co sublattice is faster than that of the Tb sublattice. The total magnetization of the system, including the transient ferromagnetic state created during the next picosecond is mediated by the antiferromagnetic exchange interaction between the Tb and Co spins. Due to this interaction and for a small window of time (\(t=\) 7.8–8.2 ps), the magnetic moments of the Co sublattice grow and align parallel with respect to the magnetic moments of the Tb sublattice. Finally, Tb magnetic moments will also cross zero and the remagnetization process (\(t>\) 8.2 ps) will reorient the total magnetization in an antiparallel magnetic polarity with respect to the initial state. The transient ferromagnetic state, is characteristic of a thermally induced magnetization switching as was already discussed by Radu et al.25 and now also confirmed by our simulations. As can be seen in the right part of plots in Fig. 3b, the electronic temperature increases up to 1600 K and rapidly decreases close to 500 K before t = 8 ps. The temporal evolution of the electronic temperature during the simulations indicates that the maximum temperature (1600 K) of the electronic bath is higher than the Curie temperature of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] multilayer (bottom inset of Fig. 3c).
(a) Normalized magnetization of Tb and Co as a function of the laser fluence. \(M_z/M_s\) ratios were extracted from the magnetization curves of two different simulations of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] system. (b) Time evolution of the magnetization of Tb and Co sublattices, \(\text {M}^\text {Co}\) (blue) and \(\text {M}^\text {Tb}\) (red), after the application of a single laser pulse at t = 7 ps in the system [\(\text {Tb}_1\)/\(\text {Co}_2\)] for two different fluences: 40 mJ/\(\text {cm}^2\) (top) and 42 mJ/\(\text {cm}^2\) (bottom). The two different events for each fluence value are represented by the solid and dashed lines. Inset: \(M_Z\)/\(M_S\) during 7.6 \(<t<\) 8.4 ps. Evolution of the electronic temperature after the single laser pulse interacts with the system at t = 7 ps (pink background). The single laser pulse duration is 50 fs. (c) Magnetization of Tb and Co sublattices in a magnetization phase map. Upper inset: Pathways followed by the magnetizations in a vicinity of zero. The magnetic moment of the Co sublattice always crosses zero before the magnetic moment of the Tb sublattice. Lower inset: Magnetization of the sublattices as a function of the temperature. The effective Curie temperature of the system is close to 750 K.
In order to visualize the evolution of the magnetization during the ferromagnetic transient regime, we plot the magnetic moment of the Tb and Co sublattices in a magnetization phase diagram (Fig. 3c) based on the model proposed by Davies et al.26,27. The top inset of Fig. 3c zooms in the trajectories followed by the Tb and Co magnetic moments during the demagnetization and reorientation processes within the third quadrant of the diagram phase (corresponding to the ferromagnetic transient state). In accordance with the duration of the transient ferromagnetic state extracted from the inset of Fig. 3b (0.4 ps), the pathways that follow the magnetic moments of the sublattices are contained within the values 0 < \(\text {M}^\text {Co}_z\)/\(\text {M}_s<\) 0.05 and 0 < \(\text {M}^\text {Tb}_z\)/\(\text {M}_s<\) 0.05 for a very short period of time. This also indicates that the Tb sublattices only posses a small magnetic moment value to transfer to the Co sublattice once the Co magnetic moment reaches zero. On the other hand, comparing the duration of the transient ferromagnetic state of our [\(\text {Tb}_1\)/\(\text {Co}_2\)]: 0.4 ps, with the one experimentally observed in GdFeCo layers that is around 1 ps25, or with the one extracted from numerical simulations in CoGd-based systems, also around 1 ps, it might be expected an increase of the switching probability in the [Tb/Co]-based system by increasing the duration of the transient ferromagnetic state. Experimentally, this might be achieved by altering the demagnetization rate though the tunning of the magnetic anisotropy during the growth of the Co and Tb layers.
As shown above, a fully reversal of the magnetization is not always observed for repeated attempts at a given fluence. Figure 3b showed the magnetization dynamics upon the incidence of a single laser pulse with a fluence F = 42 mJ/\(\text {cm}^2\), in which the final state of the magnetization of the Co and Tb sublattices does not show switching (dashed lines) after 20 ps. To further investigate the stochasticity of the magnetization reversal, we calculated the switching probability in the [\(\text {Tb}_1\)/\(\text {Co}_2\)] multilayer by repeating 20 times the simulation of the magnetization dynamics. Figure 4a illustrates the switching probability for the [\(\text {Tb}_1\)/\(\text {Co}_2\)] system, which has a maximum probability at 35 mJ/\(\text {cm}^2\). It can also be observed that for the fluence window between 35 and 50 mJ/\(\text {cm}^2\), indicated by the gray rectangle the switching probability is above 80%. For fluence values above this reversal window, the stochasticity increases, while the switching probability rapidly decreases for higher fluence values. We also explored the influence of the damping constant value (\(\lambda \)) in the switching probability of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] multilayer. By decreasing the damping constant from \(\lambda \) = 0.05 to \(\lambda \) = 0.03 in the Landau-Lifshitz-Gilbert equation, we were able to adjust the switching probability and the fluence window of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] system. As can be seen in Fig. 4a, the switching probability as a function of the laser fluence over 20 events increased to 100% for a fluence window between 30 and 50 mJ/\(\text {cm}^2\) for \(\lambda \) = 0.03. The decrease of the damping constant impacts the remagnetization time increasing the switching probability of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] system. Although, we have used the same damping constant value for Tb and Co atoms to determine the switching probability, as was done in similar works14,25, the exploration of element-specific damping could provide more information about the switching dynamics of the system as suggested by recent reports28. Our simulation and experiments12 were performed using fs-long pulses to evaluate the optical switching in [Tb/Co]-based systems and only [\(\text {Tb}_1\)/\(\text {Co}_2\)] and [\(\text {Tb}_1\)/\(\text {Co}_3\)] showed a clear fluence window for the observation of optical switching using 50 fs pulses (Fig. S2 in the Supplementary data file). The possibility to switch the magnetization in these systems with ps-long laser pulses using atomistic simulations remains unexplored.
(a) Switching probability as a function of the laser fluence for the system [\(\text {Tb}_1\)/\(\text {Co}_2\)] over 20 events with two different values of damping constant \(\lambda \) = 0.03 and 0.05. There is a maximum of the switching probability for F = 35 mJ/\(\text {cm}^2\) for the case of \(\lambda \) = 0.05. (b) Switching probability as a function of the laser fluence for the system [\(\text {Tb}_4\)/\(\text {Co}_8\)] and for a homogeneous TbCo alloy with 66% of Co. Black and grey solid lines are a guide to the eyes.
As expected, the switching probability significantly drops for thicker Tb and Co layers, and similarly to the case of [\(\text {Tb}_1\)/\(\text {Co}_3\)] and [\(\text {Tb}_3\)/\(\text {Co}_9\)] discussed before, the magnetization dynamics of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] and [\(\text {Tb}_4\)/\(\text {Co}_8\)] multilayers, which have the same proportion of Tb and Co thickness, also presents completely different behaviors. The switching probability is below 60% and the minimum fluence value to excite or induce magnetization reversal is around 60 mJ/\(\text {cm}^2\), as can be seen in Fig. 4b. The section S3 of the Supplementary information file shows the normalized magnetization of Tb and Co as a function of the laser fluence and also the time evolution of both magnetization after the laser heats the sample. Contrary to the switching of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] system, the results obtained in the multilayer [\(\text {Tb}_4\)/\(\text {Co}_8\)] showed that when there is no controlled magnetization switching, as a result of the noise introduced by the thermal fluctuations in the two-temperature model, the construction of the trajectories in the magnetization phase diagram does not follow the model suggested by Davies et al.26,27 (see Fig. S3c in the supplementary data file).
In order to assess the hypothesis that the antiferromagnetic exchange coupling is relevant to the TIMS process, we investigated the influence of the intermixing in the alloy structure and calculated the switching probability in a TbCo alloy system with 66% of Co. Fig. S4 of the supplementary data file shows the distribution of the Co and Tb atoms for the ideal sample [\(\text {Tb}_4\)/\(\text {Co}_8\)] as well as that of the samples affected by the intermixing at the interfaces. Figure 4b shows that the dependence of the thermally-induced magnetization switching probability as a function of the fluence is comparable to that of the [\(\text {Tb}_1\)/\(\text {Co}_2\)] multilayer system. The study performed in the alloy with 66% of Co in the TbCo alloy can also be compared with the one of 68% of Co in the TbCo alloy reported in Ref. 11, for which a narrow switching window is expected at a much lower pulsed laser fluence. On the other hand, the enhancement of the probability in the alloy case is also in agreement with other reports highlighting the role of intermixing in synthetic ferrimagnets29.
The results presented in this work indicate that not only thermal fluctuations but also the values of the damping constant and the intermixing are critical and expected to play an important role in the successful optical switching of ultra-thin [\(\text {Tb}_n\)/\(\text {Co}_m\)]. Further studies will be needed to understand the influence of these parameters in detail. On the other hand, signatures of stochastic switching have been found in preliminary electrical test of magnetic junctions with Tb/Co-based multilayer electrodes (Sec. S5 of the supplementary data file), therefore additional experiments that fully support our simulations need to be performed.
In conclusion, we applied an atomistic spin model to [Tb/Co] multilayered structures and calculated the thermal variation of its magnetization and the magnetization dynamics upon incidence of single laser pulses. After introducing the physical influence of the laser pulse into the model, the all-optical switching of [Tb/Co] multilayers was simulated. We found that [\(\text {Tb}_1\)/\(\text {Co}_2\)] and [\(\text {Tb}_1\)/\(\text {Co}_3\)] present a fluence window between 35 and 50 mJ/\(\text {cm}^2\) in which the multilayers can be switched by a single laser pulse, presenting an error rate below 5%. Stochasticity of the switching process was determined by the calculation of the probability over 20 switching attempts. We also pointed out the role of the damping constant value and the intermixing in the switching occurrence. Lower values of damping constant and intermixed interfaces enhance the switching probability in our system. These results predict for the first time the possibility of TIMS in ultra-thin [Tb/Co] multilayers by using atomistic spin simulations and suggest the occurrence of sub-picosecond magnetization reversal using single laser pulses in a [Tb/Co]-based multilayer structures. The deterministic switching and the estimation of the required fluence are in agreement with recent experimental observations of single-shot all-optical switching in [Tb/Co]\(_5\)-based multilayers11,12. The prediction of single-shot switching in [Tb/Co] multilayers based on numerical simulation results and its confirmation by experimental observations, show that this type of simulations can be used as a guiding tool to find appropriate systems exhibiting single pulse all-optical switching. The physical mechanisms underlying thermally induced magnetization switching in [\(\text {Tb}_1\)/\(\text {Co}_2\)] is in agreement with the model proposed by Davies et al.26,27 Furthermore, in thicker layers the demagnetizing processes dominate, the magnetization of the two ferromagnetic layers is divided in small magnetic domains preventing thus the complete reversal.
Stanciu, C. D. et al. All-optical magnetic recording with circularly polarized light. Phys. Rev. Lett. 99, 047601. https://doi.org/10.1103/PhysRevLett.99.047601 (2007).
Article ADS CAS PubMed Google Scholar
Kimel, A. V. Three rules of design. Nat. Mater. 13, 225. https://doi.org/10.1038/nmat3886 (2014).
Mangin, S. et al. Engineered materials for all-optical helicity-dependent magnetic switching. Nat. Mater. 13, 286–292. https://doi.org/10.1038/nmat3864 (2014).
El-Ghazaly, A. et al. Ultrafast magnetization switching in nanoscale magnetic dots. Applied Physics Letters 114, 232407. https://doi.org/10.1063/1.5098453 (2019).
Beens, M., Lalieu, M. L. M., Deenen, A. J. M., Duine, R. A. & Koopmans, B. Comparing all-optical switching in synthetic-ferrimagnetic multilayers and alloys. Phys. Rev. B 100, 220409. https://doi.org/10.1103/PhysRevB.100.220409 (2019).
Khorsand, A. R. et al. Role of magnetic circular dichroism in all-optical magnetic recording. Phys. Rev. Lett. 108, 127205. https://doi.org/10.1103/PhysRevLett.108.127205 (2012).
Kimel, A. V. et al. Ultrafast non-thermal control of magnetization by instantaneous photomagnetic pulses. Nature 435, 655–657. https://doi.org/10.1038/nature03564 (2005).
Hadri, M. S. E., Hehn, M., Malinowski, G. & Mangin, S. Materials and devices for all-optical helicity-dependent switching. J. Phys. D Appl. Phys. 50, 133002. https://doi.org/10.1088/1361-6463/aa5adf (2017).
Gorchon, J. et al. Single shot ultrafast all optical magnetization switching of ferromagnetic Co/Pt multilayers. Appl. Phys. Lett. 111, 042401. https://doi.org/10.1063/1.4994802 (2017).
Lalieu, M. L. M., Peeters, M. J. G., Haenen, S. R. R., Lavrijsen, R. & Koopmans, B. Deterministic all-optical switching of synthetic ferrimagnets using single femtosecond laser pulses. Phys. Rev. B 96, 220411. https://doi.org/10.1103/PhysRevB.96.220411 (2017).
Avilés-Félix, L. et al. Integration of Tb/Co multilayers within optically switchable perpendicular magnetic tunnel junctions. AIP Adv. 9, 125328. https://doi.org/10.1063/1.5129821 (2019).
Avilés-Félix, L. et al. Single-shot all-optical switching of magnetization in Tb/Co multilayer-based electrodes. Sci. Rep. 10, 5211. https://doi.org/10.1038/s41598-020-62104-w (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Vomir, M., Albrecht, M. & Bigot, J.-Y. Single shot all optical switching of intrinsic micron size magnetic domains of a Pt/Co/Pt ferromagnetic stack. Appl. Phys. Lett. 111, 242404. https://doi.org/10.1063/1.5010915 (2017).
Moreno, R., Ostler, T. A., Chantrell, R. W. & Chubykalo-Fesenko, O. Conditions for thermally induced all-optical switching in ferrimagnetic alloys: Modeling of TbCo. Phys. Rev. B 96, 014409. https://doi.org/10.1103/PhysRevB.96.014409 (2017).
Cheng, W., Li, X., Wang, H., Cheng, X. & Miao, X. Laser induced ultrafast magnetization reversal in TbCo film. AIP Adv. 7, 056018. https://doi.org/10.1063/1.4975659 (2017).
Alebrand, S. et al. Light-induced magnetization reversal of high-anisotropy TbCo alloy films. Appl. Phys. Lett. 101, 162408. https://doi.org/10.1063/1.4759109 (2012).
Alebrand, S. et al. Subpicosecond magnetization dynamics in TbCo alloys. Phys. Rev. B 89, 144404. https://doi.org/10.1103/PhysRevB.89.144404 (2014).
Zhang, G. P., Murakami, M., Si, M. S., Bai, Y. H. & George, T. F. Understanding all-optical spin switching: Comparison between experiment and theory. ArXiv e-prints (2018). arXiv:1808.06875.
Ostler, T. A. et al. Ultrafast heating as a sufficient stimulus for magnetization reversal in a ferrimagnet. Nat. Commun. 3, 666. https://doi.org/10.1038/ncomms1666 (2012).
Anisimov, S. I., Kapeliovich, B. L. & Perel'man, T. L. Electron emission from metal surfaces exposed to ultrashort laser pulses. Sov. Phys. JETP 39, 375 (1974).
Gridnev, V. N. Ferromagneticlike states and all-optical magnetization switching in ferrimagnets. Phys. Rev. B 98, 014427. https://doi.org/10.1103/PhysRevB.98.014427 (2018).
Gridnev, V. N. Role of compensation points in all-optical magnetization switching. Phys. Rev. B 100, 174405. https://doi.org/10.1103/PhysRevB.100.174405 (2019).
Cornelissen, T. D., Córdoba, R. & Koopmans, B. Microscopic model for all optical switching in ferromagnets. Appl. Phys. Lett. 108, 142405. https://doi.org/10.1063/1.4945660 (2016).
Xu, C., Ostler, T. A. & Chantrell, R. W. Thermally induced magnetization switching in Gd/Fe multilayers. Phys. Rev. B 93, 054302. https://doi.org/10.1103/PhysRevB.93.054302 (2016).
Radu, I. et al. Transient ferromagnetic-like state mediating ultrafast reversal of antiferromagnetically coupled spins. Nature 472, 205. https://doi.org/10.1038/nature09901 (2011).
Davies, C. et al. Pathways for single-shot all-optical switching of magnetization in ferrimagnets. Phys. Rev. Appl. 13, 024064. https://doi.org/10.1103/PhysRevApplied.13.024064 (2020).
Mentink, J. H. et al. Ultrafast spin dynamics in multisublattice magnets. Phys. Rev. Lett. 108, 057202. https://doi.org/10.1103/PhysRevLett.108.057202 (2012).
Ceballos, A. et al. Role of element-specific damping on the ultrafast, helicity-independent all-optical switching dynamics in amorphous (Gd,Tb)Co thin films. arXiv: Materials Science (2019). arXiv:1911.09803.
Beens, M., Lalieu, M. L. M., Duine, R. A. & Koopmans, B. The role of intermixing in all-optical switching of synthetic-ferrimagnetic multilayers. AIP Adv. 9, 125133. https://doi.org/10.1063/1.5129892 (2019).
This research has received funding from the European Union's Horizon 2020 research and innovation program under FET-Open Grant Agreement No. 713481 (SPICE).
L. Avilés-Félix
Present address: Centro Atómico Bariloche, Comisión Nacional de Energía Atómica (CNEA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Avenida Bustillo 9500, Bariloche, Rio Negro, Argentina
Univ. Grenoble Alpes, CEA, CNRS, Grenoble INP, IRIG-SPINTEC, 38000, Grenoble, France
L. Avilés-Félix, L. Farcis, Z. Jin, L. Álvaro-Gómez, B. Dieny, R. C. Sousa, I. L. Prejbeanu & L. D. Buda-Prejbeanu
Institute for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ, Nijmegen, The Netherlands
G. Li, K. T. Yamada, A. V. Kimel & Th. Rasing
FELIX Laboratory, Radboud University, 7 Toernooiveld, 6525 ED, Nijmegen, The Netherlands
A. Kirilyuk
L. Farcis
Z. Jin
L. Álvaro-Gómez
K. T. Yamada
A. V. Kimel
Th. Rasing
B. Dieny
R. C. Sousa
I. L. Prejbeanu
L. D. Buda-Prejbeanu
Z.J. and L.D.B.P developed the numerical solver, L.F., Z.J. and L.A.G. performed the numerical simulations; L.A.F., L.F., I.L.P., L.D.B.P., B.D. and R.C.S. discussed and analyzed the numerical results; L.A.F., G.L and K.Y performed the electrical measurements; L.A.F., G.L., K.Y., A.K., A.V.K., Th.R., I.L.P., R.C.S. and L.D.B.P. discussed and analyzed the electrical measurements, L.A.F. wrote the manuscript with contributions from I.L.P, R.C.S, Th.R, B.D. and L.D.B.P.; I.L.P., R.C.S. B.D. and L.D.B.P. supervised the project.
Correspondence to L. Avilés-Félix.
Supplementary Information.
Avilés-Félix, L., Farcis, L., Jin, Z. et al. All-optical spin switching probability in [Tb/Co] multilayers. Sci Rep 11, 6576 (2021). https://doi.org/10.1038/s41598-021-86065-w
Dynamic complex opto-magnetic holography
M. Makowski
J. Bomba
A. Stupakiewicz
About Scientific Reports
Guide to referees