The magnetic and electronic properties of Oxyselenides - influence of transition metal ions and lanthanides

C. Stock1 and E. E. McCabe2 1School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3FD, UK 2School of Physical Sciences, University of Kent, Canterbury, CT2 7NH, UK (Dated: October 15, 2018)

Magnetic oxyselenides have been a topic of research for several decades, firstly in the context of photoconductivity and thermoelectricity owing to their intrinsic semiconduct- ing properties and ability to tune the energy gap through metal ion substitution. More recently, interest in the oxyselenides has experienced a resurgence owing to the possible relation to strongly correlated phenomena given the fact that many oxyselenides share a similar structure to unconventional superconducting pnictides and chalcogenides. The two dimensional nature of many oxyselenide systems also draws an analogy to cuprate physics where a strong interplay between unconventional electronic phases and localised magnetism has been studied for several decades. It is therefore timely to review the physics of the oxyselenides in the context of the broader field of strongly correlated magnetism and electronic phenomena. Here we review the current status and progress in this area of research with the focus on the influence of lanthanides and transition metal ions on the intertwined magnetic and electronic properties of oxyselenides. The emphasis of the review is on the magnetic properties and comparisons are made with iron based pnictide and chalcogenide systems.

CONTENTS 8. Summary of the magnetic structure variation with transition metal ion 21 I. Introduction 1 E. Magnetic interactions 22 1. La2O2Fe2OSe2 and Ce2O2FeSe2 22 2+ 2+ II. Experimental techniques 2 2. Ln2O2M 2OSe2 for M=Mn or Co 25 A. Synthesis routes 2 B. Magnetic neutron scattering 3 V. Electronic properties 25 C. X-ray and optical spectroscopy 5 A. ZrCuSiAs structures and related phases 26 D. Transport and thermodynamic measurements 5 B. Ln2O2M 2OSe2 structures and related phases 27 C. Superconductivity and the Oxyselenides 29 III. Oxyselenide structural families 6 A. Ln – O – Se phases 6 VI. Conclusion 30 B. Ln –O– M – Se phases 7 References 30 1. ZrCuSiAs structures and related cation-ordered phases 7 2. β-La2O2MSe2 and related polymorphs 8 3. ZrCuSiAs-modified structures 8 I. INTRODUCTION 3+ 4. LnMOSe2 - M ions coordinated by selenide 8 5. Ln O M OSe materials 9 2 2 2 2 Materials based upon transition metal ions have con- IV. Magnetic Properties 10 sistently been a source of interest owing to novel elec- A. Local magnetism 10 arXiv:1704.01441v1 [cond-mat.str-el] 5 Apr 2017 tronic, magnetic, and structural properties they possess. 1. Rare Earth local magnetism 10 While the description of structural and insulating mag- 2. Transition metal ion local magnetism 11 netic transitions has been well understood with a robust B. Magnetic structures 13 C. Ln – O – Se phases 13 formalism to describe such systems, metal-insulator tran- D. Ln –O– M – Se phases 13 sitions and the new phases which exist near these critical 1. LnCrOS2 (oxysulfides) 13 points are not described by any such theory. The ex- 2. ZrCuSiAs structures 14 ploration of materials that host such transitions has led 2+ 3. Ln2O2M 2OSe2 phases with M=Fe 15 2+ to several notable discoveries including high temperature 4. Ln2O2M 2OSe2 phases with M=Co 18 2+ superconductivity in the cuprates in 1986 (Bednorz and 5. Ln2O2M 2OSe2 phases with M=Mn 19 6. Mixed Fe/Mn-based oxyselenides 20 Muller, 1986) and more recently with the discovery of 7. Critical scattering 20 superconductivity in LaFeAsO1−xFx in 2008 (Kamihara 2 et al., 2008). This report of iron-based superconductivity condensed matter physics. Materials near such transition lead to the discovery of many more iron-based supercon- provide the possibility for the discovery of new phases of ducting materials which have challenged both theories matter as illustrated by the discovery of high temperature and experiments, including iron arsenides, iron selenides superconductivity in the cuprates and more recently in and mixed-anion iron oxyarsenides. In this review, we iron based pnictides and chalcogenide systems (Johnston, discuss a class of related mixed-anion materials: oxyse- 2010; Paglione and Greene, 2010; Stewart, 2011; Wen lenides. Whilst this family of materials is generally not et al., 2011). Many oxyselenides contain iron on a two superconducting, their structural similarities with sev- dimensional lattice, structurally similar to iron based su- eral classes of the iron-based superconductors, and the perconducting systems and indeed there were suggestions ability to tune their semiconducting nature, makes them that oxyselenides may represent candidate superconduct- relevant to current research in strongly correlated elec- ing phases during early work on cuprate superconduc- tron systems. tors (Otzschi et al., 1999). This analogy has led to recent Selenium (named after the Greek, selene, for moon due investigations of magnetism and electronic phenomena in to its many similarities with tellurium which was named such systems. Given the semiconducting nature and the after the Latin, tellus, for Earth) (Emsley, 2001), can ability to tune the electronic band gap through transition adopt a wide range of oxidation states from 0 to +6 (in metal ion substitution, the oxyselenides provide a unique selenates) and +4 (in selenites) to -2 in selenides which series of compounds to search for new electronic phases are most commonly formed with the more electropositive and critical properties. group 1 and group 2 elements and lanthanides (Green- This review is divided into six sections including this wood and Earnshaw, 2005). In this review, we focus on introduction. We first present some of the experimen- transition metal oxyselenides and explore the role of the tal techniques used to investigate the magnetic and elec- lanthanide and transition metal ions in the crystal struc- tronic properties of the oxyselenides. We then outline ture and magnetic and electronic properties of these sys- the oxyselenide structural families discussed in this re- tems. view followed by a discussion of the magnetic structures The semiconducting properties of the oxyselenides and and properties. We then review the electronic proper- the ability to tune these properties has resulted in these ties of the oxyselenides and finish with a summary and materials being the topic of research for several decades. conclusion. Early interest in the oxyselenides arose in the area of photoconductivity which arguably led to the discovery of LaCuOS and related materials (Minami et al., 1974). II. EXPERIMENTAL TECHNIQUES More recently, work has focussed heavily on thermo- electrics (Snyder and Toberer, 2008) and, again, the un- In this section, we provide an outline of the various derlying tunable semiconducting properties resulted in experimental techniques used to synthesise and to probe many of these systems being identified as good thermo- the structural, magnetic and electronic properties of the electrics (Zhao et al., 2010). oxyselenides which are the focus of this review. Although this review focuses on a class of selenides, it is worth mentioning selenite materials such as Cu2OSeO3 which have received significant attention as hosts (Ruff A. Synthesis routes et al., 2015; Seki et al., 2012) for Skyrmions, vortex-like topological spin structures (Skyrme, 1962) originally dis- Unlike the synthesis of oxides which can often be car- covered in MnSi (Muhlbauer et al., 2009; Seki et al., 2012; ried out in air, preparing oxyselenides requires more Yu et al., 2010). The Cu2OSeO3 crystal structure is com- control to avoid oxidation of selenium or other selenide posed of edge- and corner-linked CuO5 polyhedra linked reagents to higher selenium oxidation states. Many of the 2− 4+ by SeO3 selenite groups, with the role of the Se lone synthetic routes reviewed by Sefat to prepare iron-based pair evident in its coordination (Meunier and Bertraud, superconductors (Sefat, 2013) are relevant to oxyselenide 1976). While these materials constitute an important synthesis and we give only a brief overview here. and rapidly evolving area of physics and materials sci- Polycrystalline samples of oxyselenides can be pre- ence, we will not be addressing this topic and the materi- pared by reacting stoichiometric quantities of reagents als in this review which is confined to the electronic and in evacuated, sealed quartz ampoules. Quartz begins magnetic properties of selenides relevant for the broad to soften for temperatures > 1150◦C and so for higher theme of the review paper discussing interplay between temperature reactions, a small pressure of inert gas may metallic and magnetic properties. be required to prevent the reaction tube imploding (Alt- Perhaps one of the greatest areas of interest in the mannshofer and Johrendt, 2008; Hawai et al., 2015; He search of novel electronic and magnetic properties is to et al., 2011). Care must be taken to ensure all reagents study materials that are close to metal-insulator transi- are thoroughly dry and that unwanted gaseous phases tions. This remains one of the least understood areas of will not be formed during the reaction (to avoid reac- 3 tion tubes exploding). We note that selenium is toxic, perature “metathesis” reactions have proven very suc- and is volatile above ≈ 640◦C and so low heating rates cessful for accessing low-temperature polymorphs and (and if necessary, a low temperature dwell) should be metastable phases (Martinolich et al., 2016, 2015; Mar- used to prevent a build-up of excess Se pressure and the tinolich and Neilson, 2014) tube exploding. Some reagents have been found to react Characterisation of single crystal samples of the ar- with quartz and this can often be minimised by plac- senide, oxyarsenide, and chalcogenide families of iron- ing the reagents in alumina or graphite crucibles inside based superconductors has given a much deeper under- the quartz tubes. For stoichiometry, an “oxygen getter” standing of the unusual (and often anisotropic) behaviour is often used (which can give more flexibility in choice of these materials than could have been gained with of starting reagents). Typical “getters” include Al or Zr polycrystalline samples (Nitsche et al., 2010; Wen et al., powders but the getter should be chosen based on the rel- 2011). Single crystals of LnFeAsO materials were first ative thermodynamic stability of various oxides, for ex- grown from a NaCl/KCl flux at high pressure, but the ample, see the Ellingham diagram for the relevant sample low solubility of the flux at the reaction temperature led cations and synthesis temperature. The metal powder re- to very slow crystal growth (Karipinski et al., 2009; Pro- acts with oxygen to form the oxide, but the reaction may zorov et al., 2009; Zhigadlo et al., 2008). Similar chal- be predominantly a surface reaction and so some excess lenges are faced when preparing single crystals of oxyse- oxygen-getter may be necessary. Some cations with two lenides. Most oxyselenide single crystals have been grown oxidation states close in energy can cause synthesis prob- from a flux, usually KCl, (Ijjaali et al., 2003)or a eutec- lems, as highlighted by Pitcher et al. (Pitcher et al., 2009) tic mix of NaI/KI, (Nitsche et al., 2014; Peschke et al., for CeCuOS; some earlier reports on “CeCuOS” had been 2015) but CsI (Tuxworth et al., 2015) I2 (Meerschaut + on samples with Cu vacancies and with some oxidation et al., 1998, 2001a) and Na2Sex (Park et al., 1993) fluxes of Ce3+ to Ce4+. This is likely to occur in oxyselenide have also been used. systems also (for example, slightly different unit cell pa- rameters have been reported for Ce2O2Fe2OSe2 (Mc- Cabe et al., 2014a; Ni et al., 2011)), and quenching B. Magnetic neutron scattering Ce2O2Fe2OSe2 and Ce2O2FeSe2 was found to minimise any phase separation on cooling (McCabe et al., 2014a). Neutrons are sensitive to both structural and magnetic Ainsworth et al. made a careful study of the effect of ex- properties making them ideal for studying the proper- cess getter and reaction temperature in the sample purity ties of magnetic oxyselenides. Given the erratic variation physical properties in Ce2O2ZnSe2 and were able to tune of the nuclear cross section with atomic number, neu- the Ce oxidation state and in turn, the band gap (and trons provide complementary structural information to colour) and unit cell volume (Ainsworth et al., 2015b). that of x-ray scattering. We focus here on recent devel- Post-synthesis reactions and treatments can also be opments in neutron inelastic scattering. While reviews of used to tune and optimise properties. For example, the technique have been provided elsewhere in the con- topotactic reactions could be carried out on many of text of strongly correlated electronic systems (Birgeneau these systems, particularly the layered materials. Hyett et al., 2006; Collins, 1989; Kastner et al., 1998; Shirane et al. (Hyett et al., 2007) were able to tune the an- et al., 2001), we focus here on recent developments in tiferromagnetic ordering temperature and the size of neutron inelastic scattering. Reviews on neutron diffrac- the ordered moment in Sr4Mn3O7.5Cu2Q2 by topotac- tionm including magnetic neutron diffraction, are pro- tic (oxidative) fluorination or (reducing) oxygen de- vided in Ref. Bacon, 1975; Izyumov and Ozerov, 1970; intercalation reactions and the same group have used and Willis and Carlile, 2009. lithium exchange and reductive Li+ insertion reactions In all neutron experiments, a neutron with a fixed inci- to tune magnetic and electrochemical properties (In- dent energy is either elastically scattered off a sample or dris et al., 2006; Rutt et al., 2006). Ammonia inter- inelastically scattered either by fluctuations in the lattice calation into layered FeSe has also been used to tune (phonons) or through magnetic interactions (magnons). its superconducting behaviour (Sedlmaier et al., 2014). In the case of unpolarised neutron scattering the scat- Other post-synthesis treatments to optimise microstruc- tering, cross section is a function of momentum transfer ~ ~ ture can have a significant effect of properties, for ex- Q~ = ki − kf and energy transfer ¯hω = Ei − Ef . Due to ample, textured samples of the thermoelectric BiCuOSe instrumentation and source qualities, neutron scattering (prepared by hot-forging) doubled the carrier mobility has historically been most successful at the “thermal” or and led to a dramatic improvement of thermoelectric be- lower energy range with typical energy transfers on the haviour (Sui et al., 2013). Lower-temperature solution order of ∼ meV. This energy scale is well matched for the methods have also been used to prepare oxyselenides, study of spin interactions and low energy lattice vibra- including nanoplates of Ln4O4Se3, (Gu et al., 2013) as tions, however the technique and selection rules associ- well as the single-step hydrothermal synthesis of poly- ated with it are limited for studying higher energy scales crystalline BiCuOSe (Stampler et al., 2008). Lower tem- on the order of ∼ eV. Given that magnetism provides 4 a probe of the underlying electronic ground state, we spin state to be derived. The elastic ordered magnetic will focus on the magnetic contribution to the scattering moment derived from neutron scattering is equal to gS cross section here. We first outline the cross sections for where g is the Lande factor and S is the spin of the the study of spin-spin correlations which is important in magnetic species. On integrating the differential neutron extracting interactions and the magnetic structure. We scattering cross section over all momentum transfers and then outline the cross section for single-ion excitations energy, the following relation is obtained. and the selection rules associated with these transitions. Spin Correlations: The differential neutron scattering Z Z cross section per element of solid angle dΩ, for wavevector d(¯hω) d3QS(Q,~ ω) = S(S + 1) (5) transfer Q~ and energy transferhω ¯ is, While the zero moment sum rule provides information 2 2 on the spin value which can be related to the local crystal- d σ (γr0) kf = S(Q,~ ω), (1) lline electric field, the first moment sum rule can provide dΩdω 4 k i information regarding local spin interactions. In partic- (γr )2 ular, as shown in Ref. Hohenberg and Brinkman, 1974, where 0 is 73 mbarns sr−1. In the dipole approxi- 4 in the case of isotropic Heisenberg exchange interactions, mation which is valid at small momentum transfers, the energy integrated S(Q) takes the form,

~ 2 2 X ˆ ˆ αβ ~ S(Q, ω) = g f (Q) (δαβ − QαQβ)S (Q, ω), (2) 2 1 X S(Q~ ) = − J hS~ · S~ i[1 − cos(Q~ · d~)] (6) αβ 3 ~ d 0 d~ (Q) i where the Sαβ(Q,~ ω) is related to the space and time where the sum (indicated by the subscript i) is over near- Fourier transform of the spin-correlation function, est neighbours and the bond distances d~. This sum rule becomes particularly practical in the simplifying case 1 X Z of where a single mode dominates the inelastic spec- Sαβ(Q,~ ω) = exp(iQ~ · (R~ − R~ )) dt... (3) 2π i j trum, termed the single mode approximation. In this ij case the full scattering structure factor can be written iωt α β e hSi (t)Sj (t)i. as S(Q,~ ¯hω) = S(Q~ )δ(¯hω − (Q~ )), where the δ-function forces energy conservation. This sum rule has in the past Neutron scattering therefore provides a direct probe of been applied to spin-chains but also to low-dimensional spin-spin correlations and hence the coupling between magnetic systems to extract dimensionality of the inter- ˆ ˆ them. The geometric term in the sum (δαβ − QαQβ) actions (Hammar et al., 1998; Hong et al., 2006; Stock provides a selection rule that neutron scattering is only et al., 2009; Stone et al., 2001). sensitive to the component of the magnetic moment per- Information on the intrinsic spin state can be deter- ~ pendicular to the momentum transfer (Q). The elastic mined through an integral of the total intensity, as dis- scattering (¯hω=0) component provides information on cussed below in the context of the local crystal field en- the magnetic structure through the correlation functions vironment. Also, from the first moment sum rule in- above while the inelastic component (¯hω 6=0) gives infor- formation on the dimensionality of the spin interactions mation on the coupling and energy terms of the magnetic can be obtained through a means which is model inde- Hamiltonian. pendent. This is particular important when measuring The scattering intensities measured from the differen- powders where directional information of the crystal is tial cross section can also be related to susceptibility us- heavily smeared by the powder averaging. ing thermodynamic techniques. In particular, the struc- In summary, neutron diffraction provides a means of ~ ture factor S(Q, ω) can be related to the imaginary part deriving the magnetic structure and sum rules obtained 00 of the susceptibility (χ ) using the fluctuation dissipation from the dynamic structure factor (measured through theorem which states, the inelastic neutron cross section) provide information about the local environment and also the exchange pa- αβ ~ rameters. αβ ~ −1 Im(χ (Q, ω)) S (Q, ω) = h  i. (4) Single-ion excitations: Whilst the above discussion in- π 1 − exp − hω¯ kB T troduced information that can be determined from neu- tron scattering cross sections, measurements of single- The usual static susceptibility can be found by taking the ion excitations can also give insight into the local crystal limit as Q~ and ω tend to zero. field environment. Given the energy scale of neutron in- Magnetic neutron scattering is subject to strict sum strumentation, these types of measurements are currently rules which allows direct information about the local best suited for rare earth or lanthanide elements where 5 spin-orbit coupling is large and the magnetism can be un- C. X-ray and optical spectroscopy derstood in terms of assigning a total angular momentum 2 J~ to the ion with eigenvalues of the operators J and Jz While efforts have been made to extend neutron scat- being equal to j(j +1) and m. The eigenstates can there- tering to higher energies approaching the ∼ eV energy fore be written in terms of these two operators as |j, mi. range, (Cowley et al., 2013; Kim et al., 2011; Stock The energy scale associated with changing j, termed in- et al., 2010b) studying single-ion transitions in d tran- termultiplet transitions, for the lanthanides is typically sition metal ions with neutrons is limited (owing to the quite large and of the order of ∼ eV while the energy energy scale and also current instrumentation) and opti- scale associated with changing m is much smaller, on or- cal and x-ray techniques are required. der ∼ meV (see for example the case of Pr in Refs. Taylor One important technique for studying high en- et al., 1988 and Turberfield et al., 1971). We therefore ergy transitions is Resonant Inelastic X-ray Scattering confine the discussion here to transitions where j is fixed (RIXS). If we consider the case where the incident pho- and only the eigenvalue m changes. ton energy is close to or above the core electron exci- In the dipole approximation for localised magnetic mo- tation threshold, the RIXS intensity can be written as ments, the neutron scattering cross section at small mo- follows, (Kotani and Shin, 2001) mentum transfers can be written as follows, X X hj|T |iihi|T |gi F (Ω, ω) = | |2... (8) Eg − Ω − Ei − iΓi 2 2 j i d σ (γr0) kf X = f 2(Q) ρ |hn|J |mi|2... (7) dΩdω 4 k n ⊥ ×δ(Eg + Ω − Ej − ω), i n,m

