Short- and long-range magnetic order in LaMnAsO

Short- and long-range magnetic order in LaMnAsO


The magnetic properties of the layered oxypnictide LaMnAsO have been revisited using neutron scattering and magnetization measurements. The present measurements identify the Néel temperature  = 360(1) K. Below the critical exponent describing the magnetic order parameter is  = 0.330.35, consistent with a three dimensional Heisenberg model. Above this temperature, diffuse magnetic scattering indicative of short-range magnetic order is observed, and this scattering persists up to  = 650(10) K. The magnetic susceptibility shows a weak anomaly at and no anomaly at . Analysis of the diffuse scattering data using a reverse Monte Carlo algorithm indicates that above nearly two- dimensional, short-range magnetic order is present with a correlation length of 9.3(3) Å within the Mn layers at 400 K. The inelastic scattering data reveal a spin-gap of 3.5 meV in the long-range ordered state, and strong, low-energy (quasi-elastic) magnetic excitations emerging in the short-range ordered state. Comparison with other related compounds correlates the distortion of the Mn coordination tetrahedra to the sign of the magnetic exchange along the layer-stacking direction, and suggests that short-range order above is a common feature in the magnetic behavior of layered Mn-based pnictides and oxypnictides.

I Introduction

The study of layered iron pnictide and chalcogenide superconductors has led to the discovery of other interesting properties in isostructural analogues. The structure types adopted by these materials are amenable to a variety of chemical substitutions, which allow access to a wide range of physical behaviors, including multiple types of superconducting, magnetic, and electronic states. This is especially true of the 1111-type materials with the ZrCuSiAs structure type Pöttgen and Johrendt (2008) adopted by the prototype iron-based, high-temperature superconductor parent phase LaFeAsO Quebe et al. (2000).

Soon after the discovery of superconductivity in F-doped LaFeAsO Kamihara et al. (2008), substitution of La with other rare-earth ions (Ln) was seen to significantly increase the superconducting transition temperature Chen et al. (2008); Ren et al. (2008), and to induce superconductivity in the case of partial substitution of Th Wang et al. (2008). Arsenic can be replaced by antimony in the case of LnMnSbO and LnZnSbO Schellenberg et al. (2008), or by phosphorus, producing low temperature superconductors LnFePO with anisotropic gap structures Kamihara et al. (2006); Fletcher et al. (2009). Oxygen can be fully replaced by fluorine to form compounds like AFeAsF (A = Ca, Sr), whose properties are quite similar to those of LnFeAsO, including coupled structural and antiferromagnetic phase transitions and superconductivity when doped, due to the presence of isoelectronic FeAs layers Tegel et al. (2008); Matsuishi et al. (2008a, b); Han et al. (2008). Doping other transition metals on the Fe site can induce superconductivity Sefat et al. (2008), while full replacement of Fe produces diverse and interesting phenomena. These include giant magnetoresistance and rare-earth magnetism in LnMnAsO Marcinkova et al. (2010a); Emery et al. (2010); Tsukamoto et al. (2011); Wildman et al. (2012), itinerant ferromagnetism in LnCoAsO Yanagi et al. (2008); Ohta and Yoshimura (2009); McGuire et al. (2010); Marcinkova et al. (2010b), low temperature superconductivity and Kondo physics in LnNiAsO Watanabe et al. (2008); Li et al. (2008); Matsuishi et al. (2012); Luo et al. (2011), and glassy magnetism in LnRuAsO McGuire et al. (2014).

Of particular interest is the series of Mn-based 1111 materials LnMnAsO, which are antiferromagnetic (AFM) insulators exhibiting giant magnetoresistance Emery et al. (2010); Wildman et al. (2012) and, when doped, metal-insulator transitions Sun et al. (2012); Hanna et al. (2013) and low temperature crystallographic distortions Wildman et al. (2015). Electronic properties of LaMnAsO thin films indicate p-type conduction and a band gap of 1.01.4 eV Kayanuma et al. (2009), and calculations suggest that LaMnAsOLaZnAsO alloys may provide a source of spin-polarized carriers for spintronic applications Li et al. (2013).

CeMnAsO, PrMnAsOF, and NdMnAsO have been shown via neutron diffraction to have antiferromagnetic ordering of Mn moments below = 347, 340 and 359 K, respectively Zhang et al. (2015); Wildman et al. (2015); Marcinkova et al. (2010a). Powder neutron diffraction measurements on LaMnAsO have shown AFM order at 300 K, and the absence of magnetic Bragg reflections at 400 K Emery et al. (2011). Based on magnetization measurements, values of 317325 K have been reported for the of LaMnAsO Emery et al. (2010); Sun et al. (2012).

Like LnMnAsO, the analogous 122 material BaMnAs is also an antiferromagnetic Mott insulator An et al. (2009); Singh et al. (2009a, b); Johnston et al. (2011), and dilute magnetic semiconductor behavior has been observed in materials produced by simultaneously alloying K and Zn on the the Ba and Mn sites, respectively Zhao et al. (2013). Magnetization measurements Johnston et al. (2011) and neutron diffraction Singh et al. (2009b) have determined the Néel temperature () of BaMnAs is near 620 K, significantly higher than for LnMnAsO. The magnetic anomaly at in BaMnAs is subtle, and the susceptibility does not obey a Curie Weiss law above this temperature. These observations suggest magnetic correlations exist above the long-range ordering temperature Johnston et al. (2011).

