Tailoring non-collinear magnetism by misfit dislocation lines

Tailoring non-collinear magnetism by misfit dislocation lines

Aurore Finco    Pin-Jui Hsu    André Kubetzka    Kirsten von Bergmann    Roland Wiesendanger Department of Physics, University of Hamburg, D-20355 Hamburg, Germany
August 30, 2019

The large epitaxial stress induced by the misfit between a triple atomic layer Fe film and an Ir(111) substrate is relieved by the formation of a dense dislocation line network. Spin-polarized scanning tunneling microscopy (SP-STM) investigations show that the strain is locally varying within the Fe film and that this variation affects the magnetic state of the system. Two types of dislocation line regions can be distinguished and both exhibit spin spirals with strain-dependent periods (ranging from to ). Using a simple micromagnetic model, we attribute the changes of the period of the spin spirals to variations of the effective exchange coupling in the magnetic film. This assumption is supported by the observed dependence of the saturation magnetic field on the spin spiral period. Moreover, magnetic skyrmions appear in an external magnetic field only in one type of dislocation line area, which we impute to the different pinning properties of the dislocation lines.


I Introduction

Strain-induced control of complex magnetic states such as spin spirals or skyrmions von Bergmann et al. (2014); Nagaosa and Tokura (2013) is a multifaceted approach towards manipulation of spin structures. Mechanical or piezoelectrical Cui et al. (2013) setups as well as direct influence on the sample growth are used to investigate such effects with multiple experimental techniques on various systems. The observed phenomena are ascribed to changes of the effective anisotropy: uniaxial mechanical compression tunes the stability of the skyrmion lattice phase, e.g. in the helimagnet MnSi Nii et al. (2015), where hydrostatic pressure also decreases the helix period and the critical temperature Fåk et al. (2005). On the other side, epitaxial strain created by growing multiferroic BiFeO thin films on different substrates can destroy the cycloidal spin spiral and stabilize an antiferromagnetic state Sando et al. (2013), again via effective anisotropy. The non-collinear magnetic structures mentioned above are stabilized by the Dzyaloshinskii-Moriya interaction (DMI), which can also be affected by strain. This mechanism was found in FeGe, where uniaxial tensile strain dramatically distorts the skyrmion lattice Shibata et al. (2015).

Although previous studies were mostly focused on the effects of spatially uniform strain, we investigate here the influence of strain variations on a local scale in an ultrathin magnetic film using spin-polarized scanning tunneling microscopy (SP-STM) Wiesendanger (2009). Our system consists of three atomic layers of Fe deposited on Ir(111). Because of the misfit between Ir (fcc, lattice constant ) and Fe (bcc, lattice constant ), the triple layer Fe film exhibits a dislocation line network. We observe that the epitaxial strain relief is not uniform within the ultrathin film, resulting in different regions exhibiting spin spirals with various periods and vanishing at different magnetic field intensities. We attribute these differences to spatial variations in the strength of the effective exchange coupling. Depending on the atom arrangement in the magnetic film, these spirals may or may not split up into single magnetic skyrmions Hsu et al. .

Ii Experimental details

The experiments were performed in a ultrahigh vacuum system with a base pressure below . Different chambers are used for substrate cleaning, Fe deposition and STM measurements. The Ir single crystal substrate was prepared by cycles of Ar-ion sputtering at and annealing up to for . The Fe film was then evaporated onto the clean substrate at elevated temperature (about ) at a deposition rate around 0.2 atomic layer per minute. Four different scanning tunneling microscopes were used for this work, three low temperature microscopes with He bath cryostats operating at respectively , and as well as a variable temperature system equipped with a He flow cryostat. Out-of-plane external magnetic fields were applied during the low temperature measurements using superconducting coil magnets. We used antiferromagnetic bulk Cr tips for the spin-resolved measurements. All the STM topography images were measured in constant-current mode, where the stabilization current is kept constant by a closed feedback loop. The differential conductance maps were simultaneously recorded at a fixed sample bias voltage using lock-in technique.