δ (En − Em − ¯hω) , where the operator T represents the radiative transition and Γi represents spectral broadening. The operator T is often calculated considering the dipole approximation where |ni, |mi are states belonging to a given J multiplet. and therefore subject to similar selection rules stated The operator J⊥ is the component of the total angular above for neutron scattering. The above expression for momentum operator perpendicular to the scattering vec- F (Ω, ω) illustrates that RIXS is the coherent second- ~ tor Q. The δ function enforces energy conservation. The order process consisting of the x-ray absorption from |gi effect of the form factor f(Q) is to decrease the magnetic- to |ii and the x-ray emission from |ii to |ji. Given the dipole transition intensities as the momentum transfer high energy scale of x-rays, this technique is particularly Q is increased. At higher momentum transitions, mag- useful for studying electronic excitations which are typi- netic octupole and higher-order transitions are possible. cally on the energy scale of ∼ eV and provides a comple- These will not be discussed in this review. Single-ion mentary technique to comparatively lower energy scale of transitions corresponding to intramultiplet transitions of neutrons which is primarily useful for investigating col- the form |ni → |mi are distinguished between spin-spin lective excitations. correlations discussed above. We note that single-ion ex- citations typically lack a strong momentum dependence while spin-spin excitations typical vary rapidly with Q. D. Transport and thermodynamic measurements Dispersing single ion excitations can exist and need to be treated in terms of a multi-level spin-wave analysis. We Electrical resistivity and thermodynamic measure- will not discuss these cases in this paper. ments are key tools used to characterise oxyselenides, Neutron scattering is a powerful and continuously which often have semiconducting properties that can be evolving technique for the study of magnetism and elec- tuned towards metallic or insulating extremes. Metallic tronic phenomena. Historically, it has been confined to materials show a vanishing resistivity as the tempera- lower energy transitions for the spin correlations and in- ture is lowered and insulators show a diverging resistiv- tramultiplet transitions outlined above. One reason for ity. Semiconductors show behaviour between these and this originates from kinematics and the fact that the neu- the simplest model for the electrical conductivity is an tron has mass meaning that high energy transfers usually activated behaviour with ρ = ρa exp(Ea/kBT ) with Ea correspond to large momentum transfers. While the en- being an approximate measure of the band gap. Resis- ergy scales are typically relevant for the study of transi- tivity therefore provides a fairly simple means of charac- tions of interest in condensed matter physics, one limita- terising the electronic properties of materials. tion is the study of local environment around d transition To study the localised properties associated with mag- metal ions where the crystal field excitations are close to netic moments, magnetic susceptibility provides a means ∼ eV. In this context x-ray and optical techniques have of studying the real part of the magnetic suscepti- played a pioneering role allowing high energy transitions bility χ (defined above). The magnetic susceptibil- to studied. ity can be related to neutron scattering techniques via 6

units. A previous review of layered oxychalcogenides and oxypnictides describes the structures of many transition metal oxyselenides (Clarke et al., 2008) and allows us to highlight some structural features common to most oxy- selenides:

• The different sizes and characters of 1st row oxide O2− and 3rd row selenide Se2− anions usually give rise to anion-ordering.

– The hard, polarising O2− anions tend to be co- ordinated by harder cations (often Ln3+ ions, as in Ln10OSe14 phases (Weber et al., 2012)) FIG. 1 The structures of La2O2S, La2O2Te and La4O4Se3. and it can be helpful to consider the parts of La=green, O=red, S/Te/Se=yellow, respectively. This figure the structure in terms of oxide-centred tetra- is reproduced from Ref. Tuxworth et al., 2015. Fluorite-like 2+ hedra (Krivovichev et al., 2013). [La2O2] layers are highlighted in red. – “Softer” transition metals are usually coordi- nated by the more covalent Se2− ions. the fluctuation-dissipation theorem stated above or also through the Curie-Weiss constant derived through high • This anion-ordering often gives layered structures temperature measurements. Neutrons are sensitive to with quite different properties associated with the the local spin-spin interactions from which a set of ex- more ionic oxide and more covalent selenide layers change constants can be derived based on a model mag- which may be electronically isolated from one an- netic Hamiltonian. Assuming Heisenberg exchange, the other. This anisotropy is often key to understand- ing the properties of the whole material. Curie constant (ΘCW ) derived from susceptibility can be related to the exchange constants, Compared with oxides, mixed-anion systems, such as the oxyselenides discussed here, allow us to prepare mate- 1 X rials with unusual cation coordination environments (of- Θ = S(S + 1) J (9) CW 3 n ten in low oxidation states). This, combined with their n layered nature, gives them often unique properties. The with the sum being over coupled neighbours. The ex- structures of several oxyselenides and the relationships change constants derived from scattering experiments between them have been described in Ref. Clarke et al., can therefore be cross checked against susceptibility and 2008 and we discuss below only some of the dominant in particular magnetisation. or more recent oxyselenide structure types to allow us In the case of itinerant magnetism, the resistivity to consider their electronic structures and properties in from the spin fluctuations can be calculated from the more depth later in this review. neutron scattering cross section using the following for- mula, (Moriya et al., 1990), A. Ln – O – Se phases

Z ∞ E E  eE/T Z The Ln – O – Se phases with no transition metal pro- ρ(T ) ∝ T d d3qχ00(~q, E). E/T 2 vide illustrations of the structural role of the two anions −∞ T T (e − 1) of different characters. The Ln4O4Se3 (Ln = La - Nd, Here, χ00 is the spin susceptibility related to the mea- Eu - Er, Yb and Y) family adopts crystal structures com- ~ ~ 1 4+ sured neutron intensity I(Q, E) ∝ S(Q, E) = π [n(E) + posed of fluorite-like [Ln4O4] oxide sheets built from 1]χ00(Q,~ E) with [n(E) + 1] being the bose factor. The edge-linked Ln4O tetrahedra. These sheets are separated 2− 2− temperature dependence of the resistivity can therefore by layers containing both Se and diselenide Se2 an- be related to the spin fluctuation in itinerant magnets. ions (Strobel et al., 2008; Tuxworth et al., 2015). Whilst 4+ This is discussed below in the context of comparing iron the [Ln4O4] sheets change little with Ln, the arrange- 2− 2− based chalcogenides and oxyselenides. ment of Se and Se2 anions in the interlayers varies with Ln3+ ionic radii (Tuxworth et al., 2015). Long- range antiferromagnetic order occurs for Ln = Gd, Tb III. OXYSELENIDE STRUCTURAL FAMILIES and Dy phases but has not been observed down to 1.8 K for other analogues. (Strobel et al., 2008; Tuxworth It is helpful to begin by considering the structural et al., 2015) Strobel et al. highlight the magnetic frus- chemistry of oxyselenides and the common structural tration within the Ln − O network as a result of arrang- 7

ƌƵ^ŝƐͲƐƚƌƵĐƚƵƌĞ ƌƵ^ŝƐʹƌĞůĂƚĞĚƐƚƌƵĐƚƵƌĞƐ >ĂLJĞƌĞĚƐƚƌƵĐƚƵƌĞƐǁŝƚŚĨůƵŽƌŝƚĞ Ğ͘Ő͘>ĂƵK^Ğ ďͿ ĐͿ &Ğ ŽƌĂŶƚŝĨůƵŽƌŝƚĞůĂLJĞƌƐ ĂͿ ď ^Ğ >Ŷ ĨͿ K Ă Ğ ^Ğ ^Ğ Ƶ Ƶ Ă K K Đ K Ğ ď Đ >Ŷ ď Đ ^Ğ Ă Ě ŝ Ă ď ď Ă ĞϮKϮ&Ğ^ĞϮ >ĂϮKϮĚ^ĞϮ

ŝϮ>ŶKϰƵϮ^ĞϮ ŚͿ ŐͿ ^ď >Ŷ DŶ K ^ƌ ^Ğ D K ĚͿ DŶ DŶ ^Ğ ĞͿ Đ DŶ^ĞϮŝŶ;Ğ͕>ĂͿϮKϮDŶ^ĞϮ ď Ŷ ^Ğ Ă

^ƌϮDŶKϮDŶϮ^ďϮ >ŶϮKϮDϮK^ĞϮ Ŷ^ĞϮŝŶ>ĂϮKϮŶ^ĞϮ 

FIG. 2 (a) ZrCuSiAs structure adopted by LaCuOSe shown in centre; right hand side shows ZrCuSiAs-related structures with (b) checkerboard-ordered La2O2CdSe2,(c) stripe-ordered Ce2O2FeSe2 and MSe2 layers in (d) (CeLa)2O2MnSe2 and (e) La2O2ZnSe2; left hand side of figures shows layered structure types including (f) Bi2LnO4Cu2Se2,(g) Sr2MnO2Mn2Sb2, and (h) Ln2O2M 2OSe2. La/Ce/Sr, Bi and M (M=Cu, Cd, Fe, Mn, Zn), O and Se shown in green, purple, blue, red, and yellow, 2+ respectively. Fluorite-like [Ln2O2] layers are highlighted in red whilst antifluorite-like [MxSe2] layers are highlighted in blue to emphasise the relationships between structure types. ing magnetic Ln3+ ions on a tetrahedral motif. (Stro- –O– M – Se phases is the ZrCuSiAs structure, (Johnson bel et al., 2008) Four compositions have been reported and Jeitschko, 1974) also referred to as the 1111 structure for this family of materials including A10OSe14 (A=La- adopted by LnFeAsO parent phases to the iron-based su- Nd) (Weber and Schleid, 2001; Weber et al., 2012; Wu perconductors, (Kamihara et al., 2008) and this structure and Huang, 2007), A2OSe2 (A=Pr, Gd) (Tougait and type and variations upon it dominate oxyselenide chem- Ibers, 2000; Weber and Schleid, 2001), A4O4Se3 (A=La- istry. Stoichiometric oxyselenides in this family include Nd, Sm) (Dugue et al., 1970; Strobel et al., 2008; Weber LnCuOSe (Ln = La - Sm) (Hiramatsu et al., 2004a,b, and Schleid, 2001), and A2O2Se (A=La,Pr,Nd, Sm, Gd, 2008; Kamihara et al., 2008; Llanos et al., 2006; Llanos Er, Ho, Yb) (Eick, 1960; Weber and Schleid, 2001). The and Pena, 2005; Mizoguchi and Hosono, 2011; Ueda et al., structures are illustrated in Fig. 1. 2000, 2003) containing monovalent Cu+ ions. These ma- terials are particularly well known for their wide band gap, optically transparent and p-type semiconducting B. Ln – O – M – Se phases behaviour (Hiramatsu et al., 2008; Ueda et al., 2000, 2003). Their tetragonal crystal structure is built from 1. ZrCuSiAs structures and related cation-ordered phases fluorite-like oxide layers of edge-linked OLn4 tetrahedra separated by anti-fluorite-like layers of edge-linked CuSe4 + The fluorite-like layers of edge-linked LnO4 tetrahedra tetrahedra (Fig. 2 (a)). Replacing monovalent Cu are widespread in oxyselenide structural chemistry and ions with divalent M 2+ ions (e.g. M = Mn2+, Fe2+, 2+ 2+ most Ln –O– M – Se phases contain this motif or varia- Zn , Cd ) leads to half-occupancy of M Se4 tetrahe- tions of it. Perhaps the simplest structure adopted by Ln dra. These M 2+ sites can be occupied in a disordered 8

ĂͿ ďͿ leaves the M (2) site in roughly tetrahedral (M (2)Se4) 2− 'Ě sites, while M (1) ions are coordinated by both O and dŝ Se2− anions in pseudo-octahedral coordination (McCabe et al., 2010). This structure is very closely related to the D;ϮͿ ^Ğ 4+ >Ă Ln4O4TiSe4 structure with Ti cations occupying only K K the M (1) sites (Fig. 3 (b)) (Meerschaut et al., 2001b; Tuxworth and Evans, 2014). D;ϭͿ Nitsche et al., 2014 revealed the polymorphism of the ^Ğ Ln2O2FeSe2 (Ln = La, Ce) systems and were able to tune the iron coordination environment with synthe- sis temperature. At high temperatures, stripe-ordered

ĐͿ >Ă ZrCuSiAs-related phases were formed with tetrahedrally- 2+ K coordinated Fe ions; at intermediate temperatures the β-phases were formed with both tetragonal and pseudo- ^Ğ octahedral Fe2+ coordination. At low temperatures a &Ğ new phase, with only pseudo-octahedral coordination of Fe2+ was formed (Fig. 3 (c)).