We recently synthesized LaMnAsO in our laboratory Idr () and found no magnetic anomaly near the reported Néel temperature, indicating that the presence of ferromagnetic MnAs ( = 318 K Guillaud (1951); SUZUKI and IDO (1982); Wada and Tanabe (2001)) as an secondary phase in some samples is likely responsible for the large magnetic anomaly usually observed near this temperature Sun et al. (2012). We performed high temperature magnetization measurements and observed a subtle anomaly near 650 K, reminiscent of the feature observed at = 618 K in Ref. 38 for BaMnAs. However, as noted above, a previous neutron diffraction study showed that of LaMnAsO is below 400 K Emery et al. (2011).

To resolve the uncertainty regarding the true Néel temperature of LaMnAsO and understand the nature of the magnetic anomaly that we observe at high temperature, we complemented our magnetization study with temperature dependent powder neutron diffraction and inelastic scattering measurements. Here we report the findings of these investigations, which include long-range magnetic order below = 360(1) K and short-range magnetic order above persisting up to about 650 K. We also observe intense quasi-elastic scattering in the short range ordered state, and strongly dispersing magnetic excitations in both the short- and long-range ordered states. A spin gap of about 3.5 meV is measured at room temperature. These findings provide a better understanding of the magnetic nature of LaMnAsO that appears to extend to other compounds with similar [MnAs] layers.

Ii Experimental Details

LaMnAsO powder samples were made by first reacting freshly filed La powder with As pieces to form LaAs. This was combined in a 1:1 molar ration with MnO and pressed into pellets. The pellets were heated at 1000 C for 48 h, ground and re-pelletized, and heated again at 1000 C for 48 h. All reactions were done inside evacuated and flame sealed silica-glass ampoules.

Powder x-ray diffraction (PXRD) data were collected using a PANalytical X’Pert Pro MPD diffractometer with monochromatic Cu K radiation. An Oxford PheniX cryostat was used for PXRD below room temperature. Rietveld refinement of room temperature PXRD data indicated the sample used in this study was 94% pure with 3% each of MnO and LaO by weight. Energy dispersive spectroscopic analysis carried out using a Bruker Quantax 70 x-ray spectrometer in a Hitachi TM3000 tabletop scanning electron microscope gave a composition of LaMnAsO determined by averaging multiple measurements on different crystallites in the sample, consistent with the expected 1111 stoichiometry. Magnetization measurement were performed with a Quantum Design Magnetic Property Measurement System, with the sample sealed inside a silica tube with approximately 1/3 of an atmosphere of ultra-high purity argon gas.

Neutron diffraction data was collected from 5 g of LaMnAsO powder at the High Flux Isotope Reactor at Oak Ridge National Laboratory (ORNL) using the Neutron Powder Diffractometer (beamline HB-2A). The sample was held in a cylindrical vanadium container, loaded in a closed cycle refrigerator capable of reaching temperatures from 6 K to 800 K. Measurements were performed using the wavelengths = 1.541 and 2.414 provided by a vertically focusing Ge monochromator. For data collection, a detector array consisting of 44 He tubes was scanned to cover the total 2 range of 7150 in steps of 0.05. Inelastic neutron scattering (INS) data was collected using the Hybrid Spectrometer (HYSPEC) at the Spallation Neutron Source at ORNL. HYSPEC is a direct geometry spectrometer that combines time-of-flight spectroscopy with focusing Bragg optics. The incident neutron beam is monochromated using a Fermi chopper and is then vertically focused by Bragg scattering from a monochromator onto the sample position. HYSPEC employs He linear position sensitive tube detectors that are assembled into 20 sets of 8-packs that cover an angular range of 60 in the horizontal scattering plane and a vertical acceptance of 15 Winn et al. (2015). Data on LaMnAsO was collected using polarized neutrons with incident energies of = 15 meV selected by a Heusler monochromator. A Mezei flipper was used to reverse the polarization direction of the incident beam, and the polarization analysis was done using a wide-angle supermirror bender. The polarization vector was oriented within the horizontal scattering plane by using a guide field generated from a set of orthogonal XYZ coils situated around the sample position. For the INS study the powder was loaded in aluminium cans and placed in a vacuum furnace.

Rietveld refinement of x-ray and neutron diffraction data was performed using Fullprof Ful (). Diffuse magnetic scattering was modeled using the reverse Monte Carlo algorithm employed by the program SPINVERT Paddison et al. (2013).

Iii Results and Discussion

iii.1 Crystal Structure and Thermal Expansion

Crystallographic details determined at 20, 300, and 730 K from x-ray and neutron diffraction are collected in Table 1. The unit cell parameters and atomic positions agree well with previous reports Nientiedt et al. (1997); Emery et al. (2010). Upon cooling from 300 to 20 K, a and c both contract by 0.22%, indicating remarkably isotropic thermal expansion in this tetragonal compound. Figure 1 confirms this, showing the lattice constants, normalized to 300 K values, determined from 20 to 730 K using x-ray and neutron diffraction.