Iii Morphology of the Fe film

The STM topography map presented in Fig. 1(a) shows a typical triple layer Fe film on an Ir(111) substrate. The surface of the film exhibits a dense network of dislocation lines. These lines follow the three equivalent high-symmetry directions of the fcc(111) surface Hsu et al. as already reported for the double layer Fe film Hsu et al. (2016).

Figure 1: (a) Constant-current STM topography map of the triple layer Fe film on Ir(111). Double and single lines are indicated by the arrows. (b),(c) Zoom-in on double line areas. At positive bias, one can see that bright and dark lines alternate and show a double line feature, whereas at negative bias, the lines look all very similar. (d),(e) STM topography of single line regions. The single lines have the same appearance at any sample bias voltage. (f) Constant-current STM topography map of a triple layer Fe island on top of a double layer Fe film. The numbers in green circles indicate the local thickness of the film. The color scale was adjusted separately for the two terraces in order to highlight the matching of the dislocation lines. The areas between the bright lines on the double layer are bcc. (g),(h) Proposed structure models for the triple layer Fe film from the experimental observations presented in (a)-(f). Measurement parameters: and (a) , , ; (b) , ; (c) , ; (d) , ; (e) , ; (f) , , .

Two types of dislocation lines can be distinguished from the STM image of Fig. 1(a): at positive sample bias voltage, the dislocation lines in the region spanned by the white dashed arrow show an alternating bright and dark contrast as well as double line features (see detail in Fig. 1(b)). These areas will thus be named double lines in the following. The spacing for double lines is ranging from to , which corresponds to half the structural period. On the other hand, in the region marked by the gray arrow, the lines are denser (spacing between and ) and they all have the same appearance at any bias voltage (see Fig. 1(d) and Fig. 1(e)). Hence these lines will be designated as single lines. At negative sample bias voltage, the distinction between single and double lines from the sole topography becomes challenging as illustrated in Fig. 1(c) and Fig. 1(e).

Figures 1(g) and 1(h) show proposed atomic structure models for the triple layer Fe film. In both cases, the first layer Fe is pseudomorphic with respect to the Ir(111) substrate lattice Heinze et al. (2011). The double lines are located exactly on top of the dislocation lines of the double layer Fe film as shown in Fig. 1(f). The double layer lines seem to get closer when they approach the island: the strain relief is larger in the triple than in the double layer. In the atomic model, the second layer is uniaxially compressed with respect to a pseudomorphic layer and the arrangement of the Fe atoms is getting close to bcc(110). The lines correspond to the positions where the atoms are located on the fcc or hcp sites and the triple layer follows a bcc stacking on top of the double layer. In the case of the single lines, also the second layer grows pseudomorphically but in bcc stacking (see Fig. 1(h)) and only the third layer on top is uniaxially compressed. All the lines are identical and induced by lateral shifts of the atoms. As will be shown later by looking at the magnetic structure, the double line regions have mirror planes along the dislocation lines, whereas two mirror symmetric domains exist for the single line regions, in agreement with the structure model.

Iv Zero field non-collinear magnetic structure

The difference between single and double line areas becomes more striking when we look at the magnetic structure. Figures 2(a) and 2(b) show respectively a spin-resolved constant-current map and a simultaneously spin-resolved differential conductance map measured with an out-of-plane spin-sensitive Cr bulk tip at low temperature. These SP-STM measurements reveal that spin spirals propagate along both types of dislocation lines. The nature of the spin spirals is cycloidal: this has been shown already for the double layer Fe spirals Hsu et al. (2016) as well as for the double line regions in the triple layer Fe Hsu et al. and is assumed also to be the case for the single line areas since the Fe-Ir interface-induced DMI is large Heinze et al. (2011). However, the appearance of the spin spirals varies in the different regions of the Fe film. The wavelength of the spirals is smaller in the double line regions ( to ) than in the single line regions (between and ).

