Dichotomy Between Orbital and Magnetic Nematic Instabilities in BaFe2S3

Nematic orders emerge nearly universally in iron-based superconductors, but elucidating their origins is challenging because of intimate couplings between orbital and magnetic fluctuations. The iron-based ladder material BaFe2S3, which superconducts under pressure, exhibits antiferromagnetic order below TN ~ 117K and a weak resistivity anomaly at T* ~ 180K, whose nature remains elusive. Here we report angle-resolved magnetoresistance (MR) and elastoresistance (ER) measurements in BaFe2S3, which reveal distinct changes at T*. We find that MR anisotropy and ER nematic response are both suppressed near T*, implying that an orbital order promoting isotropic electronic states is stabilized at T*. Such an isotropic state below T* competes with the antiferromagnetic order, which is evidenced by the nonmonotonic temperature dependence of nematic fluctuations. In contrast to the cooperative nematic orders in spin and orbital channels in iron pnictides, the present competing orders can provide a new platform to identify the separate roles of orbital and magnetic fluctuations.


I. INTRODUCTION
The discovery of low-dimensional superconductivity in iron-based ladder compounds provides a new point of view in the study of iron-based superconductors [1]. Reduced dimensionality changes electronic structures dramatically, and the ground states of iron-based ladder materials show insulating properties in stark contrast to bad metal behaviors in the typical quasi-two dimensional BaFe 2 As 2 (122) system. In spite of the totally different ground states, the ladder materials still have a stripe-type antiferromagnetic order similar to that of 122system, suggesting the common physics between these two systems. More intriguingly, pressure-induced superconductivity emerges with the suppression of antiferromagnetism [1,2]. These common features suggest unconventional pairing mechanisms of iron-based superconductivity robust against dimensionality.
The insulating nature of the ladder materials implies strong correlation effects due to the reduction of their dimensionality. Heavily hole-doped 122 system also has strong correlation effects since the nominal electron filling 3d 5.5 in iron atoms, which is closer to half-filled states than 3d 6 of non-doped 122-system, pushes towards the putative Mott insulating phase [3]. In contrast to the single-orbital case in e.g. cuprates, where only Hubbard interaction U is dominant, in the multiorbital systems, Hund's coupling J H plays an important role for increasing orbital-dependent correlation effects, leading to orbital selective Mott states [3,4]. An incoherent bad-metal conduction in these systems can be considered as a precursor phenomenon in the vicinity of Mott phases, and * These authors contributed equally to this work.
indeed some spectroscopy measurements reveal strongly orbital-dependent renormalized bands [5][6][7], although the system locates still far away from the parent Mott insulating state. The ladder materials show insulating behaviors, and thus the dimensionality may also serve as another promising parameter to tune the system to the Mott regime, which gives a new route to the study of orbital selective Mott phases. Superconductivity in BaFe 2 S 3 under pressure occurs without significant crystal structure changes [1], indicating that electronic state at ambient pressure smoothly connects to the superconductivity. Therefore, understanding the electronic states of BaFe 2 S 3 at ambient pressure is fundamentally important to study the effects of dimensionality in iron-based materials. It has been established that BaFe 2 S 3 exhibits the antiferromagnetic transition at T N ∼ 117 K with ferromagnetically aligned spins along rung directions [1]. One of the unsolved mysteries in BaFe 2 S 3 is the so-called T * anomaly characterized by a weak bump-like features found in the temperature dependence of resistivity at ∼ 180 K. The origin of the T * anomaly possibly derives from an orbital-involved phase transition, as magnetization, neutron scattering, and muon spin relaxation measurements do not detect any magnetic signatures of phase transitions [8][9][10]. Orbital ordering has been intensively discussed in terms of electronic nematic orders in iron-based superconductors, as one of the promising candidates producing large inplane anisotropy in the electronic state [11]. Therefore, it is effective to evaluate anisotropy of the electronic state in order to verify the orbital order. Polarization dependence of X-ray absorption spectroscopy (XAS) measurements have shown anisotropic electronic states, suggesting the existence of orbital order [12]. However, the observed polarization differences preserve up to 300 K, and there is no characteristic change in spectrum around T * .  Three experimental geometries represent out-of-ladder, in-ladder, and as-grown plane anisotropy measurements as illustrated in the upper insets of (a-c) including our definitions of x, y, z axes in this study and field angle θ with field rotation plane in each geometry. Angle dependence of MR can be fitted to Eq. (1) (solid lines are the fitting results). Temperature dependence of two-fold oscillation term is shown in (d) and (e) (broken lines are guides for the eyes) with error bar based on the standard deviation obtained from the fittings. The guide line for T * , is determined from ER data as shown in Fig. 2(d).
One possibility is that the T * anomaly involves change of low-energy states outside the scope of the XAS measurements which generally probes valence bonds whose energy level is deep below the Fermi level. It is desirable to measure the anisotropy of transport properties that is very sensitive to low-energy quasiparticle excitations.
Here we report anisotropy of electronic states by measuring both field-angle resolved magnetoresistance (MR) and elastoresistance (ER) at zero field. From MR measurements, we can estimate the degree of anisotropy in the order parameter, and in contrast, ER couples to anisotropic electronic instability which corresponds to the nematic fluctuations. Measurements of these two physical values can provide complementary information on the rotational symmetry of the electronic states. The electronic properties of BaFe 2 S 3 are known to be quite sensitive to the growth conditions [13]. Here, by employing the method in Ref. 8, we prepare single crystals without iron deficiency showing T N ∼ 117 K and T * ∼ 180 K, which guarantee to acquire intrinsic information on the orbital order. Our measurements find the anisotropic state above T * indicating the formation of leg-directed orbital ordering at high temperatures, consistent with XAS measurements. At around T * , more isotropic electronic state forms, implying the formation of another orbital ordering (most likely d x 2 −y 2 orbitals). Approaching T N far below T * , we find the enhancement of nematic fluctuations that are distinctly different from the high-temperature fluctuations, suggesting the presence of competing orbital and magnetic orders in this system.

