Manipulating surface magnetic order in iron telluride

See allHide authors and affiliations

Science Advances  01 Mar 2019:
Vol. 5, no. 3, eaav3478
DOI: 10.1126/sciadv.aav3478


Control of emergent magnetic orders in correlated electron materials promises new opportunities for applications in spintronics. For their technological exploitation, it is important to understand the role of surfaces and interfaces to other materials and their impact on the emergent magnetic orders. Here, we demonstrate for iron telluride, the nonsuperconducting parent compound of the iron chalcogenide superconductors, determination and manipulation of the surface magnetic structure by low-temperature spin-polarized scanning tunneling microscopy. Iron telluride exhibits a complex structural and magnetic phase diagram as a function of interstitial iron concentration. Several theories have been put forward to explain the different magnetic orders observed in the phase diagram, which ascribe a dominant role either to interactions mediated by itinerant electrons or to local moment interactions. Through the controlled removal of surface excess iron, we can separate the influence of the excess iron from that of the change in the lattice structure.


Multiple competing interactions in strongly correlated electron materials lead to a plethora of emergent phases, which are highly sensitive to external stimuli and offer tremendous potential for applications. Exploiting these requires interfacing them to the outside world, yet relatively little is known about the influence of reduced symmetries and the interface itself. Iron telluride (Fe1+xTe) is such a strongly correlated electron material with a complex magnetic phase diagram as a function of excess iron concentration x. At low excess iron concentration (x < 0.11), a bicolinear antiferromagnetic (AFM) order with an ordering wave vector qAFM = (0.5, 0, 0.5) (in units of the reciprocal lattice) is observed in a crystal structure with monoclinic distortion (13). With increasing x, the crystal structure becomes orthorhombic (x > 0.11), with a reduction in the difference in the lattice constants in the a and b directions (24). This change is accompanied by the development of patches with incommensurate magnetic order, which coexists globally with the bicolinear order for 0.11 < x < 0.14 (5). For x > 0.14, the order becomes fully incommensurate, and a helimagnetic spin spiral develops (1, 2).

The bicolinear magnetic structure at low excess iron concentrations (x < 0.11) is well reproduced by density functional theory (DFT) calculations (6, 7), whereas accounting for the influence of interstitial iron has been more challenging. The incommensurate spiral structure can be reproduced by assuming that the interstitial Fe atoms lead to electron doping (8) or by considering additional nearest-neighbor coupling due to the randomly distributed interstitial Fe atoms (9). Even for low interstitial iron concentrations, multi-q plaquette order has been predicted as a result of magnetic frustration and quantum fluctuations (1012).

It is only very recently that real space imaging of the surface magnetic order in iron telluride has been demonstrated by spin-polarized scanning tunneling microscopy (SP-STM) (1315). In this work, we use atomic-scale imaging by low-temperature SP-STM (1618) to determine and manipulate the surface magnetic order in iron telluride at different excess iron concentrations. Our results enable us to assess the impact of the structural distortion and excess iron concentration on the magnetic structure. By manipulating the excess iron concentration of the surface layer, we discover a double-q magnetic order, which is stabilized as the bulk crystal structure becomes orthorhombic.

Low-temperature SP-STM measurements are performed at temperature T = 2K to study samples with excess Fe concentration ranging between 0.04 and 0.2. Ferromagnetic tips for SP-STM are prepared by picking up excess iron atoms from the sample surface (13, 14). This ability to manipulate the excess iron concentration in the surface layer offers the opportunity to control the surface magnetic structure and discriminate between the various scenarios for the influence of interstitial iron. The lattice constants of the surface layer remain locked to the crystal structure of the bulk material independent of the excess iron concentration in the surface layer. By removing the excess Fe atoms from the surface, we are therefore decoupling the effects of the presence of the excess Fe atoms and the orthorhombic crystal distortion on the surface magnetic order. Effectively, this allows for the determination of how the magnetic order of stoichiometric Fe1+xTe would change if the crystal were strained from a monoclinic to an orthorhombic structure. This method of removing Fe atoms from the surface layer of orthorhombic Fe1+xTe thereby provides an alternative way to engineer an effective strain in the surface layer as opposed to using externally applied strain (19).