3. ZrCuSiAs-modified structures FIG. 3 (a) β-La2O2MSe2 structure adopted for M = Mn, Fe, (b) Ln O TiSe structure and (c) low temperature mono- 4 4 4 Clarke et al., 2008 describe how the ZrCuSiAs clinic polymorph of La O FeSe with only pseudo-octahedral 2 2 2 structure of LnCuOSe can be modified to accom- FeSe4O2 coordination. Ln, M (M= Mn, Fe, Ti) O and Se ions are shown in green, blue, red and yellow, respectively. modate thicker fluorite- or antifluorite-like layers, or even additional layers. One example of the latter − is Bi2LnO4Cu2Se2 in which [Bi2LnO4] blocks (with 3+ fashion (as originally reported for CeMn0.5OSe) (Ijjaali Ln in square-prismatic coordination sandwiched be- et al., 2003) but this partial occupancy often gives rise to tween two Bi - O fluorite-like layers) alternate with 2+ + ordering of the M ions. This cation ordering might fol- [Cu2Se2] sheets (Fig. 2 (f)), stabilising the mixed- low a checkerboard arrangement (as in La2O2CdSe2, (Hi- valent Cu ions and giving metallic conductivity (Chou ramatsu et al., 2004a,b) Fig. 2 (b)), a stripe arrangement et al., 2015; Evans et al., 2002). Another possibility in- (as in Ce2O2FeSe2, (McCabe et al., 2011, 2014a) Fig. 2 volves swapping the fluorite-like oxide layers for oxygen- (c)) or a combination of these. Intermediate structures, deficient perovskite-like A2MnO2 oxide layers to give the containing stripe- and checkerboard-ordered regions were Sr2MnO2Mn2Sb2-type structure (Fig. 2 (g)), (Brechtel reported for Ln2O2ZnSe2 (Ainsworth et al., 2015b; Tux- et al., 1979) variations of which are adopted by several worth et al., 2013) (Fig. 2 (e) ) and (La,Ce)2O2MnSe2 oxychalcogenides (Herkelrath et al., 2008; Jin et al., 2012; (Fig. 2 (d)) (Peschke et al., 2015) Ainsworth, Wang Otzschi et al., 1999; Tan et al., 2014; Zhou et al., 2014; et al. (Ainsworth et al., 2015a; Wang et al., 2015) Zhu, 1997; Zhu et al., 1997). These can also be modi- have shown recently that the “infinitely adaptive” or- fied to include thicker antifluorite-like Cu2S2 layers in the dering in the La2−xCexO2M Se2 solutions can be tuned series Sr2MnO2Cu2m−δSm+1 (Barrier and Clarke, 2003; between the stripe and checkerboard extremes (via Gal et al., 2006). incommensurately-modulated structures) by the Ln3+ ionic radius and Peschke and Johrendt (Peschke et al., 3+ 2015) have extended this to the La2−xLnxO2MnSe2 (Ln 4. LnMOSe2 - M ions coordinated by selenide = Pr, Nd) solid solutions and confirmed the role of Ln3+ ionic radius. Relatively few Ln -O- M - Se phases with M 3+ coor- dinated by selenide are known (although M 3+ oxysulfide chemistry is more extensive (Ceolin and Rodier, 1976;

2. β-La2O2MSe2 and related polymorphs Dugue et al., 1980a,b; Jaulmes, 1978; Jaulmes et al., 1982; Luo et al., 2013; Ogisu et al., 2008; Winterberg Initial attempts to prepare ZrCuSiAs-related et al., 1989) presumably due to the redox chemistry and 2− La2O2M Se2 (M = Mn, Fe) phases lead to the dis- the challenge of maintaining the reduced Se anions 3+ covery of a new structural family, β-La2O2M Se2. In in the presence of the more oxidising M ions. This 2+ this β-structure, the fluorite-like [Ln2O2] layers are challenge has been overcome in oxyselenide systems for 3+ 3+ 3+ broken into ribbons and arranged in a herringbone-like M = Ga and Cr . LnGaOSe2 is again built from 2+ fashion, separating M Se2 layers (Fig. 3 (a)). This fluorite-like [La2O2] oxide layers but these are sepa- 9

ĂͿ ďͿ   
>Ŷ  ^ K
 ƌ;ϮͿ >Ă
 Y ƌ;ϭͿ  b K ƌ
 a
 ď c 

Ă  Đ



FIG. 4 Structure of (a) LaCrOQ2,(b) LnCrOS2 (Ln = Pr, 


  
 

!  
 
  
 Nd). (c) Ln, Cr, O and Q ions are shown in green, blue, red 


  
! °  



  
 °  
 and yellow, respectively. 

FIG. 5 Structure of Ln2O2M 2OSe2 showing (a) tetragonal unit cell, (b) pseudo-octahedral FeSe4O2 coordination polyhe- 2− rated by [Ga2Se4] double layers with the relatively dra, and (c) Fe-Se-O layers illustrating exchange interactions. 3+ small Ga ions in quite distorted GaSe4 tetrahedra (Be- Ln, M (M = Mn, Fe, Co) O and Se ions are shown in green, nazeth et al., 1984). blue, red and yellow, respectively. The crystal structures of the LnCrOSe2 phases are quite different and vary with Ln3+ (Fig. 4). The oxyse- lenide and oxysulfide LaCrOQ 2 (Q = S, Se) are isostruc- BiX2 and fluorite-like AO layers in the LnOBiS2 super- tural with ribbons of edge-linked pseudo-octahedral conducting family, whilst the A2O2B2Se3 family consists CrQ5O extending along c, linked by fluorite-like chains of quasi-one-dimensional ribbons of edge-linked BiSe4O (Fig, 3 (a)) (Dugue et al., 1980a; Wintenberger et al., square pyramids, linked to fluorite-like SrO fragments by 1987). The crystal structure of CeCrOSe2 is again differ- the apical Bi - O bond. This appreciation of the square ent (and analogous to that adopted by several LnCrOS2 pyramidal units and their connectivity with fluorite-like phases (Wintenberger et al., 1987; Winterberg et al., units opens the possibility to design and prepare further 1989)) with the Cr-Q ribbons of LaCrOQ 2 broken to oxychalcogenides (of varying dimensionality) containing form chains of edge-linked CrSe4O2 pseudo-octahedra ex- “inert pair” ions (Panella et al., 2016). tending along c and linked along a by fluorite-like oxide chains. Unlike in LaCrOQ 2, the mixed Ce - O - Cr - Se layers are also linked along b by edge-linked CrSe6 5. Ln2O2M 2OSe2 materials octahedra (Wintenberger et al., 1987). The structural chemistry and physical properties of these LnCrOSe2 An important family of oxyselenides adopt an “anti” phases have not been fully explored but the related oxy- form of the Sr2MnO2Mn2Sb2-type structure with cation sulfides show interesting magnetic behaviour with strong and anion sites swapped: in the Ln2O2M 2OQ 2 (M = 3+ 3+ 2+ coupling between the Ln and Cr magnetic sublat- Mn, Fe, Co; Q = S, Se), fluorite-like [Ln2O2] lay- 2+ tices (Takano et al., 1999, 2002; Wintenberger et al., ers and [Fe2O] sheets (in an anti-CuO2 type arrange- 1987; Winterberg et al., 1989). ment (Park et al., 1993)) are separated along c by lay- 2− Recently a series of insulating A2O2B2Se3 (A = Sr, ers of Se anions (Figs. 2 (h) and 5). This tetrago- Ba; B = Bi, Sb) oxyselenides, consisting of double- nal structure was first reported by Mayer et al. (Mayer chains of edge-linked BiSe4O square pyramids, have et al., 1992) for La2O2Fe2OQ 2 (Q = S, Se) and has since been reported (Panella et al., 2016). These materials been extended to include a range of Ln3+ and M 2+ ions are structurally-related to the LnOBiS2 superconduct- (Ln = La, Ce, Pr, Nd, Sm; M = Mn, Fe, Co) (Free ing family (Lei et al., 2013; Lin et al., 2013; Mizoguchi et al., 2011; Fuwa et al., 2010a,b; Kabbour et al., 2008; and Hosono, 2011; Yazici et al., 2013) which also contain Lei et al., 2012; Liu et al., 2011, 2015b; Ni et al., 2010, Bi-centred square-based pyramids, compatible with the 2011; Popovic et al., 2014). 3+ 3+ 2+ “inert pair” Bi and Sb ions. The structures adopted The M environment in the Ln2O2M2OQ 2 structure by these two families can be considered to be built from is unusual with coordination by two O2− ions within fluorite-related [AO] units (layers, ribbons or discrete the plane and by four Se2− ions above and below the units) separating [B2X3] layers: the superconducting bis- plane, giving M O2Q 4 pseudo-octahedral coordination muth sulfides can be written [LnO]2[Bi2S3]S and the (similar to that of the M (1) sites in the β-La2O2M Se2 A2O2B2Se3 family as [AO]2[B2X3], (reflecting the ad- and Ln4O4TiSe4 structures described above). Closely 2+ ditional anion needed for charge balance with the triva- related Na2OFe2S2 has the same [Fe2O] sheets but 3+ lent Ln ions in the LnOBiS2 superconducting fam- here, the fluorite-like oxide layers are replaced by lay- + ily) (Panella et al., 2016). This results in two-dimensional ers of Na ions (He et al., 2011) and Na2−xCu2Se2Cu2O 10

2− also has analogous [Cu2O] sheets as well as [Cu2Se2] layers (Park et al., 1993). All these systems order anti- ferromagnetically on cooling with N´eeltemperatures that can be tuned with Ln ionic radius and will be discussed in more detail below. Free et al. (Free et al., 2011) considered the role on Ln3+ size on the range of first row transition metals that could be accommodated on the M sites in terms of size- 2+ 2+ mismatch between [Ln2O2] and [M 2O] layers, with M=Fe2+ found to be compatible with the widest range of lanthanides. Similar compounds were discussed by Ni et al. (Ni et al., 2011) for M=Fe2+. The redox chemistry of the M ions (ie balancing the oxidising ability of the M3+ ions in the presence of selenide ions) is also likely to play a role. 2+ The [Ln2O2] fluorite-like layers are analogous to those in the ZrCuSiAs structure type and are subject to similar structural distortions induced by Ln3+ ions. While extra structural peaks were reported in Ref. Free et al., 2011 in temperature dependent x-ray diffraction data for M=Mn2+ variants, these could not be in- dexed by a commensurate unit cell. Besides this, no FIG. 6 The local coordination site taken from Ref. Rossat- Mignod et al., 1974 which studied the local crystalline electric observable low temperature structural transitions have 3+ field in Tb2O2S. Nearest neighbours of a R ion in oxy- been reported for these compounds with the exception sulides and oxyselenides (Tb,Dy,Ho) O (S,Se). The empty 3+ 3+ 2 2 2 of Ln=Pr (Free et al., 2011). Pr (4f , a non- and cross-hatched circles are oxygen and sulfur/selenium, re- Kramers ion) on the high symmetry (C 4v or 4mm site spectively. in the I 4/mmm tetragonal crystal structure) tends to drive an orthorhombic distortion, lowering the symme- 3+ try of the Pr coordination environment (e.g. to C 2v scattering. or 2mm in the low temperature Immm phase). This has been observed in the mixed-anion ZrCuSiAs-structure phases PrMnSbO (Kimber et al., 2010) and PrMAsO A. Local magnetism (M = Mn, F e), (Kimber et al., 2008; Wildman et al., 2015) as well as in the oxyselenides Pr2O2Mn2OSe2 and Many oxyselenides are based upon two magnetic sites Pr2O2Fe2OSe2 (Free et al., 2011; Oogarah, 2016). For with one being a 3d transition metal ion and the other a Pr2O2Mn2OSe2 a structural transition to an orthorhom- rare earth lanthanide site. The local magnetism on both bic unit cell is found for temperatures below ∼ 50 K (Free sites is treated differently with the rare-earth magnetism et al., 2011) with a=4.08616(3) A˚ and b=4.09417(4) treated in terms of j − j coupling and the 3d site being A˚ (Free et al., 2011). understood in terms of L − S coupling. Here we discuss the results for the two sites with a discussion of the rare earth local magnetism followed by a discussion of the IV. MAGNETIC PROPERTIES local magnetism on the 3d transition metal ion site.

In this section we provide a review of the magnetic properties of the oxyselenides. This section is divided 1. Rare Earth local magnetism into three parts discussing the local magnetism, magnetic structures, and then magnetic interactions. We first pro- Oxyselenides containing only rare earth magnetic ions vide a discussion of the local magnetic environment in were a topic of interest in the 1970’s. Rare earth lo- terms of the crystalline electric field which is central to cal magnetism is treated in terms of j −j coupling where understanding the magnetic ground state in any material. the spin-orbit coupling is much larger than the local crys- We then discuss the various magnetic structures reported talline electric field. Typical energy scales for the crys- in the oxyselenides. With the recent interest in strongly talline electric field are ∼ meV while the spin-orbit cou- correlated electron systems, the studies have primarily pling is on order of ∼ 1 eV. In the j − j coupling scheme, focussed on two dimensional variants. We then finish the spin-orbit Hamiltonian is therefore diagonalised first with a discussion of the available work on the magnetic with the crystalline electric field treated as a perturba- interactions primarily probed through neutron inelastic tion on the ground state. Following the formalism in 11

FIG. 7 Crystal field scheme determined using Stevens opera- FIG. 8 Schematic illustrating the d orbital splitting for Co2+ tors and also measured with optical spectroscopy. The figure with the crystal field being octahedral (a) to square planar is taken from Ref. Rossat-Mignod et al., 1974 and shows a (d) via a slightly tetragonally elongated octahedron (b) and comparison between the observed and calculated levels of 7F J the highly tetragonally elongated octahedron (b). The fig- manifold. ure is taken from Ref. Smura et al., 2011 in the context of (Sr,Ba)2CoO2Cu2S2. terms of Steven’s operators, the single ion Hamiltonian (H ) can be written as, CEF oxypnictides and also two dimensional oxyselenides that these systems display well defined crystal field excitations X m m and levels from the rare earth site. Generally these ex- HCEF = Vl θlOl (10) l,m citations do not disperse substantially indicating weak coupling between the rare earth sites. The local nature m where Vl are the adjustable crystal field parameters, of these excitations has been used to probe crystalline m Ol the Stevens operators (which are functions of the electric field symmetry in rare earth oxypnictides. (Gore- total angular momentum operator J, with J 2|j, mi = mychkin et al., 2011; Xiao et al., 2013) The well defined j(j+1)|j, mi and Jz|j, mi = m|j, mi), and θl are the mul- nature of the crystal field excitations is also indicative tiplicative factors that depend on the free ion level and of their localised nature in contrast to expectations for are tabulated (Hutchings, 1964). An important result in itinerant systems which would not show well defined ex- understanding crystalline electric fields is Kramer’s the- citations in energy. orem which states that the crystal electric field scheme of ions with half-integer j is made up of doublets. This is always the case for such ions and the degeneracy can 2. Transition metal ion local magnetism only be broken with a field which breaks time reversal symmetry, such as a magnetic field but not an electric. The basic building block of the magnetic properties of Some of the first rare earth oxyselenides studied were the oxyselenides is the local crystalline electric field envi- (Tb,Dy,Ho)2O2(S,Se) and the local coordination geome- ronment surrounding the transition metal ion (examples 2+ 2+ try is shown in Fig. 6 for Tb2O2S (Abbas et al., 1973, include Fe or Co ). As outlined elsewhere (Abragam 1974). The case of Tb2O2(S,Se) (Abbas et al., 1974; and Bleaney, 1986; Ballhausen, 1962; Griffiths, 1964; Mc- Rossat-Mignod et al., 1974) is particularly interesting as Clure, 1959), there are two competing interactions in the magnetic structure involves a large canting of the understanding the local single-ion magnetic properties - Tb3+ moment which is not expected given the symmetry the crystal field splitting which splits the degeneracy of of the lattice. It was therefore concluded in this sys- the d orbital levels and the Hund’s energy scale which tem that the ground state could not be considered as a characterises the energy barrier for allowing double oc- doublet, but as a set of doublets with the energy scale cupancy of the orbitals. The competition between the between the low-energy crystal fields to be small and of two energy scales is described through Tanabe-Sugano order ≤1-2 meV. These were studied using optical spec- diagrams which show the energy of electronic states rela- troscopy and compared against predictions from the sus- tive to the ground electronic state as a function of crystal ceptibility and the magnetic structure with reasonable field strength. agreement being obtained (shown in Fig. 7). The mag- The competition between these two energy scales can netic structure of Dy2O2(S,Se) was found to be more un- be seen to cause an uncertainty in the value of S for cer- axial consistent with expectations from the crystal sym- tain d transition metal ions and hence an ambiguity in metry. the sum rules discussed above in the context of neutron These early studies of rare earth oxyselenides and sul- scattering and also the magnetic structure. For example, fides illustrate a common theme reflected in studies on an octahedral coordination environment will split the d 12

(a) 3 (b) Co 3d xz 2 xy 22 x−y22 3z −r 1 yz 0 -2 -1.5 -1 -0.5 0 0.5 1 1.5 1 (c) Co 3d yz La 1 xz Se O 1 Co xy 1 x z x−y2 2 1 y 3z2 −r2 1 Se 4p

Density of States1 (states/eV) O2p 0 -8 -6 -4 -2 0 2 4 6 Energy (eV)

FIG. 9 The crystal field scheme and local crystalline electric field taken from Ref. Wu, 2010. (a) shows the crystal struc- ture of La2O2Co2OSe2 having Co2Se2O layers. (b) Orbitally resolved Co 3d density of states for the nonmagnetic state taken from the first principle calculations in Ref. Wu, 2010. ~ 2 2 The low-lying t2g-like orbitals (x − y , yz, and xz) and the 2 2 higher eg (xy and 3z − r ) are split. (c) Density of states of the high-spin state.

o orbitals into a lower energy triplet, referred to as the |t2gi orbitals and a higher energy doublet, the |e2gi orbitals. FIG. 10 The zeroeth moment sum rule applied to powders Conversely, a tetrahedral coordination environment will of La2O2Fe2OSe2 from McCabe et al. (McCabe et al., 2014b) give a much smaller splitting and a lower energy doublet a) illustrates the low temperature powder averaged spectrum (|ei orbitals) and a higher energy triplet (|t2i orbitals). taken using the MARI spectrometer with Ei=100 meV. b) il- For a tetrahedrally coordinated Fe2+, there are two sce- lustrates the averaged integrated spectral weight as a function narios for populating the five 3d orbitals with 6 elec- of Q which can be related to the zeroeth moment sum rule. The dashed line is predicted value for S=2. trons: for weak or intermediate crystal fields, the Hund’s energy scale dominates and the “high spin” configura- 3 3 tion is found, with egt2 (S = 2); for larger crystal fields, the crystal field splitting dominates and the “low spin” ple calculations in panel (b) also resolve the different or- configuration is found, with e4t2 (S = 1). These scenar- bital contributions with the relative energies consistent g 2 2+ ios give different S and also different possible magnetic with the local D2h (mmm) Co site symmetry. The Hamiltonians with the later introducing an orbital degree deviation from a perfect octahedral environment has di- 2+ of freedom. This distinction between weak/intermediate rect implications for the case of Fe oxyselenides and and strong crystal field limits has been suggested to give their magnetic excitation spectrum. As noted in Dai and 6 rise to a spin state transitions or additional orbital terms Whangbo, 2005, high spin d ions in an undistorted octa- in the magnetic Hamiltonian (Gretarsson et al., 2013; hedral environment cannot have uniaxial magnetic prop- Kruger et al., 2009; Stock et al., 2014). erties. The above discussion does not solve the ambiguity In Figs. 9 and 8, we consider the example of over which energy scale dominates and hence - whether 2+ we should expect high spin or low spin Fe2+ ions. To ad- Co in the two dimensional La2O2Co2OSe2 with its dress this we consider evidence from magnetic diffraction pseudo-octahedral CoO2Se4 coordination (taken from Refs. Smura et al., 2011 and Wu, 2010). Figure 9 shows and sum rules from neutron inelastic scattering for the 2+ the crystal structure highlighting the local environment case of Fe oxyselenides. around the Co2+ site. Figure 8 illustrates the crystal Fig. 10 shows an analysis of the zeroeth moment sum field scheme, assuming a point charge distribution, given rule in a powder sample of La2O2Fe2OSe2 from inelastic different environments with (a) illustrating a perfect oc- neutron scattering experiments (from the supplementary tahedron and a highly distorted case in (d). First princi- information in Ref. McCabe et al., 2014b). Panel a) 13 shows the low temperature powder averaged spectrum wide range of materials containing both a series of lan- and panel b) shows a plot of the average spectral weight thanides and also various d transition metal ions (par- as a function of Q, ticularly Mn2+, Fe2+, and Co2+). Because of the strong cross section for magnetic moments and the ability to study spin correlations, magnetic neutron scattering has R +∞ R Q 3 dE d qS(~q, E) played a central role in these studies. Given the interest I˜(Q) = −∞ 0 = S(S + 1). (11) R Q d3q in iron based systems and the structural similarities with 0 iron based superconductors, there has been a number The dashed line in Fig. 10 (b) is the theoretical value for of studies of two-dimensional variants iron based oxyse- S=2. The good agreement between the powder average lenides. We first briefly outline currently available results spectral weight and the zeroeth sum rule value for S=2 in oxyselenides that only host a lanthanide magnetic ion implies that Fe2+ has 4 unpaired electrons, consistent and then discuss oxyselenides where both d transition with a high spin configuration for the pseudo-octahedral metal ions and lanthanides are present. The structure of this section largely follows the outline presented above environment in La2O2Fe2OSe2. Evidence supporting the high spin state S = 2 (weak- classifying the structural types. intermediate crystalline electric field) is also found from neutron diffraction. Table I shows the refined or- C. Ln – O – Se phases dered magnetic moments for the oxyselenides where mag- netic diffraction and complete refinement have been per- The magnetic structural properties of Ln-O-Se phases formed. From diffraction theory, the ordered moment where no d transition metal ion is present have not been should be equal to gS and it can be seen that all of the as fully investigated as transition metal oxyselenide sys- ordered moments are ∼ 3 µB. Taking the Lande factor tems. Early work on Tb2O2(S,Se) and Dy2O2(S,Se) de- g=2, these results are inconsistent with low spin state scribed above motivated crystal field work to understand (strong crystal field limit) yet clearly closer to S = 2. the unusual canting measured in the Tb variant (Ab- The values are very close to the reported value of 3.3 2+ bas et al., 1973, 1974; Rossat-Mignod et al., 1974). For µB for FeO further corroborating the fact that Fe is A4O4Se3 compounds long-range antiferromagnetic order in a high spin state in these compounds (Roth, 1958). occurs for Ln = Gd, Tb and Dy phases but has not Given the total moment sum rule applied to the neutron been observed down to 1.8 K for other analogues (Strobel inelastic scattering spectrum and a summary of the or- et al., 2008; Tuxworth et al., 2015). Interestingly, these dered moments from neutron diffraction and magnetic later compounds have been suggested to show geometric refinement, we conclude the local crystalline electric field frustration based on the lack of obvious magnetic order environment surround the transition metal ion in the oxy- 10+ and the fact that spins on the [OM4] tetrahedron can- selenides is in a weak-intermediate crystal field limit. not all minimise the magnetic Hamiltonian in analogy to The high spin S=2 nature of the local magnetism in rare earth pyrochlore lattices (Gardner et al., 2010). iron oxyselenides is a distinguishing point over their iron arsenide, pnictide, and chalcogenide counterparts (Yin et al., 2011). For example, in Fe1+xTe the momentum D. Ln – O – M – Se phases and energy integrated spectral weight over the spectrum up to ∼ 150 meV only yields a value consistent with a 1. LnCrOS2 (oxysulfides) value of S slightly larger than S=1, and well below the value expected for high spin S=2 (Stock et al., 2014). Although most work on LnCrOQ2 systems has been Consistent with this general statement, magnetic diffrac- carried out on the oxysulfides, the magnetic interactions tion studies of FeAs only find an ordered moment 0.5 ± present and the resulting magnetic structures give some 0.05 µB (Rodriguez et al., 2011b), a value much less than insight into the structurally related oxyselenide materials the ordered moment expected for S=2. Indeed, as tab- discussed below. For completeness we therefore provide ulated in Table 10 in Ref. Johnston, 2010, the magnetic an overview of the results in this section and then re- ordered moment in the iron based pnictides and chalco- turn to our review of oxyselenides. As described above, genides is universally less than ∼ 1 µB. the crystal structure of LaCrOQ2 is composed of dou- ble chains of edge-linked CrQ5O pseudo-octahedra with intrachain Cr - Cr distances of ∼ 3.4 A˚ and ∼ 3.7 A,˚ B. Magnetic structures and Cr - Q - Cr angles of ∼ 90◦ and ∼ 100◦. Sus- ceptibility measurements indicate that Cr3+ ions order Having discussed the local magnetic properties of oxy- ferromagnetically (TC = 35 K (LaCrOS2); TC = 51 K selenides, we now discuss the magnetic structure which is (LaCrOSe2)) (Wintenberger et al., 1987) and recent mag- sensitive to exchange interactions. The magnetic struc- netization and specific heat studies suggest that LaCrOS2 tures of the oxyselenides have been investigated for a can be described in terms of Ising chains with relatively 14

TABLE I Ordered Magnetic moments in several Fe2+ oxyselenides and selenides (Note that T refers to a tetrahedral coordi- nation environment; ps-O refers to pseudo-octahedral coordination environment.)

Compound µ (µB ) Coordination environment Site symmetry Reference

BaFe2Se3 2.80(8) T C1 (1) (Caron et al., 2011)

Sr2F2Fe2OSe2 3.3(1) ps-O D2h (mmm) (Zhao et al., 2013)

La2O2Fe2OSe2 3.50(5) ps-O D2h (mmm) (McCabe et al., 2014b)

Ce2O2Fe2OSe2 3.33(3) ps-O D2h (mmm) (McCabe et al., 2011)

Ce2O2FeSe2 3.14(8) T D2 (222) (McCabe et al., 2014b)

CaOFeS 2.59(3) T C3ν (3m) (Jin et al., 2015)

weak ferromagnetic interchain coupling (Takano et al., a) b) 1999, 2002). For smaller lanthanides Ln3+ = Pr, Nd, a slightly different structure is formed with single chains of edge- linked CrS6 and CrO2S4 pseudo-octahedra with only the ˚ Fe Fe longer intrachain Cr - Cr distance of ∼ 3.7 A. LnCrOS2 Se Se (Ln = Pr, Nd) phases order antiferromagnetically (TN = Oe Oe 83 K for PrCrOS2,TN = 72 K for NdCrOS2) with ferro- Ce Ce ◦ magnetic coupling within chains (presumably via ∼ 100 c a Cr - S - Cr exchange), but antiferromagnetic coupling be- b c ◦ a tween chains (via ∼ 125 Cr - S - Cr exchange). Moments b are close to the [010] direction. Ln3+ moments order at c) d) fairly high temperatures in both systems suggesting that Fe there is significant interaction between the magnetic or- 3+ 3+ dering on the Cr sublattice and Ln moments (Win- c tenberger et al., 1987; Winterberg et al., 1989). To the J1 best of our knowledge, the magnetic behaviour of analo- b gous oxyselenide systems has not been fully investigated. a a

b 2. ZrCuSiAs structures FIG. 11 The magnetic structure taken from Ref. McCabe et al., 2014a. (a) and (b) show the orthorhombic and mono- Ce2O2FeSe2- Ce2O2FeSe2 orders antiferromagneti- clinic magnetic unit cells of Ce2O2FeSe2 (Ce-green, Fe=blue, cally at TN=171 K as seen by a local maximum in sus- O=red, and Se=yellow spheres). (c) and (d) displays isolated ceptibility (McCabe et al., 2011, 2014a). Neutron pow- sheets of the FeSe4 tetrahedra. der diffraction experiments indicate that the magnetic structure of Ce2O2FeSe2 (Fig. 11) consists of ferro- magnetic chains of edge-linked FeSe4 tetrahedra (Fe - Ce2O2FeSe2 (McCabe et al., 2014a; Wang et al., 2015) Fe nearest neighbour distance ∼ 2.84 A,˚ Fe - Se - Fe on cooling, similar to observations made for similar ◦ angle ∼ 71 ) with antiferromagnetic coupling between experiment in LnCrOS2 (see above), consistent with chains. These exchange interactions inferred from the two components contributing to the low-temperature magnetic structure are consistent with Goodenough rules magnetically-ordered state. An ordered moment of and bond angle analysis compiled experimentally from 1.14(1) µB on the cerium site was also found to be re- the cuprates (Mizunoa et al., 1998; Shimizu et al., 2003). quired in the magnetic powder refinement. Interestingly, Fe2+ moments are consistent with high spin d6 Fe2+ ions the onset of magnetic order for the rare earth cerium was (3.14(8) µB at 4 K) and are oriented along [010] (McCabe found in Ref. McCabe et al., 2014a to be coincident with et al., 2011). The Fe2+ magnetic structure is illustrated the iron ordering and onset at high temperatures. The in Fig. 11 (c) and (d). A discussion of the magnitude of coupling between the rare earth cerium site and the Fe2+ the exchange interaction is discussed below in the context lattice will be discussed below in the context of the lo- of neutron spectroscopy measurements. calised crystal field excitations of cerium. The combined A change in relative intensity of magnetic Bragg iron and cerium magnetic structures are shown in Fig. reflections in neutron powder data is observed for 11 (b). 15