II. FIELD-ANGLE RESOLVED MAGNETORESISTANCE
MR was measured up to 17.5 T and field angle dependence is resolved by using rotating probes in a superconducting magnet with three geometric configurations labeled by G1, G2 and G2', where magnetic field H is rotated within out-of-ladder plane ((001) or xy plane), in-ladder plane ((100) or yz plane), and as-grown plane ((110)-leg plane), respectively. Hereafter, we define z and y axes along leg and rung directions, respectively, to discuss the electron d orbitals [upper insets in Figs. 1(a-c)]. From the view of space group Cmcm of BaFe 2 S 3 [15,16], our introduced x, y, z directions correspond to rung (a = y), plane (b = x), and leg (c = z), where a, b, c are the crystal axes. In all cases, the direction of current I is along the leg ([001] or z) direction. Magnetic field angle dependence of resistivity can be fitted by where θ is the angle from the rung axis as shown in Figs. 1(a-c) and its angle dependence of each term is determined through the restriction of orthorhombic symmetry [17]. Here, ρ H represents the contribution from the Hall effect possibly caused by its strongly anisotropic structure. The ρ M term is the amplitude of angle dependence of MR and can be unambiguously separated from the Hall contribution because of the parity difference to magnetic field. The ratio of ρ M /ρ 0 represents the degree of field-angle dependence reflecting the anisotropy of the electronic states in the field-rotation planes. The field-angle dependence of MR in G1 geometry starts to grow just above T N and MR shows maximum (minimum) when magnetic field H is applied perpendicular to the ladder plane (parallel to the rung direction). Magnetic moments aligned ferromagnetically along rung directions [1] can cause such angle dependence. The temperature dependence of the ratio ρ M /ρ 0 shows a sharp peak at T N as shown in Fig. 1(d). This suggests that antiferromagnetic fluctuations evolve towards T N , and such fluctuations start to develop below T * .
To study in-plane anisotropy, G2 and G2' geometries are used, in which field rotation plane is parallel to and slightly tilted from the ladder plane, respectively, but both measurements reproduce essentially the same behavior as shown in Fig. 1(e). Below T N , MR shows strong angle dependence similar to that of G1 geometry and takes a minimum when the H is along the rung direction. At higher temperatures above T * , however, ρ M /ρ 0 has the opposite sign and the magnitude of MR decreases with decreasing temperature and becomes very small around T * . These results indicate that BaFe 2 S 3 has an anisotropic electronic structure above T * , which is distinct from that in antiferromagnetic state below T N , and exhibits a phase transition around T * leading to more isotropic electronic state.