We will first demonstrate that SP-STM characterization of Fe1+xTe as a function of excess iron concentration reproduces the magnetic phase diagram found by bulk characterization techniques and then show how the magnetic order changes by manipulating the excess iron concentration of the surface layer.


Surface magnetic order of Fe1+xTe

Figure 1A shows a topographic STM image of the surface of Fe1.06Te obtained immediately after first approaching the sample and was recorded with a nonmagnetic PtIr tip prepared by field emission on a gold target. In the image, the square lattice corresponds to the Te atoms in the top surface layer. From large-scale STM imaging of the surface, the bright protrusions are found to have the same concentration as the excess Fe in any particular Fe1+xTe sample studied, and hence, they are identified as the excess Fe atoms. Imaging the surface with a magnetic tip, prepared by collecting Fe atoms from the sample surface (13, 14), leads to the appearance of an additional stripe-like modulation (Fig. 1B) running along the crystal b axis. This modulation has a periodicity along the crystal a axis twice that of the Te lattice and results in the development of an additional peak in the Fourier transform at a wave vector of qAFM = (0.5, 0) [relative to the reciprocal lattice vectors with qTe = (1, 0), (0, 1) that arise from the atomic lattice; compare the Fourier image insets of Fig. 1, A and B]. This additional modulation arises from the bicolinear AFM order in Fe1+xTe with low levels of excess Fe concentration (x < 0.11) (1315). The AFM magnetic order can be seen in the STM images due to tunneling magnetoresistance: For a fixed tip-sample distance, the current depends on the relative magnetization of the tip and the sample. The images shown here are constant current images, resulting in the magnetic contrast being observed in the topographic SP-STM image. The strength of this contrast depends on the relative magnetizations of tip and sample and on the spin polarization of the charge carriers in the tip and the sample at the Fermi energy. Imaging the magnetic structure of Fe1.06Te with the same ferromagnetic tip but with its magnetization aligned along two opposite in-plane directions results in phase reversal in the appearance of the stripe-like modulation in SP-STM images (Fig. 1, C and D) (1315). Subtraction of these SP-STM images produces a real-space image of the magnetic structure (Fig. 1E). The height difference Δz is proportional to the spin polarization of the tunneling current and provides information on the local magnetic order projected onto the magnetization direction of the tip. The inset in Fig. 1E shows for comparison a spin-polarized image simulated on the basis of a DFT calculation, showing excellent agreement with the experimental data (see section S1).

Fig. 1 Spin-polarized STM of Fe1.06Te.

(A) Topographic STM image taken with a nonmagnetic tip (14.5 by 14.5 nm2). Protrusions are excess Fe atoms. Inset: Fourier transform (FT) image of (A). Peaks that arise from the Te lattice are highlighted with solid circles. (B) Topographic SP-STM image taken at the same position as (A) with a magnetic tip. Stripes arise from the AFM order. Inset: FT image of (B) showing additional peaks due to the AFM order. (C and D) Topographic SP-STM images taken at the same position with the tip polarized along two opposite in-plane directions (1.9 by 2.6 nm2). Tunneling parameters for (A) to (D): V = 100 mV, I = 50 pA. (E) Difference image of (C) and (D). The height difference is proportional to the spin polarization of the tunneling current. Inset at the bottom left: Structural model of the Fe1+xTe surface, showing the spin order in the Fe lattice (red). Inset at top right: DFT calculation of the magnetic contrast due to the spin polarization at the Fermi energy (see section S1).