FIG. 12 The magnetic structure of (Ce0.78La0.22)2O2MnSe2 at 30 K. (a) shows the view along the c axis and (b) illustrates the structure along the a axis. The figure is taken from Ref. Wang et al., 2015.

(Ce, La)2O2MnSe2- (Ce,La)2O2MnSe2 orders anti- axis of the electric field gradient tensor was measured ferromagnetically at TN=150 K and the magnetic struc- to be parallel to the Fe-O bond and the change in ture has been investigated using neutron powder diffrac- the electric field quadrupole splitting was found to be tion in Ref. Wang et al., 2015. (Ce,La)2O2MnSe2 has temperature independent indicating antiferromagnetic both corner- and edge-linked MnSe4 tetrahedra and, in order. The principle axis of the electric field gradient contrast to the case of Ce2O2FeSe2, both nearest neigh- is shown in Fig. 13 (a) taken from Ref. Fuwa et al., bour and next nearest neighbour exchange interactions 2010b. As noted in Fuwa et al., 2010b, high-spin Fe2+ is are antiferromagnetic (Mn - Mn nearest neighbour dis- expected to have an anisotropy in this crystalline electric tance ∼ 3.18 A˚ and Mn - Se - Mn angles of ∼ 75◦ and field environment due to the possibility of unquenched ∼ 100◦ for nearest and next nearest neighbour interac- orbital magnetic moments and also the deviation from tions). Based on magnetic neutron diffraction, the Mn2+ a perfect octahedral environment described above. moments (4.12(1) µB at 30 K) are oriented perpendic- Because of this magnetic anisotropy, combined with the ular to the layers as shown in Fig. 12. The ordered tetragonal structural symmetry, it was suggested that moment direction is also perpendicular to the Fe2+ mo- the directions between nearest neighbour Fe2+ moments ment direction in R2O2Fe2OSe2 discussed above. Above along the Fe-O bond axes are perpendicular to each TN, momentum broadened peaks in the neutron response other with the moments oriented in the a − b plane (see indicate the presence of short-range correlations and are Fig. 13 (c) and (d) for examples). possibly the reason for the transition being less obvious The exchange paths are illustrated in Fig. 5 (c) where in the susceptibility than the Fe2+ counterpart discussed three different exchange interactions are defined. Given above. Such short range correlations were not observed that antiferromagnetic nearest neighbour interactions are in Ce2O2FeSe2. frustrated, the next nearest neighbour interactions are A similar variation in relative intensity of mag- expected to be key in determining the magnetic struc- netic Bragg reflections on cooling is also observed for ture. There are two such interactions with one mediated 3+ ◦ (Ce,La)2O2MnSe2 and the onset of Ce ordering is by a 97 bond through selenium and a second through an ∼100 K. Ce3+ moments are oriented in-plane, parallel oxygen. Based on Goodenough rules the 180◦ Fe-O-Fe in- to the nearest Mn - Mn vector with ordered moments of teraction is expected to be antiferromagnetic. Therefore, 0.85(1) µB (Wang et al., 2015) at 30 K and the magnetic the next nearest neighbour interaction through selenium structure is illustrated in Fig. 12. Unlike the case of is expected to be central and the two possible magnetic the pnictide CeOMnAs (Corkett et al., 2014; Tsukamoto structures are displayed in Fig. 13 for (c) antiferromag- et al., 2011; Zhang et al., 2015), no evidence of a reorien- netic interactions and (d) ferromagnetic interactions. tation of the Mn2+ is observed at low temperatures when These two possibilities have very different conse- cerium magnetic order is present. quences for neutron diffraction studies of the magnetic structure. For antiferromagnetic alignment through the Fe-Se-Fe (Fig. 13 (c)), a single propagation wave vector 2+ 1 1 3. Ln2O2M 2OSe2 phases with M=Fe of ( 2 , 2 ) would be needed while for ferromagnetic Fe-Se- 1 Fe (Fig. 13 (d)), two propagation vectors of (0, 2 ) and 1 Nd2O2Fe2OSe2- The magnetic structure in ( 2 , 0) would be required. This is now discussed in the the 2+ Nd2O2Fe2OSe2 was initially investigated using Moss- other Fe compounds where magnetic neutron diffrac- bauer spectroscopy and compared with first principles tion studies have been completed. calculations in Ref. Fuwa et al., 2010b. Magnetic (Sr, Ba)2F2Fe2O(Se, S)2 and La2O2Fe2OSe2- ordering was observed at TN=90 K and the principal Susceptibility data collected for the oxide-fluoride- 16

FIG. 13 Proposed magnetic structures for Nd2O2Fe2OSe2 taken from Fuwa et al. (Fuwa et al., 2010b) (a) shows the principle axis of the electric field gradient(EFG) tensor. (b) illustrates the magnetic exchange pathways. (c) − (d) show the magnetic structures are proposed based upon an analysis of the electric field gradient for Fe-O-Fe antiferromagnetic and Fe-O-Fe ferromagnetic respectively. Note that in (c) the next nearest neighbour interaction mediated by selenium is anti- ferromagnetic while in (d) the interaction is ferromagnetic.

chalcogenides (Sr,Ba)2F2Fe2O(Se,S)2 (Kabbour et al., 2008) indicate similar magnetic behaviour to the oxy- selenides Ln2O2Fe2OSe2 (including Nd2O2Fe2OSe2 described above) with TN ∼ 90 K and significant deviation from Curie-Weiss behaviour at temperatures above TN, consistent with short-range magnetic corre- lations. (Mayer et al., 1992) Kabbour et al. were the FIG. 14 Thermodynamic data and diffraction results on first to investigate the magnetic order using neutron Ba2F2Fe2OSe2 and Sr2F2Fe2OS2 taken from Ref. Kabbour powder diffraction in this class of materials. Their et al., 2008. (a) is the magnetic susceptibility and (b) is the work, using reactor neutron powder diffraction data for specific heat for Ba2F2Fe2OSe2 with (c) showing a logarithmic plotting showing the divergence near the Neel temperature. Ba2F2Fe2OSe2, suggested that the magnetic structure was incommensurate with modulation wavevector ~q = (d) is the evolution of the magnetic Bragg peak intensity with temperature taken from neutron diffraction. (e) example neu- (0.42, 0.00, 0.00), see Fig. 14. (Kabbour et al., 2008) tron diffraction taken on Ba2F2Fe2OSe2 collected above and High resolution neutron powder diffraction data were below the magnetic transition. Magnetic Bragg peaks are ob- ◦ collected for La2O2Fe2OSe2 by Free and Evans and anal- servable particularly at scattering angles 2θ ∼ 25 . ysis of these data first suggested a single-k vector model with Fe2+ moments oriented in the a − b plane with collinear spins (Fig. 15 (a)). (Free and Evans, 2010) powder diffraction data, although there may be some 1 disorder of the O(2) sites within the Fe O layers. In This collinear model is consistent with the ~q = ( 2 , 0, 2 1 La O Fe OSe , this collinear model has J , J and J 0 2 ) propagation vector observed from these neutron pow- 2 2 2 2 1 2 2 der diffraction data and is very similar to the magnetic interactions all partially frustrated and so is difficult to structure adopted by Fe1+xTe (with the same modula- justify on energy grounds and seems surprising given the tion vector). However, in Fe1+xTe for small values of symmetry of the nuclear structure. x, this magnetic transition is accompanied by a tetrag- Zhao et al. found that high resolution neutron pow- onal - monoclinic structural transition (Bao et al., 2009; der diffraction data for Sr2F2Fe2OS2 were consistent with Koz et al., 2013; Rodriguez et al., 2013, 2011b; Turner a commensurate magnetic structure (Zhao et al., 2013) et al., 2009). There is no evidence for such a distor- (in contrast to earlier work on Ba2F2Fe2OSe2 mentioned tion in La2O2Fe2OSe2 from these high resolution neutron above (Kabbour et al., 2008)). Zhao et al. were able 17 to show that the 2 − k model proposed by Fuwa et the magnetic ordering in La2O2Fe2OSe2 using the lo- al. (with ferromagnetic Fe - Se - Fe exchange, Fig. cal probes of Mossbauer spectroscopy and muon spin 13) and Free and Evans’ collinear model (Fig. 15 (a)) rotation (Gunther et al., 2014). These Mossbauer are indistinguishable using neutron powder diffraction studies are consistent with those described above for 2+ data and that the 2 − k model is more appropriate for Nd2O2Fe2OSe2, (Fuwa et al., 2010b) with the Fe mo- Sr2F2Fe2OS2 (Zhao et al., 2013). They observed signifi- ments directed along the Fe - O bond axes (Gunther cant anisotropic broadening of magnetic Bragg reflections et al., 2014). Their recent muon spin rotation experi- (which was fitted by a Warren-like lineshape), consistent ments confirmed the static, long-range magnetic order with shorter range correlations along c. The in-plane on Fe2+ sites (consistent with the stacking faults sug- correlations were found to be resolution-limited, indicat- gested by neutron powder diffraction work described ing an in-plane correlation length > 300 A˚ and out-of- above (McCabe et al., 2014b; Zhao et al., 2013), rather plane correlations with ξc= 17(3) A˚ gave a good fit to than slow magnetic dynamics contributing to the unusual the data (Zhao et al., 2013). peak shapes) but also revealed a dynamic component At about the same time, the magnetic structure of due to muons affected by a large hyperfine coupling con- La2O2Fe2OSe2 was investigated using high-flux (lower stant (Gunther et al., 2014). resolution) neutron powder diffraction data. As ob- R2O2Fe2OSe2 (R = Ce, P r, Nd, and Sm)- Free served for Sr2F2Fe2OS2, the collinear and 2 − k mod- els gave equivalent fits to the low temperature neutron et al. explored the compositional flexibility of the Ln2O2M 2OSe2 structure and found that M = Fe is com- powder diffraction data for La2O2Fe2OSe2, and simi- 3+ lar anisotropic broadening of magnetic Bragg reflections patible with the widest range of Ln radii in the fluorite- (Fig. 16) (McCabe et al., 2014b) was also observed. like layers from Ln = La (eight-coordinate ionic radius 1.16 A)˚ to Ln = Sm (1.079 A)˚ (Free et al., 2011). This These were fitted using a model to simulate stacking 3+ faults (Her et al., 2007) in the magnetic structure. The series allows us to investigate the role of the Ln ion in high-resolution neutron powder diffraction data and x- the magnetic behaviour of these materials. ray powder diffraction data collected for La2O2Fe2OSe2 The first point to note is that there is no change revealed no defects in the nuclear crystal structure, indi- in the ordered magnetic structure as the Ln3+ ion is cating that these stacking faults exist only in the mag- changed (Free et al., 2011; McCabe et al., 2014a; Ni et al., netic ordering and that La2O2Fe2OSe2 and Sr2F2Fe2OS2 2011). This suggests that despite the decrease in unit cell have similar magnetic microstructures (McCabe et al., volume (and particularly the contraction within the ab 2014b). Neutron powder diffraction data were also col- plane) as Ln3+ radius decreases, the 2−k magnetic order lected for La2O2Fe2OSe2 on cooling through the mag- is robust in terms of the effects of chemical pressure (Mc- netic phase transition and showed that the onset of mag- Cabe et al., 2014a). With decreasing Ln3+ radius, there netic order had 2D-Ising like character (discussed fur- is a slight increase in TN (Table II). This is presumably ther below, Fig. 19) (McCabe et al., 2014b), consistent due to increased overlap of orbitals involved in magnetic with the results from Mossbauer spectroscopy studies exchange interactions as the unit cell (and therefore Fe - of Nd2O2Fe2OSe2 (Fuwa et al., 2010b) described above O bond lengths) decrease. and with a predominantly two-dimensional character to the magnetism (rather than resulting from nearer one- Ordering of Nd3+ moments has not been observed dimensional chains). A Warren-like peak was observed for Nd2O2Fe2OSe2 (McCabe et al., 2014a). However, 3+ in neutron powder diffraction data for La2O2Fe2OSe2 Ce moments are thought to order below 16 K in in a narrow temperature range (∼13 K) above TN, in- Ce2O2Fe2OSe2, with a similar in-plane arrangement to dicating short-range magnetic ordering within the Fe2O the Fe moments but with some out-of-plane canting, per- layers immediately above the three-dimensional ordering haps indicating some coupling between Fe2+ and Ce3+ temperature. This is on contrast to La2O2Mn2OSe2 for sublattices (McCabe et al., 2014a). The behaviour of which a similar Warren peak is observed up to 140 K the Pr2O2M 2OSe2 analogues is unusual and the low above TN (Ni et al., 2010). This indicates a low degree temperature tetragonal - orthorhombic distortion has of frustration in the La2O2Fe2OSe2 magnetic structure, been discussed above. Analogous distortions in related consistent with the 2 − k model and not the collinear PrMnSbO (Kimber et al., 2010), PrMnAsO (Wildman model. This conclusion was consistent with neutron in- et al., 2015) and PrFeAsO (Kimber et al., 2008) are ac- elastic scattering data: simulated inelastic spectra for companied by long-range ordering of Pr3+ moments. Ni various models were obtained and were comparable with et al. observed a peak in heat capacity data at ∼ 23 experimentally observed spectra (Figure 25) but only the K for Pr2O2Fe2OSe2 and slight changes in neutron pow- exchange constants determined for the 2 − k model are der diffraction data at low temperature suggesting that consistent with the magnetism observed. (Fuwa et al., Pr3+ moments may order but this is not yet fully under- 2010b; Kabbour et al., 2008; McCabe et al., 2014b) stood (Ni et al., 2011). Ordering of Sm3+ moments is Most recently, Gunther et al. have investigated also thought to occur below 6 K (Ni et al., 2011). 18

ĂͿ ďͿ K K &Ğ &Ğ

^Ğ ^Ğ

ď ď

Ă Ă ĐͿ

FIG. 15 Collinear model showing in-plane spin arrange of the (a) full collinear magnetic structure; (b) 2-k model and (c) illutrating the in-plane spin arrangement viewed as two interpenetrating square lattices.

TABLE II Summary of key structural and magnetic ordering behaviour in Fe2O materials; unit cell parameters and Fe-O bond lengths are from references listed and TN is from diffraction results except for Na2Fe2OSe2 and Sm2O2Fe2OSe2 for which TN ∗∗ is from mangetic susceptibility measurements. indicates TN extracted from heat capacity measurements.

3+ Compound Ln ionic radius (A)˚ a (A)˚ (295-300 K) dF e−O (A)˚ TN (K) µ (µB ) Ref.

Na2Fe2OSe2 - 4.107(8) 2.054(8) 75 - He et al., 2011

La2O2Fe2OSe2 1.16 4.084466(9) 2.042233(9) 89.50(3) 3.50(5) Free and Evans, 2010 McCabe et al., 2014b

Ce2O2Fe2OSe2 1.143 4.06134(5) 2.03067(5) 92.3(2) 3.32(1) McCabe et al., 2014a ∗∗ Pr2O2Fe2OSe2 1.126 4.0447(1) 2.0224(1) 92.09(2); 88.6 3.36(2) Oogarah, 2016 Ni et al., 2011

Sm2O2Fe2OSe2 1.079 3.9976(1) 1.9988(1) 85.3 - Ni et al., 2011

2+ 4. Ln2O2M 2OSe2 phases with M=Co possible structures noted in Free et al., 2011 and Fuwa et al., 2010a are both consistent with the next nearest neighbour Co-Se-Co (bond angle ∼ 99◦) and Co-O-Co ◦ La2O2Co2OSe2- The magnetic structure for (bond angle 180 ) interactions being antiferromagnetic. La2O2Co2OSe2 has been reported in Ref. Fuwa et al., This differs from La2O2Fe2OSe2 where, as noted above, 2010a with order parameter shown in Fig. 17 (TN the Fe-Se-Fe interaction is ferromagnetic. =217 K). The magnetic structure for La2O2Co2OSe2 is different from Fe2+ analogue discussed above in that While momentum broadened scattering in the diffrac- it is determined by a single propagation wave vector tion data was observed in the range of 225-250 K, sharp 1 1 of ( 2 , 2 , 0). As noted in Ref. Free et al., 2011, an Bragg peaks were present at lower temperatures. The ambiguity exists in the magnetic structure from powder high temperature momentum broadened scattering was diffraction data with magnetic moments directed either considered to originate from magnetic diffuse scattering. along or perpendicular to the Co-O bonds. The two The momentum broadened peaks at high temperatures 19

La2O2Fe2OSe2 84.20 % 80 R = 6.976% Structure 15.80 % a) wp ;&ƵůůƉĂƩĞƌŶͿ 60 28 independent parameters

40

20

0

10 15 20 25 30 35 40 45 50 55 60 65 2ɽ / °

La2O2Fe2OSe2 84.20 % b) Structure 15.80 % 8 Rwp = 6.976%;EŽĚŽŵĂŝŶƐͿ 28 independent parameters 6

4

Intensity 2

(arbitrary units) 0

La2O2Fe2OSe2 84.20 % c) magnetic 15.80 % 8 Rwp = 5.375%;ŶƟƉŚĂƐĞĚŽŵĂŝŶƐͿ 29 independent parameters 6

4

2

0

504846444240383634323028262422201816 2ɽ / °

FIG. 16 Neutron diffraction data on La2O2Fe2OSe2 taken from Ref. McCabe et al., 2014b. (a) shows Rietveld refin- ments with the 2−k model showing both nuclear (blue arrows) and magnetic (black tick marks) phases. (b) shows the refine- ment with same peak shape for both nuclear and magnetic phases. (c) shows refinement including antiphase boundaries FIG. 17 Neutron diffraction and magnetic structure for in the magnetic phase. Note the increase in the quality of the La2O2Co2OSe2 reported in Ref. Fuwa et al., 2010a. (a) neu- fit with the inclusion of just one additional free parameter. tron diffraction patter taken at 10 K. The upper and lower Note that the tick marks do not include a refined zero offset ◦ sets of vertical markers in the pattern are the calculated nu- of ∼ 0.4 . clear and magnetic peak positions, respectively. (b) shows possible magnetic structure model for the cobalt moments. (c) is the temperature dependence of the lattice constants. occur over a similar temperature range where strong de- (d) is a plot of the ordered magnetic moment as a function of viation from Curie-Weiss behaviour is observed in suscep- temperature. tibility (Fuwa et al., 2010c). Below TN, no anisotropic lineshape to the magnetic Bragg peaks was reported in contrast to analogous La2O2Fe2OSe2 and Sr2F2Fe2OS2. 2+ similarities to that reported in PrOMnSb (Kimber et al., The magnetic moment for Co was measured to be 2010). 3.53 ± 0.01 µ in Ref. Fuwa et al., 2010a and 3.29(3) µ B B The magnetic moment was reported to be 4.147± 0.028 in Ref. Free et al., 2011. It was noted in Ref. Fuwa et al., µ (Ref. Ni et al., 2010) and 4.5 ± 0.3 µ (Ref. Free 2010a that this was much larger than the theoretical first B B et al., 2011). This value is in excellent agreement with the principle calculations that predicted 2.70 µ (Wu, 2010) B magnetic moments reported in BaMn P (4.2(1) µ ) and leading to the suggestion of a possible orbital contribu- 2 2 B BaMn As (3.88(4) µ ) (Brock et al., 1994; Singh et al., tion to the magnetic moment. 2 2 B 2009). They are also in agreement with MnO (4.892 µB at 10 K) and Mn2SiSe4 (4.36 µB at 2 K) (Bodenan et al., 2+ 1996; Nonfante et al., 1972). 5. Ln2O2M 2OSe2 phases with M=Mn Similar to the case of the iron based oxyselenides La2O2Mn2OSe2- The magnetic and structural prop- discussed above, there are three magnetic interac- 2+ erties of La2O2Mn2OSe2 were investigated and reported tions between Mn ions that need to be considered. in Free et al., 2011 and Ni et al., 2010. The nuclear and These include the superexchange interaction through magnetic structure (from Ref. Ni et al., 2010) refined the 180◦ Mn-O-Mn pathway, the superexchange inter- from neutron powder data is shown in Fig. 18 with long- action through the ∼ 95◦ Mn-Se-Mn channel, and the 2+ range magnetic order being reported below TN=163 K interaction between the nearest Mn ions. These are (Ref. Ni et al., 2010) and 168.1 K (Ref. Free et al., denoted as J1, J2, and J20 in Fig. 18. Based on 2011) with magnetic propagation vector ~q=(0,0,0). The Goodenough rules, the 180◦ Mn-O-Mn exchange is ex- ◦ magnetic structure of La2O2Mn2OSe2 is very different pected to be antiferromagnetic while the ∼ 95 Mn-Se- to those of the Fe2+ and the Co2+ analogues: the Mn2+ Mn coupling is expected to be ferromagnetic. Unlike spins are oriented perpendicular to the Mn2O (along c) La2O2Fe2OSe2 discussed above, the nearest-neighbour planes with different relative spin arrangements. The antiferromagnetic J1 exhange interactions dominate in magnetic structure of La2O2Mn2OSe2 does have strong La2O2Mn2OSe2, resulting in a G-type antiferromagnetic 20