Figure 2: (a),(b) Spin-resolved constant-current map and spin-resolved differential tunneling conductance map of a triple layer Fe film on Ir(111) measured simultaneously using a Cr bulk tip with out-of-plane spin sensitivity. The magnetic sensitivity of the tip is simply deduced from the magnetic contrast on the spirals, which is identical for all the propagation directions. The magnetic contrast is visible on both the constant-current and differential conductance maps. The ellipses on image (b) mark a zone where the period of the spiral is changing: this can be correlated with the spacing between the dislocation lines indicated in (a). (c) Detail of the spin-resolved differential tunneling conductance map of a double line region. The spin spiral propagates along the lines and its wavefront exhibits a zigzag shape. (d),(f) Line profiles taken in the double and single line regions, respectively, at the positions indicated in (b). The displayed profiles are the mean values within the rectangles, each data line was laterally shifted to straighten the wavefront and offsets were applied to the resulting profiles for clarity. The red dashed lines are the results of fits with sine functions, showing that the spirals are homogeneous regardless of the considered region. (e) Spin-resolved differential conductance map of an area with single lines. Spin spirals with various wavelengths are visible and their wavefronts are tilted by about with respect to the lines. Measurement parameters: no external magnetic field, (a),(b) , , ; (c) , , ; (e) , , .

Furthermore, the wavefront of the double line spirals has a zigzag shape Hsu et al. (2016, ), whereas that of the spirals in the single line regions is straight but tilted with respect to the lines as can be seen from the details of spin-resolved differential conductance maps in Fig. 2(c) and Fig. 2(e). The proposed structure models can explain these shapes for the wavefronts: the wavevector prefers to follow the bcc[001]-like rows of atoms as observed for the double layer Fe on Ir(111) Hsu et al. (2016), Cu(111) Phark et al. (2014) and W(110) Meckler et al. (2009). For the double lines, the direction of the rows alternate and this creates the magnetic zigzag structure, whereas the direction does not change for the single lines, as shown by the bcc(110)-like unit cells marked in red in Fig. 1(g) and Fig. 1(h). However, for the single line regions, two mirror symmetric structural domains are possible and found in the SP-STM data.

For a strict coupling of the wavevector to the bcc [001] direction, the expected angle between the wavefront and the dislocation lines (as defined in Fig. 2(c) and Fig. 2(e)) can be computed from the structure models:


where is the line spacing and the in-plane interatomic distance for Ir(111). This expression is the same for both structures. This angle is below and decreases with the line spacing. Yet the angle obtained from the data is slightly larger (up to ) and seems to change randomly. This deviation which is going towards a straighter wavefront for the double line regions was previously observed for the double layer spirals Hsu et al. (2016). The guiding of the wavevector might compete here with the presence of energetically disadvantageous kinks in the wavefront. For the single line areas, the angle is also larger, meaning that the wavefront is more perpendicular to the dislocation lines than expected. This effect could be attributed to boundary effects as well as to domain wall-like structures preferring to be as short as possible. Line profiles of spirals (Fig. 2(d) and Fig. 2(f)) taken in various regions can be fitted with sine functions. This sinusoidal shape indicates that the spirals are homogeneous, i.e. no direction of the magnetic moments is more favorable than the others. Yet epitaxial strain in ultrathin films creates an effective anisotropy via magnetoelastic coupling Bruno and Renard (1989) and since the strain relief is not uniform in the triple layer Fe film, the total effective magnetic anisotropy is expected to vary between the different regions. This should result in more or less pronounced distortions of the spiral profiles depending on the dislocation line spacing. We therefore conclude from the homogeneity of all the observed spiral profiles that the effective magnetic anisotropy and its variations in the triple layer Fe film are small enough to be neglected in the following.

