Influence of interstitial Mn on magnetism in room-temperature ferromagnet MnSb
We report elastic and inelastic neutron scattering measurements of the high- ferromagnet MnSb. Measurements were performed on a large, K, single crystal with interstitial Mn content of . The neutron diffraction results reveal that the interstitial Mn has a magnetic moment, and that it is aligned antiparallel to the main Mn moment. We perform density functional theory calculations including the interstitial Mn, and find the interstitial to be magnetic in agreement with the diffraction data. The inelastic neutron scattering measurements reveal two features in the magnetic dynamics: i) a spin-wave-like dispersion emanating from ferromagnetic Bragg positions (HK2), and ii) a broad, non-dispersive signal centered at forbidden Bragg positions (HK2+1). The inelastic spectrum cannot be modeled by simple linear spin-wave theory calculations, and appears to be significantly altered by the presence of the interstitial Mn ions. The results show that the influence of the interstitial Mn on the magnetic state in this system is more important than previously understood.
pacs:75.50.Cc, 75.25.-j, 71.15.Mb, 78.70.Nx
MnSb is a high Curie temperature (), highly anisotropic, metallic ferromagnet. The observation of a ferromagnetic state in MnSb is unusual in its class of materials. Most Mn-alloys are antiferromagnetic, and other 3 transition metal mono-antimonides are not ferromagnetic — CrSb, FeSb and CoSb are antiferromagnetic, TiSb is paramagnetic and NiSb is diamagnetic Motizuki et al. (2009). The ferromagnetic state in MnSb is highly sensitive to substitutions; for example, Cr-doping quickly tunes the system towards an antiferromagnetic state Takei et al. (1963); Kohnke et al. (1996). In-fact substitutions of the cation, the anion, or presence of an interstitial can each alter , the magnetic anisotropy, or the type of magnetic order Kohnke et al. (1995, 1996); Takei et al. (1963); Chen et al. (1978); Reimers et al. (1982); Nakamura et al. (1992); Okita and Makino (1968); Teramoto and Van Run (1968). This potential for tuning the properties in MnSb, and closely related MnBi, has attracted considerable attention because the materials show promise as alternatives to rare-earth containing permanent magnets, and as magneto-optic mediums Anand et al. (2014); Carey et al. (1990); Coehoorn et al. (1985); Continenza et al. (2001); Di et al. (1992); Hong et al. (2013); Kharel et al. (2011); Zarkevich et al. (2014); Ashizawa et al. (2002).
The spontaneous magnetization in MnSb is along the crystallographic -axis at high temperatures. However, the anisotropy decreases on cooling and passes through zero at the spin reorientation temperature, , so that the magnetization is in-plane at low temperatures. In nominally polycrystalline MnSb the is as high as Guillaud (1949). However, single crystal studies have been unable to produce the material without the inclusion of interstitial Mn ions Okita and Makino (1968). These interstitial Mn ions, Mn2, are present in addition to the fully-occupied main site Mn ions, Mn1. The Mn2 ions are hosted within the hexagonal NiAs crystal structure as shown in Fig. 1(c). Compositions in the range of have been found to be stable.
The presence of interstitials significantly alters the properties of MnSb. As increases, the lattice parameter and unit cell volume increase, but the lattice parameter, , and the total magnetization decrease Okita and Makino (1968); Chen et al. (1978); Teramoto and Van Run (1968). Studies combining magnetization measurements with chemical analysis Okita and Makino (1968); Yamaguchi et al. (1976) found that the interstitial Mn results in a reduction of the ferromagnetic ordering temperature as captured in a simple relationship between and :
The electronic structure and magnetic state of MnSb have been investigated by a number of authors via electronic structure calculations Coehoorn et al. (1985); Continenza et al. (2001); Katoh et al. (1986); Podloucky (1984); Ravindran et al. (1999); Vast et al. (1992). They found that MnSb is described as a metallic system with localized magnetic moments, and show that the ferromagnetic state is stabilized due to the significant hybridization between the Mn 3 and Sb 5 orbitals. All these calculations, however, have failed to describe key experimental observations such as the behavior of the magnetic anisotropy Ravindran et al. (1999), and the effect of the interstitial Mn on the magnetic properties Coehoorn et al. (1985). Coehoorn et al. Coehoorn et al. (1985) simply discuss Mn2 as an electron donor, and were unable to explain its influence on the magnetic properties.
Given the importance of the interstitials in tuning the properties of MnSb, it is surprising that detailed theoretical investigations of the effect of Mn2 are lacking. Experimental results are also limited, and have focused on the empirical determination of the dependencies of cell parameters, and . This is probably due to a prevailing view in the literature that there is no magnetic moment associated with Mn2. This perspective originated from the analysis of polarized neutron diffraction data Yamaguchi et al. (1978); Reimers et al. (1983), using a model with a complex aspherical magnetic form factor for Mn1. A large asymmetry in the form factor allows the Mn1 site to contribute magnetic intensity to Bragg reflections that would otherwise be forbidden for Mn1 scattering. In these analyses, low- Bragg peaks had to be excluded to model the data. Analysis of powder neutron diffraction data Bouwma et al. (1971); Yamaguchi and Watanabe (1983), x-ray Compton profile measurements Nakamura et al. (1995), and electronic structure calculations Coehoorn et al. (1985) raised questions as to how robust the aspherical magnetic form factor model is. Additionally, a density functional theory (DFT) calculation on the influence of Pt doping in MnBi did include Mn at the interstitial site, created via a vacancy on a main Mn site, and predicted a magnetic moment of order associated with the interstitial site Kharel et al. (2011). The absence of Mn2 from most calculations may be attributable to the inherent difficulty in including a partially occupied, disordered arrangement of interstitial ions into calculations. The role of the interstitial Mn in determining the magnetic properties of MnSb, therefore, is an open question which could provide key insights into the material properties and tunability.
Here we show that the presence of the interstitial Mn ion has a strong influence on the magnetic state of MnSb, beyond the role of a simple electron donor. We use neutron diffraction to measure a large number of Bragg peaks from a single crystal sample at temperatures above and below both and . We find that the data is described across the entire -range by a simple model including a magnetic moment on the Mn2 site, aligned antiparallel to the moment on the Mn1 site. We perform DFT calculations including an interstitial Mn, and find these to carry a sizable magnetic moment in agreement with our diffraction results. The DFT results are also consistent with the observed reduction in due to the presence of Mn2. We also investigate the magnetic dynamics of MnSb via inelastic neutron scattering (INS) measurements performed on a large single crystal. We identify a spin-wave-like signal associated with the ferromagnetic Mn1 ions, and are able to partially reproduce its dispersion using a localized-moment Heisenberg model. In addition to the spin-wave-like scattering, we also identify a broad, intense, magnetic response in the inelastic spectrum. This signal is not observed in MnBi, and therefore is likely attributable to the presence of the interstitials in MnSb. These results highlight the strong influence that the presence of the interstitial Mn has on the magnetic properties of MnSb.
The single crystal of MnSb used in this investigation is the same cm crystal that was studied in Ref. Radhakrishna and Cable (1996). Pieces taken from the large crystal were used for diffraction and magnetization measurements.
To determine the spin reorientation temperature, a piece of the crystal was aligned at room temperature to within 15 of the easy axis (-axis) using a permanent magnet. The magnetization (M) was then measured in a Quantum Design Magnetic Property Measurement System, with applied fields H either parallel or perpendicular to this axis. This measurement clearly reveals a spin reorientation temperature of K, as shown in Fig. 2(a). The same crystal was then used to obtain the Curie temperature of 434 K, see Fig. 2(b). This high-temperature measurement employed the Sample Space Oven from Quantum Design, which utilizes thin quartz holders that precluded alignment of the crystal. Here, we have defined as the intercept of the greatest tangent to the M/H versus data, with the greatest tangent being observed at 423 K. According to Eq. 1, K gives an interstitial content of 16%, consistent with the value predicted from the lattice parameters in Ref. Radhakrishna and Cable (1996). With an applied field of 6 T, the moment is essentially saturated at 2 K, and for H we measure a saturated moment of 3.20/f.u., where f.u. is a formula unit MnSb. The high-temperature data shows a minor onset of magnetization magnetization at K, see Fig. 2(b). This is likely associated with a small close-to-stoichiometric MnSb impurity, or a small amount of an unknown phase.
Neutron diffraction measurements were performed on the HB-3A four-circle diffractometer at the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL), using a neutron wavelength of 1.003. Full datasets of more than 190 reflections were collected at temperatures of 10 K (), 200 K () and 450 K (), in a closed-cycle refrigerator (CCR). 450 K data were counted for 2 s per point and 200 K and 10 K data were counted 5 s per point. In addition, a set of 10 peaks, which are predicted by model 1 (see Results section) to be highly sensitive to the direction of the interstitial moment with respect to the Mn1 site moment, were counted for up to 60 s per point at all temperatures.
INS measurements were performed on the SEQUOIA spectrometer at the Spallation Neutron Source (SNS), ORNL with the crystal aligned with (H0L) in the scattering plane. Measurements were performed with an incident energy () of 150meV with chopper frequency of Hz, and with meV and Hz. This gave energy resolutions at the elastic line of 11 meV and 4 meV, respectively. A CCR was used to reach sample temperatures between 10K and 443K. Full data sets were collected at 10K, 350K and 443K with meV over at least a 130 range of sample rotation angles with a step size of 0.5. Limited angular range data sets with 1 step size were collected at 40K, 100K, 125K, 150K, 170K, 200K, 250K and 300K. Measurements were performed with meV at 10K and 350K over an 80 range with a 2 step size. The angular step data were combined, a slight misorientation of the crystal was corrected for, and cuts through the data were performed using the Horace software package hor (). Additional INS measurements were performed on the US/Japan Cold Neutron Triple-Axis Spectrometer (CTAX) at HFIR, ORNL, using guide-open-80-open collimation and fixed final neutron energy of 5meV.
Results from the neutron diffraction measurements are shown in Fig. 1. The temperature dependencies of and Bragg reflections across are shown in Fig. 1(a). These reflections have structure factor contributions from all three atoms in the unit cell. The intensities of these Bragg peaks respond to , indicating that they are magnetic Bragg reflections, and confirming the identification of from magnetization. Figure 1(b) shows the Bragg reflection at 443 K and 10 K. The Mn1 site in the NiAs crystal structure has the reflection condition , where is an integer. Therefore, the peak has zero nuclear contribution to its structure factor from Mn1, and also zero magnetic contribution if the magnetic form factor of Mn1 is spherical. Both Mn2 and Sb ions contribute scattering intensity to the reflection.
To determine the origin of the increased intensity at low temperature of , Bragg peaks, we investigated two possible models for the magnetic state in MnSb for all reflections collected. Both models use the NiAs structure, hexagonal space group (No. 194), which has two formula units of MnSb per unit cell, with Mn1 on the 2 Wyckoff site (0 0 0), Sb on the 2 Wyckoff site , and Mn2 on the 2 Wyckoff site . We use standard Miller indices (h k l) to index the reflections. Note that because of the hexagonal symmetry the Miller-Bravais notation can be used to illustrate equivalent reflections, given by (h k i l) where i=-(h+k) and permutations of h k i give equivalent reflections Frank (1965). Recent x-ray diffraction results from the closely related compound MnBi identify a slight distortion from hexagonal symmetry below , however the resolution of our neutron diffraction measurement is not sufficient to detect a distortion of this size McGuire et al. (2014). The difference between the two models we investigate concerns the magnetic component of scattering below .
The first model comprises a ferromagnetic arrangement of Mn1 moments, with Mn2 moments aligned antiparallel to Mn1. A depiction of this model is shown in Fig. 1(c) for , i.e. with spins aligned within the - plane. For the Mn1 and Mn2 moments are aligned along the -axis, but still antiparallel to each other. The 3-5 hybridization is expected to induce a small moment on the Sb site antiparallel to Mn1 Coehoorn et al. (1985). However, considering the expected moment size and the steep magnetic form factor for Sb 5 electrons Yamaguchi et al. (1978), the overall contribution of Sb to the magnetic diffraction pattern is expected to be small. Therefore we do not include magnetic intensity from Sb in the fitting for either model 1 or model 2, consistent with previous reports Bouwma et al. (1971); Yamaguchi et al. (1978); Reimers et al. (1983). As Sb and Mn2 have the same reflection conditions, this may have the effect of slightly increasing the moment size assigned to the interstitial Mn.
The model was fit to the data using Rietveld refinement in the FullProf software suite Rodríguez-Carvajal (1993). For a random distribution of interstitial Mn ions, the average structure can be modeled by assuming a uniform distribution of Mn2 on every 2 Wyckoff site and scaling the scattering intensity by the occupancy . The atomic displacement parameters for Sb and Mn2 were constrained to be equal, as they contribute to the same reflections. The initial fit was performed against the K (i.e. ) data set, allowing us to determine the interstitial content of the sample without the influence of magnetic Bragg scattering. The standard empirical correction for extinction used in FullProf was utilised in this refinement Rodríguez-Carvajal (1993), resulting in an extinction coefficient of 580(60) which gives a maximum extinction correction of 68 % on the strongest Bragg peak. Using this fitting procedure the value of the interstitial content was determined to be , which is reasonably close to the value estimated from comparison of the magnetization measurement (see Fig. 2) to Eq. 1. For subsequent fits of the magnetic model against 10 K and 200 K data sets the interstitial content and the extinction parameter were kept fixed. The magnetic form factor of Mn was used Reimers et al. (1983); Dianoux and Lander (2003). For the 10 K data set, the model included three domains in equal proportions with spins along , , and [-1 -1 0] directions respectively. The results of these fits are summarized in Fig. 3(a), (b) and (c) and Table 1. Excellent agreement between data and model 1 is found.
|Temperature (K)||()||()||Vol. ()||()||()||ext||scale|
To highlight the magnetic component of the scattering, Fig. 3(e) shows a comparison of the 10K data minus the calculated structural component of the scattering, with the calculated magnetic contribution to the scattering. The intensity of the magnetic contribution is seen to decrease with for both data and calculation, as expected due to the magnetic form factor. The effects of the structure factor and polarization factor, which selects scattering from only moments perpendicular to , in the neutron scattering cross section are present in this data, which is why there is not an exact form factor dependence. In addition, some forbidden Bragg positions were measured and appear at zero integrated intensity. Attempts to fit an interstitial Mn moment aligned parallel to Mn1, instead of antiparallel, were unsuccessful. The model with an antiparallel moment on the Mn2 site describes the data very well.
The second model is based on that proposed by Haneda et al. Haneda et al. (1977) for MnAs and discussed in detail for MnSb by Yamaguchi et al. Yamaguchi et al. (1978), in which Mn2 does not have a magnetic moment associated with it. In this case, they explain the observed magnetic scattering at and other odd Bragg positions Bouwma et al. (1971); Yamaguchi et al. (1978); Reimers et al. (1983) using a highly aspherical magnetic form factor for the main site Mn. Having an aspherical magnetization density breaks the symmetry conditions that normally result in systematic absences for reflections from the 2 Wyckoff site, allowing intensity at , positions. The model includes the combined scattering intensity for structural contributions from all ions, plus the magnetic intensity from the Mn1 ions for all reflections. We have added to model 2 the same extinction correction as used in model 1. The overall structure factor resulting from the Yamaguchi model, with adjusted interstitial content and electron occupancies to match our data, is compared to the 10 K data in Fig. 3(d). For this model we find and , compared with and for model 1 at 10 K (Table 1). The purely magnetic component of model 2 is shown in Fig. 3(f), along with the difference between the data and calculated structural component. In this case there is significant deviation between the model and the data over a large range of , with the intensities of only some magnetic peaks being replicated. Therefore model 2 is not as effective as model 1 in describing the results of our neutron diffraction experiment. This will be examined further in the Discussion Section.
To gain a better understanding of the influence of the Mn interstitials on MnSb we performed DFT calculations including an interstitial in the unit cell (see dft () for technical details). First we readdress the question of whether the interstitial Mn has a magnetic moment, and if so how it is directed with respect to the ferromagnetically aligned moments of the host Mn atoms. To this end we considered the 221 supercell MnSb depicted in Fig. 4, which contains a single interstitial Mn atom. To account for the influence of Coulomb interactions among the Mn-3 electrons we use the PBE+U approximation, and vary the Hubbard parameter between 0 and 8 eV. For all cases, we find that it is not possible to stabilize a configuration in which the interstitial Mn is non-magnetic. The tables in Fig. 4 show that for , 2 and 4 eV the configuration with Mn2 aligned parallel to Mn1 is higher in energy than the configuration in which it is aligned antiparallel, consistent with the results from neutron diffraction. For and 8 eV, however, the configuration with the parallel alignment of the interstitial Mn is energetically favored. It is reasonable to expect that the Coulomb - interactions are at the lower end of the scale, given that the very large wave functions of the Sb 5 electrons can effectively screen the transition metal interactions. Related to this, we find for all cases that the system remains metallic, even for a Hubbard of 8 eV.
Next we examine the effect of the Mn interstitials on the ferromagnetic configuration of the host Mn moments. To investigate the stability of the ferromagnetic ground state (FM) we compare its energy with that of two higher energy configurations in which the atomic positions remain unchanged, but the host Mn moments assume an antiferromagnetic configuration either along the -axis (AFM-c) or the -axis (AFM-a) as shown in Fig. 5. Comparing Fig. 5(a) and (c), the presence of the Mn interstitial lowers the energy difference between the FM ground state and the high energy AFM configurations. The energy difference between the FM and the AFM-c configuration reduces quite significantly from 286 to 212 meV and the energy difference between the FM and the AFM-a configuration decreases from 326 to 306 meV. The results in Fig. 5 are for a eV, but for a of 8 eV we find that this qualitative conclusion remains unchanged U8e ().
It is important to note that the simulations in Fig. 5(a) and (c) include the indirect influence of Mn2 via their induced changes in the lattice constants. From comparing the lattice parameters between the undoped and doped system in Fig. 5(a) and (c) we see that, just like in the experimental observations Teramoto and Van Run (1968), the Mn interstitials tend to increase the lattice parameter while decreasing the lattice parameter. This leads to the question of whether the change in lattice parameters is solely responsible for the damaging influence of the Mn interstitials on the FM state. To address that question, we present in Fig. 5(b) the stability of the FM state again, for a system which has the strained lattice parameters of the doped system but does not include the Mn interstitials. The energy differences of this undoped strained system shown in Fig. 5(b) differ significantly from those in the doped system shown in Fig. 5(c). This illustrates that Mn2 affects the ferromagnetic state not just via the change of lattice parameters, suggesting that the magnetic exchange interactions between the host and the interstitial Mn moments and/or doping effects play an important role. Previous attempts to explain the magnetic properties of MnSb just via doping effects were not successful Coehoorn et al. (1985).
Inelastic neutron scattering
To investigate the magnetic state in MnSb further, we performed INS experiments to probe the magnetic dynamics. An overview of the results from SEQUOIA is given in Figs. 6(a) and 7, which show HL-plane and -plane color maps, respectively, as well as constant-energy cuts through the data in Fig. 7(d). These inelastic scattering spectra contain contributions from both magnetic and phonon scattering. However, the phonon scattering contribution is effectively negligible in the results we present. Mn and Sb have nuclear scattering cross sections of opposite sign, therefore their scattering contributions largely cancel-out and as a result the acoustic phonon scattering from the sample is weak. Additionally, the phonon intensity follows a dependence and is therefore weak in the low- zones we investigate. In the low temperature data the phonon scattering is suppressed by the Bose-population factor in the energy range investigated, whereas at higher temperature phonons (attributable to both sample and background contributions) simply result in an overall increase in the background scattering, particularly at low energies, as evident in Fig. 6(a). Therefore the results presented are attributable to the magnetic dynamics of MnSb.
In Fig. 6(a) rings of scattering in the HL plane are seen dispersing with energy out of the ferromagnetic positions such as and ; compare to Fig. 6(c) which depicts the reciprocal lattice plane. Only L = even positions are magnetic zone centers because the Mn1 magnetic lattice is half the size of the structural lattice along the -axis, see Fig. 1(c). The rings persist above both and . The intensity of the rings is rapidly suppressed with increasing at all temperatures, and it persists up to energies of meV (see Fig. 7(a)), well above the phonon cutoff of the sample meV, indicating the magnetic origin of the signal. We observe this signal down to low energies in the data from our measurement on CTAX, Fig. 8. This shows that the dispersion of the magnetic signal along H can already be observed out of with an energy transfer of 0.5 meV, see Fig. 8. This dispersive signal is reasonably reproduced by a spin-wave Heisenberg model for Mn1 ions, presented in Fig 9, discussed further below, and a qualitatively similar signal is observed in MnBi wil (), which contains no interstitial Mn, see Fig. 6(b).
In the inelastic spectra, in addition to the spin-wave type signal, a broad signal is evident between the rings in the HL plane at all temperatures, see Fig. 6(a). This broad signal was not observed in MnBi, corroborating that it is a distinct feature of the excitation spectrum of MnSb, Fig. 6(b). The broad signal is centered on and equivalent positions. is not an allowed Bragg reflection in this system (see Fig. 6(c)), and no Bragg peak was observed at this position in any of the measurements we performed. Figure 7(c) highlights that this signal is broad but localized in -space, and seemingly non-dispersive. Again, the intensity suppression with and the energy range of the signal at all temperatures both indicate that it is magnetic in origin. We observe this signal down to the lowest energies we were able to resolve in our experimental set-up on SEQUOIA, meV, but are unable to determine if there is an associated elastic diffuse signal. This signal was previously observed at in the triple-axis-spectroscopy work of Radhakrishna and Cable Radhakrishna and Cable (1996), but we have now been able to map its full and dependence and confirm the signal’s magnetic origin, in addition to observing it at .
Further investigation of the scattering centered on is presented in Fig. 10. Cuts taken along and directions through datasets collected at 10 K, 350 K and 443 K are shown in Fig. 10(a) and (b). As the direction is the zone boundary for the spin wave signal (see Fig. 6(c)) the cuts in Fig. 10(b) and the color map in Fig. 7(c) essentially show the broad signal in isolation. The direction, however, passes through the zone center for both the spin wave and the broad signal, therefore both are observed in Figs. 10(a) and 7(d). We show the temperature dependence of the signal in Fig. 10(c). The integrated intensity was determined from cuts made along the direction through the limited-angular-range data measured at meV for all temperatures. The data were averaged over 10–20 meV and folded along H and L to improve statistics. Two Gaussian functions on a flat background were fit to the data, with widths constrained to be equal, but amplitudes varying independently, and centers fixed at H=0 and H=1. The resulting area of the Gaussian centered at H=0 was corrected for the Bose population factor, , at each temperature and the result gives the integrated intensity plotted in Fig. 10(c).
These results show that the broad, -centered scattering is influenced by , Fig. 10(c). There is a dramatic increase in intensity of the signal on warming past , which is expected for scattering from a transverse magnetic fluctuation as the spins reorient from -plane to -axis alignment, because the magnetic cross section depends only on the component of the magnetic moment perpendicular to . This is further supported by the observation that the (0 0 1) scattering is more sensitive to than the scattering at (1 0 1), compare the cuts in Fig. 10(b). These observations are a further confirmation of the magnetic origin of the signal, and indicate coupling between this scattering and the ferromagnetic state of the system.
We attempted to model the inelastic magnetic response of the system using linear spin-wave theory. Isotropic exchange interactions including up to sixth nearest neighbors were the minimum required to model the spin-wave frequencies. The six nearest neighbor exchange parameters, s, are illustrated in Fig. 9(c). Single-ion anisotropy was not included because there was no evidence for a gap in the spin-wave spectrum down to at least 0.5 meV (<1 % of the bandwidth), see Fig. 8. The inelastic neutron cross section for undamped spin waves was calculated using the 1/S formalism outlined in Ref. Haraldsen and Fishman, 2009 and appendix A of Ref. Fishman et al., 2013. For comparison with experimental intensities, the effects of the magnetic form factor and an approximation for the instrumental resolution were included in the calculation. We used the magnetic form factor for Mn from Ref. Dianoux and Lander, 2003. The resolution function was approximated as a Gaussian in energy with a full width at half-maximum of 7.5 meV and the results averaged over the same volume of as the experimental data (Fig. 7) and averaged over six domains. To determine the exchange parameters, the model was fit to the dispersions along (1 0 0), (0 0 1) and (1 1 0) directions, extracted from the experimental data via a series of constant-energy cuts and fits like those in Fig. 7(d). We determined the energy at the (H 0 2) zone boundary, 70(4) meV, by fitting a Gaussian to a constant- cut at (1.5 0 2) from which a similar cut at (2 0 2) had been subtracted as a background.
The results of this modeling are shown in Fig. 9(a) and (b), which show reasonable correspondence to the spin-wave data presented in Fig. 7(a) and (b). However, the bandwidth of the dispersion does not appear to be well reproduced along the H direction, see Fig.. 9(a). Only can be uniquely determined from our results. Instead, the linear combinations
were fit. The values best fitting these measurements are = 20.0 meV, = -0.8 meV, = 4.8 meV, and = 2.9 meV. Fig. 9 shows the calculated frequencies and intensities along two of the measured directions, for comparison to Fig. 7(a) and (b). Previously Refs. Yamaguchi et al. (1992) and Radhakrishna and Cable (1996) reported triple-axis spectroscopy measurements on MnSb which covered limited ranges of and , not reaching the zone boundary. In Ref Yamaguchi et al. (1992) they fit a Heisenberg model up to fourth nearest neighbors (i.e. for comparison to our model) to the low energy data measured on a MnSb crystal. They found meV, meV and meV. These results extrapolate to a bandwidth of 130 meV along (H 0 2), which is much larger than the 70(4) meV we observe, see Fig. 7.
The results of our neutron diffraction experiment on MnSb are well described by the magnetic structure shown in Fig. 1(c), with ferromagnetic alignment of Mn1 ions, and Mn2 ions aligned antiparallel to Mn1. There appears to be no need to include the model of an aspherical form factor for Mn1 to account for magnetic scattering at , Bragg peak positions for our data set. This is consistent with early neutron powder diffraction results on MnSb Bouwma et al. (1971), but in contrast to the interpretations of Yamaguchi et al. Yamaguchi et al. (1978) and Reimers et al. Reimers et al. (1983) from their polarized-neutron, single-crystal experiments.
The single crystals investigated in Refs. Yamaguchi et al. (1978) and Reimers et al. (1983) were of composition MnSb and MnSb respectively, lower interstitial contents than our crystal, and in both cases they excluded some low data from their analysis. This may explain the seeming discrepancy between their interpretation and our results. In the former case, the scattering from Mn2 sites would have been weaker, and therefore harder to interpret. The higher data that they included in their analysis is less sensitive to magnetic scattering from the interstitial site, but more sensitive to the asymmetric form factor they describe. In the case of our data, only model 1, including an antiparallel moment on Mn2, can account for all peaks observed below across all , see Fig. 3(c)–(f). Our data is less sensitive to any asymmetry in the Mn1 form factor because it was not a polarized neutron measurement and the magnetic contribution at high is small. Our results do not preclude some asphericity in the magnetic form factor of Mn1. It is likely, however, that the asphericity found in the previous models is too large, as some of the scattering intensity should have been accounted for by the moment on Mn2. This helps to explain the discrepancies between the Yamaguchi model Yamaguchi et al. (1978) and results from x-ray Compton profile analysis Nakamura et al. (1995), and electronic structure calculations Coehoorn et al. (1985), as well as the failure of the model at low Yamaguchi et al. (1976, 1978); Reimers et al. (1983).
The size of the magnetic moments found from our FullProf refinement of model 1 are consistent with magnetization measurements. The saturated moment per formula unit was found to be 3.20 at 2 K from magnetization. The results from HB-3A give a total moment per formula unit of 3.2(1) at 10 K, taking into account the Mn2 occupancy of , see Table 1. The determination of a magnetic moment associated with the Mn2 site is further supported by the results from our DFT calculation. Previous modeling of MnSb did not include an interstitial Mn, and those that attempted to account for it merely adjusted the electron count in stoichiometric MnSb models Coehoorn and Groot (1985); Continenza et al. (2001); Katoh et al. (1986); Podloucky (1984); Ravindran et al. (1999); Sandratskii et al. (1981); Vast et al. (1992); Chen et al. (1978). We find a sizable moment on Mn2 in the calculation, see Fig. 4.
Given the presence of a magnetic moment on the interstitial site, in a simple local Heisenberg picture we would expect the antiparallel alignment of the interstitial to enhance the overall ferromagnetic state. However, it has clearly been observed that the presence of the interstitial acts to reduce Chen et al. (1978); Okita and Makino (1968); Teramoto and Van Run (1968). Our DFT calculation shows that the interstitial Mn degrades the stability of the ferromagnetic state beyond the simple influence of changing lattice parameters (see Fig. 5), in agreement with the experimental observations. This indicates that a purely localized moment model may not be sufficient to describe the magnetic state in MnSb.
Our inelastic neutron scattering measurements of MnSb serve to highlight the influence the interstitial Mn has on the magnetic state. The spin-wave-type excitations appear to be less well defined than in the non-interstitial containing sister compound MnBi (see Fig. 6), and are modified such that a simple Heisenberg model does not describe the observed dispersion. For a local moment system, linear spin-wave theory is expected to capture the form of the dispersion by inclusion of a small number of nearest neighbor exchange parameters. Although the observed dispersion (e.g. Fig. 7 (a)) appears to be associated with the ferromagnetic wavevector, and a similar feature is observed in interstitial-free MnBi (Fig. 6(b) wil ()), we have shown that the measurement is not fully described by the Heisenberg model for localized Mn1 spins, Fig. 9.
For the spin wave calculation along (H 0 2), attempting to match both the bandwidth and the gradient of the dispersion at lower energies requires a large interaction (2.9 meV), which probably reflects the inadequacy of the model. The shape at the zone boundary remains quite different, resulting in the 70(4) meV bandwidth not being reproduced, see Fig. 9(a). The calculation does not produce a signal emanating from (0 0 1) as is observed in the data, see Fig. 7(c). As (H 0 1) is a zone boundary for the Mn1 magnetic lattice, as illustrated in Fig. 6(c), the spin-wave intensity at (0 0 1) is concentrated around the zone boundary energy (see Fig. 9(b)). Our crude attempts to extend the model to include an interstitial Mn2 moment failed to better reproduce the spin wave data, nor did they produce a signal like the observed scattering along . These attempts are limited at the present time by the absence of general theoretical and computational tools to include the effects of a disordered component in a system on the inelastic magnetic response. Including a small moment at one interstitial site in our model calculation artificially introduces a symmetry element which is not present in the real system. Extending the calculation to include thousands of unit cells, which can then be randomly assigned interstitials, is currently difficult to implement and extremely computationally intensive.
In addition to the spin wave, we identified a second magnetic signal in the inelastic spectrum which is broad in , and centered on an unexpected position in reciprocal space (see Figs. 6, 7 and 10). Like the spin-wave scattering, this signal responds to , but unlike the spin-wave scattering, no similar signal is observed in interstitial-free MnBi. is not an allowed nuclear or magnetic reflection (for either magnetic structure model discussed above). The presence of a Bragg peak would require the correspondence between the upper and lower halves of the unit cell shown in Fig. 1(c) to be broken, for example, if Mn1 was antiferromagnetically aligned along , or if Mn2 was only present in the lower half of the cell. We do not observe a Bragg peak at in any of the neutron data sets collected, implying that no such long-range order is present. However, the periodicity of the signal in reciprocal space implies that underlying long-range correlations are present.
Our results show that the dynamic correlations observed at are not associated with a long-range order in the system. They could, however, be associated with a short-range structural or magnetic order in the system. This proposition may be supported by the broad appearance of the signal in reciprocal space (see Fig. 6), which suggests a short correlation length in real space. We made a rough estimate for the correlation length by fitting Gaussians to cuts similar to those in Fig. 10(a) and (b), integrated over 10–15 meV, from the meV, K data. We find correlation lengths of ~16 in the - plane, and in the range 10–18 along (where it is harder to define, due to the overlap of the signal with the spin-wave-type scattering), i.e. ~4 unit cells in plane, and ~3 units cells along . We were unable to determine whether there is an elastic, diffuse scattering signal at the position due to the neutron instrumentation used in this investigation. We do observe the signal down to the lowest energies that were probed on SEQUOIA, meV. A short-range order could be present in the system if the interstitial Mn ions are not randomly distributed, but instead are structurally ordered over a few unit cells. An observation of elastic, diffuse scattering at (0 0 1) would establish that short-range order is present. Alternatively, the presence of interstitials could modify the interactions of neighboring main-site Mn1 ions, with the associated magnetic excitations unable to propagate over large length scales due to disruption from a random distribution of magnetic Mn2 ions. Further studies, including diffuse neutron scattering experiments, are highly desirable to investigate the potential short-range order in this system.
The spectral weight of the scattering, i.e. the total intensity of the signal integrated over space, can in principal be used to indicate the strength of the magnetic excitations and therefore the size of the fluctuating magnetic moment, which could give an indication of the origin of the scattering. Unfortunately, due to the overlap of the two magnetic excitations in space, we cannot quantitatively compare the spectral weights of the two signals. Qualitatively, however, by inspecting a series of slices and cuts through the data, such as those shown in Figs. 6, 7 and 10, we see that the strength of the signal emanating from is comparable to the spin-wave signal. This implies that the (0 0 1) signal is associated with a relatively large magnetic moment per formula unit in the sample, of the order of the Mn1 moment size. This makes it unlikely to be purely the result of the Mn2 moments behaving independently of Mn1. Identification of the origin of this signal would be extremely significant for understanding the magnetic state in MnSb and determining the possibilities of tuning the state with interstitial ions of Mn or other 3 transition metals.
We have shown that a magnetic moment on the interstitial Mn in MnSb is required to describe the neutron diffraction data from a single crystal of MnSb. This magnetic moment is aligned antiparallel to the main site Mn magnetic moment. We performed DFT calculations which find that the interstitial Mn is magnetic, consistent with the diffraction data. The DFT results are also consistent with previous experimental evidence that the interstitial reduces and lattice parameter, and increases in MnSb. We find that the presence of the magnetic interstitial Mn also has a substantial effect on the magnetic excitation spectrum of MnSb, measured by inelastic neutron scattering. It results in the appearance of an intense, broad signal in -space, that is likely associated with short-range order. While a Heisenberg Hamiltonian model calculation captures the essential -space dependence of the spin-wave type signal observed, the additional broad signal cannot be explained by straightforward attempts to include an interstitial ion in this model.
We thank J. Q. Yan, D. Mandrus, M. A. McGuire and B. C. Sales for useful discussions. The research at ORNL’s Spallation Neutron Source and High Flux Isotope Reactor was supported by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy (DOE). A.F.M. and R.F. were supported by the U.S. DOE, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division. T.B. and T.J.W. are supported as Wigner Fellows at ORNL. Work by TB was performed at the Center for Nanophase Materials Sciences, a DOE Office of Science user facility. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. DOE under Contract No. DE-AC02-05CH11231. S.E.H. acknowledges support by the Laboratory’s Director’s fund, ORNL.
This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).
- Motizuki et al. (2009) K. Motizuki, H. Ido, T. Itoh, and M. Morifuji, Electronic Structure and Magnetism of 3d-Transition Metal Pnictides (Springer Science & Business Media, 2009).
- Takei et al. (1963) W. J. Takei, D. E. Cox, and G. Shirane, Phys. Rev. 129, 2008 (1963).
- Kohnke et al. (1996) H. Kohnke, C. Kleeberg, V. Dankelmann, K. Bärner, and J. Schünemann, J. Alloy Compd 239, 150 (1996).
- Kohnke et al. (1995) H. J. Kohnke, V. Dankelmann, C. Kleeberg, J. W. Schünemann, K. Bärner, A. Vetcher, and G. A. Govor, Phys. Stat. Sol. (b) 191, 511 (1995).
- Chen et al. (1978) T. Chen, D. Rogowski, and R. M. White, J. Appl. Phys. 49, 1425 (1978).
- Reimers et al. (1982) W. Reimers, E. Hellner, W. Treutmann, and G. Heger, J. Phys. C: Solid State Phys. 15, 3597 (1982).
- Nakamura et al. (1992) J. Nakamura, Y. Kobayashi, F. Ambe, K. Asai, and N. Yamada, Hyperfine Interact. 68, 311 (1992).
- Okita and Makino (1968) T. Okita and Y. Makino, J. Phys. Soc. Jpn. 25, 120 (1968).
- Teramoto and Van Run (1968) I. Teramoto and A. Van Run, J. Phys. Chem. Solids 29, 347 (1968).
- Anand et al. (2014) K. Anand, J. Pulikkotil, and S. Auluck, J. Magn. Magn. Mater. 363, 18 (2014).
- Carey et al. (1990) R. Carey, B. Thomas, and G. Bains, J. Magn. Magn. Mater. 83, 21 (1990).
- Coehoorn et al. (1985) R. Coehoorn, C. Haas, and R. A. de Groot, Phys. Rev. B 31, 1980 (1985).
- Continenza et al. (2001) A. Continenza, S. Picozzi, W. T. Geng, and A. J. Freeman, Phys. Rev. B 64, 085204 (2001).
- Di et al. (1992) G. Di, S. Iwata, S. Tsunashima, and S. Uchiyama, J. Magn. Magn. Mater. 104-107, 1023 (1992).
- Hong et al. (2013) Y.-K. Hong, J. Park, O. N. Mryasov, S.-G. Kim, S. Kim, J. Lee, G. S. Abo, C.-J. Choi, and J. Lee, AIP Advances 3, 052137 (2013).
- Kharel et al. (2011) P. Kharel, R. Skomski, P. Lukashev, R. Sabirianov, and D. J. Sellmyer, Phys. Rev. B 84, 014431 (2011).
- Zarkevich et al. (2014) N. A. Zarkevich, L.-L. Wang, and D. D. Johnson, APL Materials 2, 032103 (2014).
- Ashizawa et al. (2002) Y. Ashizawa, S. Saito, and M. Takahashi, J. Appl. Phys. 91, 8240 (2002).
- Guillaud (1949) C. Guillaud, Ann. Phys-Paris 4, 671 (1949).
- Yamaguchi et al. (1976) Y. Yamaguchi, H. Watanabe, and T. Suzuki, J. Phys. Soc. Jpn. 41, 703 (1976).
- Katoh et al. (1986) K. Katoh, A. Yanase, and K. Motizuki, J. Magn. Magn. Mater. 54-57, 959 (1986).
- Podloucky (1984) R. Podloucky, Solid State Commun. 50, 763 (1984).
- Ravindran et al. (1999) P. Ravindran, A. Delin, P. James, B. Johansson, J. M. Wills, R. Ahuja, and O. Eriksson, Phys. Rev. B 59, 15680 (1999).
- Vast et al. (1992) N. Vast, B. Siberchicot, and P. G. Zerah, J. Phys.: Condens. Matter 4, 10469 (1992).
- Yamaguchi et al. (1978) Y. Yamaguchi, H. Watanabe, and T. Suzuki, J. Phys. Soc. Jpn. 45, 846 (1978).
- Reimers et al. (1983) W. Reimers, E. Hellner, W. Treutmann, and P. Brown, J. Phys. Chem. Solids 44, 195 (1983).
- Bouwma et al. (1971) J. Bouwma, C. F. van Bruggen, C. Haas, and B. van Laar, Le Journal de Physique Colloques 32, C1 (1971).
- Yamaguchi and Watanabe (1983) Y. Yamaguchi and H. Watanabe, J. Magn. Magn. Mater. 31-34, 619 (1983).
- Nakamura et al. (1995) J. Nakamura, T. Takeda, K. Asai, N. Yamada, Y. Tanaka, N. Sakai, M. Ito, A. Koizumi, and H. Kawata, J. Phys. Soc. Jpn. 64, 1385 (1995).
- Radhakrishna and Cable (1996) P. Radhakrishna and J. W. Cable, Phys. Rev. B 54, 11940 (1996).
- (31) Horace: Software for visualising and manipulating S(,) measured in all four dimensions. T. G. Perring et al. http://horace.isis.rl.ac.uk.
- Frank (1965) F. C. Frank, Acta Crystallographica 18, 862 (1965).
- McGuire et al. (2014) M. A. McGuire, H. Cao, B. C. Chakoumakos, and B. C. Sales, Phys. Rev. B 90, 174425 (2014).
- Rodríguez-Carvajal (1993) J. Rodríguez-Carvajal, Physica B 192, 55 (1993).
- Dianoux and Lander (2003) A.-J. Dianoux and G. Lander, Neutron Data Booklet, 2nd ed. (Institut Laue-Langevin, 2003).
- Haneda et al. (1977) S. Haneda, N. Kazama, Y. Yamaguchi, and H. Watanabe, J. Phys. Soc. Jpn. 42, 1201 (1977).
- (37) We employ the plane wave projector augmented wave (PAW) method Blöchl (1994) as implemented in the VASP code Kresse and Hafner (1993); Kresse and Furthmüller (1996); Kresse and Joubert (1999). Atomic relaxations are converged down to 1 meV/Å. The number of plane waves is determined by an energy cut-off of 650 eV. For the MnSb, MnSb and MnSb unit cells a 8x8x5, 4x8x5 and 4x4x5 k-mesh is employed respectively. For the PBE+U calculations the double counting correction scheme by Dudarev et al. Dudarev et al. (1998) is used. Atomic images were produced with the VESTA program Momma and Izumi (2011).
- (38) For eV we find that the energy difference between the FM and the AFM-c configuration reduces from 212 to 156 meV and the energy difference between the between the FM and the AFM-a configuration lowers from 219 to 185 meV.
- (39) T. J. Williams et. al, unpublished.
- Haraldsen and Fishman (2009) J. T. Haraldsen and R. S. Fishman, J. Phys.: Condens. Matter 21, 216001 (2009).
- Fishman et al. (2013) R. S. Fishman, J. T. Haraldsen, N. Furukawa, and S. Miyahara, Phys. Rev. B 87, 134416 (2013).
- Yamaguchi et al. (1992) Y. Yamaguchi, S. Tomiyoshi, M. Harada, and G. Shirane, J. Magn. Magn. Mater. 103, 50 (1992).
- Coehoorn and Groot (1985) R. Coehoorn and R. A. d. Groot, J. Phys. F: Met. Phys. 15, 2135 (1985).
- Sandratskii et al. (1981) L. M. Sandratskii, R. F. Egorov, and A. A. Berdyshev, Phys. Stat. Sol. (b) 103, 511 (1981).
- Blöchl (1994) P. E. Blöchl, Phys. Rev. B 50, 17953 (1994).
- Kresse and Hafner (1993) G. Kresse and J. Hafner, Phys. Rev. B 48, 13115 (1993).
- Kresse and Furthmüller (1996) G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996).
- Kresse and Joubert (1999) G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999).
- Dudarev et al. (1998) S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Phys. Rev. B 57, 1505 (1998).
- Momma and Izumi (2011) K. Momma and F. Izumi, J. Appl. Crystallogr. 44, 1272 (2011).