Using the above approach and aligning the magnetization of the tip along three orthogonal directions, it becomes possible to reconstruct the surface magnetic structure of a sample (20)—provided that the sample magnetic structure remains unaltered while changing the magnetization of the tip. In the case of Fe1+xTe, because of the high exchange interaction between the Fe atoms (15) and the considerable magnetic anisotropy energy of each Fe atom (13), the magnetic order remains unaffected by the applied fields of up to 5 T used here (14). For Fe1.06Te, we find that the magnetic moments have an out-of-plane component, pointing into a direction 28 ± 3° away from the surface plane (see section S2 and fig. S1). The magnetic ordering wave vector we find at low excess iron concentrations is in excellent agreement with previous SP-STM (1315) and neutron scattering (1, 2) studies, as well as with calculations (6, 7). The out-of-plane angle of the magnetic order differs from that observed in neutron scattering, where the magnetization of the iron atoms is parallel to the b direction but is fully consistent with previous SP-STM studies (15).

At high excess Fe concentrations x > 0.11, the structure of Fe1+xTe transforms from monoclinic to orthorhombic (2, 3, 21). Figure 2A shows a topographic image of the surface of Fe1.16Te, with three images of the magnetic order obtained with the tip magnetized in three different angles in the b-c plane (Fig. 2, B to D). As the magnetization of the tip rotates, the magnetic contrast in the difference image translates along the a direction, similar to what has been observed in previous SP-STM studies of other systems with spin spiral magnetic orders (22). This demonstrates the presence of a unidirectional spin spiral propagating along the a axis with spins rotating in the b-c plane. Analysis of the line cuts taken from the difference images (Fig. 2E) reveals the spin spiral with a wave vector of q = (0.43, 0), slightly incommensurate with the crystal lattice (see section S3 and fig. S2). The spin spiral found here is in full agreement with that detected in bulk samples at high excess iron concentrations x > 0.12 by neutron scattering.

Fig. 2 Spin spiral in Fe1.16Te.

(A) Topographic SP-STM image (8.2 by 4.6 nm2, V = 50 mV, I = 200 pA). (B to D) Magnetic images taken at out-of-plane angles θ = −90°, −30°, and 30°. Close inspection reveals that the positions of the maxima of the magnetic order shift as a function of out-of-plane angle θ. (E) Line cuts through magnetic images as shown in (B) to (D) along a for different out-of-plane field angles. The line cuts show the shift of the maxima of the stripes. (F) Plot of the phase of the stripes shown in (E) as a function of field angle θ. The phase has been extracted using the maximum marked by an arrow in (E). Measurements were taken at an in-plane angle φ = 120° from the crystal a axis.

The key result from this section is that SP-STM at the surface of Fe1+xTe yields consistent result with neutron scattering. Differences appear merely in details, such as the out-of-plane component of the ordered moment.

Manipulation of the excess iron concentration

In addition to imaging the surface magnetic structure at the atomic scale, STM also allows us to manipulate the surface composition, which for the case of Fe1+xTe can be achieved by fast, high tunneling current scanning (~2 nA) at a slow feedback loop response time. The resulting change in composition is illustrated in Fig. 3: The interstitial iron atoms originally present on the surface layer of Fe1+xTe are removed, leaving the remaining Fe interstitials at the lower part of the surface layer intact (for details, see section S4). As a consequence, this provides an opportunity to study the magnetic order in the surface layer with only half of the excess Fe concentration present, but with the same lattice structure as that of the bulk. Because the magnetic coupling between the layers is rather weak, as determined from inelastic neutron measurements of the spin wave excitations (23), a study of the surface magnetic structure provides us with information about how the bulk of the material would respond to similar conditions.

Fig. 3 Manipulation of surface excess iron concentration.