V Tuning of the spiral period with strain relief

In Fig. 2(b), a single line area in which the wavelength of the spiral is changing is marked with ellipses. This variation occurs because the dislocation lines become more distant in the green zone than in the surrounding red and blue ones. The distance between the lines in each ellipse is given in the constant-current map of Fig. 2(a). The closer the lines, i.e. the larger the compression in the third layer Fe, the larger the spiral period.

Figure 3: Dependence of the spiral period on the spacing of the dislocation lines for the triple layer Fe on Ir(111), both double and single line regions. The SP-STM measurements were performed at low temperature ( or ) on several different samples and the spiral periods were extracted using 2D Fourier transformation and fits of the real space data to sine functions, keeping only points with an error bar below 15%. For the double line regions, the zigzag shape of the wavefront was not considered and the period was measured along the lines. The actual wavelength along the direction of the wavevector might be 10% smaller. On the single line regions, the angle of the wavefront was taken into account in the determination of the period. It is evident that if the distance between the dislocation lines decreases, the spiral period increases, as observed already in Fig. 2. The insets show the data (both spin-resolved constant-current map and differential conductance map as well as their Fourier transforms) used to determine the point marked with the red square. The axis on the right indicates that the effective exchange parameter should vary between and to create spirals with wavelengths between and from the simple model described in the text.

The graph in Fig. 3 presents a systematic investigation of the effect of the strain relief on the wavelength of the spirals, for both double and single line regions. The spiral period and the spacing between dislocation lines were measured using Fourier transforms and fits of the real space data to sine functions in several different regions on several samples. The trend observed in Fig. 2 is confirmed, the spiral period decreases when the line spacing increases. The line spacing is linked to the amplitude of the strain relief by the structure models: the mean value of the distance between the Fe atoms in the fcc [11̄0] direction in the top layer decreases with the line spacing. To give an order of magnitude, this distance is 5% smaller for a than for a line spacing.

When the distance between the dislocation lines changes, all the interatomic distances in the Fe film are modified, thus it is expected that the strengths of the magnetic interactions will be affected. The relevant interactions for this system are the exchange couplings and the interface-induced DMI Fert and Levy (1980); von Bergmann et al. (2014). Based on the data shown in Fig. 2(d) and Fig. 2(f), we assume that the effective anisotropy is negligible. The DMI originates mainly from the Fe-Ir interface and since the first layer Fe layer is pseudomorphic regardless of the considered region, it should not be significantly influenced (contrary to bulk systems like FeGe Shibata et al. (2015)). We therefore infer that the strain relief is acting mostly on the exchange couplings to modify the period of the spirals.

Ab initio calculations for a free-standing Fe bcc(110) layer Shimada et al. (2012) also assigned the period variation of the spirals as well as their stability under in-plane strain to modifications of the exchange couplings. However, the spin spiral state is only stable in the free layer when a compressive strain is applied and the compression reduces the period of the spiral. The calculated evolution of the spirals with the strain is thus opposite to what we measured but this discrepancy could result from the presence of the substrate and the complicated atomic structure of the Fe film.

In order to estimate the magnitude of the strain effect, we consider a simplified one-dimensional micromagnetic model derived from the one proposed by Bogdanov and Hubert Bogdanov and Hubert (1994) in which only the effective isotropic exchange coupling and the DMI are kept. We completely ignore here that the film is not spatially uniform, resulting in spatially inhomogeneous magnetic coupling constants. Only average values are taken into account. The energy density is thus:


where is the effective exchange stiffness constant, the effective DMI constant and the reduced dimensionless magnetization. In this particular case, there is always a stable cycloidal spin spiral state in the system and its period is simply:


For the monolayer Fe on Ir(111), the value of the DMI constant from ab initio calculations is   Heinze et al. (2011) which gives the micromagnetic parameter for a three-layer-thick Fe film. As shown in Fig. 3, the exchange stiffness parameter should then vary between and to obtain magnetic periods ranging from to . These values are similar to found for the PdFe bilayer on Ir(111)  Romming et al. (2015) with a spiral period of about . Since the DMI originates from the interface, its effect should decrease when the thickness of the film increases, i.e. the typical length scale of the magnetic structure is expected to become larger. Indeed, the period of the nanoskyrmion lattice on the monolayer Fe on Ir(111) is  Heinze et al. (2011), the wavelength of the double layer spiral is about  Hsu et al. (2016) and for patches of a size below (larger ones were not found on the samples), the quadruple layer appears ferromagnetic (as, for example, in Fig. 5).

Vi Skyrmions and domain walls in magnetic field

In external out-of-plane magnetic fields between and , the spin spirals in the double line regions split up into distorted magnetic skyrmions Hsu et al. , as illustrated in Fig. 4. In contrast, no skyrmion is created on the single line regions. There, the spirals become inhomogeneous and behave like an assembly of independent domain walls which disappear one by one when the field is increased. Adding a Zeeman term to the model proposed before (equation (2)) does not allow to understand this difference. For a isotropic film without effective anisotropy, there is a stable skyrmionic phase for any couple under the appropriate external magnetic field. The crucial point for the triple layer Fe film on Ir(111) are the dislocation lines which break the rotational symmetry. Thus one cannot expect the existence of the typical circular shaped skyrmions Hagemeister et al. (2016) as observed, for example, in the PdFe bilayer on Ir(111) Romming et al. (2013) and a priori, even the presence of skyrmions is not obvious. Nevertheless, distorted skyrmions appear in the double line regions. Their bean-like shape is induced by the local arrangement of the Fe atoms similarly to the zigzag shape of the spin spiral wavefront. They are always located on top of three dislocation lines, two identical ones at the ends and a third different one in the center. This particular preferred position results from the local variation of the magnetic interaction within the layer. The skyrmions are hence pinned to the lines and naturally aligned on “tracks” (see Fig. 4 and Fig. 5(g)). The atom arrangement is different for the single lines and it appears that the skyrmion pinning effect is absent and only domain walls are observed in magnetic field.

Figure 4: Spin-resolved differential tunneling conductance maps of a region with single lines surrounded by double lines (left) and of a region with double lines surrounded by single lines (right) in increasing out-of-plane magnetic field, measured with an out-of-plane spin-sensitive Cr bulk tip. When the external field increases, the spin spirals in the double line regions transform into distorted skyrmions whereas those in the areas with single lines become inhomogeneous and can be described as an assembly of independent domain walls. When the width of these walls is similar to the width of the skyrmions in the adjacent areas, they couple to them. Measurement parameters: Left: , , ; Right: , , .

Vii Metastability and transition fields

A more quantitative investigation of the influence of an external magnetic field is shown in Fig. 5. For this measurement series, the magnetic field was increased up to in steps. At , the sample has fully reached the ferromagnetic state. Then, the magnetic field was decreased again in steps to zero. Every region in the film behaves differently because of its structure and its interactions with adjacent areas (see Fig. 4). Comparisons between pictures taken at the same field value during the up-sweep and the down-sweep reveal, however, that almost all the areas show hysteresis. The skyrmions and domain walls vanish at a much higher field than the one needed to produce them again. This indicates that the magnetic states are metastable.

Figure 5: (a)-(h) Spin-polarized differential tunneling conductance maps of an ultrathin Fe film on Ir(111). The numbers in green circles in (a) indicate the local Fe coverage. An external out-of-plane magnetic field was applied, increased in steps of up to (at this field the sample is completely polarized, i.e. ferromagnetic) and then decreased again to zero. Comparison between scans taken at the same field value shows a hysteretic behavior for the triple layer Fe, both in the regions with the double and single lines. (i) Plot of the magnetic state during the field sweep of the four areas marked in blue in (a). Areas I and II are double line regions, III and IV single line regions. Measurement parameters: , out-of-plane spin-sensitive Cr bulk tip.