Calculation of the linear coefficient of thermal expansion using the data between 200 and 300 K gives 1.05 K along a and 1.03 K along c. Nearly isotropic thermal expansion near room temperature is also seen in NdMnAsO Emery et al. (2011), where c changes with temperature about 1.2 times as fast as a, and PrMnAsOF Wildman et al. (2015), where this ratio is 1.3. This can be compared to isostructural LaFeAsO McGuire et al. (2008) and NdCoAsO McGuire et al. (2010), in which the ratios of linear thermal expansion along the different crystallographic directions c and a are 2.0 and 5.4, respectively, indicting that this anisotropy likely arises from the transition metal 3d-orbitals.

Table 2 lists the MnMn interatomic distances determined near room temperature in LaMnAsO along with several related 1111 and 122 compounds that also contain layers of Mn in edge sharing tetrahedral coordination. Also shown is the Mn bond angle formed by to two neighboring pnictogen (Pn) atoms located above a Mn layer. The angle , formed by two Pn on opposite sides of the Mn layer, is related to by cos() = -[1+cos()]. For ideal tetrahedral geometry, both angles are 109.5 degrees. Long range antiferromagnetic ordering temperatures (), magnetic moments on Mn determined at low temperature (), and the magnetic structure defined by how the layers stack along the direction are also listed. Correlations between crystal structure and magnetic behavior observed in these data will be discussed below.

(K) 20 300 300 730
radiation x-ray x-ray neutron neutron
(Å) 1.540598 1.540598 1.541 1.541
a (Å) 4.1109(1) 4.1200(1) 4.1235(2) 4.1473(1)
c (Å) 9.0259(3) 9.0462(2) 9.0532(5) 9.1040(3)
0.6328(3) 0.6323(3) 0.6325(4) 0.6319(3)
0.1680(5) 0.1684(4) 0.1685(6) 0.1687(3)
) 0.45(8) 0.79(6)
) 0.58(14) 1.34(12)
) 0.65(10) 1.29(9)
) 0.62(10) 1.01(9)
1.89 1.65 1.86 1.08
Table 1: Crystallographic parameters for LaMnAsO determined by Rietveld refinement of powder x-ray and neutron diffraction data using space group with La at (, , z), Mn at (, , 0), As at (, , z), and O at (, , ).
Figure 1: Temperature dependence of the lattice parameters of LaMnAsO. Neutron and x-ray diffraction results are shown for and , respectively. Each data set was separately normalized to 300 K values. See also Table 1.
compound (Mn) stacking Refs.
Å Å deg. K
LaMnAsO 2.913 9.046 107.0 360(1) 3.34(2) FM this work, 34
CeMnAsO 2.889 8.956 106.7 347(1) 3.32(4) FM 47; 33
PrMnAsOF 2.877 8.921 105.6 340 3.7 FM 30
NdMnAsO 2.864 8.905 105.0 359 3.54(4) FM 34
NdMnAsOF 2.863 8.897 105.2 356(2) 3.83(2) FM 18
LaMnPO 2.869 8.843 111.2 375 3.28(5) FM 49
PrMnSbO 2.961 9.472 230 3.69(3) FM 50; 7
BaMnAsF 3.022 9.588 110.3 338(1) 3.65(5) AFM 51
BaMnSbF 3.167 9.830 107.2 272(1) 3.66(3) AFM 51
BaMnAs 2.948 6.734 108.7 625 3.88(4) AFM 37
BaMnP 2.854 6.526 111.8 4.2(1) AFM 52
Table 2: Comparison of MnMn distances near room temperature and the angle () formed by two coplanar pnictogen (Pn) atoms and their shared Mn atom, long range magnetic ordering temperatures, ordered magnetic moments at low temperature, and stacking pattern of checkerboard AFM layers in some manganese pnictides containing layers of Mn in edge-sharing tetrahedral coordination.

iii.2 Magnetization

Results from magnetization measurements are shown in Figure 2. Isothermal magnetization curves at 5, 200, and 380 K (Fig. 2a) show no indication of significant ferromagnetic components. The H = 0 intercept extrapolated from the high field portion of the 200 K M vs H data shows that the sample contains less than 0.001% MnAs by weight Ido (1985). The cusp in magnetic susceptibility () near 120 K (Fig. 2b) arises from the small MnO impurity observed by diffraction, which orders antiferromagnetically at this temperature Tyler (1933). The low temperature upturn upon cooling is expected to arise from a small concentrations of paramagnetic moments associated with defects or secondary phases. The Curie constant determined by fitting the low temperature tail corresponds to 0.03 spin 1/2 moments per formula unit.

No significant anomaly is observed in (Fig. 2b) between 300 and 400 K, the temperature range in which is expected to lie Emery et al. (2011). A change in the temperature dependence of is observed at 650 K. This similar to the anomaly observed at 618 K in the closely related material BaMnAs Johnston et al. (2011), which was shown to coincide with the onset of long-range antiferromagnetic ordering of Mn moments Singh et al. (2009b). However, neutron diffraction results presented below will show that in LaMnAsO this high temperature feature arises from the development of short-range magnetic order, and long-range order occurs only upon cooling below 360 K, with no associated anomaly in the magnetic susceptibility.

Figure 2: Results of magnetization measurements on LaMnAsO. (a) Isothermal magnetization curves showing the absence of any significant ferromagnetic signal. (b) Temperature dependence of the magnetic susceptibility, with the anomaly at 650 K indicated on the figure and the inset. Note the small cusp associated with the MnO impurity, and the absence of any magnetic anomaly near room temperature.

iii.3 Neutron scattering: magnetic order