Model of Fe1+xTe (A) before and (B) after removal of excess iron. Red (yellow) spheres are Fe (Te) atoms. Dashed open circles mark the interstitial/excess Fe atoms at the surface layer that are removed during surface manipulation, leading to a 50% reduction in the concentration of Fe interstitials in the surface layer (to x/2) compared to that of the bulk (x). (C) SP-STM image of Fe1.12Te (65.3 by 28.2 nm2, V = 150 mV, I = 50 pA) showing an area where surface excess iron has been removed by the tip (blue) next to one where the excess iron has been left untouched (green). (D) Fourier transform of (C) showing magnetic peaks due to the bicolinear order (blue arrow) and at an incommensurate position (green arrow, see fig. S5 for Fourier transforms of regions with high and low excess iron concentrations at the surface). (E) Line cut from the Fourier transform in (D) taken along the aTe direction. Peaks corresponding to the bicolinear order at q = (0.5, 0) and the incommensurate order q = (0.39, 0) are highlighted by a blue arrow and a green arrow, respectively. (F and G) Maps of the intensity of the magnetic order at the wave vector of the bicolinear order [q = (0.5, 0)] and of the incommensurate order [q = (0.39, 0)]. Both have been obtained through Fourier filtering at the respective wave vector and then low pass filtering of the modulus. The maps show that the bicolinear order predominantly exists in regions that have been cleaned of Fe, while the incommensurate order is dominant in regions where Fe is still present. a.u., arbitrary units.

Figure 3C shows an SP-STM image of the surface of Fe1.12Te, on one half of which the majority of the surface excess Fe atoms have been removed, whereas they are still present in the other half. The Fourier transform of the area shown in Fig. 3C plotted in Fig. 3D reveals two sets of magnetic ordering peaks: one at the bicolinear ordering vector q = (0.5, 0) and one at an incommensurate vector q = (0.39, 0). Fourier filtering the topographic image (Fig. 3C) at q = (0.5, 0) (Fig. 3F) and at q = (0.39, 0) (Fig. 3G), respectively, reveals that the bicolinear order is concentrated to the areas where the excess Fe atoms at the surface have been mostly removed while the incommensurate order is confined to areas where the initial surface excess Fe concentration has been left untouched. This incommensurability has also been observed in neutron scattering conducted at similar excess Fe concentrations. The STM image in Fig. 3C demonstrates that while the interstitial iron concentration of the surface layer can be manipulated using STM, the lattice structure remains commensurate with the bulk as no additional superstructure is seen, which would arise from a structural incommensurability of the surface layer with the bulk.

Impact of excess iron concentration on surface magnetic order

A detailed investigation of the magnetic structure of the area cleaned of excess iron in Fe1.12Te (compare Fig. 3) reveals a complex picture: Fig. 4 (A and B) shows images of the magnetic order in the surface layer projected onto two different magnetization directions of the tip, in-plane and out-of-plane. While the image obtained with an in-plane magnetization of the tip (Fig. 4A) exhibits only the unidirectional bicolinear order, measurements with an out-of-plane magnetized tip (Fig. 4B) reveal domains of checkerboard-like double-q order. In the Fourier transformation (Fig. 4C), the component of the magnetic order along the a axis is characterized by a sharp peak at q = (0.5, 0), whereas the new component along b manifests itself as a broad peak at q ~ (0, 0.5), reflecting its localized nature. Unlike the single-q magnetic order found at low excess Fe concentration, which has the magnetization of the iron atoms oriented in opposite directions in the b-c plane, in the double-q order, the magnetization is also modulated in the a direction, leading to the formation of a spin spiral–like order. A model that reproduces this behavior is described in section S6 and equation S1. In the orthorhombic phase at higher excess iron concentration, we find an even stronger change of the magnetic structure after removing the surface interstitial iron atoms: Fig. 5A, obtained from a surface layer of Fe1.1Te on an Fe1.2Te sample, shows a topographic, image together with three images of the magnetic order (Fig. 5, B to D) for the tip magnetized along two in-plane (Fig. 5, B and C) and an out-of-plane direction (Fig. 5D). All three magnetic images reveal a double-q magnetic order, of which both components, along the a and b directions, are characterized by commensurate wave vectors of q = (0.5, 0) and (0, 0.5). As revealed by their Fourier images (insets of Fig. 5, B to D), the strength of both components varies with the magnetization of the tip. To resolve the magnetic structure, we have rotated the magnetization of the tip in the surface plane in steps of 18° through 180° and recorded a magnetic image for each angle (for details, see section S5). The results are summarized in Fig. 5E. The intensities of both components [I (qa), I (qb)] vary as |cos φ| and reach their maxima at ~30° from the b axis. Both exhibit an almost identical angular dependence under the in-plane rotation. However, when the magnetization of the tip is rotated out of the surface plane, the strength of the component of the magnetic order along a varies strongly, while that along b remains unchanged (Fig. 5F). This can be accounted for by a model magnetic structure consisting of two spin spirals along the [110] and Embedded Image directions (see section S7, eq. S2, and fig. S7). The spin spirals alternate between clockwise and counterclockwise winding on alternate rows of Fe atoms.