netic structures of Mn2+ and Co2+ materials (Free et al., 2011; Fuwa et al., 2010a,c; Ni et al., 2010), compared with the Fe2+ analogues (McCabe et al., 2014a,b; Zhao et al., 2013). R2O2Mn2OSe2 (R = Ce, P r, Nd, and Sm) - The magnetic properties of two dimensional Mn2+ based oxy- selenides have been investigated for a series of lanthanides and reported in Ref. Free et al., 2011. There is no re-
N 3+ M
ported changed in the magnetic structure with Ln ion and the results reported in Table III show no observable
change of the ordered magnetic moment with rare earth N@ 2+ substitution. However, a decrease in Mn TN in pro- portion to the rare earth radius is observed. To allow a direct comparison with the Fe2+ materials discussed above, bond distance information is reproduced from Ref. Free et al., 2011 in Table. III. The manganese ordering temperature for two dimensional oxyslenides in- creases with decreasing rare-earth radius and has been as- FIG. 18 The magnetic structure of La2O2Mn2OSe2 taken from Ref. Ni et al., 2010. (a) illustrates the chemical sociated with increased overlap of the orbitals involved in and magnetic structure and (b) shows a view of a sin- the magnetic coupling. Finally, no evidence for rare earth 2− gle [Mn2OSe2] with the three exchange interactions high- ordering in either Pr2O2Mn2OSe2 or Ce2O2Mn2OSe2 has lighted. The labelling of the exchange interactions have been been reported. edited to be consistent with the text of this review.

6. Mixed Fe/Mn-based oxyselenides structure with nearest-neighbour Mn2+ spins antiparal- lel and next-nearest-neighbour spins parallel. (Ni et al., Nd O (Fe Mn ) OSe - Several groups have ◦ 2 2 1−x x 2 2 2010) This leaves the 180 Mn-O-Mn J20 interactions investigated the magnetic behaviour of compositions frustrated. within the Ln2O2Fe2−xMnxOSe2 solid solution and mag- Susceptibility measurements on La2O2Mn2OSe2 sug- netic susceptibility measurements suggest some ferro- gest (Free et al., 2011; Ni et al., 2010). Indications of magnetic behaviour on cooling (Landsgesell et al., 2013; high-temperature short-range correlations are indicated Lei et al., 2012; Liu et al., 2015a). We note that by a deviation from Curie Weiss behaviour. At low tem- this has also been observed for powder samples of peratures the published data shows several features and La2O2Mn2OSe2 (Free et al., 2011; Ni et al., 2010) a marked difference between field and zero field cooled whilst single crystal studies suggested only antiferromag- responses. However, as noted in Ref. Free et al., 2011, netic ordering on cooling (Liu et al., 2011). Lands- it is likely that the features are due to trace quantities gesell et al. have investigated the magnetic ordering of Mn3O4 which has an ordering temperature of TN=42 in La2O2MnFeOSe2 with a disordered arrangement of 2+ 2+ K (Regmi et al., 2009; Seo et al., 2004). Features near 140 Mn and Fe ions within the M2O layers. Their re- K in the magnetic susceptibility can also be explained sults from neutron powder diffraction experiments indi- ~ by trance amounts of LaMnO3 (Ghivelder et al., 1999; cate that La2O2MnFeOSe2 orders with k = (0 0 0) (as for Huang et al., 1997). La2O2Mn2OSe2) but with moments within the ab plane, Evidence for short range, two-dimensional magnetic although the exact spin arrangement has yet to be con- order in La2O2Mn2OSe2 was observed from neutron pow- firmed (Landsgesell et al., 2013). der diffraction, where a Warren peak, indicative of two dimensional ordering (Ni et al., 2010), was observed over a wide temperature range. This peak diminished on cool- 7. Critical scattering ing and was unobservable at 6 K. At 100 K, the out-of- plane correlation length was extracted to be 60 A˚ (Ni The temperature dependence of the magnetic order et al., 2010). Broadening of magnetic Bragg reflections parameter can provide helpful information on the uni- (as observed for Ln2O2Fe2OSe2 (Ln = La, Ce, Nd) (Mc- versality class and also the dimensionality of the mag- Cabe et al., 2014a,b) and Sr2F2Fe2OS2 phases (Zhao netism (Collins, 1989). The critical properties classify et al., 2013)) is not observed for La2O2Mn2OSe2 (Free the phase transition and allow commonalities to be es- et al., 2011; Ni et al., 2010). This suggests more simi- tablished with other systems. To this end, the magnetic lar magnetic correlation lengths within, and perpendic- order parameter has been studied in a number of mate- ular to the M2O planes in the three-dimensional mag- rials and relies on the fitting of a critical exponent to the 21

TABLE III Summary of key structural and magnetic ordering behaviour in Mn2O materials; unit cell parameters and Mn-O bond lengths are from references listed and TN is from diffraction results. Data taken from Ref. Free et al., 2011.

3+ Compound Ln ionic radius (A)˚ a (A)˚ (295-300 K) dMn−O (A)˚ TN (K) µ (µB )

La2O2Mn2OSe2 1.16 4.138921(4) 2.06435(1) 168.1(4) 4.5(2)

Ce2O2Mn2OSe2 1.143 4.11304(2) 2.05124(1) 174.1(2) 4.8(3)

Pr2O2Mn2OSe2 1.126 4.09739(2) 2.04308(1) 180.3(4) 4.5(1)

) oxyselenides have exponents ∼ 0.1. The critical expo- nents for the 2D Ising universality class is β= 0.125 and 2D XY is 0.13. 3D Ising has an exponent of 0.326 and 3D Heisenberg is 0.36 (Collins, 1989). The two dimen- sional Fe2+ based oxyselenides clearly display 2D charac- ter and the exponents are similar to ReFeAsO materials studied in Ref. Wilson et al., 2010 with β ∼ 0.125, close to the ideal 2D Ising universality class. The exponents are also similar to Fe1+xTe which, for large values of in- terstitial iron, displays β=0.15 (Rodriguez et al., 2013) and FeAs with β=0.16(2) (Rodriguez et al., 2011b). The manganese compounds and structurally one dimensional Ce2O2FeSe2 are closer to the 3D limit and are similar to to the critical exponents in (Ba,Sr)Fe2As2 (Christianson et al., 2009; Wilson et al., 2009) and also more recently FIG. 19 The magnetic order parameters for La2O2Mn2OSe2 extracted in SrMn2As2 (Das et al., 2016). (blue) and La2O2Fe2OSe2 (red) from Ref. McCabe et al., Chen et al., 2009 has considered the dimensionality of 2014b. The data was obtained from the magnetic intensities the order parameters in two dimensional pnictide iron derived from magnetic neutron powder diffraction. based systems in terms of coupling to an orbital degree of freedom. This orbital order parameter possess 2D- T β Ising character and results from a spin-orbital Hamilto- magnetisation M(T ) = M0(1 − ) . A summary of re- TN sults obtained from several groups using neutron diffrac- nian. This idea is broadly consistent with the response discussed above in the case of the magnetic oxyselenides. tion to measure the magnetisation is presented in Table 2+ IV. Fe based materials with 6 d electrons potentially have an orbital degree of freedom as discussed above while Mn2+ has only 5 d electrons and therefore no orbital TABLE IV Critical exponents for the magnetic order param- component in the weak/intermediate crystal field limit. eter extracted from magnetic neutron diffraction Based on the orbital model proposed in Ref. Chen et al., 2+ Compound β Reference 2009, it is therefore expected that Fe would display a stronger 2D-Ising character as illustrated in the data in La O Mn OSe 0.24(7) Free et al., 2011 2 2 2 2 Table IV. Ce2O2Mn2OSe2 0.29(7) Free et al., 2011

Pr2O2Mn2OSe2 0.27(4) Free et al., 2011 8. Summary of the magnetic structure variation with transition Ce2O2FeSe2 0.28(1) McCabe et al., 2014a metal ion

La2O2Fe2OSe2 0.122(1) McCabe et al., 2014b While the magnetic structures of the Ln2O2M2OSe2

Ce2O2Fe2OSe2 0.11(1) McCabe et al., 2014b do not vary drastically with changing Ln ion, there are dramatic changes with transition metal ion (M) substitu- Ba2F2Fe2OSe2 0.118 Kabbour et al., 2008 tion. While we have discussed these structures in depth Sr2F2Fe2OS2 0.15 Kabbour et al., 2008 above, here we summarise the three magnetic structures observed for M=Fe2+, Co2+, and Mn2+ in Fig. 20. It should also be noted that there is a large difference in the The critical exponents can be seen to fall into two critical properties between M=Fe2+ and Mn2+ as noted broad categories with the manganese variants having ex- above. ponents ∼ 0.2-0.3 and the iron based two dimensional The magnetic structures adopted by the 22

perconductors. In this section we discuss the magnetic interactions with first a review of neutron inelastic work  on two dimensional Fe-based oxyselenides and then a
 short summary of the results found for other materials 
 
   discussed above in the context of the static structure. A summary of the magnetic interactions is presented in Fig. 5 (c). The nearest neighbour M - M dis- tance within the M2O planes ranges from ∼ 3.02 A˚ (in Ba2F2OMn2Se2) to ∼ 2.83 A(in˚ Sm2O2Fe2OSe2). The nearest neighbour J1 exchange could be direct, or could proceed via 90◦ M -O- M superexchange or ∼ 64◦ M - Se - M superexchange, and is expected to be antiferro- 
 
 magnetic for M = Mn, Fe and Co (Kabbour et al., 2008; Koo and Whangbo, 2012; Wang et al., 2010; Wu, 2010; Zhao et al., 2013; Zhu et al., 2010). The next nearest ◦ neighbour J20 exchange interaction (180 superexchange via M - O - M) is also expected to be antiferromagnetic from Goodenough-Kanamori rules for M=Mn, Fe and ◦ Co. The next nearest neighbour J2 exchange (∼ 97 M - Se - M superexchange) is expected to be ferromagnetic for M = Fe (Kabbour et al., 2008; Zhao et al., 2013; 
  Zhu et al., 2010) and Co (Wang et al., 2010; Wu, 2010), despite the apparent antiferromagnetic interaction from the magnetic structure. The case of M=Mn is different for the J2 exchange (∼ 97◦ M - Se - M superexchange) with density functional calculations differing in terms of the sign of J2 (Koo and Whangbo, 2012; Liu et al., 2011). Antiferromagnetic or- FIG. 20 A summary of the reported magnetic structure for dering of (Ce,La)2O2MnSe2 described above (with M - ◦ ◦ the Ln2O2M2OSe2 series of compounds with the transition Se - M angles of ∼ 75 and ∼ 100 ) might suggest an- 2+ 2+ metal ion M= (a) Fe ,(b) Co (note the ambiguity of mag- tiferromagnetic J2 in La2O2Mn2OSe2. A key factor for netic structures and the two possible results) and (c) Mn2+. the relative strength of the three exchange interactions in Fig. 5 (c) with M is the electronegativity: for M = Mn, nearest neighbour J1 interactions dominate and J20 Ln2O2M2OSe2 materials are governed by the ex- interactions are frustrated while for the more electroneg- change interactions within the M 2O planes. We now ative M = Co, J20 interactions dominate at the expense discuss the magnetic interactions governing these of nearest neighbour J1 interactions. structures discussed above.

1. La2O2Fe2OSe2 and Ce2O2FeSe2 E. Magnetic interactions Given the semiconducting nature of the oxyselenides, The sign and relative strength of the magnetic inter- a localised model for the spin interactions is appropriate actions can be postulated based upon magnetic diffrac- and therefore the dominant term in the magnetic Hamil- P ~ ~ tion data establishing the magnetic structure of the ma- tonian that needs to be considered is H = J i,j Si · Sj, terial. These can also be compared with Goodenough where the sum is performed over nearest neighbours. rules and expectations on systematically characterised This model is much more applicable to the oxyselenides systems like the cuprates. A definitive measure of the over the cuprates or iron based superconductors which coupling strength is obtained through spectroscopy and derive from metallic ground states and hence display neutron inelastic spectroscopy is ideal given the sensitiv- strong evidence of coupling between electronic and mag- ity to magnetic moments and the energy scale that is netic moments. This is particularly evident in the high sensitive to. However, few low-energy spectroscopy mea- energy neutron scattering response in the cuprates (Stock surements have been reported with the focus on the iron et al., 2007, 2010a) and iron based systems (Stock et al., based two dimensional oxyselenides given their possible 2014). interesting Mott insulating behaviour and also the close One of the key questions that arise from the mag- analogy to the iron based pnictide and chalcogenide su- netic diffraction data is how to stabilise the orthogonal 23

2 − k magnetic structure reported for La2O2Fe2OSe2. As noted in Ref. McCabe et al., 2014b, second order terms in the spin involving either antisymmetric (such as Dzyaloshinskii-Moriya interactions ∼ D~ · (S~i × S~j)) or symmetric (such as Heisenberg) interactions are not able to stabilise the 2−k structure for these tetragonal crystal structures. Other terms that may be relevant to the magnetic Hamiltonian and discussed in the literature include the biquadratic spin-spin interactions. These terms have the P ~ ~ 2 form H1 = −K i,j(Si · Sj) and are required to under- stand (Stanek et al., 2011; Wysocki et al., 2011; Yu et al., 2012) spin excitations in the pncitides near the Bril- louin zone boundary (Harriger et al., 2011; Zhao et al., 2009). Without consideration of this term in the mag- netic Hamiltonian, anisotropic exchange terms need to be considered which are difficult to reconcile given the tetragonal nuclear structure. Studies of the magnetic structure are not able to uniquely determine which term FIG. 21 Neutron inelastic data and calculations using the is present in the Hamiltonian and therefore neutron in- first-moment sum rule combined with the single mode analysis elastic scattering is required to obtain a better under- on Ce2O2FeSe2 . The figure is taken from Ref. McCabe et al., standing of the magnetic interactions. 2014a. (a − c) shows temperature dependent data illustrating Due to lack of large single crystals, complete neutron a cerium crystal field peak near ∼ 12 meV and low-energy scattering data is currently quite sparse for the oxyse- ferromagnetic fluctuations near Q=0. (d) shows an estimate lenide materials. While powders provide limited informa- of the iron contribution subtracting off the cerium crystal field tion and, most importantly, are not able to determine the using high-angle data. (e − h) show calculations from which is concluded that the Fe-Fe exchange is ferromagnetic and ∼ dispersion near the zone boundary, they can provide help- 25 meV in magnitude. Note the sign convention that negative ful information on the integrated intensity (a rough mea- (-) is antiferromagnetic and positive (+) is ferromagnetic. sure of the overall exchange constant), anisotropy gap (directly related to the local crystalline electric field en- vironment), and also the integrated spectral weight which helps in understanding the spin state. The coupling between the rare earth site and the iron Inelastic data on powders has been obtained on pow- site were also studied by investigating the response of the 3+ ders of La O Fe OSe and Ce O FeSe allowing the Ce crystal field excitations to Neel ordering on the iron 2 2 2 2 2 2 2 3+ 5 magnetic interactions to be investigated with increas- site. Ce can be assigned a J = 2 and the crystal field ing levels of structural complexity. We first outline the scheme consists of three doublets. Kramers theorem en- sures that the degeneracy of these doublets is not split results for the structurally one-dimensional Ce2O2FeSe2 and then discuss this in the context of two-dimensional unless there is a field which breaks time reversal symme- try such as a magnetic field. Crystalline electric fields La2O2Fe2OSe2. The powder averaged magnetic excitations for a sam- will not split the doublet degeneracy alone. Since no splitting of the Ce3+ crystal field doublets was observed ple of Ce2O2FeSe2 is shown in Fig. 21 with data taken from the MARI direct geometry spectrometer at ISIS. at low temperatures, it was concluded in Ref. McCabe 3+ 2+ et al., 2014a that the coupling between the iron and rare Ce2O2FeSe2 has two magnetic sites (Ce and Fe ) which complicates the neutron excitation spectrum as it earth sites is weak. This is in contrast to rare earth sub- consists of both crystal field excitations from the Ce3+ stituted pnictides where a strong coupling is observed be- sites and also collective excitations for S=2 Fe2+ mo- tween the iron and rare earth sites as demonstrated by a ments. However, the crystal field excitations are mo- splitting of the crystal field doublets at temperatures be- mentum independent while the collective excitations of low the iron ordering (Chi et al., 2008). High resolution Fe2+ moments are comparatively localised in momentum. neutron spectroscopy studies even observe a dispersion This difference was used in Ref. McCabe et al., 2014a to of the crystal field excitations implying coupling between subtract the crystal field contribution from the powder the rare-earth sites (Li et al., 2014). No such effects have been reported in rare earth substituted oxyselenides. average neutron inelastic spectrum in Ce2O2FeSe2. The results are shown in Fig. 21 where the remaining spectral An estimate of the Fe-Se-Fe exchange constant was weight after subtraction is concentrated near Q=0, indi- obtained by modelling the powder average neutron spec- cating ferromagnetic interactions consistent with neutron trum using the first moment sum rule combined with the magnetic diffraction. single mode approximation outlined above. Given the 24

It is interesting to compare the magnetic anisotropy gap and its temperature dependence to Fe1+xTe. The o magnetic anisotropy is very similar to that measured in single crystals of Fe1+xTe where commensurate ~q = 1 1 ( 2 , 0, 2 ) is observed (Stock et al., 2011). As shown in Ref. Rodriguez et al., 2013, Fe1+xTe for small values of x undergoes a metal to “semi metal” transition char- acterised by a sharp response and change in slope in the resistivity. This change in slope is coincident with a gapping of the magnetic fluctuations and it was pos- tulated that this temperature dependence in the gapped spin fluctuations was responsible for the “metallic-like” behavior in Fe1+xTe (Rodriguez et al., 2013). As noted above, the resistivity from spin fluctuations can be cal- culated from S(Q,~ ω) and the electronic scattering from these fluctuations and the removal of low-energy decay channels was found to explain the change in resistivity. While spin excitation gap magnitude, and also the order- ing wavevector, is similar in Fe1+xTe and La2O2Fe2OSe2, o no such metal-“semimetal” transition has been reported in La2O2Fe2OSe2 despite the spin fluctuations showing a qualitatively similar temperature dependence across TN. While the magnetic excitations are gapped in the mag-

o o netically ordered Neel state at low temperatures, they are gapless at high temperatures. This may explain the

FIG. 22 (a)−(c) Powder averaged spectra for La2O2Fe2OSe2 origin of the observation of momentum broadened dif- measured on DCS. (d) Momentum-integrated energy scan at fuse scattering in La2O2Co2OSe2 given neutron diffrac- 2 K (upper) and 150 K (lower); the curves are calculations tion measurements are typically done in two-axis mode using a single-mode analysis with a 1D model, a 2D model, and hence energy integrating. We emphasise, though, and a 3D model. (f) − (h) Plots of the powder-averaged that no inelastic data has been reported for this particu- temperature spectra taken on the MARI spectrometer. The figure is taken from McCabe et al. (McCabe et al., 2014b) lar Co compound and this observation is speculative with currently available data. The intensity distribution at the edge of the low tem- perature excitation gap is sensitive to the dimensionality kinematics of neutron spectroscopy which masks lower of the magnetic interactions. The dimensionality of the momentum transfers at higher energies, the errorbar on interactions is established in Fig. 22 (d) where a first- this analysis is large however an estimate can be ob- moment sum rule analysis suggests that the interactions tained. As shown in Fig. 21 J=-25 meV provides a in La2O2Fe2OSe2 are two dimensional. The figure shows reasonable description of the results. Note in Ref. Mc- the momentum integrated data compared against calcu- Cabe et al., 2014a, positive (+) exchange indicates fer- lations based on the single-mode approximation for an romagnetic exchange while usually (-) is often taken and isotropic dispersion in one-dimensional (1D) chain, 2D is the convention used here. plane, or 3D structure. The 2D model gives the best de- Having analysed the chain compound and established scription consistent with the 2D-Ising critical properties the ferromagnetic exchange, we now discuss work done on discussed above. two dimensional La2O2Fe2OSe2. A summary of the mag- Scans that probe larger energy transfers are shown in netic interactions and the bond angles governing them Figs. 22 (f) − (h) where it is shown that the magnetic in this compound is illustrated in Fig. 5 (c). The pow- excitations extend up to energy transfers of ∼ 25 meV. der average inelastic spectrum for La2O2Fe2OSe2 is illus- This small band width accounts for all of the expected trated in Fig. 22 at several temperatures. Panels (a)-(c) spectral weight, confirmed by integrating the intensity illustrate the low-energy part of the magnetic structure in momentum and energy and comparing against the ze- showing a gap of ∼ 6 meV which softens with increasing roth sum rule discussed above. The total integral of both temperature. The presence of a gap is also confirmed by elastic and inelastic spectral weight was found to give an NMR studies in Ref. Gunther et al., 2014 where the re- integral of 5.9 (4) which is close to the S=2 value of 6. 2 −∆/T laxation rate was fitted to (1/T1) ∝ T e and a gap To give an estimate of the exchange coupling constants, value of ∆=55 K was extracted which is close to that Ref. McCabe et al., 2014b compared the data against measured with neutron inelastic scattering. calculations with a large single-ion anisotropy to fix the 25