III. ELASTORESISTANCE
ER, resistance change induced by uniaxial strain, is a powerful probe of nematic fluctuations [18,19]. We apply this technique to BaFe 2 S 3 to study putative orbital ordering at T * . The ER measurements require ideal geometry of thin bar shapes to control strain through the direct attachment of samples to piezoelectric device. Despite their fragile nature of single crystals, these requirements are achieved by cutting the samples with a wire-saw and polishing them carefully. In this study, both strain and current are applied along the [001] leg direction [ Fig. 2(c)]. We measured ER under two geometries with different conduction planes: the as-grown (110) and ladder (010) planes, the latter of which is a more appropriate configuration to probe the nature of ladder nematicity. In both cases, strain along the leg (z) direction most effectively changes electronic states in ladder plane, because strain along other directions would significantly alter the separation between ladders, which makes it more difficult to extract the ER response in the ladder. ER in tetragonal materials linearly couples to nematic fluctuations, but here the relationships between ER and nematicity fluctuations should be modified by the orthorhombic crystal structure of BaFe 2 S 3 . Thus the Curie-Weiss analysis, which widely works for typical twodimensional iron-based superconductors [18][19][20], may not be valid in this case. However, as the leg distortions work as a conjugate field to ladder nematicity, strain-induced changes in resistance still couple in some form to nematic fluctuations even with distorted ladder structure.
Figures 2(a) and (b) show the strain response of resistance at low and high temperatures, respectively. The ER response is dominated by a linear slope with negative sign, in contrast to the metalization induced by hydrostatic pressure, which supports the successful control of anisotropic strain under little effect of symmetric strain. At high temperatures above T * , the magnitude of ER slope decreases with decreasing temperatures. Below T * , ER signal turns to grow and is enhanced towards T N . The temperature dependence of ER slope shown in Fig. 2(d) exhibits two anomalies; a kink around T N ∼ 117 K and a broad minimum around T * ∼ 180 K. Here, three measured samples including different experimental geometries show similar behaviors, implying that our ER data successfully capture the essential features of nematic fluctuations. The observed clear anomaly of ER provides strong evidence for the existence of electronic phase transition at T * . In iron-based superconductors with Fe square lattice, ER nematic signal increases from high temperature side towards the nematic transition, but in the present case it shows an opposite decreasing trend. The behavior taking a minimum in ER can be explained by assuming multiple components of nematic fluctuations and competitive relationship between them.

IV. DISCUSSION
Our measurements of field-angle dependent MR and ER consistently indicate a clear electronic change to a more isotropic state below T * . Both measurements indicate that the origin of T * is distinguished from antiferromagnetism below T N and most likely involves orbital degrees of freedoms as previously speculated [8,12]. Although orbital order in iron-based superconductors has been discussed in the view of nematic order which induces in-plane anisotropy, the T * phase in BaFe 2 S 3 rather reduces the anisotropy of electronic structure.
Next we discuss possible orbital states above and below T * . Strong field-angle dependence of MR and large ER coefficient at ∼ 270 K is likely originating from legdirected orbital ordering. From the view of low-energy states, effective two-band model with the d x 2 −y 2 orbitals and d xz hybridized with d xy orbitals have been proposed [21,22]. Since transport coefficients are generally sensitive to low-energy quasiparticle excitations, our transport anisotropy should be related to the imbalance between these two orbitals. Namely, the origin of the anisotropic state along the leg direction in high temperature region should come from d xz -like orbital order in this model. Furthermore, the reduction of ER towards T * from room temperature indicates a suppression of this leg-directed orbital fluctuations. This accomplishment of an isotropic electronic state can naturally be explained by the cancellation of anisotropy by the developments of another orbital d x 2 −y 2 , which is elongated along rung direction and competes with the d xz -like orbital. Here, the evolution of rung-directed order is gradual, so that it just compensates the anisotropy derived from leg-directed orbital order in range of T N < T < T * , which is summarized in a schematic phase diagram in Fig. 3. We note that if we assume d 3z 2 −r 2 orbital ordering in the place of d xz -like orbitals as previously suggested from XAS [12], discussion on gradual orbital switching to d x 2 −y 2 at T * is also applicable. These orbital decoupling behaviors emphasize the importance of the orbital selective Mott physics with Hund's coupling.
On the other hand, strong increase of ER below T * can be understood by the contribution from magnetic nematicity associated with stripe-type antiferromagnetism, which also coincides with the evolution of MR. It is contrasting that orbital order becomes isotropic around T * while the anisotropic antiferromagnetic order evolves. This is a quite unusual case because in most iron-based superconductors d xz orbitals and stripe-type antiferromagnetism give a cooperative formation of nematic orders with the help of spin-orbit coupling. Another exam-ple that is distinct from the present case can be found in a related iron-based ladder compound BaFe 2 Se 3 , which exhibits a structural distortion that supports the magnetic ordering [23,24]. However, d x 2 −y 2 orbital order hybridized along rung directions may support antiferromagnetic correlations along rungs, which is opposite to actual magnetic structure that magnetic moments align ferromagnetically along rungs [1]. We should also note that in BaFe 2 S 3 magnetic ordering temperature T N increases under pressure with a suppression of T * [14], which is consistent with the competing nature of these two orders. Such an unusual competition of two orders in BaFe 2 S 3 may provide hints to clarify the independent roles of spin and orbital degrees of freedoms.
The T * anomaly is originally characterized by the slight improvement of conductivity [8]. In fact, this can also be explained by the present orbital competition model discussed above: d xz -like bands are strongly correlated and the evolution of d x 2 −y 2 orbitals, which are more weakly correlated, leads to a decrease in resistivity. This resembles with incoherent-coherent crossover observed in hole-doped iron-based superconductors originating from the coexistence of localized and itinerant electrons due to orbital selective correlation effects [25]. The density matrix renormalization group analysis on BaFe 2 S 3 shows the orbital selective features in the perturbation response to hole doping: d xz -like orbitals have a tendency to localize while d x 2 −y 2 orbitals become coherent [22]. Experimentally, the coexistence of local and itinerant electrons is revealed by photoemission spectroscopy [26]. The orbital switching from d xz -like to d x 2 −y 2 orbital proposed here based on the electronic anisotropy is thus consistent with the incoherent-coherent crossover associated with orbital selective Mottness.