Fig. 4 Magnetic order in a surface layer of Fe1.06Te on Fe1.12Te.

(A) Image of the magnetic order (20.5 by 20.5 nm2) projected onto an in-plane direction of the magnetization as indicated by the arrows. (B) As in (A), taken with the same tip with an out-of-plane direction of the magnetization (V = 100 mV, I = 50 pA). (C) FT image of (B), with intensity at the center suppressed for clarity. Peaks due to the AFM order along a and b are marked with pink and cyan circles, respectively. (D) Normalized line cuts taken from the origin along a (red) and b (blue) in (C). Dashed lines indicate the positions qTe and qAFM along both a and b directions.

Fig. 5 Plaquette order in on Fe1.1Te on Fe1.2Te.

(A) Nonmagnetic image (for details, see section S5 and fig. S6). Inset: FT image showing peaks due to the Te lattice. (B) Image of the magnetic structure for φ = −27°. AFM order can be seen along both lattice directions. Inset: FT image of (B), with peaks due to AFM order in the a (b) direction marked with pink (blue) circles. (C and D) As in (B), with (C) φ = 27° and (D) out-of-plane magnetic field (θ = 90°). Insets show the corresponding FT images. Images (A) to (D) were recorded in the same area (24.1 by 24.1 nm2), V = −40 mV, I = 100 pA. (E) Integrated intensities of the magnetic peaks in the FT image as a function of in-plane angle φ. The horizontal dashed line indicates the integrated intensity of a point away from the magnetic peaks. Blue (red) markers are the intensities of the qb (qa) peak. Solid lines are numerical fits of I = I0|cos (φ − φ0)| + c to the data. Vertical dashed lines indicate the in-plane field directions parallel to the a and b axes, as well as that of maximum intensity (φ = 117°). (F) As in (E), plotted as a function of out-of-plane angle θ at in-plane angle of maximum intensity of the magnetic order, φ = 117°. Data shown here were recorded in the same location and with the same tip.


By combining spin-polarized imaging with the ability to manipulate the surface Fe atoms, our results enable us to extract a comprehensive picture of the magnetic phase diagram of iron telluride and assess the impact of interstitial iron on the magnetism in this material.

In Fig. 6A, we show models of the magnetic structure deduced from our study of the surface magnetic structure after removal of the interstitial iron, effectively for a surface layer of Fe1+x/2Te on Fe1+xTe. In the monoclinic phase, we observe the same magnetic order after removal of interstitial iron as in the bulk. With reduced asymmetry of the lattice constants in the a and b directions, we see an increasing deviation of the spins from the b-c plane. At the transition of the bulk from the monoclinic to the orthorhombic structure at x = 0.12, the surface layers exhibit patches of double-q order, with apparent competition between the bicolinear order in the a direction and a developing spin density wave (SDW) order in the b direction. At even higher excess iron concentration in the bulk, the cleaned surface layer adopts a double-q magnetic order. As opposed to the bulk, the surface magnetic order remains commensurate after removing the excess iron atoms.