Neutron powder diffraction patterns and Rietveld refinement results at 300 and 730 K are shown in Figure 3a and 3b, respectively. At 730 K the data are fully described by the nuclear structures with no indication of magnetic scattering. Additional sharp Bragg reflections indicating long-range magnetic order (LRO) are apparent at 300 K, and were modeled well using the magnetic structure first reported by Emery et al. in Refs. 16 and 34. The magnetic structure is shown in the inset of Figure 3a. Moments on Mn form a checkerboard pattern in the ab-plane and these planes stack ferromagnetically along the c-direction. Such a spin arrangement is described by the magnetic space group (with Mn1 at (3/4,1/4,1/2 0,0,) and Mn2 at (1/4,3/4,1/2 0,0,)). This same magnetic structure is adopted at room temperature by the other LnMnAsO compounds with Ln = Ce, Pr, and Nd Marcinkova et al. (2010a); Zhang et al. (2015); Wildman et al. (2015). The fit shown in Figure 3 indicates a magnetic moment of 2.27(3)  per Mn at 300 K. As expected, this is slightly smaller than the value of 2.43(1)  per Mn reported at 290 K Emery et al. (2010).

Structural and magnetic properties of 1111 and 122 Mn pnictide compounds are collected in Table 2. In all of the compounds the ordered moment at low temperature saturates to values significantly smaller than the full moment of , with values in the range of 3.24.2 . This moment reduction has been attributed to quantum fluctuations Johnston et al. (2011) or MnAs hybridization An et al. (2009). The largest magnetic moments are seen in 122 compounds. The long range ordering temperature is highest for BaMnAs, which could be attributed to the shortest interlayer MnMn distance; this is expected to be the direction of weakest magnetic exchange, and therefore should determine the onset of LRO. Within the 1111 type compounds there is no clear overall trend in with MnMn distance, although the lowest values, for PrMnSbO and BaMnSbF, occur with two of the largest intralayer MnMn separation. The most interesting correlation seen in the table is between this bond angle and the stacking sequence of the checkerboard-ordered Mn layers along the c axis. With one exception (LaMnPO), all compounds with have layers stacked FM resulting in a C-type AFM structure, while those with stack AFM giving a G-type AFM structure. This shows that the magnetic exchange along the c-axis is determined by the details of the crystal field splitting of the Mn 3d orbitals. This suggests that increasing the angle with applied or chemical pressure may switch the magnetic order from C-type to G-type in some 1111 Mn compounds.

Data collected at intermediate temperatures was used to determine the lattice parameters shown in Figure 1, and analyzed for indications of magnetic scattering. As expected based on Ref. 34, the sharp magnetic Bragg peaks were absent at 400 K; however additional diffuse scattering centered near the magnetic 100 and 101 Bragg peaks was observed. This is shown Figure 3c, which compares diffraction patterns collected at 300, 400, and 730 K. The diffuse scattering is not detectable at 730 K, and nearly completely suppressed in the long-range ordered state at 300 K. This is identified as a signature of short-range magnetic order (SRO) above the Néel temperature. An estimate of the correlation length () of the SRO was determined by fitting the Q dependence of the intensity to a Lorentzian lineshape ], where . The fit is shown in Figure 3d, which gives = 8.0(3) Å. This SRO is examined in more detail below.

Order parameters were determined as a function of temperature for the LRO and SRO by counting with the detector at the center of the sharp and diffuse peaks in Figure 3d, respectively. The results are shown in Figure 4. Upon heating from room temperature, the LRO is observed to vanish near 360 K. Fits of the order parameter (defined as the square root of the measured intensity) above 330 K to the functional form give  = 361.4(5) K and  = 0.35(2) using the 101 reflection and  = 360.5(3) K and  = 0.33(1) using the 100 reflection using either the 100 or 101 reflection. Thus, an Néel temperature of 360(1) K is determined. The values determined for the critical exponent are close to those expected for a three dimensional Heisenberg systems. For BaMnAs a similar value of 0.35(2) has been reported Singh et al. (2009b). Among Mn-based 1111 materials, reports include a significantly higher  = 0.47 for CeMnAsO Zhang et al. (2015) and a significantly lower  = 0.27(1) for NdMnAsO Marcinkova et al. (2010a). The fitting range used in Ref. Marcinkova et al. (2010a) for NdMnAsO extends to T = 0, which may account for that discrepancy. For CeMnAsO, the authors of Ref. Zhang et al. (2015) note that large values of can arise from strongly anisotropic exchange interactions, though it is not clear why CeMnAsO would be more anisotropic than LaMnAsO. Values similar to those found here are reported for BaMnPF,  = 0.31(1), and BaMnAsF,  = 0.36(3) Saparov et al. (2013).

The SRO decreases upon cooling below this , and is nearly absent at 300 K. Upon heating above the SRO decreases gradually, and vanishes near the temperature at which the magnetic susceptibility anomaly is observed (Fig. 2b). This is used to identify the onset of magnetic correlations at  = 650(10) K.