parameters derived for the collinear model. Therefore, through a combination of neutron diffraction and spec- troscopy, Ref. McCabe et al., 2014b concluded that the magnetic structure of La2O2Fe2OSe2 is the 2 − k struc- ture. It is interesting to note that while consistency is obtained in the sign of the exchange constants be- tween Ce2O2FeSe2 and La2O2Fe2OSe2, the single mode analysis suggests that the ferromagnetic exchange in Ce2O2FeSe2 is much larger than La2O2Fe2OSe2. This can be attributed to the different local bond environ- 2+ ment in both materials. In Ce2O2FeSe2, the Fe ion is in a local tetrahedral environment while La2O2Fe2OSe2 is more square planar or pseudo-octahedral. In this con- 2+ text the Fe site La2O2Fe2OSe2 is quite different than iron based pnictide and chalcogenide systems where the iron is in a tetrahedral framework. The small exchange constants in La2O2Fe2OSe2 de- rived from this heuristic model and the bandwidth of the magnetic excitations are remarkable in the context of observations in the cuprates and also iron based pnic- tides. Mott insulating La2CuO4 (Coldea et al., 2001) and YBa2Cu3O6+x (Hayden et al., 1996) both have band- widths of over 300 meV and the parent phases of the pnictides have magnetic excitations that extend up to

o about ∼ 100 meV for BaFe2As2 and ∼ 150 meV in CaFe2As2 (Dai, 2015; Zhao et al., 2009). The excita- FIG. 23 Low temperature neutron inelastic scattering for tions in Fe1+xTe chalcogenides extend up to ∼ 150-200 La2O2Fe2OSe2 compared against calculations based on a meV and the high energy excitations account for a large heuristic spin-wave model. The results are taken from Ref. fraction of the total spectral weight (Stock et al., 2014). McCabe et al., 2014b. (a) shows a inelastic neutron scat- Therefore, while it is tantalising to make a connection be- tering data compared against models supporting (b) 2 − k tween the two dimensional oxyselenides to the cuprates structure and (c) collinear structure. The effect of weak an- and iron based superconductors owing to the qualita- tiferromagnetic and ferromagnetic values on the spectra are tively similar electronic phenomena, the magnetic excita- shown in panels (d) and (e). tions are very different with the oxyselenides displaying significantly smaller coupling. moment direction and considering only Heisenberg spin 2+ 2+ exchange. This model is somewhat artificial as it uses 2. Ln2O2M 2OSe2 for M=Mn or Co anisotropy terms in the magnetic Hamiltonian to fix the moment direction to allow consistency with possible mag- At the time of writing this review, there has been no netic structures found from magnetic diffraction data. As reports of neutron inelastic scattering data on Mn2+ or shown in Fig. 23, the experimental spectrum can be re- Co2+ analogues of the two dimensional oxyselenides dis- produced reasonably well for the 2 − k ground state with cussed above. Future work study the fluctuation spec- J1=0.75 meV, J2=-0.10 meV, and J20 =1.0 meV. Con- trum in these materials will be useful in comparison to sistent models could be obtained for a collinear model, the work presented above in the context of the iron based however these required an antfierromagnetic J2 which oxyselenides. is inconsistent with first principles calculations and also Goodenough rules. Perhaps more conclusively from an experimental perspective, the requirement of antiferro- V. ELECTRONIC PROPERTIES magnetic J2 in the collinear model is surprising given the neutron spectroscopy data reviewed above on the We now discuss the electronic properties of oxyse- chain compound Ce2O2FeSe2. Also, comparisons with lenides in terms of resistivity, optical measurements, x- magnetic high temperature susceptibility data find bet- ray spectroscopy, and calculations. The ZrCuSiAs ma- ter agreement for the Weiss temperature with the param- terials have been investigated in the context of thermo- eters derived from 2 − k model than the corresponding electric properties while interest in compounds based on 26

values). These calculations indicate that the hole carri- 2+ ers are confined to the [Cu2Se2] layers by the insulat- 2+ ing [La2O2] layers, giving significant two-dimensional character, consistent with features at the absorption edge of the material (Ueda et al., 2004a). More recently, den- sity functional theory calculations have investigated the origin of the p-type semiconductivity and indicate that although this can be induced by aliovalent doping, Cu+ vacancies are easily formed and are likely to be the domi- nant acceptor defect in samples (Hiramatsu et al., 2010b; Scanlon et al., 2014). With the band gap dominated by 2+ the [Cu2Se2] layers, similar optical and electronic prop- erties are observed for LnCuOS (Ln = La, Pr, Nd), with a slight decrease in band gap with decreasing Ln3+ ionic radius (Ueda et al., 2003).