Fig. 6 Magnetic structures and phase diagram of a surface layer of Fe1+x/2Te.

(A) Model structures of surface layers of Fe1.03Te/Fe1.06Te, Fe1.06Te/Fe1.12Te, and Fe1.1Te/Fe1.2Te that are consistent with the SP-STM results. The magnetic unit cell in each case is highlighted (for details of the model, see section S7). Arrows indicate Fe spins and are colored blue if they have a positive component along b and red if they have a negative component. Green and gray spheres represent the upper and lower Te atoms. (B) Phase diagram of the magnetic order in the surface layer after removal of excess iron (with concentration Fe1+x/2Te) as a function of bulk excess Fe concentration x. Red dots represent the ratio of the intensities I(qb)/I(qa) of the magnetic order along the lattice directions a and b taken from the Fourier transforms of the STM data. The blue diamonds show the wave vector of the magnetic order in terms of the lattice spacing from neutron scattering from (2); the black dotted lines depict the “mixed spin density wave (SDW)” phase defined there. The gray background highlights how the lattice parameters change with excess Fe concentration and is defined by Embedded Image, using values for a and b given in (2) (see also table S1). (C to E) SP-STM images (7 by 7 nm2) of out-of-plane polarization of a surface layer of (C) Fe1.03Te in the monoclinic phase and (D) Fe1.06Te and (E) Fe1.1Te in the orthorhombic phase.

This result is summarized in the phase diagram in Fig. 6B, where we compare the incommensurability seen in neutron scattering with the appearance of magnetic order at a second q-vector along b in our SP-STM measurements and the anisotropy of the crystal structure. To this end, we plot the ratio I(qb)/I(qa) of the intensity I(q) of the magnetic order in the b direction at qb = (0, 0.5) to that in the a direction, at qa = (0.5, 0), as a function of excess iron concentration x. The doping dependence of this ratio shows that the magnetic order at qb appears once the bulk has undergone the transition from a monoclinic crystal structure to an orthorhombic one.

The magnetic structure of the surface iron telluride layer adopts a staggered magnetic order in the a-b plane, with the spins of the iron atoms alternating between fixed angles pointing away from and parallel to the b axis direction, while the component along the c axis alternately points out of (or into) the a-b plane. This model of the magnetic order consists of a pair of coexisting spin spirals along the Fe-Fe [110] and Embedded Image directions, where the spirals alternate between right and left handedness for every other row of Fe atoms along [110]. The order remains commensurate and is similar to the plaquette order predicted theoretically in scenarios where interstitial iron is neglected and the lattice structure approaches tetragonal symmetry (10, 11) or when biquadratic exchange interactions are included (24). Djaloshinskii-Moriya interactions, which may occur for the surface layer, would promote formation of noncolinear order in the surface layer but is not expected to lead to an out-of-plane angle.

Our results suggest that a description of the interstitial iron atoms as modifying the local couplings, and thus the magnetic order, is more appropriate than there being an overall charge doping, which changes the nesting of the band structure. The latter would imply that removal of surface excess iron recovers the magnetic order found at lower excess iron concentrations, continuously changing the ordering wave vector, which, however, is not what we observe.

Our findings show some parallels to the C4 magnetic order found in the iron pnictide superconductors when the magnetostructural phase transition is suppressed (2527). In both cases, the reduction of the structural asymmetry results in formation of a double-q magnetic order, although based on a different order in the undoped compound.


Our measurements demonstrate how atomic manipulation combined with SP-STM can be used to understand the influence of defects on emergent orders (28) in strongly correlated electron materials and provide a new pathway to control emergent phases by atomic manipulation. They show that in iron telluride, it is through the monoclinic distortion that the single-q magnetic order is favored over a double-q order and that once the distortion is reduced, the magnetic structure is closer to a plaquette type of magnetic order. Our results, thus, show new opportunities to use STM-based manipulation to understand emergent magnetic orders in strongly correlated electron materials.