Figure 3: Neutron powder diffraction data from LaMnAsO collected with a neutron wavelength of 2.414 Å. (a) Rietveld refinement of data collected at 300 K. (b) Rietveld refinement of data collected at 730 K. Upper set of ticks locate nuclear reflections, lowest set of ticks at 300 K locate magnetic reflections, other sets indicate reflections from MnO and LaO impurities. (c) Comparison of diffraction patterns showing sharp peaks at 300 K indicating long-range magnetic order, a broad peak at 400 K indicating short-range order, and no magnetic scattering at 730 K. (d) A fit to the short-range order peak using a Lorentzian function. The two sharp reflections that appear in the Q range of interest are modeled in the same way.
Figure 4: Temperature dependence of magnetic scattering from LaMnAsO. (a) Counts measured at the center of the 101, 100, and short-range order (SRO) peaks. (b) The behavior of the data shown in (a) near the Néel temperature.

To allow further analysis of the SRO in LaMnAsO, the diffuse scattering intensity was extracted from the neutron powder diffraction pattern collected at 400 K. This was accomplished by using the background determined at 300 K together with a fit of the nuclear scattering profile at 400 K. The resulting difference curve reflects the broad diffuse scattering present at 400 K, as well as sharp features corresponding to errors in the fitted intensities of the Bragg peaks. The sharp features were excluded and the resulting data is shown in Figure 5a.

The diffuse scattering data was analyzed using SPINVERT Paddison et al. (2013). This program uses a reverse Monte Carlo algorithm to fit non-periodic spin configurations to measured neutron powder diffraction data. A supercell of the known nuclear structure is generated and a randomly oriented magnetic moment is assigned to each site occupied by a magnetic ion, in this case Mn. A 101010 supercell was used here. The scattering from the spin configuration is calculated and compared with the experimental data. The orientation of a spin is then changed and a determination is made of whether the change improved the fit or not. If it is improved, the new configuration is accepted and the process repeats until a convergence criterion is met or a maximum number of “moves” are performed. This is typically done multiple times with different starting configurations to ensure the robustness of the results. Here five runs were performed with a maximum of 300 moves per spin in each run. For more details see Ref. 46 and references therein.

Figure 5: (a) Measured diffuse magnetic scattering (circles) and calculated scattering (lines) from the spin configurations determined by reverse Monte Carlo using the program SPINVERT Paddison et al. (2013). (b) The radial dependence of the spin-pair correlation function corresponding to the fits shown in (a). (c) Data from (b) grouped into sets corresponding spins separated along to the indicated crystallographic directions. (d) Absolute value of the data in (b) for pairs of spins within a common layer (closed symbols) and in neighboring layers (open symbols). The vertical error bars on the spin-pair correlation values in panels b,c,d are approximately the size of the data markers.

The diffuse scattering intensity calculated from the resulting spin configurations is shown in Figure 5a as a solid line for an isotropic model and as a dashed line for an anisotropic model with spins directed along the c-axis. The results are nearly identical for the two models, and reproduce the measured intensity well. The spin-pair correlation function is shown for both models in Figure 5b, where each data point corresponds to a distance between two Mn sites within the LaMnAsO crystal structure. Larger absolute values of this function indicate stronger preference for collinear arrangement of the spins separated by the distance r, and the sign indicates whether the spins are parallel (+) or antiparallel (-). Again, there is little difference observed between the results of the isotropic and anisotropic spins models, and further analysis will focus only on the isotropic results.

The spin-pair correlations shown in Figure 5b are strongly diminished as r increases up to 20 Å, and include both positive and negative values. This is indicative of short-range antiferromagnetic correlations. Based on the interatomic distance, each data point can be identified with a specific pair of Mn atoms within the LaMnAsO structure. For reference, a single Mn square is shown in the inset of Figure 5c. The shortest MnMn distance is along the edge of the Mn square lattice lying in the ab-plane. Thus, since the Mn square lattice is rotated by 45 degrees with respect to the square face of the tetragonal crystal lattice, the nearest MnMn distance is along the [110] direction and equal to  = 2.92 Å, and the next nearest MnMn distance, which is the diagonal of the Mn square lattice, is in the [100] direction and equal to a = 4.12 Å. By grouping together points in Figure 5b that are at multiples of these distances, spin-pair correlations along [100] and [110] can be separately isolated. The results are shown in Figure 5c. It can be seen from this Figure that ferromagnetic alignment of the spins is preferred along the [100] direction, and antiferromagnetic alignment of the spins is preferred along the [110] direction. Distances along the c axis, or [001], are also shown. The correlations for these pairs of spins are positive and significantly smaller than the in-plane correlations. The signs of magnetic interactions along these three crystallographic directions are the same as those observed in the LRO state below (Fig. 3a). The reduced magnitude of along [001] relative to [100] and [110] shows that the magnetic interactions are relatively weak between the Mn layers.

The correlation length associated with the SRO can be defined by how rapidly the decrease with r. For this purpose, the absolute value of the spin-pair correlation function is examined. Figure 5d displays the data for pairs that lie within a Mn layer (intraplanar, solid symbols) and pairs which lie within two neighboring Mn layers, separated by one unit cell along the c axis (interplanar, open symbols). These sets form two well-separated branches on the plot. Each is fitted with an exponential decay toward zero . The resulting correlation lengths are shown on the figure. As expected, the correlation length within the layers is larger than that between the layers. The values are consistent with that approximated by the width of the strongest diffuse scattering feature in Figure 3d. In the ab plane, corresponds to approximately two unit cells. For pairs of spins in neighboring layers, the magnitude of is likely less reliable, since it is smaller than the smallest distance used to determine it. However, it does indicate weaker magnetic correlations along the c axis, as may be expected. In the related compound BaMnAs, which is structurally less anisotropic than LaMnAsO, exchange interactions along the layer stacking direction have been estimated to be about 10 times weaker than in-plane interactions Johnston et al. (2011).