V. CONCLUSION
In summary, we performed field-angle resolved MR, and ER measurements in iron-ladder material BaFe 2 S 3 . Our results reveal strong anisotropy above T * indicating leg-directed orbital ordering such as d xz orbitals hybridized with d xy , and the existence of an electronic phase transition leading to an isotropic state at T * . A plausible microscopic origin of this phase transition is the orbital switching from d xz -like orbitals to d x 2 −y 2 coincided with incoherent-coherent crossover related to orbital selective Mott physics. Our proposed orbital order below T * is an exotic state different from nematic orders in iron-based superconductors. One of the intriguing features is the dichotomy between orbital and magnetic nematic instabilities, which provides a new avenue for studying the roles of orbital and magnetic fluctuations.

ACKNOLEDGEMENTS
We thank K. Takubo, T. Mizokawa, J. C. Palmstrom, T. Worasaran, and K. Izawa for helpful comments or fruitful discussions. We also thank C. W. Hicks for technical supports on designing the wire-saw, and N. Abe, Y. Tokunaga, and T. Arima for experimental supports on Xray diffraction measurements. This work was supported by Grants-in-Aid for Scientific Research (KAKENHI) (Grant Nos. JP15H02106, JP16H4007, JP16H04019, JP16K17732, JP17H05473, JP17H05474, JP17H06137, JP17J11382, JP18H01159, JP18H04302, JP18H05227, JP19H00649, and JP19H04683) and on Innovative Areas "Quantum Liquid Crystals" (Nos. JP19H05823, JP19H05824, and JP20H05162) from Japan Society for the Promotion of Science (JSPS). This work was partly performed at the High Field Laboratory for Superconducting Materials, Institute for Materials Research, Tohoku University (Project No.17H0202, 18H0204). X-ray diffraction measurements were partly performed using facilities of the Institute for Solid State Physics, the University of Tokyo.

A. APPENDIX:SAMPLE PREPERATION
High-quality single crystals of BaFe 2 S 3 were grown by the melt-growth method. The starting materials are elemental Ba shots, Fe powders, and S powders with the ratio of 1 : 2.1 : 3. The excess of iron is necessary to prevent the iron deficiency in the products as reported in Ref. [8]. The starting materials in a carbon crucible were sealed into an evacuated quartz ampoule. The ampoule was slowly heated up to 1373 K, kept for 48 hours, and slowly cooled to room temperature. We characterized samples by measuring the temperature dependence of resistivity and magnetic susceptibility. According to Hirata et al. [8], the sample with the best quality has the highest antiferromagnetic transition temperature (T N ) and the clearest anomaly in the resistivity at the possible orbital-order transition temperature (T * ). Our electrical resistivity and magnetic susceptibility data shown in Fig. 4 exhibit clear anomalies at T N = 117 K and T * = 180 K. These features are consistent with those reported in Ref. [8]. We can therefore conclude that the sample used in this study is appropriate to acquire the intrinsic properties of BaFe 2 S 3 . In the main text, we deal with the magnetoresistance under the rotating magnetic field. We here show the temperature dependence of magnetoresistance at fixed magnetic field direction along the rung direction (Fig.  5). The data were obtained from the temperature dependence of resistivity (ρ) of in a constant magnetic field of 0 T and 17.5 T. One can clearly see the suppressed magnetoresistance with a bump-like structure at around T * and enhanced negative magnetoresistance at T N , although the T * position estimated from ER data is slightly deviated from the top of the bump structure. The former and the latter correspond to the formation of the isotropic electronic states and the enhancement of magnetic fluctuations, respectively. The details of each transitions are discussed in the main text.