STM experiments

The STM experiments were performed at 2 K using a low-temperature STM equipped with a vector magnet that enables application of magnetic fields of up to 5 T in any direction with respect to the tip-sample geometry (29). To obtain a pristine, impurity-free surface for imaging, Fe1+xTe samples were prepared by in situ cleaving at a temperature of ~20 K. Pt-Ir tips were conditioned by field emission on an Au(111) sample. Ferromagnetic tips used for SP-STM measurements were prepared by picking up the interstitial Fe atoms from the Fe1+xTe sample in the STM to create a ferromagnetic cluster of Fe atoms at the tip apex (13, 14). The magnetization of the tip was controlled through the applied magnetic field. The influence of the magnetic field on the magnetic structure of Fe1+xTe was assumed to be negligible due to substantial magnetocrystalline anisotropy (13). The magnetization direction of the magnetic tip is denoted by (θ, φ), where θ represents the out-of-plane angle and φ the in-plane angle measured from the crystal axes of the FeTe sample. All STM images were obtained at 2 K.

Sample growth

Single crystals of Fe1+xTe were grown by the self-flux method (30, 31). The excess iron concentrations reported here were determined using energy-dispersive x-ray (EDX) analysis. Throughout the main text, the excess iron concentration of bulk samples (i.e., before removal of surface excess iron) refers to the off-stoichiometric part x of the composition of the material as extracted in EDX, which, in principle, can originate either from interstitial iron or a tellurium deficiency. Interstitial or excess iron refers to iron between the FeTe layers or at the surface. Characterization of the crystals indicates that excess iron concentration and interstitial iron concentration are identical within the errors of our measurements.


Supplementary material for this article is available at

Section S1. DFT calculation of magnetic contrast

Section S2. SP-STM study of the magnetic structure of Fe1.06Te

Section S3. Incommensurate order in Fe1.16Te

Section S4. Manipulating the surface excess Fe concentration

Section S5. Alternative method to determine sample spin polarization

Section S6. Model for the magnetic structure at x = 0.12

Section S7. Model for the magnetic structure at x = 0.2

Fig. S1. Spin-polarized imaging at low excess iron concentrations x < 0.12.

Fig. S2. Spin-polarized imaging at high excess iron concentrations x > 0.12.

Fig. S3. Manipulation of surface excess iron with aggressive tunneling parameters.

Fig. S4. Manipulation of surface excess iron with moderate tunneling parameters.

Fig. S5. Manipulating surface magnetic order.

Fig. S6. Extracting surface spin polarization.

Fig. S7. Simulated SP-STM images for x = 0.12.

Table S1. Crystal structure of Fe1+xTe at different excess iron concentrations x.

References (3237)

This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Acknowledgments: Funding: C.T., C.M.Y., and P.W. acknowledge funding from EPSRC through EP/L505079/1, EP/I031014/1, and EP/R031924/1, and C.S. through EP/M01052X/1. V.T. and A.L. acknowledge funding from the Deutsche Forschungsgemeinschaft (DFG) via the Transregional Collaborative Research Center TRR 80 (Augsburg, Munich, Stuttgart). C.H. acknowledges support from the Austrian Science Fund (FWF) project no. J3806-N36 and the Vienna Science Cluster. F.G. acknowledges support from the Leverhulme Trust (grant RL-2012-001) and the U.K. EPSRC Research Council (EP/M020517/1). Author contributions: P.W. designed the experiment. C.T. and C.M.Y. performed the STM measurement and analyzed the data. C.H. and F.G. performed density function theory calculations. D.C., V.T., A.L., E.E.R., and C.S. grew Fe1+xTe single crystals. C.T. and P.W. wrote the manuscript. All authors were involved in the discussion of the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Underpinning data will be made available at the University of St Andrews Research Portal ( Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article