This analysis of the diffuse magnetic scattering reveals that the magnetic order in LaMnAsO is short-range and nearly two dimensional between 650 and 360 K. It is interesting to consider whether similar magnetic correlations exist above the long-range ordering temperature in related materials, including other Mn based 1111 and 122 compounds. Several examples can be found in the literature. Wildman et al. noted some diffuse intensity attributed to SRO at 350 K in PrMnAsOF ( = 340 K) Wildman et al. (2015). Similarly, in CeMnAsO ( = 347 K) a broad magnetic peak was observed up to about 420 K Zhang et al. (2015). In addition, neutron diffraction data collected at 600 K from BaKMnAs ( = 575 K) Lamsal et al. (2013) and at 650 K from BaMnAs ( = 625 K) Singh et al. (2009b) show the presence of diffuse scattering. Thus, it appears that SRO above is a common feature of the magnetism in Mn-based 1111 and 122 type compounds, and may have been overlooked in previous studies of the other materials listed in Table 2 that focused primarily on the long range magnetic order.

iii.4 Neutron scattering: magnetic excitations

Figure 6: (a-c) Energy versus momentum transfer (Q) slices of the polarized inelastic neutron scattering measured in spin-flip configuration at 300 K, 370 K and 385 K. Highly dispersive spin-wave excitation are located at momentum transfer Q=1.6Å near the 100 and 101 magnetic peak positions. The excitation spectrum at 300 K displays an anisotropy gap that vanishes with the long-range order. (d) Constant-Q cuts near Q = 1.55 Å revealing a spin-gap of approximately 3.5 meV in the LRO state, at 300 K, and the gapless nature of the excitation spectrum in the SRO state at 370 K.

Polarized inelastic neutron scattering experiments were performed to examine the magnetic excitation spectrum of LaMnAsO in both the LRO and SRO states. Results of measurements using the spin-flip configuration are shown in Figure 6. Measurements at 300 K, below , reveal the presence of well-defined spin-wave excitations located near a momentum transfer  = 1.6 . The feature was not observed in the non-spin-flip configuration, demonstrating that it is magnetic in nature. The momentum transfer at which the strong inelastic scattering is centered is close to that where the strongest magnetic Bragg peaks and the maximum in the diffuse magnetic scattering are located (Fig. 5d). At 300K, a gap is observed in the magnetic excitation spectrum, which arises from magnetic anisotropy in the LRO state. The gap is estimated to be 3.5 meV from the present data. As noted above, the spectrum is very narrow in Q, indicating the spin waves producing this feature are very sharply dispersing. This scattering is quite similar to the polycrstalline averaged spectrum produced from spin-wave calculations on BaMnAs Johnston et al. (2011). Non-polarized measurements using higher incident energy neutrons showed that the magnetic spectrum in LaMnAsO extends to beyond 30 meV, and evidence for magnetic excitations extending up to 100 meV is seen in BaMnAs Johnston et al. (2011). The stiff spin wave spectra in these materials indicate strong magnetic exchange interactions and are consistent with the observed high magnetic ordering temperatures.

As the temperature is raised above  = 360 K, the inelastic magnetic scattering (Fig. 5b,c) remains narrow in Q; however, the spin-gap vanishes with the long-range order. Just above , at 370 K, very intense quasielastic scattering associated with the short range magnetic order is observed at 1.6 . This quasielastic scattering is suppressed strongly as temperature is increased further, as seen by comparison of the data at 370 K (Fig. 5b) and 385 K (Fig. 5c). Fixed-Q cuts through the spectra are collected in Figure 5d, which shows the spin-gap at 300 K, strong, low-energy scattering at 370 K, and weaker magnetic scattering at 385 K that still extends up beyond 12 meV. The scaling of the energy gap with the ordered moment is expected as the the thermal energy overcomes the magnetic anisotropy energy, leading to a gradual decay of the spin correlations and thus to a short range order. However, the strong exchange interaction makes the magnetic excitations persist at temperatures well above the LRO regime.

Iv Summary and Conclusions

The long range antiferromagnetic order present at room temperature in LaMnAsO vanishes above  = 360 K. No corresponding anomaly is detected in the magnetic susceptibility, because short range magnetic order is present above . This short range order extends up to  = 650 K. A change is observed in the temperature dependence of the magnetic susceptibility at . Analysis of the diffuse magnetic scattering shows that the magnetic correlations associated with the short range order are anisotropic, with longer correlation length ( 10Å) in the ab-plane than along the c-axis ( 5Å). The magnetic exchange is similar to that in the long range ordered state: antiferromagnetic along the edge of the Mn square lattice and ferromagnetic between layers. Comparison of the long range ordered magnetic structures of related compounds shows that the magnetic exchange along the stacking direction is determined by the Mn d-orbital crystal field splitting, with a strong correlation between Mn bond angle and the layer stacking sequence. This suggests chemical substitutions or applied pressure may be used to switch the interlayer magnetic coupling. A gap of about 3.5 eV in the magnetic excitation spectrum of LaMnAsO is observed at room temperature (below ), arising from magnetic anisotropy. Steeply dispersing spin waves extend beyond 30 meV and indicate strong magnetic exchange interactions consistent with high magnetic ordering temperatures. Very strong quasielastic scattering associated with the short range magnetic order is seen just above . Comparison with neutron powder diffraction data from other LnMnAsO and BaMnAs based materials shows that short ranged order above may be a common fundamental behavior of these quasi-two dimensional magnetic materials.