TABLE V Estimated band gaps from optical measurements taken from Ref. Hiramatsu et al., 2008. Compound ∆(eV ) LaCuOS 3.1 BiCuOS 1.1 FIG. 24 Schematic band structures of (a) BiCuOS and (b) LaCuOS reproduced from Ref. Hiramatsu et al., 2008. BiCuOSe 0.8 LaCuOSe 2.8 LaCuOTe 2.4 magnetic Fe2+ have been pursued in the context of un- conventional electronic properties in relation to super- conductivity. Bi3+ ions can also be accommodated in the fluorite- like oxide layers and give a dramatic change in proper- ties. BiCuOSe has a much higher electron conductiv- A. ZrCuSiAs structures and related phases ity than LaCuOSe and a smaller band gap (∼0.8 eV with absorption in the near-infrared region) and simi- The optical transparency and p-type semiconducting lar behaviour is observed for other BiCuOQ phases (Hi- properties of LaCuOQ materials have prompted several ramatsu et al., 2008). Although electron conductivity experimental and theoretical studies to understand their was found to be different for BiCuOSe and LaCuOSe, electronic structure. Diffuse reflectance spectra measured they have similar hole conductivities, implying that dif- for LaCuOSe indicate a band-gap of ∼2.8 eV, (Ueda ference in the band structure can be attributed to the and Hosono, 2002) although optical absorption measure- conduction band. While the Bi3+ 6s states are 2 - 5 eV ments reveal sub-band gap absorptions (Hiramatsu et al., below the Fermi energy, the 6p states form the bottom 2010a). LaCuOQ (Q = S, Se) phases show high con- of the conduction band, deepening the conduction band ductivities and Seebeck and Hall measurements confirm and decreasing the band gap (Fig. 24) (Hiramatsu et al., that they behave as p-type semiconductors (Ueda and 2008; Zou et al., 2013). The low thermal conductivity Hosono, 2002; Ueda et al., 2003). Hole-doping (for ex- of BiCuOQ phases, coupled with their semiconducting 2+ 2+ ample La1−xAxCuOSe; A = Sr , Mg ) can further behaviour makes them promising thermoelectric materi- enhance this conductivity (Hiramatsu et al., 2007; Ueda als (Zhao et al., 2014). and Hosono, 2002). Related cation-ordered La2O2CdSe2 has an even larger Energy band calculations show that the wide band gap band gap ( 3.3 eV from diffuse reflectance measurements) 2+ arises primarily from the [Cu2Se2] layers with the va- and high electrical resistivity, and attempts to induce lence band maximum composed of antibonding Cu 3d semiconducting behaviour by aliovalent doping were un- and Se 4p states, whilst Cu 4s states make up the conduc- successful (Hiramatsu et al., 2004a). Density functional tion band minimum (see Fig. 24) (Ueda et al., 2004a,b). theory calculations suggest that the valence band is pre- The connection of the band gap with Se states is further dominantly composed of Se 4p states, similar to that of confirmed by investigations as a function of substitut- LaCuOSe, whilst the conduction band is composed of ing S for Se which showed a large change in the band Cd 5s states and is much narrower than that of LaCu- gap with doping (see Table V below for pure compound OSe. This narrower conduction band dispersion may 27

measurements on La2O2MnSe2 indicate an optical band gap of 2.31 eV (Peschke et al., 2015). Again, the Ce ana- logue Ce2O2MnSe2 has slightly different properties, with polycrystalline samples being purple (single crystals or- ange (Wang et al., 2015)), an activation energy for elec- tronic conduction of 0.41(1) eV and room temperature conductivity of ∼ 9 × 10−6 Ω−1 cm−1, presumably due to the influence of Ce 4f states near EF (Peschke et al., 2015). The β- and monoclinic polymorphs of La2O2FeSe2 allow us to consider the effect of the Fe coordination en- vironment on the electronic structure. β–La2O2FeSe2 (with Fe(1)Se4O2 and Fe(2)Se4 sites) is a black semicon- ductor with a band gap of ∼ 0.7 eV and room tempera- ture resistivity ∼ 102 Ω cm (McCabe et al., 2011). The FIG. 25 Electronic structure of (a) Ce2O2FeSe2 and (b) monoclinic polymorph of La2O2FeSe2 with only FeSe4O2 (BaFe2Se3 in the nonmagnetic state. The figure is taken from sites is also a semiconductor but with a smaller band gap Ref. Li et al., 2014. (∼0.3 eV) and slightly lower room temperature resistiv- ity (Nitsche et al., 2014). account for the difficulty in doping La2O2CdSe2 (Hi- ramatsu et al., 2004b). La O ZnSe behaves similarly 2 2 2 B. Ln2O2M 2OSe2 structures and related phases (with high electrical resistivity, a band gap of 3.4(2) eV from diffuse reflectance measurements, and difficul- Due to the connection with superconducting cuprates ties with aliovalent doping) and density functional theory and iron based compounds, the Ln2O2M2OSe2 series of calculations indicate that the conduction band is mainly materials have been studied in depth. These compounds composed of La states, leading to a larger band gap than have typically smaller band gaps than the ZrCuSiAs com- in LaCuOSe (Tuxworth et al., 2013). Unlike the Ln = pounds discussed above. La systems, Ln = Ce analogues often have much smaller The oxyselenides are all, nearly universally, semicon- band gaps and higher conductivities due to the Ce 4f and ductors and sometimes described as “bad-metals” with 5d bands near the band gap (Ainsworth et al., 2015b; the resistivity increasing with decreasing temperature. Pitcher et al., 2009; Ueda et al., 2003). This occurs, for example in La2O2Fe2OSe2, even though The cation-ordered ZrCuSiAs-related phases magnetic fluctuations are gapped as described above. Ln2O2FeSe2 tend to have smaller band gaps and This differs from the case of Fe1+xTe where gapped mag- semiconducting behaviour. Polycrystalline samples netic excitations were found to coincide with a transition of Ce2O2FeSe2 are black and exhibit semiconducting from a “semi metallic/bad metal” state to a metallic re- behaviour (with electronic band gap ∼0.64 eV) and sistivity. The connection between low-energy spin fluc- room temperature resistivity of ∼20 Ω cm (McCabe tuations and the resistivity using the formula outlined et al., 2011). Electronic structure calculations were con- above in the experimental section was made in Ref. Ro- sistent with this and indicate itinerant Fe 3d states and driguez et al., 2013. A summary of the activation ener- Ea/kB T suggest that, despite the “stripe” ordering, the electronic gies extracted by fitting the resistivity to ρ = ρ◦e structure is far from being pseudo-one dimensional (Li is summarized in Table VI for a series of oxyselenides. et al., 2014). This is illustrated in Fig. 25 which shows The data on these oxyselenides is difficult to interpret, band structure presented in Ref. Li et al., 2014. The but it generally seems to follow a trend that for decreas- oxyselenide Ce2O2FeSe2 is shown in Fig. 25 (a) where ing rare earth radius the activation energy decreases. it is seen that three bands cross the Fermi level. Ref. Li This might imply a correlation between lattice constant et al., 2014 calculated the possible magnetic structures and cell volume and activation energy. This trend is con- and found an ordered moment of 3.12 µB consistent with firmed in Fig. 26 which plots the resistivity fit to a “small experimental work of 3.33 µB (McCabe et al., 2011). polaron hopping” model ρ(T ) = AT exp(Ea/kBT ). With Density functional calculations suggest that, like the increasing Mn2+ doping the lattice constants also in- copper systems, the valence band and the conduc- crease and therefore this plot confirms the trend which tion band are composed mainly of Fe 3d and Se 4p was suggested from the data in Table VI that the acti- states (McCabe et al., 2014b). By contrast, the man- vation energy scales with the lattice constant. We note ganese analogues Ln2O2MnSe2 generally have wider that the trend of decreasing cell volume corresponding band gaps and more insulating behaviour: polycrys- to lower activation energies is reflected in comparing talline Ln2O2MnSe2 (Ln = La, Pr, Nd) samples tend data on Sr2OBi2Se3 (Eg=0.0092(1) eV) and Ba2OBi2Se3 to be orange-brick red in colour and diffuse reflectance (Eg=0.11(1) eV) (Panella et al., 2016) and indeed may 28

7 TABLE VI Activation energies extracted from resistivity Fe1 3d Ea/kB T 6 (a) data (ρ = ρ◦e ). Se 4p 5 O1 total Compound Ea Reference 4 O2 2p Fe1(S) 3d La2O2Fe2OSe2 0.19 Zhu et al., 2010 3 Ce2O2Fe2OSe2 0.26 Ni et al., 2011 2

Pr2O2Fe2OSe2 0.15 Ni et al., 2011 1 Partial DOS (States/eV) Nd2O2Fe2OSe2 0.15 Ni et al., 2011 6 Sm2O2Fe2OSe2 0.18 Ni et al., 2011 5 (b) Fe 3d As 4p La O Fe OS 0.24 Zhu et al., 2010 2 2 2 2 4 O 2p Ce2O2FeSe2 0.32 McCabe et al., 2011 3 BaFe2Se2O 0.29 Lei et al., 2012 2 La2O2Mn2OSe2 0.24 Free et al., 2011 1 La2O2Co2OSe2 0.35 Free et al., 2011 Partial DOS (States/eV) 0 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 Energy (eV)

FIG. 27 The partial density of states for paramagnetic La2O2Fe2OSe2 (a) and LaOFeAs (b) taken from Ref. Zhu et al., 2010. The curve fits are to the “small polaron hopping” model described in the main text. An increase of activation energy with increasing lattice constant is observed.

FIG. 26 The resistivity and activation energies for Nd2O2(Fe1−xMnx)2OSe2 taken from Ref. Liu et al., 2015a. be a general feature across the oxyselenides. Electronic structure calculations have been performed for a number of Ln2O2M2O(Se,S)2 oxyselenides. In par- FIG. 28 Resistivity as a function of temperature and pressure ticular for La2O2Fe2O(Se,S)2 the electronic structure was for La2O2Co2OSe2. The figure is taken from Ref. Wang et al., calculated using density functional theory and is reported 2010. in Ref. Zhu et al., 2010. Fig. 27 shows the projected density of states for both compounds were it is shown that the 3d electrons on the iron site contribute strongly ing that the behavior is insulating in contrast to iron to the density of states at the Fermi energy. The other based pnictides and chalcogenides which are typically point noted in the calculation is that the d-electron den- metallic. Resistivity data reported in Ref. Zhu et al., sity of states for iron is mostly confined between -2 and 2010 find activation energy gaps of 0.19 and 0.24 eV for 1.2 eV. This represents a considerable narrowing of the La2O2Fe2O(Se,S)2 respectively. Based on the electronic iron d-electron band compared with pnictides such as the behaviour combined with the antiferromagnetic ordering ZrCuSiAs-related as LaFeAsO (Fig. 27 panel (b)) where at low temperatures, these materials were classified as the band occurs between -2.2 and 2 eV and also FeTe Mott insulators (Zhu et al., 2010). or FeSe based chalcogenides (Subedi et al., 2008). The A similar picture for the electronic structure is found narrower iron 3d electronic bands point to enhanced cor- for La2O2Co2O(Se,S)2 and reported in Ref. Wu, 2010. relations. This is substantiated by resistivity data show- Density functional calculations find an even narrower 29 electronic bandwidth than for the iron analogue with Co 3d bandwidth of ∼ 2.8 eV (Wu, 2010). Resistivity data was reported in Ref. Wang et al., 2010 with an activation energy of 0.35 eV which led to the conclusion that the 3d electrons associated with Co2+ were localized. In the context of the discussion in relation to Mott insulating behavior, Ref. Wang et al., 2010 investigated the pres- sure dependence of the resistivity where a sharp (several orders of magnitude) drop in resistivity was reported. Based on this, and a comparison to first principle calcu- lations, the authors concluded that La2O2Co2OSe2 was a “marginal” Mott insulator and the pressure dependence (Fig. 28) suggests that La2O2Co2OSe2 is proximate to a metallic state.

While the first principle calculations above have pointed towards Mott insulating behaviour where in- sulating, or semiconducting, properties are the result of electron correlations, recent x-ray inelastic scattering data has come to a slightly different conclusion (Freelon et al., 2015). By combining resonant inelastic x-ray spec- troscopy with first principle calculations, Ref. Freelon et al., 2015 suggested that the electronic properties of La2O2Fe2OSe2 are more reminiscent of a Kondo insulator where a gap opens due to hybridisation of orbitals. This was established through density functional calculations of the orbitally resolved self-energies. It was proposed that La2O2Fe2OSe2 was a “Mott-Kondo” insulator.

The suggestion of a combination of electronic cor- relations (termed Mott insulators) and orbital effects (Kondo insulators) mimics recent proposals for a new FIG. 29 The resistivity of LaO0.5F0.5BiSe2 as a function of temperature and magnetic field taken from Ref. Krzton- type of metallic state term “Hunds metals” (Georges Maziopa et al., 2014. et al., 2013). The idea of Hunds metals has evolved from a proposal based on local density approximation calcula- tions in LaOFeAs where it was noted that the splitting of the crystal fields due to a tetragonal distortion are C. Superconductivity and the Oxyselenides comparable to the overall crystal splitting between the |ei and |ti states. In this case, it was noted that even Iron oxychalcogenides FeQ (Q = S, Se, Te) have been a small Hunds coupling would result in a spin transi- shown to display unconventional superconductivity with tion from S=2 to S=1 (Haule and Kotliar, 2009). Such properties tuned by intercalation chemistry between the a framework could explain the low ordered moments in antifluorite-like FeQ layers (Borg et al., 2016; Dagotto, the pnictides (for example gS=0.5 ± 0.05 µB in FeAs 2013; Lynn et al., 2015; Rodriguez et al., 2011a; Tay- from neutron diffraction) where much larger moments are lor et al., 2013; Vivanco and Rodriguez, 2016). To the clearly expected in the case of weak Hunds coupling or best of our knowledge, the only oxychalcogenide super- in the intermediate crystal field description. The model conductors are of the LnOBiS2 family (Tanaka et al., also provides a means of explaining a strongly correlated 2014) including LaO0.5F0.5BiSe2 with superconducting metal and has been applied to LaO1−xFxFeAs (Haule Tc = 2.6 K (Fig. 29) (Krzton-Maziopa et al., 2014; Mi- et al., 2008). The tuning from a Hunds metal to a Mott zoguchi and Hosono, 2011), LaO0.5F0.5BiSeS with a Tc insulator has been proposed to be sensitive to the Fe-Fe = 3.8 K (Wang et al., 2015), and related rare earth sub- distance (Yin et al., 2010, 2012) which is interesting in stituted compounds (Demura et al., 2015, 2013). the context of the resistivity measurements under pres- LaO0.5F0.5BiSe2 displays metallic conductivity (rather sure noted above. than semiconducting or bad metallic conductivity de- 30 scribed above for other oxyselenides). The low super- widths) and magnetic structures. conducting transition temperatures have been attributed The final concept to highlight is the local environ- to the significant distortion of the BiSe2 layers (Tanaka ment and particularly crystal field effects. While re- et al., 2014), with transition temperatures suppressed sults from magnetic neutron diffraction are consistent further with pressure (Kotegawa et al., 2012). Electronic with weak/intermediate crystal fields for transition metal studies combining photoemission spectroscopy (Nagira sites, recent suggestions imply that coupling between et al., 2014; Ye et al., 2014; Zeng et al., 2014) and photo- electronic and orbital properties may occur. (Chen et al., electron spectroscopy (Saini et al., 2014) on these and re- 2009) Crystal field effects have been shown to have a lated oxyselenides give good agreement with band struc- role in the magnetic and structural behaviour of other ture calculations, suggesting that correlation effects may mixed-anion systems (Kimber et al., 2008, 2010; Wild- not be important. However, it has been suggested that man et al., 2015) and may also be relevant to understand- LaO0.54F0.46BiS2 is close to a topology change in the ing the low-temperature behaviour of Pr2O2M2OSe2 (M Fermi surface (Terashima et al., 2014). This conclu- = Mn, Fe) (Free et al., 2011; Ni et al., 2011). Oxyse- sion and their more conductive nature clearly set these lenides provide a diverse series of materials in which novel compounds apart from the other oxyselenides discussed magnetic and electronic phenomena can be studied. An above. These LnOBiQ2-related materials represent a example of this is the orthogonal, 2−k magnetic structure new development in the field of oxyselenide research. adopted by La2O2Fe2OSe2 and Sr2F2Fe2OS2, which, to the best of our knowledge, is unique among magneti- cally ordered systems (McCabe et al., 2014a; Zhao et al., VI. CONCLUSION 2013) This structure results from coupling between or- bital and electronic properties and from competition be- There are several concepts that underlie the magnetic tween anisotropy and competing exchange interactions, and electronic properties of the oxyselenides. The first is which will continue to challenge theory. dimensionality, usually as a result of the anion-ordering, While a number of studies on oxyselenides have been giving layered crystal structures. This influences the performed, the lack of large single crystals is hampering electronic structures, highlighted by the band narrow- efforts to fully understand the magnetism of these sys- ing in La2O2Fe2OQ2 predicted by Zhu et al. (Zhu et al., tems and synthetic efforts to produce sufficiently large 2010) and consistent with the small values for magnetic crystals would be enhance the field. The lack of un- exchange interactions observed experimentally (McCabe conventional superconductivity among magnetic oxyse- et al., 2014a), and also by the confinement of holes in lenides (in contrast to the selenides and to oxypnictides) 2− LaCuOSe materials confined to the [Cu2Se2] layers by is interesting, especially given the strong electron cor- 2+ the insulating [La2O2] layers (Ueda et al., 2004a). This relation effects. The connections between cuprates and dimensionality also influences the magnetic ordering with iron-based superconductors will remain a point of future magnetic stacking faults and longer magnetic correlation study. lengths within layers found for several Ln2O2M2OSe2 materials (McCabe et al., 2014a,a; Ni et al., 2010; Zhao et al., 2013). REFERENCES The second key concept is connectivity which influ- Abbas, Y., J. Rossat-Mignod, G. Quezel, and C. Vettier ences the magnetic and electronic structures: the mag- (1973), Solid State Commun. 12, 985. netic frustration resulting from the tetrahedral arrange- Abbas, Y., J. Rossat-Mignod, G. Quezel, and C. Vettier 3+ ment of magnetic Ln ions in Ln4O4Se3 materials dis- (1974), Solid State Commun. 14, 1115. cussed in Section 3a is a good example of this, as well as Abragam, A., and B. Bleaney (1986), Electron paramag- the strong interplay between Cr3+ and Ln3+ magnetism netic resonance of transition ions (Dover Publications, New 3+ 3+ York). in LnCrOS2 materials with Cr and Ln coordina- tion polyhedra linked via oxide anions (Takano et al., Ainsworth, C. M., C. H. Wang, , H. Johnston, E. E. McCabe, M. G. Tucker, and H. B. J. S. O. Evans (2015a), Inorg. 1999, 2002; Wintenberger et al., 1987). The Ce2O2FeSe2 Chem. (in press). oxyselenide with one-dimensional chains of FeSe4 tetra- Ainsworth, C. M., C. H. Wang, M. G. Tucker, and J. S. O. hedra (McCabe et al., 2011, 2014a) is surprising in this Evans (2015b), Inorg. Chem. 54, 1563. respect: in contrast to its one-dimensional connectivity, Altmannshofer, S., and D. Johrendt (2008), Z. Anorg. Allg. its electronic structure is far from one-dimensional (Li Chem. 634, 1361. et al., 2014). The preparation of several polymorphs of Bacon, G. E. (1975), Neutron Diffraction (Monographs on Ln O FeSe (Ln = La, Ce) built up from FeSe tetrahe- the physics and chemistry of materials) (Oxford University 2 2 2 4 Press, Oxford). dra, FeSe4O2 pseudo-octahedra and from combinations Ballhausen, C. J. (1962), Ligand Field Theory (McGraw-Hill, of these (Nitsche et al., 2014) coordination environments New York). will provide an ideal means to investigate the role of Bao, W., Y. Qiu, Q. Huang, M. A. Green, P. Zajdel, M. R. the coordinating anion on the electronic structure (band Fitzsimmons, M. Zhernenkov, S. Chang, M. Fang, B. Qian, 31

E. K. Vehstedt, J. Yang, H. M. Pham, L. Spinu, and Z. Q. Dugue, J., T. Vovan, and J. Villiers (1980b), Acta Cryst. B Mao (2009), Phys. Rev. Lett. 102, 247001. 36, 1294. Barrier, N., and S. J. Clarke (2003), Chem. Commun. 1, 164. Eick, H. A. (1960), Acta Cryst. 14, 161. Bednorz, J. G., and K. A. Muller (1986), Z. Phys. B 64, 189. Emsley, J. (2001), Nature’s building blocks: an A-Z guide to Benazeth, S., M. Guittard, and P. Laruelle (1984), Acta the elements (Oxford University Press, Oxford). Cryst. C 40, 345. Evans, J. S. O., E. B. Brogden, A. L. Thompson, and R. L. Birgeneau, R. J., C.Stock, J. M. Tranquada, and K. Yamada Cordiner (2002), Chem. Commun. 8, 912. (2006), J. Phys. Soc. Jpn. 75, 111003. Free, D. G., and J. S. O. Evans (2010), Phys. Rev. B 81, Bodenan, F., V. B. Cajipe, G. Ouvard, and G. J. Andrew 214433. (1996), J. Magn. Magn. Mater. 164, 233. Free, D. G., N. D. Withers, P. J. Hickey, and J. S. O. Evans Borg, C. K. H., Z. Ahou, C. Eclberg, D. J. Campbell, S. R. (2011), Chem. Mater. 23, 1625. Saha, J. Paglione, and E. E. Rodriguez (2016), Phys. Rev. Freelon, B., Y. H. Liu, J. L. Chen, L. Craco, M. S. Laad, B 93, 094522. S. Leoni, J. Chen, L. Tao, H. Wang, R. Flauca, Z. Yamani, Brechtel, E., G. Cordier, and H. Schaefer (1979), Naturforsch. M. Fang, C. Chang, J. H. Guo, and Z. Hussain (2015), B: Anorg. Chem. Org. Chem. 34, 777. Phys. Rev. B 92, 155139. Brock, S. L., J. E. Greedan, and S. M. Kauzlarich (1994), J. Fuwa, Y., T. Endo, M. Wakeshima, Y. Hinatsu, and Solid State Chem. 113, 303. K. Ohoyama (2010a), J. Am. Chem. Soc. 132, 18020. Caron, J. M., J. R. Neilson, D. C. Miller, A. Llobet, and Fuwa, Y., M. Wakeshima, and Y. Hinatsu (2010b), J. Phys.: T. M. McQueen (2011), Phys. Rev. B 84, 180409(R). Condens. Matter 22, 346003. Ceolin, R., and B. Rodier (1976), Acta Cryst. B 32, 1476. Fuwa, Y., M. Wakeshima, and Y. Hinatsu (2010c), Solid Chen, C. C., B. Moritz, J. van den Brink, T. P. Devereaux, State Commun. 150, 1698. and R. R. P. Singh (2009), Phys. Rev. B 80, 180418(R). Gal, Z. A., O. J. Rutt, C. F. Smura, T. P. Overton, N. Barrier, Chi, S., D. T. Adroja, T. Guidi, R. Bewley, S. Li, J. Zhao, S. J. Clarke, and J. Hadermann (2006), J. Am. Chem. Soc. J. W. Lynn, C. M. Brown, Y. Qiu, G. F. Chen, J. L. Lou, 128, 8530. N. L. Wang, and P. Dai (2008), Phys. Rev. Lett. 101, Gardner, J. S., M. J. Gingras, and J. E. Greedan (2010), 217002. Rev. Mod. Phys. 82, 53. Chou, T. L., G. C. Tewari, T. S. Chan, Y. Y. Hsu, J. M. Georges, A., L. de’Medici, and J. Mravlje (2013), Annual Chen, H. Yamauchi, and M. Karppinen (2015), Eur. J. Review of Condensed Matter Physics 4, 137. Inorg. Chem. 15, 2574. Ghivelder, L., I. A. Castillo, M. A. Gusmao, J. A. Alonso, Christianson, A. D., M. D. Lumsden, S. E. Nagler, G. F. Mac- and L. F. Cohen (1999), Phys. Rev. B 60, 121284. Dougall, M. A. McGuire, A. S. Sefat, R. Jin, B. C. Sales, Goremychkin, E. A., R. Osborn, C. H. Wang, M. D. Lumsden, and D. Mandrus (2009), Phys. Rev. Lett. 103, 087002. M. A. McGuire, A. S. Sefat, B. C. Sales, D. Mandrus, H. M. Clarke, S. J., P. Adamson, S. J. C. Herkelrath, O. J. Rutt, Ronnow, Y. Su, and A. D. Christianson (2011), Phys. Rev. D. R. Parjer, M. J. Pitcher, and C. F. Smura (2008), Inorg. B 83, 212505. Chem. 47, 8473. Greenwood, N., and A. Earnshaw (2005), Chemistry of Coldea, R., S. M. Hayden, G. Aeppli, T. G. Perring, C. D. the Elements, Vol. 2nd edition (Elsevier Butterworth- Frost, T. E. Mason, S. W. Cheong, and Z. Fisk (2001), Heinemann). Phys. Rev. Lett. 86, 5377. Gretarsson, H., S. R. Saha, T. Dyre, J. Paglione, J. Kim, Collins, M. F. (1989), Magnetic Critical Scattering (Oxford D. Casa, T. Gog, W. Wu, S. R. Julian, and Y. J. Kim University Press, Oxford). (2013), Phys. Rev. Lett. 110, 047003. Corkett, A. J., D. G. Free, and S. Clarke (2014), J. Inorg. Griffiths, J. (1964), The Theory of Transition-Metal Ions Chem. 54, 1178. (Cambridge University Press, London). Cowley, R. A., W. J. L. Buyers, C. Stock, Z. Yamani, C. Frost, Gu, J., Z. Q. Zhao, Y. Ding, H. L. Chen, Y. W. Zhang, and J. W. Taylor, and D. Prbhakaran (2013), Phys. Rev. B 88, C. H. Yan (2013), J. Am. Chem. Soc. 135, 8363. 205117. Gunther, M., S. Kamusella, R. Sarker, T. Goltz, H. Luetkens, Dagotto, E. (2013), Rev. Mod. Phys. 85, 849. G. Pascua, S. H. Do, K. Y. Choi, H. D. Zhou, C. G. F. Dai, D., and M. H. Whangbo (2005), Inorg. Chem. 44, 4407. Blum, S. Wurmehl, B. Buchner, and H. H. Klauss (2014), Dai, P. (2015), Rev. Mod. Phys. 87, 855. Phys. Rev. B 90, 184408. Das, P., N. S. Sangeetha, A. Pandey, Z. A. Benson, T. W. Hammar, P. R., D. H. Reich, C. Broholm, and F. Trouw Heitmann, D. C. Johnston, A. I. Goldman, and A. Kreyssig (1998), Phys. Rev. B. 57, 7846. (2016), unpublished (arXiv:1605.02004v1). Harriger, L. W., H. Q. Luo, M. S. Liu, C. Frost, J. P. Hu, Demura, S., K. Deguchi, Y. Mizuguchi, K. Sato, R. Honjyo, M. R. Norman, and P. Dai (2011), Phys. Rev. B 84, A. Yamashita, T. Yamaki, H. Hara, T. Watanabe, S. J. 054544. Denholme, M. Fujioka, H. Okazaki, T. Ozaki, O. Miura, Haule, K., and G. Kotliar (2009), New J. Phys. 11, 025021. T. Yamaguchi, and H. T. ad Y. Takano (2015), J. Phys. Haule, K., J. H. Shim, and G. Kotliar (2008), Phys. Rev. Soc. Jpn. 84, 024709. Lett. 100, 226402. Demura, S., Y. Mizuguchi, K. Deguchi, H. Okazaki, H. Hara, Hawai, T., Y. Nambu, K. Ohgushi, F. Du, Y. Hirata, T. Watanabe, S. J. Denholme, M. Fujioka, T. Ozaki, H. Fu- M. Avdeev, Y. Uwatoko, Y. Sekine, H. Fukazawa, and juhisa, Y. Gotoh, O. Miura, T. Yamaguchi, H. Takeya, and J. Ma (2015), Phys. Rev. B 91, 184416. Y. Takano (2013), J. Phys. Soc. Jpn. 82, 033708. Hayden, S. M., G. Aeppli, T. G. Perring, H. A. Mook, and Dugue, J., C. Adolphe, and P. Khodadad (1970), Acta Cryst. F. Dogan (1996), Phys. Rev. B 54, R6905R. B: Struct. Crystallogr. Cryst. Chem. 26, 1627. He, J. B., D. M. Wang, H. L. Shi, H. X. Yang, and G. F. Dugue, J., T. Vovan, and J. Villiers (1980a), Acta Cryst. B Chen (2011), Phys. Rev. B 84, 205212. 36, 1291. Her, J. H., P. W. Stephens, Y. Gao, G. L. Soloveichil, J. Ri- 32

jssenbeek, M. Andrus, and J. C. Zhao (2007), Acta Cryst. C. Hess, G. Berh, A. Kondrat, and B. Buchner (2008), B 63, 561. Phys. Rev. B 78, 140503(R). Herkelrath, S. J. C., I. Saratovsky, J. Hadermann, and S. J. Kimber, S. A. J., A. H. Hill, Y. Zhong-Zhong, H. O. Clarke (2008), J. Am. Chem. Soc. 130, 14426. Jeschke, R. Valenti, C. Ritter, I. Schellenberg, W. Hermes, Hiramatsu, H., T. Kamiya, T. Tohei, E. Ikenaga, T. Mizogu- R. Pottgen, and D. N. Argyriou (2010), Phys. Rev. B 82, ichi, Y. Ikuhara, K. Kobayashi, and H. Hosono (2010a), J. 100412. Am. Chem. Soc. 132, 15060. Koo, H. J., and M. H. Whangbo (2012), J. Magn. Magn. Hiramatsu, H., T. Kamiya, K. Ueda, M. Hirano, and Mater. 324, 3859. H. Hosono (2010b), Phys. Status Solidi A 207, 1636. Kotani, A., and S. Shin (2001), Rev. Mod. Phys. 73, 203. Hiramatsu, H., K. Ueda, T. Kamiya, H. Ohta, M. Hirano, Kotegawa, H., Y. Tomita, H. Tou, H. Izawa, Y. Mizuguchi, and H. Hosono (2004a), J. Mater. Chem. 14, 2946. O. Miura, S. Demura, K. Deguchi, and Y. Takano (2012), Hiramatsu, H., K. Ueda, T. Kamiya, H. Ohta, M. Hirano, J. Phys. Soc. Jpn. 81, 103702. and H. Hosono (2004b), J. Phys. Chem. B 108, 17344. Koz, C., S. Robler, A. A. Tsirlin, S. Wirth, and U. Scharz Hiramatsu, H., K. Ueda, H. Ohta, M. Hirano, M. Kituchi, (2013), Phys. Rev. B 88, 094509. H. Yanagi, T. Kamiya, and H. Hosono (2007), Appl. Phys. Krivovichev, S. V., O. Mentre, O. I. Siidra, M. Colmont, and Lett. 91, 012104. S. K. Filatov (2013), Chem. Rev. 113, 6459. Hiramatsu, H., H. Yanagi, T. Kamiya, K. Ueda, M. Hirano, Kruger, F., S. Kumar, J. Zaanen, and J. van den Brink and H. Hosono (2008), Chem. Mater. 20, 326. (2009), Phys. Rev. B 79, 054504. Hohenberg, P. C., and W. F. Brinkman (1974), Phys. Rev. Krzton-Maziopa, A., Z. Guguchia, E. Pomjakushina, V. Pom- B 10, 128. jakushin, R. Khasanov, H. Luetkens, P. K. Biswas, A. Am- Hong, T., M. Kenzelmann, M. M. Turnbull, C. P. Landee, ato, H. Keller, and K. Conder (2014), J. Phys.: Condens. B. D. Lewis, K. P. Schmidt, G. S. Uhrig, Y. Qiu, C. Bro- Matter 26, 215702. holm, and D. Reich (2006), Phys. Rev. B. 74, 094434. Landsgesell, S., E. Blumenrother, and K. Prokes (2013), J. Huang, Q., A. Santoro, J. W. Lynn, R. W. Erwin, J. A. Phys. Condens. Matter 25, 086004. Borchers, J. L. Peng, and R. L. Greene (1997), Phys. Rev. Lei, H., E. S. Bozin, A. Llobet, V. Ivanovski, V. Koteski, B 55, 14987. J. Belosevic-Cavor, B. Cekic, and C. Petrovic (2012), Phys. Hutchings, M. T. (1964), Solid State Phys. 16, 227. Rev. B 86, 125122. Hyett, G., N. Barrier, S. J. Clarke, and J. Hadermann (2007), Lei, H., K. Wang, M. Abeykoon, E. S. Bozin, and C. Petrovic J. Am. Chem. Soc. 129, 11192. (2013), Inorg. Chem. 52, 10685. Ijjaali, I., K. Mitchell, C. L. Haynes, A. D. McFarland, Li, W., C. Setty, X. H. Chen, and J. Hu (2014), Front. Phys. R. P. V. Duyne, and J. A. Ibers (2003), J. Solid State 9, 465. Chem. 176, 170. Lin, X., X. Ni, B. Chen, X. Xu, X. Yang, J. Dai, Y. Li, Indris, S., J. Cabana, O. J. Rutt, S. J. Clarke, and C. P. X. Yang, Y. Luo, Q. Tao, G. Cao, and Z. Xu (2013), Phys. Grey (2006), J. Am. Chem. Soc. 128, 13354. Rev. B 87, 020504. Izyumov, Y. A., and R. P. Ozerov (1970), Magnetic Neutron Liu, R. H., J. S. Zhang, P. Chen, X. G. Luo, J. J. Ying, Y. J. Diffraction (Springer, New York). Yan, M. Zhang, A. F. Wang, Z. J. Xiang, G. F. Ye, and Jaulmes, S. (1978), Acta Cryst. B 34, 2610. X. H. Chen (2011), Phys. Rev. B 83, 174450. Jaulmes, S., E. Godlewski, M. Palazzi, and J. Etienne (1982), Liu, Y., S. B. Zhang, W. J. Lu, L. J. Li, S. G. Tan, B. Yuan, Acta Cryst. B 38, 1707. J. Chen, and Y. P. Sun (2015a), J. Alloys Comp. 618, 263. Jin, S., X. Chen, J. Guo, J. Lin, J. Xi, W. Wang, and Liu, Y., S. B. Zhang, S. G. Tan, B. Yuan, X. C. Kan, L. Zu, W. Wang (2012), Inorg. Chem. 51, 10185. and Y. P. Sun (2015b), J. Solid State Chem. 221, 272. Jin, S. F., Q. Huang, Z. P. Lin, Z. L. Li, X. Z. Qu, T. P. Llanos, J., R. COrtes, T. Guizouarn, and O. Pena (2006), Ying, G. Wang, and X. L. Chen (2015), Phys. Rev. B 91, Mat. Res. Bull. 41, 1266. 094420. Llanos, J., and O. Pena (2005), J. Solid State Chem. 178, Johnson, C., and W. Jeitschko (1974), J. Solid State Chem. 957. 11, 161. Luo, Z. Z., C. S. Lin, W. D. Cheng, Y. B. Li, H. Zhang, W. L. Johnston, D. C. (2010), Adv. Phys. 59, 803. Whang, and Z. Z. He (2013), Dalton Trans. 42, 9938. Kabbour, H., E. Janod, B. Corraze, M. Danot, C. Lee, M. H. Lynn, J. W., X. Zhou, C. K. H. Borg, S. R. Saha, J. Paglione, Whangbo, and L. Cario (2008), J. Am. Chem. Soc. 130, and E. E. Rodriguez (2015), Phys. Rev. B 92, 060510(R). 8261. Martinolich, A. J., R. F. Higgins, M. P. Shores, and J. R. Kamihara, Y., T. Watanabe, M. Hirano, and H. Hosono Neilson (2016), Chem. Mater. 28, 1854. (2008), J. Am. Chem. Soc. 130, 3296. Martinolich, A. J., J. A. Kurzman, and J. R. Neilson (2015), Karipinski, J., N. D. Zhigadlo, S. Katrych, Z. Bukowski, J. Am. Chem. Soc. 137, 3827. P. Moll, S. Weyeneth, H. K. ad R. Puzniak, M. Tortello, Martinolich, A. J., and J. R. Neilson (2014), J. Am. Chem. D. Daghero, R. Gonnelli, I. Maggio-Aprile, Y. Fasno, Soc. 136, 15654. O. Fischer, K. Rogacki, and B. Batlogg (2009), Physica Mayer, J. M., L. F. Schneemeyer, T. Siegrist, J. V. Waszczak, C 469, 370. and B. V. Dover (1992), Angew. Chem. Int. Ed. Engl. 31, Kastner, M. A., R. J. Birgeneau, G. Shirane, and Y. Endoh 1645. (1998), Rev. Mod. Phys. 70, 897. McCabe, E. E., D. G. Free, and J. S. O. Evans (2011), Chem. Kim, Y. J., A. P. Sorini, C. Stock, T. G. Perring, J. van den Commun. 47, 1261. Brink, and T. P. Devereaux (2011), Phys. Rev. B 84, McCabe, E. E., D. G. Free, B. G. Mendis, J. S. Higgins, and 085132. J. S. O. Evans (2010), Chem. Mater. 22, 6171. Kimber, S. A. J., D. N. Argyriou, F. Yokaichiya, K. Habicht, McCabe, E. E., C. Stock, J. L. Bettis, M. H. Whangbo, and S. Gerischer, T. Hansen, T. Chatterji, R. Klingeler, J. S. O. Evans (2014a), Phys. Rev. B 90, 235115. 33

McCabe, E. E., C. Stock, E. E. Rodriguez, A. S. Wills, J. W. J. Paglione, and M. A. Green (2011a), Chem. Sci. 2, 1782. Taylor, and J. S. O. Evans (2014b), Phys. Rev. B 89, Rodriguez, E. E., C. Stock, K. L. Krycka, C. F. Majkrzak, 100402(R). P. Zajdel, K. Kirshenbaum, N. P. Butch, R. S. Saha, McClure, D. (1959), Sol. State. Phys. 9, 399. J. Paglione, and M. A. Green (2011b), Phys. Rev. B 83, Meerschaut, A., C. Boyer, A. Lafond, L. Cario, and J. Rouxel 134438. (1998), J. Solid State Chem 137, 122. Rossat-Mignod, J., J. C. Souillat, and G. Quezel (1974), Meerschaut, A., A. Lafond, V. Meignen, and C. Deudon Phys. Stat. Sol. (b) 62, 223. (2001a), J. Solid State Chem 162, 182. Roth, W. L. (1958), Phys. Rev. 110, 1333. Meerschaut, A., A. Lafond, V. Meignen, and C. Deudon Ruff, E., P. Lunkenheimer, A. Loidl, H. Berger, and (2001b), J. Solid State Chem. 162, 182. S. Krohns (2015), Scientific Reports 5, 15025. Meunier, P. G., and M. Bertraud (1976), Appl. Cryst. 9, 364. Rutt, O. J., G. R. Williams, and S. J. Clarke (2006), Chem. Minami, T., M. Hibino, and M. Tanaka (1974), J. Non-Cryst. Commun. 27, 2869. Solids 15, 141. Saini, N. L., D. Ootsuki, E. Paris, B. Joseph, A. Barinov, Mizoguchi, H., and H. Hosono (2011), J. Am. Chem. Soc. M. Tanaka, Y. Takano, and T. Mizokawa (2014), Phys. 133, 2394. Rev. B 90, 214517. Mizunoa, Y., T. Tohyama, S. Maekawa, T. Osafune, N. Mo- Scanlon, D. O., J. Buckeridge, R. A. Catlo, and G. W. Wat- toyama, H. Eisaki, and S. Uchida (1998), Phys. Rev. B 57, son (2014), J. Mater. Chem C 2, 3429. 5326. Sedlmaier, S. J., S. J. Cassidy, R. G. Morris, M. Drakopou- Moriya, T., Y. Takahashi, and K. Euda (1990), J. Phys. Soc. los, C. Reinhard, S. J. Moorhouse, D. O’Hare, P. Manuel, Jpn. 59, 2905. D. Khalyavin, and S. J. Clarke (2014), J. Am. Chem. Soc. Muhlbauer, S., B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, 136, 630. A. Neubauer, R. Georgii, and P. Boni (2009), Science 323, Sefat, A. S. (2013), Curr. Opin. Solid St. M. 17, 59. 915. Seki, S., X. Z. Yu, S. Ishiwata, and Y. Tokura (2012), Science Nagira, S., J. Sonoyama, T. Wakita, M. Sunagawa, Y. Izumi, 336, 198. T. Muro, H. Kumigashira, M. Oshima, K. Deguchi, Seo, W. S., H. H. Jo, K. Lee, B. Kim, S. J. Oh, and J. T. H. Okazaki, Y. Takano, O. Miura, Y. Mizuguchi, K. Suzuki, Park (2004), Angew. Chem. Int. Ed. 43, 1115. H. Usui, K. Kuroki, K. Okada, Y. Muraoka, and T. Yokoya Shimizu, T., T. Matsumoto, A. Goto, T. V. C. Rao, (2014), J. Phys. Soc. Jpn. 83, 033703. K. Yoshimura, and K. Kosuge (2003), Phys. Rev. B 68, Ni, N., E. Climent-Pascual, S. Jia, Q. Huang, and R. J. Cava 224433. (2010), Phys. Rev. B 82, 214419. Shirane, G., S. M. Shapiro, and J. M. Tranquada (2001), Neu- Ni, N., S. Jia, Q. Huang, E. Climent-Pascual, and R. J. Cava tron scattering with a triple-axis spectrometer (Cambridge (2011), Phys. Rev. B 83, 224403. University Press, Cambridge). Nitsche, F., A. Jesche, E. Hieckmann, T. Doert, and M. Ruck Singh, Y., A. Ellern, and D. C. Johnston (2009), Phys. Rev. (2010), Phys. Rev. B 82, 134514. B 79, 094519. Nitsche, F., R. Niklaus, and D. Johrendt (2014), Z. Anorg. Skyrme, T. H. R. (1962), Nucl. Phys. 31, 556. Allg. Chem. 640, 2897. Smura, C. F., D. R. Parker, M. Zbiri, M. R. Johnson, Z. A. Nonfante, M., F. Moussa, B. Hennion, and G. Pepy (1972), Gal, and S. J. Clarke (2011), J. Am. Chem. Soc. 133, 2691. Solid State Commun. 10, 553. Snyder, G. J., and E. S. Toberer (2008), Nature Mater. 7, Ogisu, K., A. Ishikawa, Y. Shimodaira, T. Takata, 105. H. Kobayashi, and K. Domen (2008), J. Phys. Chem. C Stampler, E. S., W. C. Sheets, M. I. Beroni, W. Prellier, T. O. 112, 11978. Mason, and K. R. Poeppelmeier (2008), Inorg. Chem. 47, Oogarah, R. K. (2016), private communication. 134514. Otzschi, K., H. Ogino, J. Shimoyama, and K. Kishio (1999), Stanek, D., O. P. Sushkov, and G. S. Uhrig (2011), Phys. J. Low Temp. Phys. 117, 729. Rev. B 84, 064505. Paglione, J., and R. L. Greene (2010), Nat. Phys. 6, 645. Stewart, G. R. (2011), Rev. Mod. Phys. 83, 1589. Panella, J. R., J. Chamorro, and T. M. McQueen (2016), Stock, C., L. C. Chapon, O. Adamopoulos, A. Lappas, Chem Mater. 28, 890. M. Giot, J. W. Taylor, M. A. Green, C. M. Brown, and Park, Y., D. C. Degroot, J. L. Schindler, C. R. Kannewurf, P. G. Radaelli (2009), Phys. Rev. Lett. 70, 077202. and M. Kanatzidis (1993), Chem. Mater. 5, 8. Stock, C., R. A. Cowley, W. J. L. Buyers, R. Coldea, C. Bro- Peschke, S., F. Nitsche, and D. Johrendt (2015), Z. Anorg. holm, C. D. Frost, R. J. Birgeneau, R. Liang, D. Bonn, Allg. Chem. 641. and W. N. Hardy (2007), Phys. Rev. B 75, 172510. Pitcher, M. J., C. F. Smura, and S. J. Clarke (2009), Inorg. Stock, C., R. A. Cowley, W. J. L. Buyers, C. D. Frost, J. W. Chem. 48, 9054. Taylor, D. Peets, R. Liang, D. Bonn, and W. N. Hardy Popovic, Z. V., M. Scepanovic, N. Lazarevic, M. M. Radonjic, (2010a), Phys. Rev. B 82, 174505. D. Tanaskovic, H. Lei, and C. Petrovic (2014), Phys. Rev. Stock, C., R. A. Cowley, J. W. Taylor, and S. M. Bennington B 89, 014301. (2010b), Phys. Rev. B 81, 024303. Prozorov, R., M. E. Tillman, E. D. Mun, and P. C. Canfield Stock, C., E. E. Rodriguez, M. A. Green, P. Zavalij, and J. A. (2009), New J. Phys. 11, 035004. Rodriguez-Rivera (2011), Phys. Rev. B 84, 045124. Regmi, R., R. Tackett, and G. J. Lawes (2009), J. Magn. Stock, C., E. E. Rodriguez, O. Sobolev, J. A. Rodriguez- Magn. Mater. 321, 2296. Rivera, R. A. Ewings, J. W. Taylor, A. D. Christianson, Rodriguez, E. E., D. A. Sokolov, C. Stock, M. A. Green, and M. A. Green (2014), Phys. Rev. B 90, 12113(R). O. Sobolev, J. A. Rodriguez-Rivera, H. Cao, and Stone, M. B., I. Zaliznyak, D. H. Reich, and C. Broholm A. Dauoud-Aladine (2013), Phys. Rev. B 88, 165110. (2001), Phys. Rev. B. 64, 144405. Rodriguez, E. E., C. Stock, P. Y. Hsieh, N. P. Butch, Strobel, S., A. Choudhury, P. K. Dorhout, C. Lipp, and 34

T. Schleid (2008), Inorg. Chem. 47, 4936. Wildman, E. J., F. Sher, and A. C. McLaughlin (2015), Inorg. Subedi, A., L. Zhang, D. J. Singh, and M. H. Du (2008), Chem. 54, 2536. Phys. Rev. B 78, 134514. Willis, B. T. M., and C. J. Carlile (2009), Experiment Neu- Sui, J. H., J. Li, J. Q. He, Y. L. Pei, D. Berardan, H. J. Wu, tron Scattering (Oxford University Press, Oxford). N. Dragoe, W. Cai, and L. D. Zhao (2013), Energy. Env. Wilson, S. D., C. R. Rotundu, Z. Yamani, P. N. Valdivia, Sci. 6, 2916. B. Freelon, E. Bourret-Courchesne, and R. J. Birgeneau Takano, T., K. Takase, G. Kitamura, C. Itoi, and K. Sekizawa (2010), Phys. Rev. B 81, 014501. (1999), J. Appl. Phys. 85, 6103. Wilson, S. D., Z. Yamani, C. R. Rotundu, B. Freelon, Takano, T., T. Tsubaki, C. Itoi, K. Takase, and K. Sekizawa E. Bourret-Courchesne, and R. J. Birgeneau (2009), Phys. (2002), Solid State Commun. 122, 661. Rev. B 79, 184519. Tan, S. G., H. C. Lei, E. J. Lu, P. Tong, L. J. Li, S. Lin, Y. N. Wintenberger, M., J. Dugue, M. Guittard, N. H. Dung, and Huang, Z. H. Huang, Y. Liu, B. C. Zhao, and Y. P. Sun V. V. Tien (1987), J. Solid State Chem. 70, 295. (2014), J. Alloys and Comp. 598, 171. Winterberg, M., V. Tien, M. Guittard, and J. Dugue (1989), Tanaka, M., M. Nagao, Y. Matsushita, M. Fujioka, S. J. Den- J. Solid State Chem. 79, 285. holme, T. Yamaguchi, H. Takeya, and Y. Takano (2014), Wu, H. (2010), Phys. Rev. B 82, 020410(R). J. Solid State Chem. 219, 168. Wu, L. B., and F. Q. Huang (2007), Z. Kristallogr. -New Taylor, A. D., R. Osborn, K. A. McEwen, W. G. Stirling, Z. A. Cryst. Struct. 222, 175. Bowden, W. G. Williams, E. Balcar, and S. W. Lovesey Wysocki, A. L., K. D. Belashchenko, and V. P. Antropov (1988), Phys. Rev. Lett. 61, 1309. (2011), Nat. Phys. 7, 485. Taylor, A. E., S. J. Sedlmaier, S. J. Cassidy, E. A. Goremy- Xiao, Y., M. Zbiri, R. A. Downie, J. W. G. Bos, T. Bruckel, chkin, R. A. Ewings, T. G. Perring, S. J. Clarke, and A. T. and T. Chatterji (2013), Phys. Rev. B 88, 214419. Boothroyd (2013), Phys. Rev. B 87, 220508(R). Yazici, Y., K. Huang, B. White, A. Chang, A. Friedman, and Terashima, K., J. Sonoyama, T. Wakita, M. Sunagawa, M. Maple (2013), Philos. Mag. 93, 673. K. Ono, H. Kumigashira, T. Muro, M. Nagao, S. Watauchi, Ye, Z. R., H. F. Yang, D. W. Shen, J. Jiang, X. H. Niu, D. L. I. Tanaka, H. Okazaki, Y. Takano, O. Miura, Y. Mizuguchi, Feng, Y. P. Du, X. G. Wan, J. Z. Liu, X. Y. Zhu, H. H. H. Usui, K. Suzuki, K. Kuroki, Y. Muraoka, and T. Yokoya Wen, and M. H. Jiang (2014), Phys. Rev. B 90, 045116. (2014), Phys. Rev. B 90, 220512(R). Yin, W. G., C. C. Lee, and W. Ku (2010), Phys. Rev. Lett. Tougait, O., and J. A. Ibers (2000), Acta Crystallogr. Sect. 105, 107004. C: Cryst. Struct. Commun. 56, 623. Yin, W. G., C. H. Lin, and W. Ku (2012), Phys. Rev. B 86, Tsukamoto, Y., Y. Okamoto, K. Matsuhira, M. H. Whangbo, 081106(R). and Z. Hiroi (2011), J. Phys. Soc. Jpn. 80, 094708. Yin, Z. P., K. Haule, and G. Kotliar (2011), Nature Mater. Turberfield, K. C., L. Passell, R. J. Birgeneau, and E. Bucher 10, 932. (1971), J. Appl. Phys. 42, 1746. Yu, R., Z. Wang, P. Goswami, A. H. Nevidomskyy, Q. Si, Turner, A. M., F. Wang, and A. Vishwanath (2009), Phys. and E. Abrahams (2012), Phys. Rev. B 86, 085148. Rev. B 80, 224504. Yu, X. Z., Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Tuxworth, A. J., and J. S. O. Evans (2014), J. Solid State Y. Matsui, N. Nagaosa, and Y. Tokura (2010), Nature Chem. 201, 188. 465, 901. Tuxworth, A. J., E. E. McCabe, D. G. Free, S. J. Clark, and Zeng, L. K., X. B. Wang, J. Ma, P. Richard, S. M. Nie, H. M. J. S. O. Evans (2013), Inorg. Chem. 52, 2078. Weng, N. L. Wang, Z. Wang, T. Qian, and H. Ding (2014), Tuxworth, A. J., C. H. Wang, and J. S. O. Evans (2015), Phys. Rev. B 90, 054512. Dalton Trans. 44, 3009. Zhang, Q., W. Tian, S. G. Peterson, K. W. Dennis, and Ueda, K., H. Hiramatsu, H. Ohta, M. Hirano, T. Kamiya, D. Vaknin (2015), Phys. Rev. B 91, 064418. and H. Hosono (2004a), Phys. Rev. B 69, 155305. Zhao, J., D. T. Adroja, D. X. Yao, R. Bewley, S. Li, X. F. Ueda, K., and H. Hosono (2002), J. Appl. Phys. 91, 4768. Wang, G. Wu, X. H. Chen, J. Hu, and P. Dai (2009), Nat. Ueda, K., H. Hosono, and N. Hamada (2004b), J. Phys. Con- Phys. 5, 555. dens. Matter 16, 5179. Zhao, L. D., D. Berardan, Y. L. Pei, C. Byl, L. Pinsard- Ueda, K., S. Inoue, S. Hirose, H. Kawazoe, and H. Hosono Gaudart, and N. Dragoe (2010), Appl. Phys. Lett. 97, (2000), Appl. Phys. Lett. 77, 2701. 092118. Ueda, K., K. Takafuji, H. Hiramatsu, H. Ohta, T. Kamiya, Zhao, L. D., J. He, D. Berardan, Y. Lin, J. F. Li, C. W. Nan, H. Hirano, and H. Hosono (2003), Chem. Mater. 15, 3692. and N. Dragoe (2014), Energ. Env. Sci. 7, 2900. Vivanco, H. K., and E. E. Rodriguez (2016), J. Solid State Zhao, L. L., S. Wu, J. K. Wang, J. P. Hodges, C. Broholm, Chem. (in press). and E. Morosan (2013), Phys. Rev. B 87, 020406(R). Wang, C., M. Q. Tan, . F. C. M, Z. F. Ma, S. Jiang, Z. A. Zhigadlo, N. D., S. Katrych, Z. Bukowski, S. Weyeneth, Xu, G. H. Cao, K. Matsubayashi, and Y. Uwatoko (2010), R. Puzniak, and J. Karpinski (2008), J. Phys.: Condens. J. Am. Chem. Soc. 132, 7069. Matter 20, 342202. Wang, C.-H., C. M. Ainsworth, D. Y. Gui, E. E. McCabe, Zhou, T., Y. Wang, S. Jin, D. Li, X. Lai, T. Ying, H. Zhang, M. G. Tucker, I. R. Evans, and J. S. O. Evans (2015), S. Shen, W. Wang, and X. Chen (2014), Inorg. Chem 53, Chem. Mater. 27, 3121. 4154. Weber, F. A., and T. Schleid (2001), Z. Anorg. Allg. Chem. Zhu, E. J. (1997), J. Solid State Chem. 130, 319. 627, 1383. Zhu, J. X., R. Yu, H. Wang, L. L. Zhao, M. D. Jones, J. Dai, Weber, F. A., C. Schurz, S. Frunder, C. F. Kuhn, and E. Abrahams, E. Morosan, M. Fang, and A. Si (2010), T. Schleid (2012), Crystals 2, 1136. Phys. Rev. Lett. 104, 216405. Wen, J. S., G. Xu, G. Gu, J. M. Tranquada, and R. J. Bir- Zhu, W. J., P. H. Hor, A. J. Jacobson, G. Crischi, T. A. geneau (2011), Rep. Prog. Phys. 74, 124503. Albright, S. H. Wang, and J. Vogt (1997), J. Am. Chem. 35

Soc. 119, 12398. Chem. A 1, 8888. Zou, D., S. Xie, Y. Liu, J. Lin, and J. Li (2013), J. Mater.