Abstract
Controlling the perpendicular magnetic anisotropy (PMA) in thin films has received considerable attention in recent years due to its technological importance. PMA based devices usually involve heavy-metal (oxide)/ferromagnetic-metal bilayers, where, thanks to interfacial spin-orbit coupling (SOC), the in-plane (IP) stability of the magnetisation is broken. Here we show that in V/MgO/Fe(001) epitaxial junctions with competing in-plane and out-of-plane (OOP) magnetic anisotropies, the SOC mediated interaction between a ferromagnet (FM) and a superconductor (SC) enhances the effective PMA below the superconducting transition. This produces a partial magnetisation reorientation without any applied field for all but the largest junctions, where the IP anisotropy is more robust; for the smallest junctions there is a reduction of the field required to induce a complete OOP transition (\(H_\text {OOP}\)) due to the stronger competition between the IP and OOP anisotropies. Our results suggest that the degree of effective PMA could be controlled by the junction lateral size in the presence of superconductivity and an applied electric field. We also discuss how the \(H_\text {OOP}\) field could be affected by the interaction between magnetic stray fields and superconducting vortices. Our experimental findings, supported by numerical modelling of the ferromagnet-superconductor interaction, open pathways to active control of magnetic anisotropy in the emerging dissipation-free superconducting spin electronics.
Similar content being viewed by others
Introduction
Control of out-of-plane (OOP) anisotropies in ultra thin ferromagnetic multilayer films have revolutionized magnetic storage and spintronics technologies by mitigating the impact of the demagnetizing energy as the bit and magnetic tunnel junction sizes diminished1,2. Tuning of perpendicular magnetic anisotropy (PMA) by careful selection of structure design3,4 and size5 has been among the main challenges of spintronics. Besides the variation of the ferromagnet thickness and interface with oxides, OOP magnetisation reorientation can be achieved by a temporary reduction of the IP-OOP barrier using, for example, heat and microwave pulses6,7 or a combination of magnetic field and low temperature8.
Recently, we demonstrated a fundamentally different route to magnetisation reorientation through the influence of superconductivity on the IP magnetisation anisotropy9. The key idea behind this effect is that the magnetisation aligns to minimize the weakening of the superconducting condensate associated with the creation of spin triplet (ST) Cooper pairs10. The spin triplet generation depends on the magnetisation direction relative to the interfacial Rashba spin-orbit field. Understanding the factors influencing this superconductivity-induced change of magnetic anisotropy is crucial for designing the next generation of cryogenic memories in the emerging field of superconducting spintronics, where control over non-volatile magnetisation states still remains a major challenge11,12,13,14.
The main underlying physical mechanisms for the transformation of ST Cooper pairs from singlet to mixed-spin and equal-spin triplet pairs are magnetic inhomogeneities15,16, two misaligned FM layers17,18 or SOC19. Previous experiments focusing on SOC-driven generation of triplets have focused on heavy metal (Pt) layers in non-epitaxial SC/FM structures20,21 and Rashba SOC in epitaxial V/MgO/Fe junctions9,22 where ST Cooper pairs are generated depending on the magnetisation orientation relative to the Rashba field.
Theoretically, it has been shown10 that a superconductor coupled to a ferromagnet by SOC could stimulate the modification not only of the IP9, but also of the OOP magnetic anisotropy below the superconducting critical temperature (\(T_C\)). Due to the stray fields, however, ferromagnetic films are expected to have a stronger interaction with the superconductor when an OOP magnetisation is present, compared to a simple IP variation23,24. Therefore, a careful consideration of the interaction of these stray field generated by the OOP magnetisation and superconducting vortices is essential to fully capture the factors influencing the effective OOP anisotropy.
Here, we investigate the superconductivity-induced OOP magnetisation reorientation in epitaxial Fe(001) films with competing IP and OOP anisotropies, both at zero field and in the presence of out-of-plane applied magnetic fields. The V/MgO/Fe(001) junctions are ideal candidates to verify the predicted modification of the effective perpendicular anisotropy in the superconducting state for several reasons10. Firstly, the Fe(001) has the required10 cubic symmetry; secondly, previous studies show that the normal state IP-OOP reorientation transition takes place at a well-defined critical field8; thirdly, the system has Rashba type SOC, which is responsible for the PMA in MgO/Fe25; fourthly, the relative contribution of the IP and OOP magnetisation anisotropies can be tuned by changing the junction lateral size, and SOC can be varied by applying an external electric field; finally, the change in magnetisation can be determined with high precision by studying the transport characteristics using a second magnetically hard Fe/Co layer which is magnetostatically decoupled from the soft Fe layer8.
For the smallest junctions, where the IP and OOP anisotropies strongly compete, we remarkably observe the full superconductivity-induced IP-OOP magnetisation reorientation predicted in Ref.10. This results in (i) a decreasing of the required field to induce the full IP-OOP transition below \(T_C\), which is not observed in bigger junctions; and (ii) a spontaneous increasing of the misalignment angle between the two FM layers below \(T_C\) in the absence of applied field, which is consistently observed in all but the largest junctions. These differences in the observed behaviour depending on the junctions dimensions are most likely due to the IP anisotropy becoming more dominant with increasing lateral size. We discard the magnetostatic interaction between supeconducting vortices and the FM layers as the main cause of the observed effects.
Results
Figure 1 shows the experimental configuration and the different types of OOP transition observed above the vanadium \(T_C\). Figure 1a shows the V(40 nm)/MgO(2 nm)/Fe(10 nm)/MgO(2 nm)/Fe(10 nm)/Co(20 nm) (N(SC)/FM1/FM2) junctions with a hard Fe/Co layer (FM2) sensing the magnetisation alignment of the 10 nm thick Fe(001) soft layer (FM1). Details about the sample growth, characterization and the experimental set-up are explained in the “Methods” section. All junctions were saturated with a 3 kOe IP magnetic field (see the alignment calibration procedure in Supplementary Material, Sect. 1) before each of the OOP magnetoresistance (TMR) measurements, in order to eliminate magnetic inhomogeneities from previous OOP measurements. All except one of the studied junctions showed OOP anisotropy below 3 kOe. On the right side of the vertical axes of Fig. 1b–d, we indicate the TMR values corresponding to the well established parallel (P), perpendicular out-of-plane (OOP) and antiparallel (AP) states for each sample, which are used to calibrate the angle between the two FM layers (\(\Delta \phi =\phi _\text {FM1}-\phi _\text {FM2}\), where \(\phi _\text {FM1}\) and \(\phi _\text {FM2}\) are the angles of each FM layer with respect to the plane of the layers, as shown in Fig. 1a) with the same procedure as described in Refs.8,9. This indicates that the IP-OOP transition also triggers a total or partial reorientation of the sensing (hard) FM2 layer, providing a resistance close to that of an AP state. Previous OOP measurements8 above \(T_C\) made on only two \(20\times 20\) \(\upmu \text {m}^2\) junctions revealed asymmetric transitions into the perpendicular alignment of the soft FM1 layer, without any subsequent transition of the sensing layer with perpendicular fields up to 3 kOe. The present study is made with a total of 16 junctions of four different lateral sizes, where about half of them also demonstrate a transition to an AP configuration when the magnetic field is further increased after the transition to the OOP state has been completed. This AP configuration could potentially be either with the two layers oriented OOP or IP, although it seems rather unlikely that both layers reorient to an IP configuration while the applied OOP field increases. We believe that the high-field-induced transition from OOP to AP alignment or, in some cases, a nearly direct P to AP transition in N/FM1/FM2 junctions could be a consequence of the enhanced antiferromagnetic coupling reported for MgO magnetic tunnel junctions with perpendicular magnetic anisotropy (see26). We cannot exclude the possibility that the AP alignment could be triggered by a partial reorientation of the hard Fe/Co layer (with only the Fe part or the atomic layers closer to the Fe/MgO interface in the hard layer orienting OOP, as shown in the sketches in Fig. 1b,c). However, since we measure the total resistance of the junctions, it is impossible to distinguish between these two cases from transport measurements alone. Therefore, we mainly focus on the influence of superconductivity on the transition between IP and OOP states and the partial OOP reorientation at zero magnetic field. Consequently, for the OOP field range reported here, we will assume that the \(\phi _\text {FM2}\) angle of the FM2 layer with respect to the in-plane configuration is fixed and close to 0.
Figure 1c shows typical OOP TMR cycles measured in two \(20\times 20\) \(\upmu \text {m}^2\) junctions, one of them switching to an AP alignment following an OOP orientation (blue) and the other one only switches to the OOP state (red). Figure 1d shows an OOP TMR for a \(30\times 30\) \(\upmu \text {m}^2\) junction where the OOP alignment of the FM1 and FM2 electrodes remains stable up to 3 kOe. Note that all junctions showed remanent OOP alignment of the soft Fe(001) electrode once the perpendicular magnetic field is removed (Fig. 1). This indicates the relatively small number of interfacial defects present in our junctions, as supported by numerical simulations analyzing the OOP configuration robustness as a function of the density of interfacial defects by studying the inverse OOP to IP transition, which are discussed in the Supplementary Material, S2.2.
The symmetry broken spin reorientation observed in the OOP TMR experiments shown in Fig. 1 b–d, has been previously explained in Ref.8 by the difference in the dislocation density present at the top and bottom surfaces of the soft Fe(001) layer due to the growth process. This differently affects the top and bottom surface anisotropies, which leads to different intensities at each interface resulting in the magnetisation being more easily reoriented into the OOP configuration for one field direction than the other. This asymmetric field behaviour might seem at odds with the Stoner–Rashba model developed in Ref.27. This model suggests that a net Rashba field related to the asymmetric top and bottom interfaces of a ferromagnetic film leads to a pseudo-dipolar contribution to the anisotropy, which would mainly favor an in-plane magnetisation and an uniaxial-like anisotropy favouring the/a perpendicular magnetisation configuration. Correspondingly, the hysteresis curve of a single magnetic (here Fe) layer is expected to be an even function with respect to the external magnetic field. However, we note that the model does not fully account for the complexities discussed below that could lead to the asymmetric hysteresis we observe in our multilayer structures.
The fact that the hysteresis curve is not an even function of the external magnetic field is simply related to the fact that the model is developed for a single ferromagnetic layer while in our complex heterostructure, we do not reverse the Fe/Co interface magnetisation. This is not unreasonable considering a large interface anisotropy. A full magnetisation reversal including interfacial magnetisation would only result in an asymmetric hysteretic response. Secondly, in our structures stray fields play a relevant/crucial/central role and importantly the stray fields seen by both interfaces are not similar. The bottom interface experiences the stray field of the Fe/Co top bilayer, while the top interface sees the contribution from the bottom Fe layer. In a macrospin model, increasing the stray fields would decrease the perpendicular anisotropy. To fully understand the complexities of the asymmetric magnetisation response, future studies such as direct OOP magnetisation measurements on the MgO/Fe/MgO structures in the absence of the sensing Fe/Co and V/MgO layers, could be performed.
It is worth mentioning a distinct feature of our junctions, having a strongly preferred IP magnetisation at room temperature8, with the OOP configuration of the soft 10 nm Fe layer only becoming non-volatile below 80 K. In the temperature range in which this study takes place (0.3 to 7 K), the magnetic field required to induce an OOP transition in the soft layer does not typically exceed 2 kOe. These relatively low values (with respect to continuous 10 nm thick Fe films) could be explained by the combined influence of a few factors. Firstly, the variation of the relation between the IP and OOP anisotropy energies could vary with temperature, possibly favouring the OOP configuration at low temperatures28. Secondly, interfacial strain has also demonstrated the potential to induce changes in the perpendicular anisotropy in thin ferromagnetic films29. Thirdly, the IP saturation in this study was carried out with a field of 3 kOe. This value was considered sufficiently high since the resistance values were stable above 1 kOe, but it could be insufficient to induce a perfect IP alignment at low temperatures. This factor could be more relevant for the smallest junctions where edge magnetic charges would have a relatively higher influence on the measured OOP switching field, qualitatively explaining the dependence of this field with the junctions lateral size, as supported by numerical simulations (see Supplementary Material Fig. S5). Finally, as mentioned before, as long as we measure the total resistance of the junction, we can’t exclude that the OOP reorientation might take place preferently in the atomic layers closer to the Fe/MgO interface (where it would be easier to reorient the magnetic moments due to the surface anisotropy). Thus, the surface OOP state (with a thickness of a few nm, close to that of the Fe magnetic exchange length30) might be realized with the aid of interface anisotropy at the Fe/MgO interface and an external OOP magnetic field. This is shown in the sketches of the spin valve configuration in Fig. 1b–d.
Superconductivity induced change of the out-of-plane anisotropy field
Figure 2a shows the temperature dependence of OOP TMRs in a \(10\times 10\) \(\upmu \text {m}^2\) junction, in a field range where the field-induced magnetisation reorientation of the FM1 layer (measured at 5 mV) takes place. A decrease of the characteristic \(H_{\text {OOP}}\) field (defined as the applied magnetic field providing a complete OOP reorientation of the FM1 layer) just below \(T_C\) can be observed upon lowering the temperature, as represented in Fig. 2c. We note that the SC-induced full IP-OOP transitions have been clearly observed in the smallest \(10\times 10\) \(\upmu \text {m}^2\) lateral size junctions. The larger junctions showed a small low field TMR increase below \(T_C\), which could be interpreted either as a partial FM1 layer reorientation or an inhomogeneous OOP alignment (Fig. 2b). For the \(20\times 20\) \(\upmu \text {m}^2\) and larger junctions, the \(H_{\text {OOP}}\) anisotropy field turned out to be nearly independent of temperature (Fig. 2c). Interestingly, our junctions also revealed spontaneous zero field TMR emerging below \(T_C\) (corresponding to a partial magnetic reorientation of the soft FM1 layer), which is more pronounced for the smaller samples and diminishes with lateral size, abruptly disappearing for the largest junctions. This is shown in Fig. 2d, where instead of the TMR, the calculated angle between the two FM layer is plotted. It is worth noting that this relative angle calculation is similar to our previous work8,9, and assumes a uniform magnetisation in the whole FM layer. However, the real scenario could be more complex (see Supplementary Material, S2.1).
Influence of electric field on the out-of-plane reorientation
The presence of the MgO barriers allows us to explore the possible influence of high electric fields on the magnetic-field-induced IP-OOP transitions above and below \(T_C\). High electric field influences the PMA anisotropy by modifying the SOC Rashba field in magnetic tunnel junctions27. Our previous study8 revealed that roughly two thirds of the voltage drop in our junctions occurs at the V/MgO/Fe barrier, resulting in a high electric field across this interface. The remaining voltage drops at the Fe/MgO/Fe interface, which is responsible for the change in the resistance providing the measured TMR depending on the relative magnetic configuration of the two FM layers.
We have therefore investigated the influence of high bias and its polarity on the IP-OOP transition in junctions with different lateral sizes. Figure 3a–d show that an applied bias of 600 mV (generating an electric field at the V/MgO/Fe interface exceeding \(2\times 10^8\) V/m) hardly affects \(H_{\text {OOP}}\) above \(T_C\), independently of the junctions size. Moreover, the application of a large electric field has also a negligible effect on the superconductivity-induced IP-OOP transition in the larger than \(30\times 30\) \(\upmu \text {m}^2\) junctions, with a dominant IP magnetisation alignment (Fig. 3c). However this changes for the smaller junctions, where IP and OOP anisotropy values are comparable leading to an entirely different behaviour. Strikingly, we observe that for \(10\times 10\) \(\upmu \text {m}^2\) and \(20\times 20\) \(\upmu \text {m}^2\) junctions, the electric field stimulates an IP-OOP transition below \(T_C\) at very small values of the applied magnetic field (below 100 Oe).
Figure 3d compares the influence of an electric bias close to 600 mV with different polarities on the magnetisation alignment below \(T_C\) (0.3 K) with an applied magnetic field of \(-\,50\) Oe, within the field range in which we observed a larger influence of the electric field on the IP-OOP transition for the smaller junctions. This field is about an order of magnitude below the first critical field of our Vanadium films, which was estimated to be close to 400 Oe22, therefore minimizing the presence of vortices in the superconducting layer. We believe that the electric field effect asymmetry could be due to the combined influence of the relatively more dominant proximity effects between the SC and FM states at the V/MgO/Fe interface in smaller junctions, and the electric-field-induced variation of the Rashba field influencing the OOP anisotropy for the non-equivalent interfaces MgO/Fe and Fe/MgO in the junctions27.
Discussion
Evaluation of magnetostatic coupling between superconducting vortices and ferromagnet
Let us start our discussion by considering different scenarios involving the possible magnetostatic coupling between the superconducting vortices and the ferromagnet31. It is tempting to consider the device edges as mainly responsible for the superconductivity-induced spin reorientation, as the edge has a more important contribution for the smallest samples, in which the minimum applied field is enough to fully reorient the magnetisation. However, a few experimental facts contradict this scenario. Firstly, the superconductivity-induced additional zero field OOP angle variation is similar for \(10\times 10\) to \(30\times 30\) \(\upmu \text {m}^2\) junctions (see inset in Fig. 2 d), which would not be the case if the change comes from the device’s edges. The superconductivity induced spin reorientation effect abruptly diminishes for the \(40\times 40\) \(\upmu \text {m}^2\) junction only (Fig. 2d). Secondly, numerical simulations show that the OOP reorientation due to magnetostatic coupling, if relevant, could potentially be triggered by the nucleation of OOP domains in the interior of the samples rather than at the edges; even if we assume the edges as the initial OOP nucleation places, the resulting vortex distribution would affect the whole FM layer (see Supplementary Material S2.1). Finally, electric field stimulates the OOP transition for relatively small junctions with competing anisotropies (see Fig. 3) which points towards the possible role of the Rashba field.
We have seen from micromagnetic simulations and as an experimental trend that, on average, the normal state \(H_{\text {OOP}}\) increases with the junctions area (Supplementary Material S2.3). This is in agreement with the gradual decrease of the partial OOP magnetisation reorientation with increasing lateral size seen in the normal state, just above the critical temperature (see Fig. 2d). Within the above picture, a lower \(H_{\text {OOP}}\) field is required to reorient the magnetisation perpendicularly in the smallest junctions, and therefore one would expect a weaker magnetostatic coupling to SC vortices.
Numerical simulations of the magnetostatic interaction of the V/MgO/Fe system during an OOP TMR experiment such as the ones shown in Fig. 1, where a varying OOP magnetic field is applied, is a complex problem which requires self-consistent treatment of the interaction between magnetic charges and stray fields of superconducting vortices32. The Supplementary Material (Sect. 2.1) introduces a simplified simulation scheme which evaluates this interaction in the presence of the Meissner effect. These results show that the vortex-mediated magnetostatic interaction might only explain a weak enhancement of \(H_{\text {OOP}}\) in the superconducting state in the largest junctions (Fig. 2c). However, we note that varying the superconducting hysteresis strength or width in the magnetostatic simulations could not explain the strong decrease of \(H_{\text {OOP}}\) below \(T_C\) which was experimentally observed in the smaller junctions. Moreover, a dominant magnetostatic coupling would contradict the observed influence of electric field on TMR below \(T_C\) for the smallest junctions (Figs. 2, 3).
Microscopic model
To explain the strong decrease in the OOP anisotropy field below \(T_{C}\) for the smallest junctions and the superconductivity-induced zero field magnetic reorientation in all except the largest ones, as well as the influence of the SOC strength through the application of an electric field, we present a microscopic model describing the observed superconductivity-assisted OOP magnetic reorientation. In heterostructures consisting of superconducting and magnetic layers, the superconducting condensate is weakened as Cooper pairs leak into the magnetic regions33. This leakage is more efficient when the spin-singlets are transformed into equal-spin triplet pairs polarized along the same axis as the magnetisation. In our system, the MgO layer boosts the Rashba SOC at the SC/FM interface allowing for a generation of equal-spin triplets that depends on the orientation of the magnetisation with respect to the interface10,20.
To show how the efficiency of the triplet leakage affects the critical field for reorienting the magnetisation OOP, we calculate the free energy of the system from a tight-binding Bogoliubov–de Gennes (BdG) Hamiltonian (see Sect. 3 in the Supplementary Material for a complete description of our method). The V/MgO/Fe structure is modelled as a cubic lattice with electron hopping between neighboring sites. We include conventional s-wave on-site superconducting pairing potential in the V layer, Rashba SOC in the MgO layer, and an exchange splitting between spins in the Fe layer. Although this model is valid in the ballistic limit, we expect similar results for diffusive materials since spin singlets are partially converted into odd-frequency s-wave triplets that are robust to impurity scattering. Moreover, the variation in the singlet population under IP to OOP reorientation of the magnetisation have previously been demonstrated both experimentally and by dirty limit calculations20.
The free energy determined from this model captures the contribution from the superconducting proximity effect, and also includes a normal-state contribution favoring an IP magnetisation. In addition, we include a normal-state anisotropy \(K_{\text {IP}}[1-\cos ^4(\Phi _{\text {FM}1})]+K_{\text {OOP}}[1-\sin ^2 (\Phi _{\text {FM}1})]\), where \(\Phi _{\text {FM}1}\) ranges from \(0^{\circ }\) (corresponding to an IP magnetisation of the soft ferromagnet) to \(90^{\circ }\) (OOP magnetisation). In total this gives a normal-state anisotropy favoring an IP magnetisation, with an additional local minimum for the OOP magnetisation direction. Here we only focus on the superconductivity-assisted deepening of these OOP quasi-minima associated with the spin singlet to spin triplet conversion. The increase in \(H_{\text {OOP}}\) below \(T_{C}\) discussed in the previous section is not covered by this theoretical framework, as it does not take into account formation of vortices or the size of the junction. The discussion here is therefore relevant to the smaller junctions where the superconductivity-assisted decrease in \(H_{\text {OOP}}\) dominates.
In Fig. 4a), we demonstrate how the local free energy minimum for an OOP magnetisation deepens as the temperature is decreased below \(T_C\). As a simple qualitative model, we calculate the external magnetic field that can be used to force the magnetisation into the OOP orientation as \(H_{\text {OOP}}=(K_{\text {anis}}+F_{\text {OOP}}-F_{\text {IP}})/\mu _0 \mu _{\text {tot}}\), where \(\mu _{\text {tot}}\) is the total magnetic moment, \(K_{\text {anis}}\) is a constant anisotropy favoring the IP orientation that includes the above mentioned parameters for the normal-state IP and OOP anisotropies, \(K_\text {IP}\) and \(K_\text {OOP}\), as well as an energy barrier associated with the reorientation; and \(F_{\text {OOP}}\) and \(F_{\text {IP}}\) are the calculated free energies in the OOP and IP states of the soft layer respectively. In Fig. 4b, we show how \(H_{\text {OOP}}\) decreases below \(T_C\) as observed for the \(10\times 10\) \(\upmu\hbox {m}^2\) junction in Fig. 2c. We have thus demonstrated that the proximity effect enables a strong decrease in \(H_{\text {OOP}}\) that cannot be explained by the coupling of the ferromagnet to superconducting vortices discussed in the previous section. Moreover, since this variation in \(H_{\text {OOP}}\) requires that SOC is present, it also explains the dependence on the electric field observed for the smaller junctions (Fig. 3). The fact that \(H_{\text {OOP}}\) decreases over a longer temperature interval than in the experiments, rather than flattening out for low T, is caused by the downscaling of the lattice that is necessary in our theoretical model. In order to scale down the superconducting coherence length so that it remains comparable to the thickness of the superconducting layer, the on-site interaction must be increased, leading to a higher \(T_C\). Since the temperature interval is larger, a smaller fraction of the temperatures exist in the low-temperature limit where the free energy is temperature independent. Keeping in mind that our measurements of \(H_{\text {OOP}}\) show a dependence on the magnitude of the Rashba SOC, we can conclude that the SOC induced change in magnetic anisotropy below \(T_C\) shown here must strongly contribute to the suppression in \(H_{\text {OOP}}\) for the \(10\times 10\) \(\upmu\hbox {m}^2\) junctions.
Conclusions
Our experiments point towards the superconductivity induced modification of the perpendicular magnetic anisotropy in the epitaxial Fe(001) films in the V(40 nm)/MgO(2 nm)/Fe(10 nm) system. The behaviour depends on the lateral dimensions of the junctions in the following way: First, for the smallest junctions, the magnetic field necessary for a full OOP magnetisation reorientation drops by an order of magnitude in the superconducting state, while for the rest of the junctions it varies only slightly. Second, in all but the largest junctions, an increase in the OOP misalignment angle between the soft Fe(10 nm) layer and the hard one is observed when the temperature is decreased below \(T_C\) without any applied field. This spontaneous reorientation is similar for \(10\times 10\) to \(30\times 30\) \(\upmu\hbox {m}^2\) junctions and disappears in the largest ones, suggesting that superconductivity could be affecting the competition between the IP and OOP anisotropies (which is more pronounced for the smaller junctions) rather than being the result of the reorientation taking place at the edges of the samples. The decreasing of \(H_\text {OOP}\) transition field in the superconducting state, which could also be stimulated by the application of electric field changing the Rashba SOC, is consistent with the theoretical prediction10 of the absolute minimum of free energy corresponding to the OOP spin direction in SC/SOC/FM hybrids with competing (IP vs OOP) anisotropies just below \(T_C\). The magnetostatic interaction between vortices and magnetic inhomogeneities could explain a weak hardening of the OOP transition in the largest junctions. A detailed theoretical analysis of the mutual interplay between the inhomogeneous magnetisation of the soft ferromagnet and the superconductor is, however, beyond the scope of this work. Our results open a route to active manipulation of perpendicular magnetic anisotropy in the expanding field of dissipation-free superconducting electronics involving spin34,35,36 or spin polarized supercurrents37.
Methods
Samples growth and characterization
The V(40 nm)/MgO(2 nm)/Fe(10 nm)/MgO(2 nm)/Fe(10 nm)/Co(20 nm) MTJ multilayer stacks have been grown by molecular beam epitaxy (MBE) in a chamber with a base pressure of \(5\times 10^{-11}\) mbar following the procedure described in Ref.38. The samples were grown on (001) MgO substrates. A 10 nm thick seed of anti-diffusion MgO underlayer is grown on the substrate to trap the C from it before the deposition of the Fe (or V). The MgO insulating layer is then epitaxially grown by e-beam evaporation up to a thickness of approximately \(\sim 2\) nm and the same process is then executed for the rest of the layers. Each layer is annealed at 450 °C for 20 min for flattening. After the MBE growth, all the MTJ multilayer stacks are patterned in micrometre-sized square junctions by UV lithography and Ar ion etching, controlled step-by-step in situ by Auger spectroscopy.
Experimental measurement methods
The measurements are performed inside a \(\hbox {JANIS}^{\tiny {\textregistered }}\) \(\hbox {He}^3\) cryostat (the minimum attainable temperature is 0.3 K). The magnetic field is varied using a 3D vector magnet consisting of one solenoid (Z axis) with \(H_\text {max}=3.5\) T and two Helmholtz coils (X and Y axis) with \(H_\text {max}=1\) T. In our system the different magnetic states can be distinguished by looking at the resistance, so the relative orientation between two electrodes can be measured. The magnetoresistance measurements are performed by first setting the magnetic field to the desired value, then applying positive and negative current up to the desired voltage (5 mV unless otherwise stated), and averaging the absolute values of the measured voltage for the positive and negative current, obtaining a mean voltage which is used to calculate the resistance at that point. The temperature is measured and controlled with a LakeShore 340 thermometer.
References
Piramanayagam, S. N. Perpendicular recording media for hard disk drives. J. Appl. Phys. 102, 011301 (2007).
Dieny, B. & Chshiev, M. Perpendicular magnetic anisotropy at transition metal/oxide interfaces and applications. Rev. Mod. Phys. 89, 025008 (2017).
Chuang, T. C., Pai, C. F. & Huang, S. Y. Cr-induced perpendicular magnetic anisotropy and field-free spin-orbit-torque switching. Phys. Rev. Appl. 11, 061005 (2019).
Yi, D. et al. Enhanced interface-driven perpendicular magnetic anisotropy by symmetry control in oxide superlattices. Phys. Rev. Appl. 15, 024001 (2021).
Sun, J. Z. Resistance-area product and size dependence of spin-torque switching efficiency in CoFeB–MgO based magnetic tunnel junctions. Phys. Rev. B 96, 064437 (2017).
Challener, W. A. et al. Heat-assisted magnetic recording by a near-field transducer with efficient optical energy transfer. Nat. Photonics 3, 220–224 (2009).
Zhu, J.-G., Zhu, X. & Tang, Y. Microwave assisted magnetic recording. IEEE Trans. Magn. 44(1), 125–131 (2008).
Martínez, I., Tiusan, C., Hehn, M., Chshiev, M. & Aliev, F. G. Symmetry broken spin reorientation transition in epitaxial MgO/Fe/MgO layers with competing anisotropies. Sci. Rep. 8, 9463 (2018).
González-Ruano, C. et al. Superconductivity-induced change in magnetic anisotropy in epitaxial ferromagnet-superconductor hybrids with spin-orbit interaction. Phys. Rev. B 102, 020405(R) (2020).
Johnsen, L. G., Banerjee, N. & Linder, J. Magnetization reorientation due to superconducting transition in heavy metal heterostructures. Phys. Rev. B 99, 134516 (2019).
Banerjee, N., Robinson, J. W. A. & Blamire, M. G. Reversible control of spin-polarized supercurrents in ferromagnetic Josephson junctions. Nat. Commun. 5, 4771 (2014).
Baek, B., Rippard, W. H., Benz, S. P., Russek, S. E. & Dresselhaus, P. D. Hybrid superconducting-magnetic memory device using competing order parameters. Nat. Commun. 5, 3888 (2014).
Gingrich, E. C. et al. Controllable 0-\(\pi \) Josephson junctions containing a ferromagnetic spin valve. Nat. Phys. 12, 564 (2016).
Satchell, N. et al. Pt and CoB trilayer Josephson \(\pi \) junctions with perpendicular magnetic anisotropy. Sci. Rep. 11, 11173 (2021).
Bergeret, F. S., Volkov, A. F. & Efetov, K. B. Long-range proximity effects in superconductor-ferromagnet structures. Phys. Rev. Lett. 86, 4096 (2001).
Keizer, R. S. et al. Spin triplet supercurrent through the half-metallic ferromagnet CrO2. Nature 439, 825 (2006).
Fominov, Ya. . V. et al. Superconducting triplet spin valve. JETP Lett. 91, 308 (2010).
Leksin, P. V. et al. Evidence for triplet superconductivity in a superconductor-ferromagnet spin valve. Phys. Rev. Lett. 109, 057005 (2012).
Bergeret, F. S. & Tokatly, I. V. Singlet-triplet conversion and the long-range proximity effect in superconductor-ferromagnet structures with generic spin dependent fields. Phys. Rev. Lett. 110, 117003 (2013).
Banerjee, N. et al. Controlling the superconducting transition by spin-orbit coupling. Phys. Rev. B 97, 184521 (2018).
Satchell, N. & Birge, N. O. Supercurrent in ferromagnetic Josephson junctions with heavy metal interlayers. Phys. Rev. B 97, 214509 (2018).
Martínez, I. et al. Interfacial spin-orbit coupling: A platform for superconducting spintronics. Phys. Rev. Appl. 13, 014030 (2020).
Dubonos, S. V., Geim, A. K., Novoselov, K. S. & Grigorieva, I. V. Spontaneous magnetization changes and nonlocal effects in mesoscopic ferromagnet-superconductor structures. Phys. Rev. B 65, 220513(R) (2002).
Milošević, M. V. & Peeters, F. M. Interaction between a superconducting vortex and an out-of-plane magnetized ferromagnetic disk: Influence of the magnet geometry. Phys. Rev. B 68, 094510 (2003).
Yang, H. X., Chshiev, M. & Dieny, B. First-principles investigation of the very large perpendicular magnetic anisotropy at Fe–MgO and Co–MgO interfaces. Phys. Rev. B 84, 054401 (2011).
Nistor, L. E. et al. Oscillatory interlayer exchange coupling in MgO tunnel junctions with perpendicular magnetic anisotropy. Phys. Rev. B 81, 220407(R) (2010).
Barnes, S. E., Ieda, J. & Maekawa, S. Rashba spin-orbit anisotropy and the electric field control of magnetism. Sci. Rep. 4, 4105 (2014).
Wang, W. et al. Out-of-plane magnetic anisotropy enhancement in \(\text{ La}_{1-x}\text{ Sr}_x\text{ CoO}_{3-\delta }\)/\(\text{ La}_{2/3}\text{ Sr}_{1/3}\text{ MnO}_3\)/\(\text{ La}_{1-x}\text{ Sr}_x\text{ CoO}_{3-\delta }\) thin films. Phys. Rev. B 101, 024406 (2020).
Wolloch, M. & Suess, D. Strain-induced control of magnetocrystalline anisotropy energy in FeCo thin films. J. Magn. Magn. Mater. 522, 167542 (2021).
Abo, G. S. et al. Definition of magnetic exchange length. IEEE Trans. Magn. 49, 8 (2013).
Fritzsche, J., Kramer, R. B. G. & Moshchalkov, V. V. Visualization of the vortex-mediated pinning of ferromagnetic domains in superconductor-ferromagnet hybrids. Phys. Rev. B 79, 132501 (2009).
Niedzielski, B. & Berakdar, J. Controlled vortex formation at nanostructured superconductor/ferromagnetic junctions. Phys. Status Solidi B 257, 7 (2019).
Eschrig, M. Spin-polarized supercurrents for spintronics: A review of current progress. Rep. Prog. Phys. 78, 104501 (2015).
Palermo, X. et al. Tailored flux pinning in superconductor-ferromagnet multilayers with engineered magnetic domain morphology from stripes to skyrmions. Phys. Rev. Appl. 13, 014043 (2020).
Shafraniuk, S. E., Nevirkovets, I. P. & Mukhanov, O. A. Modeling computer memory based on ferromagnetic/superconductor multilayers. Phys. Rev. Appl. 11, 064018 (2019).
Golod, T., Kapran, O. M. & Krasnov, V. M. Planar superconductor-ferromagnet-superconductor josephson junctions as scanning-probe sensors. Phys. Rev. Appl. 11, 014062 (2019).
Jeon, K.-R. et al. Effect of Meissner screening and trapped magnetic flux on magnetization dynamics in thick \(\text{ Nb }/\text{Ni}_{{80}}\text{ Fe}_{{20}}/\text{Nb }\) trilayers. Phys. Rev. Appl. 11, 014061 (2019).
Tiusan, C. et al. Spin tunneling phenomena in single crystal magnetic tunnel junction systems. J. Phys. Condens. Matter 19, 165201 (2007).
Acknowledgements
We acknowledge Antonio Lara and Miguel Granda for help with simulations and Yuan Lu for help in sample preparations. The work in Madrid was supported by Spanish Ministerio de Ciencia (RTI2018-095303-B-C55) and Consejería de Educación e Investigación de la Comunidad de Madrid (NANOMAGCOST-CM Ref. P2018/NMT-4321) Grants. FGA acknowledges financial support from the Spanish Ministry of Science and Innovation, through the “María de Maeztu” Program for Units of Excellence in \(R \& D\) (CEX2018-000805-M) and “Acción financiada por la Comunidad de Madrid en el marco del convenio plurianual con la Universidad Autónoma de Madrid en Línea 3: Excelencia para el Profesorado Universitario”. D.C. has been supported by Comunidad de Madrid by contract through Consejería de Ciencia, Universidades e Investigación y Fondo Social Europeo (PEJ-2018-AI/IND-10364). NB was supported by EPSRC through the New Investigator Grant EP/S016430/1. The work in Norway was supported by the Research Council of Norway through its Centres of Excellence funding scheme Grant 262633 QuSpin. C.T. acknowledges “EMERSPIN” grant ID PN-IIIP4-ID-PCE-2016-0143, No. UEFISCDI: 22/12.07.2017 and “MODESKY” Grant ID PN-III-P4-ID-PCE-2020-0230 No. UEFISCDI: 4/04.01.2021. The work in Nancy was supported by CPER MatDS and the French PIA project “Lorraine Université d’Excellence”, reference ANR-15-IDEX-04-LUE.
Author information
Authors and Affiliations
Contributions
F.G.A. and C.T. designed the junctions, and C.T. and M.H. fabricated the junctions. The measurements were performed by C.G.-R. with the participation of D.C. The modelling was developed by L.G.J. and J.L. Micromagnetic and TDGL calculations have been carried out by D.C. and C.G-R. F.G.A., C.G.-R., D.C., L.G.J., N.B. and J.L. analyzed the data and wrote the manuscript. All authors contributed to the study, discussed the results, commented and revised the manuscript. F.G.A. supervised the work.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
González-Ruano, C., Caso, D., Johnsen, L.G. et al. Superconductivity assisted change of the perpendicular magnetic anisotropy in V/MgO/Fe junctions. Sci Rep 11, 19041 (2021). https://doi.org/10.1038/s41598-021-98079-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-98079-5
This article is cited by
-
Perpendicular magnetic anisotropy, tunneling magnetoresistance and spin-transfer torque effect in magnetic tunnel junctions with Nb layers
Scientific Reports (2023)
-
Electron beam evaporation of superconductor-ferromagnet heterostructures
Scientific Reports (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.