The transition fields to the ferromagnetic (FM) state are obtained from field dependent measurements using the procedure described in Fig. 5(i): the evolution of the magnetic states for different regions is indicated for an up- and a down-sweep of the magnetic field. The transition field value corresponds to the middle of the corresponding step and the error bar to the step height.

Results for several samples are gathered in Fig. 6, where they are correlated with the spiral period in the absence of external magnetic field. Only sweeps with increasing field are considered in this figure, the field values are thus upper estimates due to the hysteresis effect. The transition field is decreasing when the spiral period increases. The trend that higher external fields are required to destroy non-collinear spin structures with smaller spatial period is consistent which the observation that the spin spiral in the double layer Fe does not change in magnetic fields up to  Hsu et al. (2016).

Figure 6: Effect of the spiral period on the magnetic fields needed to reach the ferromagnetic state during an increasing field sweep. The correlation between the spiral period and spacing of the dislocation lines in the Fe layer shown in Fig. 3 leads also to a dependence of the transition field on the strain relief. A trend that higher external fields are required to destroy non-collinear spin structures with smaller spatial period appears. The calculated transition fields were obtained from the phase diagram detailed by Bogdanov and Hubert Bogdanov and Hubert (1994) with the parameters from Fig. 3: and . The magnetic moment used is per Fe atom Heinze et al. (2011). Because of the metastability of the magnetic states revealed by the hysteresis, the experimental field values are upper estimates of the transition fields. The SP-STM measurements were performed at on three different samples and the external magnetic field was increased in steps of . The error bars correspond to the height of the steps as indicated in Fig. 5(i).

Although the simplified model (2) does not help to understand the absence of a skyrmionic phase in the single line regions, it reproduces the decrease of the transition field for larger magnetic structures once the Zeeman term is included. In the zero effective anisotropy case, the transition field can be obtained from the phase diagram provided by Bogdanov and Hubert Bogdanov and Hubert (1994). There is a threshold value for the reduced field parameter such as:


The saturation magnetization is estimated from the magnetic moment of 2.7 per atom in the monolayer Fe Heinze et al. (2011). We did not consider here a potential variation of with the strain relief. The threshold value is different for the transition from spirals to the FM state and from skyrmions to the FM state:


The transition fields are plotted as plain lines in Fig. 6, assuming again that . They are defined as the fields at which the energy of the FM state is equal to the one of the spiral or skyrmion state, respectively. As expected, the experimental values are larger than the computed ones because of the metastability of the spiral and skyrmion states in increasing magnetic field. Remarkably, the smallest field values are almost on the theoretical curve and none of them below, which indicates a rather good agreement between the model and the actual system. This behavior of the transition fields supports our assumption that the effective exchange coupling is affected by the strain relief and is responsible for the observed variation of the magnetic properties of the Fe film.

Viii Conclusion

Exploiting the strain relief in epitaxial ultrathin films is an effective way to control precisely and locally their magnetic state. Both the typical size of the spin structures and the transition fields could be tuned. Moreover, the actual uniaxial structure and pinning properties of ultrathin films exhibiting dislocation lines may allow to stabilize a skyrmion state as well as to confine skyrmions on well-defined tracks on the order of a few nanometers. Combined with the possibility to write and delete the skyrmions by a local electric field Hsu et al. , this could be of great interest in view of future racetrack-based spintronic devices Parkin et al. (2008); Fert et al. (2013); Wiesendanger (2016).

Financial support by the European Union via the Horizon 2020 research and innovation programme under grant agreement No. 665095, by the Deutsche Forschungsgemeinschaft via SFB668-A8, and by the Hamburgische Stiftung für Wissenschaften, Entwicklung und Kultur Helmut und Hannelore Greve is gratefully acknowledged.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description