V Acknowledgements

This work was supported by the U. S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division (M.A.M.) and Scientific User Facilities Division (V.O.G). The authors are grateful for helpful discussions with S. Calder, B.C. Sales, A.F. May, and J.-Q. Yan.


  1. R. Pöttgen and D. Johrendt, Z. Naturforsch. 63b, 1135 (2008).
  2. P. Quebe, L. J. Terbüchte,  and W. Jeitschko, J. Alloys Compd. 302, 70 (2000).
  3. Y. Kamihara, T. Watanabe, M. Hirano,  and H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008).
  4. G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, P. Zheng, J. L. Luo,  and N. Wang, Phys. Rev. Lett. 100, 247002 (2008).
  5. Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, Z.-C. Li, G.-C. Che, X.-L. Dong, L.-L. Sun, F. Zhou,  and Z.-X. Zhao, Chin. Phys. Lett. 25, 2215 (2008).
  6. C. Wang, L. Li, Z. Zhu, Z. Ren, Y. Li, Y. Wang, X. Lin, Y. Luo, S. Jiang, X. Xu, G. Cao,  and Z. Xu, EPL 83, 67006 (2008).
  7. I. Schellenberg, T. Nilges,  and R. Pöttgen, Z. Naturforsch. 63b, 834 (2008).
  8. Y. Kamihara, H. Hiramatsu, M. Hirano, R. Kawamura, H. Yanagi, T. Kamiya,  and H. Hosono, J. Am. Chem. Soc. 128, 10012 (2006).
  9. J. D. Fletcher, A. Serafin, L. Malone, J. G. Analytis, J.-H. Chu, A. S. Erickson, I. R. Fischer,  and A. Carrington, Phys. Rev. Lett. 102, 147001 (2009).
  10. M. Tegel, S. Johansson, V. Weiss, I. Schnellenberg, W. Hermes, R. Pöttgen,  and D. Johrendt, EPL 84, 67007 (2008).
  11. S. Matsuishi, Y. Inoue, T. Nomura, M. Hirano,  and H. Hosono, J. Phys. Soc. Jpn. 77, 113709 (2008a).
  12. S. Matsuishi, Y. Inoue, T. Nomura, H. Yanagi, M. Hirano,  and H. Hosono, J. Am. Chem. Soc. 130, 14428 (2008b).
  13. F. Han, X. Y. Zhu, G. Mu, P. Cheng,  and H. H. Wen, Phys. Rev. B 78, 180503 (2008).
  14. A. S. Sefat, A. Huq, M. A. McGuire, R. Y. Jin, B. C. Sales, D. Mandrus, L. M. D. Cranswick,  and P. W. S. adn K. H. Stone, Phys. Rev. B 78, 104505 (2008).
  15. A. Marcinkova, T. C. Hansen, C. Curfs, S. Margadonna,  and J. W. G. Bos, Phys. Rev. B 82, 174438 (2010a).
  16. N. Emery, E. J. Wildman, J. M. S. Skakle, G. Giriat, R. I. Smith,  and A. C. Mclaughlin, Chem. Commun. 46, 6777 (2010).
  17. Y. Tsukamoto, Y. Okamoto, K. Matsuhira, M.-H. Whangbo,  and Z. Hiroi, J. Phys. Soc. Jpn. 80, 094708 (2011).
  18. E. J. Wildman, J. M. S. Skakle, N. Emery,  and A. C. Mclaughlin, J. Am. Chem. Soc. 134, 8766 (2012).
  19. H. Yanagi, R. Kawamura, T. Kamiya, Y. Kamihara, M. Hirano, T. Nakamura, H. Osawa,  and H. Hosono, Phys. Rev. B 77, 224431 (2008).
  20. H. Ohta and K. Yoshimura, Phys. Rev. B 80, 184409 (2009).
  21. M. A. McGuire, D. J. Gout, V. O. Garlea, A. S. Sefat, B. C. Sales,  and D. Mandrus, Phys. Rev. B 81, 104405 (2010).
  22. A. Marcinkova, D. A. M. Grist, I. Margiolaki, T. C. Hansen, S. Margadonna,  and J.-W. G. Bos, Phys. Rev. B 81, 064511 (2010b).
  23. T. Watanabe, H. Yanagi, Y. Kamihara, T. Kamiya, M. Hirano,  and H. Hosono, J. Solid State Chem. 181, 2117 (2008).
  24. Z. Li, G. Chen, J. Dong, G. Li, W. Hu, D. Wu, S. Su, P. Zheng, T. Xiang, N. Wang,  and J. Luo, Phys. Rev. B 78, 060504(R) (2008).
  25. S. Matsuishi, A. Nakamura, Y. Muraba,  and H. Hosono, Supercond. Sci. Technol. 25, 084017 (2012).
  26. Y. Luo, H. Han, H. Tan, X. Lin, Y. Li, S. Jiang, C. Feng, J. Dai, G. Cao, Z. Xu,  and S. Li, J. Phys.: Condens. Matter 23, 175701 (2011).
  27. M. A. McGuire, V. O. Garlea, A. F. May,  and B. C. Sales, Phys. Rev. B 90, 014425 (2014).
  28. Y. L. Sun, J. K. Bao, Y. K. Luo, C. M. Feng, Z. A. Xu,  and G. H. Cao, EPL 98, 17009 (2012).
  29. T. Hanna, S. Matusishi, K. Kodama, T. Otomo, S. Shamoto,  and H. Hosono, Phys. Rev. B 87, 020401 (2013).
  30. E. J. Wildman, F. Sher,  and A. C. Mclaughlin, Inorg. Chem. 54, 2536 (2015).
  31. K. Kayanuma, H. Hiramatsu, T. Kamiya, M. Hirano,  and H. Hosono, J. Appl. Phys. 105, 3093685 (2009).
  32. X. X. Li, X. J. Wua,  and J. L. Yang, J. Mater. Chem. C 1, 7197 (2013).
  33. Q. Zhang, W. Tian, S. G. Peterson, K. W. Dennis,  and D. Vaknin, Phys. Rev. B 91, 064418 (2015).
  34. N. Emery, E. J. Wildman, J. M. S. Skakle, A. C. Mclaughlin, R. I. Smith,  and A. N. Fitch, Phys. Rev. B 83, 144429 (2011).
  35. J. M. An, A. S. Sefat, D. J. Singh,  and M. H. Du, Phys. Rev. B 79, 075120 (2009).
  36. Y. Singh, A. Ellern,  and D. C. Johnston, Phys. Rev. B 79, 79.094519 (2009a).
  37. Y. Singh, M. A. Green, Q. Huang, A. Kreyssig, R. J. McQueeney, D. C. Johnston,  and A. Goldman, Phys. Rev. B 80, 100403 (2009b).
  38. D. C. Johnston, R. J. McQueeney, B. Lake, A. Honecker, M. E. Zhitomirsky, R. Nath, Y. Furukawa, V. P. Antropov,  and Y. Singh, Phys. Rev. B 84, 094445 (2011).
  39. K. Zhao, Z. Deng, X. C. Wang, W. Han, J. L. Zhu, X. Li, Q. Q. Liu, R. C. Yu, T. Goko, B. Frandsen, L. Liu, F. L. Ning, Y. J. Uemura, H. Dabkowska, G. M. Luke, H. Luetkens, E. Morenzoni, S. R. Dunsiger, A. Senyshyn, P. Boni,  and C. Q. Jin, Nat. Commun. 4 (2013), 10.1038/ncomms2447.
  40. J. C. Idrobo et. al, unpublished (2015).
  41. C. Guillaud, J. Phys. Radium 12 (1951).
  42. T. SUZUKI and H. IDO, J. Phys. Soc. Jpn. 51, 3149 (1982).
  43. H. Wada and Y. Tanabe, Appl. Phys. Lett. 79, 3302 (2001).
  44. B. Winn, U. Filges, V. O. Garlea, M. Graves-Brook, M. Hagen, C. Jiang, M. Kenzelmann, L. Passell, S. M. Shapiro, X. Tong,  and I. Zaliznyak, EJP Web Conferences 83, 03017 (2015).
  45. J. Rodriguez-Carvajal, FullProf.2k, version 5.30, March 2012, ILL.
  46. J. A. M. Paddison, J. R. Stewart,  and A. L. Goodwin, J. Phys.: Condens. Matter 25, 454220 (2013).
  47. A. T. Nientiedt, W. Jeitschko, P. G. Pollmeier,  and M. Brylak, Zeitschrift Fur Naturforschung Section B-a Journal of Chemical Sciences 52, 560 (1997).
  48. M. A. McGuire, A. D. Christianson, A. S. Sefat, B. C. Sales, M. D. Lumsden, R. Jin, E. A. Payzant, D. Mandrus, Y. Luan, V. Keppens, V. Varadarajan, J. W. Brill, R. P. Hermann, M. T. Sougrati, F. Grandjean,  and G. J. Long, Phys. Rev. B 78, 094517 (2008).
  49. H. Yanagi, T. Watanabe, K. K. andS Iikubo, S. Shamoto,  and T. Kamiya, J. Appl. Phys. 105, 093916 (2009).
  50. S. A. J. Kimber, A. H. Hill, Y.-Z. Zhang, H. O. Jeschke, R. Valentí, C. Ritter, I. Schellenberg, W. H. adn R. Pöttgen,  and D. Argyriou, Phys. Rev. B 82, 100412(R) (2010).
  51. B. Saparov, D. J. Singh, V. O. Garlea,  and A. S. Sefat, Sci. Rep. 3, 2154 (2013).
  52. S. L. Brock, J. E. Greeden,  and S. M. Kauzlarich, J. Solid State Chem. 113, 303 (1994).
  53. H. Ido, J. Appl. Phys. 57, 3247 (1985).
  54. R. W. Tyler, Phys. Rev. 44, 776 (1933).
  55. J. Lamsal, G. S. Tucker, T. W. Heitmann, A. Kreyssig, A. Jesche, A. Pandey, W. Tian, R. J. McQueeney, D. C. Johnston,  and A. I. Goldman, Phys. Rev. B 87, 144418 (2013).
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 minumum 40 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