Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 26 February 2020
Sec. Inorganic Chemistry
This article is part of the Research Topic Flux Crystal Growth: From Fundamentals to Crystal Functionality View all 8 articles

Frustrated Magnetism in Triangular Lattice TlYbS2 Crystals Grown via Molten Flux

  • Oak Ridge National Laboratory, Materials Science and Technology Division, Oak Ridge, TN, United States

The triangular lattice compound TlYbS2 was prepared as large single crystals via a molten flux growth technique using sodium chloride. Anisotropic magnetic susceptibility measurements down to 0.4 K indicate a complete absence of long-range magnetic order. Despite this lack of long-range order, short-range antiferromagnetic interactions are evidenced through broad transitions, suggesting frustrated behavior. Variable magnetic field measurements reveal metamagnetic behavior at temperatures ≤2 K. Complex low temperature field-tunable magnetic behavior, in addition to no observable long-range order down to 0.4 K, suggest that TlYbS2 is a frustrated magnet and a possible quantum spin liquid candidate.

Introduction

The variety of basic and applied properties arising from geometrically frustrated magnets continue to motivate the study of structures with two-dimensionally layered triangular lattices prone to antiferromagnetic interactions (Chubokov and Golosov, 1991; Lee et al., 2006; Shen et al., 2016; Li et al., 2017a; Zhu et al., 2018; Bordelon et al., 2019; Ranjith et al., 2019). Such triangular lattices restrict the number of available spin degrees of freedom, resulting in quantum fluctuations that can produce degenerate ground states (Savary and Balents, 2016). Frustrated antiferromagnets with degenerate ground states have garnered significant interest for their potential as quantum spin liquid (QSL) candidates, a state characterized by dynamic entangled spins, exhibiting no long-range magnetic order, even at 0 K (Anderson, 1973; Balents, 2010; Hu et al., 2015; Starykh, 2015; Yamamoto et al., 2015; Dun et al., 2016; Savary and Balents, 2016; Shen et al., 2016; Xu et al., 2016; Li et al., 2017b; Paddison et al., 2017; Baenitz et al., 2018). While QSL candidates with 3d ions exhibit weak spin-orbit coupling (SOC) (Lee and Lee, 2005; Helton et al., 2007; Yoshida et al., 2009; Yamashita et al., 2010; Zhou et al., 2011; Shirata et al., 2012), the presence of stronger SOC in 5d/4f ions, on par with the energy scale of crystal electric field effects and the coulomb interaction U, further enhances the frustration via entangled spin and orbital degrees of freedom, and has thus shifted the search for new QSL materials to contain these heavier lanthanides (Ln) (Okamoto et al., 2007; Curnoe, 2008; Gardner et al., 2010; Onoda and Tanaka, 2010; Applegate et al., 2012; Hu et al., 2015; Li et al., 2016; Lu et al., 2017; Laconis et al., 2018; Liu et al., 2018).

Recently, Yb(III) containing compounds such as NaYbO2 (Bordelon et al., 2019; Ranjith et al., 2019), NaYbS2 (Baenitz et al., 2018), and YbMgGaO4 (Li et al., 2015, 2017a,b; Xu et al., 2016; Paddison et al., 2017) have been presented as QSL candidates, all crystallizing in a layered triangular lattice of trigonal space group R3¯m. Due to the odd number of 4f electrons and strong SOC, these materials behave as effective spin Jeff = ½. A similar family of 4f-containing delafossites, of the general formula A(I)Ln(III)Ch(II) [A = Na, Rb, K; Ch = O, S, Se, Te] (Liu et al., 2018), has also been proposed as a promising candidate. Delafossites are often free from crystallographic site-mixing, unlike the more commonly studied YbMgGaO4 that can mimic QSL behavior by eliminating long-range order through disorder (Li et al., 2017b; Zhu et al., 2017). Additionally, the modular nature of delafossite structures allows for the possibility of differences in crystal structure as a function of the ratio of ionic radii; this has been reported to result in changes in triangular lattice layer stacking, such as ABAB stacking in the hexagonal P63/mmc or ABCABC layer stacking in the trigonal R3¯m. Reports of possible inter-layer interactions on the highly sensitive magnetic ground state of such systems makes the delafossite structure advantageous to study, as such subtle interactions can be probed as a function of selective ion control (Yamamoto et al., 2015).

Despite the structural modularity and promise of QSL candidacy in such delafossite structures, the limited availability of large single crystals to study electronic and magnetic anisotropy serves as motivation for this work. Herein we discuss the flux crystal growth, structure determination, and magnetic property measurements of TlYbS2, which crystallizes in the hexagonal space group P63/mmc. This study allows for proper structure elucidation of TlYbS2 single crystals, contrary to the results of polycrystalline powders (Duczmal and Pawlak, 1994), in addition to reporting novel anisotropic magnetization results that were inaccessible on polycrystalline samples (Duczmal and Pawlak, 1994).

Experimental

Synthesis

The TlYbS2 compound was synthesized using a two-step method comprised of (1) producing the powder form via traditional solid-state synthesis, followed by (2) crystallizing the precursor powder using molten flux growth via sodium chloride. Solid pieces of Yb metal (REacton, 99.99%), S (Puratronic, 99.9995%), and Tl (REacton, 99.99%) were all stored in a glove box. Sodium chloride (Alfa Aesar, 99.999%) was dried in an oven overnight at 300°C and stored in a desiccator prior to use.

For the solid-state synthesis, 2.0 mmol of Tl and Yb were added to 4.0 mmol of S in an alumina crucible with a loose-fitting alumina frit and a second, inverted, alumina crucible on top. This second crucible was used to assist in catching any trace amounts of volatilized Tl or S. This setup of alumina crucibles was loaded and sealed inside an evacuated silica tube, with a small amount of quartz wool at the bottom to prevent cracking due to differences in thermal expansion (Figure 1A). The sealed silica tube was heated to 300°C at a rate of 10°C/h, dwelled for 24 h, ramped to 800°C at a rate of 10°C/h, dwelled at 72 h, and then the reaction was allowed to cool by shutting off the furnace.

FIGURE 1
www.frontiersin.org

Figure 1. (A) Experimental setup consists of an alumina crucible housing the starting materials with a loose-fitting frit and a second inverted alumina crucible, shown in the inset, both held within a sealed evacuated quartz tube. (B) Optical image of a typical crystal of TlYbS2, showing its red color. The c axis is out of the plane of the crystal.

Single crystals were produced by loading 0.87 mmol of polycrystalline TlYbS2 and a ten-fold excess (by mass) of NaCl (40.47 mmol) into sealed evacuated silica tubing. The reaction was heated to 850°C at a rate of 60°C/h, dwelled at 504 h (i.e. 3 weeks), and then cooled by shutting off the furnace. The resulting red crystals were mechanically separated from the remaining TlYbS2 powder and vacuum filtered using ethanol to remove any surface impurities (Figure 1B). The crystallographic c-axis is out of the plane of the paper. The purity of the polycrystalline powder (first step) and resulting crystals (second step) were both determined by powder X-ray diffraction (PXRD) using a PANalytical X'Pert Pro MPD diffractometer with Cu Kα1 radiation (λ = 1.5418 Å), shown in Figure 2.

FIGURE 2
www.frontiersin.org

Figure 2. Powder X-ray diffraction pattern of solid-state produced polycrystalline TlYbS2, shown in blue, overlaid with a diffraction pattern of ground single crystals of TlYbS2, shown in green.

Structure

The structure determination of the single crystals of TlYbS2 was performed on a Bruker Quest D8 single crystal X-ray diffractometer (SXRD). Single crystals were first sonicated in acetone to remove any surface impurities. The data were collected at room temperature utilizing a Mo Kα radiation (λ = 0.71073 Å). The crystal diffraction images were collected using Φ and ω-scans. The diffractometer was equipped with an Incoatec IμS source using the APEX III software suite for data setup, collection, and processing (Bruker, 2015). The structure was resolved using intrinsic phasing and full-matrix least square methods with refinement on F2 using the SHELXTL software suite (Sheldrick, 2008). All atoms were first refined with isotropic thermal displacement parameters and then refined anisotropically. Crystallographic information can be found in Tables 13.

TABLE 1
www.frontiersin.org

Table 1. Crystallographic data for TlYbS2 determined by single crystal X-ray diffraction.

TABLE 2
www.frontiersin.org

Table 2. Atomic coordinates and equivalent isotropic displacement parameters (Å2) for TlYbS2.

TABLE 3
www.frontiersin.org

Table 3. Bond distances (Å) and angles (°) for TlYbS2.

Energy dispersive spectroscopy (EDS) was performed using a Hitachi S-3400 scanning electron microscope equipped with an OXFORD EDX microprobe to confirm the elemental composition in the single crystal sample. The presence of Tl, Yb and S were verified, and the absence of extraneous elements such as Na and Cl were confirmed. EDS data can be found in Table 4.

TABLE 4
www.frontiersin.org

Table 4. Energy dispersive spectroscopy data.

Magnetic Susceptibility

Physical properties were measured using the Quantum Design Superconducting quantum interference device (SQUID) Magnetic Properties Measurement System (MPMS). Two thin plates of crystals were stacked with a total mass of 0.36 mg using vacuum grease, for each measurement above 2 K. The magnetization measurements made below 2 K were measured using the Quantum Design iHe3 option on four stacked single crystals (0.5 mg total) using vacuum grease. The crystals were aligned in two separate orientations: one set of measurements aligned the crystals such that the applied magnetic field was perpendicular to the c-axis, and second set where the applied field was parallel with c. The thermometer calibration was done in zero field when the applied field was parallel with c.

Discussion

Synthesis

The scarcity of known single crystal growth methods for lanthanide containing delafossites, particularly those employing the use of salt flux, has made the optimization and characterization of large single crystals of the titled composition challenging (Stowe, 1997). Fortunately, there is significant literature evidence demonstrating that the use of molten flux as a growth medium is a robust method, potentially capable of crystallizing nearly every element combination on the periodic table (Bugaris and zur Loye, 2012). The selection of an alkali halide flux was guided in part by recent reports of similar fluxes crystallizing compositions containing lanthanides and/or chalcogenides (Klepov and zur Loye, 2018; Tsujimoto et al., 2018; Usman et al., 2019a,b). Additionally, amongst the few delafossite-type structures reported as single crystals, synthesis typically involves the use of reactive alkali fluxes, such that the alkali metal in the flux incorporates into the final product, such as the use of KCl for KErSe2 (Xing et al., 2019a) or NaCl for NaYbS2 (Baenitz et al., 2018). For producing TlYbS2 crystals, the use of TlCl as a flux medium was ruled out due to the low solubility prohibiting ease of crystal separation upon completion of the reaction. We tried CaCl2 as a flux, and although single crystals were produced, they were of poor quality. Attempts to improve the quality of the crystals by introducing cooling rates also failed: a variety of polycrystalline powders and no crystals were present, indicating TlYbS2 may be a metastable kinetic phase that is “trapped” via quenching. Ultimately, the use of NaCl as a flux, in addition to quenching the reaction upon the completion of the prolonged dwelling period, resulted in high quality, large (1 mm+) single crystals that were suitable for structural and anisotropic magnetization studies.

Structure

The compound TlYbS2 was first reported as a polycrystalline powder crystallizing in the trigonal space group R3¯m, commonly referred to as the α-NaFeO2 structure, with lattice parameters a = 3.935 Å and c = 22.47 Å (Duczmal and Pawlak, 1994). This two dimensional layered structure is built from two distinct triangular lattice layers that alternate along the c axis. The first layer is built from a network of edge-shared octahedral YbS6 units that adopts an ABCABC stacking pattern (Figure 3). A second non-magnetic triangular lattice layer of edge-shared octahedral TlS6 units resides between these layers. As a result of the availability of single crystals of the titled composition, we thoroughly investigated the nuclear structure. Careful analysis revealed that the grown single crystals of TlYbS2 adopt the hexagonal P63/mmc β-RbScO2 structure type. In this hexagonal structure, Tl (Wyckoff 2c), Yb (2b), and S (4f) occupy the special positions with site symmetries of −6m2, −3m, and 3m, respectively. This is in contrast to the previously reported trigonal R3¯m structure for polycrystalline powders of TlYbS2 where Tl (3b) and Yb (3a) are in the −3m position, and S (6c) is in the 3m special position. The primary difference between the two structures is best understood by the number and stacking sequence of the triangular lattice layers. In the hexagonal β-RbScO2 structure, one unit cell is built from three YbS6 layers and two TlS6 layers, whereas the trigonal α-NaFeO2 structure unit cell is built from four YbS6 layers and three TlS6 layers. Additionally, the hexagonal β-RbScO2 structure adopts a higher symmetry ABAB triangular lattice layer stacking, in contrast to the ABCABC layer stacking in the trigonal α-NaFeO2 type structure. This result is consistent with the trend observed for smaller A-site ions in the delafossite structure crystallizing in the trigonal system, such as NaYbS2 (Baenitz et al., 2018) and larger A-site ions crystallizing in the hexagonal system, such as CsYbSe2 (Xing et al., 2019b). Comparison of powder X-ray diffraction patterns of solid-state produced polycrystalline powder, and single crystals grown via molten flux, of TlYbS2 overlay well (Figure 2) with slight differences in peak intensity and crystallinity. Although polymorphism between the polycrystals and single crystals should not be completely ruled out, slight differences in peak intensity may result from preferred orientation in such a highly anisotropic nuclear structure.

FIGURE 3
www.frontiersin.org

Figure 3. Polyhedral representation of the delafossite structure, with the hexagonal P63/mmc variant (shown left) and the trigonal R3¯m variant (shown right). The hexagonal structure is built of three distinct triangular layers of edge-shared YbS6 octahedra (Yb shown in gray, S in red) that stack in an ABAB order, meanwhile the trigonal structure is built from four distinct triangular layers of edge-shared YbS6 octahedra that stack in an ABCABC sequence. Non-magnetic TI atoms are shown as blue spheres.

The primary interest in studying this structure type is the triangular lattice that extends along the ab plane, serving as an ideal host for geometric frustration (Figure 4). The idealized hexagonal network of Yb atoms does not allow for a purely antiparallel configuration of spins, leading to enhanced quantum fluctuations that are of interest to study. Careful attention was taken to select Tl as the A site in the delafossite structure, since it is larger than Yb, both to prevent crystallographic site-mixing, and to maximize the interlayer distances, thereby minimizing inter-layer interactions via the mediating cation. The reactive-flux nature of NaCl in the growth of similar delafossite structures also guided the selection of the large cation Tl, in hopes that Na would be too small to occupy the same site. Additionally, a non-magnetic A site was selected to further simplify study of any resulting magnetic properties.

FIGURE 4
www.frontiersin.org

Figure 4. Polyhedral representation of the edge-shared octahedral network of YbS6 extending in the ab plane for the composition TlYbS2. The planar triangular lattice is an ideal host for geometric frustration for an antiparallel alignment of spins, as shown in the zoomed-in area. Yb atoms are shown in gray, and S atoms are shown in red.

Magnetic Susceptibility

The temperature dependence of the magnetization perpendicular to the applied field (in the ab plane) at 1 T indicates a broad deviation from Curie-Weiss behavior at 50 K (Figure 5). A Curie-Weiss fit of the high temperature section (200–350 K) of the inverse susceptibility taken at 1 T indicates a Weiss constant of θCW = −101.3 K and an effective moment of 4.57 μB, in good agreement with the expected J = 7/2 (4.54 μB) moment (Hashimoto et al., 2003). The large negative Weiss temperature indicates strong antiferromagnetic interactions. A fit of the low temperature (2–50 K) taken at 1 T indicates a smaller Weiss constant of θCW = −19.5 K and an effective moment of 3.41 μB. To better understand the significant difference between high temperature moment [J = 7/2 (4.54 μB)] and low temperature moment (3.41 μB, unknown J), a literature search for analogous results was conducted. Similar measurements (also taken at low temperature and field) on structurally related (α-NaFeO2-type) single crystals of NaYbS2 made in-plane (ab plane) report comparable values of θCW = −13.5 K and μeff = 3.2 μB (Baenitz et al., 2018). The slight deviation in value between the two structures may originate from differences between Na and Tl ionic radii, [1.02 and 1.50 Å, respectively] causing the structural shift from R3¯m (closer layers) to P63/mmc (greater inter-layer distances). As a result of the larger non-magnetic cation (Tl) further separating the triangular lattice layers, TlYbS2 more closely resembles a pure two-dimensional structure, compared to the Na analog. The report on NaYbS2 describes the effective magnetic moment (μeff = 3.2 μB) to originate from strong SOC and ab plane anisotropy, treating the studied material as a Jeff = ½ system according to the two-dimensional spin ½ triangular lattice Heisenberg antiferromagnet model (Anderson, 1973; Huse and Elser, 1988). This Heisenberg model has been used to understand previously reported isostructural frustrated systems in which all magnetic ions are 120° from one another, which is the case for NaYbS2, and thus is reasonable as a starting hypothesis for a possible magnetic model to explain the behavior of TlYbS2 (Liu et al., 2018). To assess the Jeff = ½ magnetic model, describing the proposed low temperature behavior for the Yb atoms in NaYbS2, Baenitz et al. conducted electron-spin-resonance (ESR) measurements on single crystals and found the g-factor to be highly anisotropic along the two crystallographic directions (gab = 3.19; gc = 0.57), an expected consequence of the triangular lattice layer (Baenitz et al., 2018). Utilizing the experimentally determined gab-factor, the Jeff = ½ magnetic model hypothesis, and considering Equation (1), an expected moment (μeff) was calculated and found to be close to their experimentally observed moment of 3.2 μB, suggesting that NaYbS2 behaves as Jeff = ½ triangular lattice Heisenberg antiferromagnet.

μeff=gJ (J+1)    (1)

Based on the success of the analysis used for NaYbS2, a similar method was employed for TlYbS2 to understand the low temperature magnetic moment. Within the triangular lattice Heisenberg AFM model, the magnetization should plateau at ~1/3 the expected saturation magnetization, according to Equation (2). In this equation, the saturation magnetization

ms =J ·g    (2)

(ms) is equal to the product of the g-factor and the total angular momentum (J). By analyzing TlYbS2 in the same way NaYbS2 was treated using Equation (1), an anisotropic g-factor of gab = 3.94 for TlYbS2 is obtained. This can be used in Equation (2) to obtain a saturation magnetization saturation (1.97 μB). Taking the ratio of the moment (0.74 μB) corresponding to the inflection point observed in field dependent magnetization measurements taken at 0.42 K (Figure 6) with the ms discussed above yields a value of 0.37, in good agreement with the expected 0.33. It should be noted that in the absence of a sufficiently strong magnetic field, a full plateau is not observed, but the inflection point is the onset of such a plateau. Considering the triangular lattice orientation of the Yb atoms in the ab plane, as shown in Figure 4, and the good agreement of 0.37 to the expected 0.33 as per the triangular lattice Heisenberg AFM model, it is reasonable to propose a magnetic model in which TlYbS2 behaves as a Jeff = ½ system in the low temperature (below 50 K) regime. This finding is further supported by previous reports of subtle changes in slope for inverse magnetic susceptibility plots, such as that observed in Figure 5, as being attributed to gradual transitions toward an isolated Kramers doublet ground state (Ranjith et al., 2019).

FIGURE 5
www.frontiersin.org

Figure 5. Temperature dependence of the magnetic susceptibility for single crystals of TlYbS2, oriented H⊥c under an applied field of 1 T. Field cooled data is shown in blue and inverse magnetic susceptibility data shown in red. Solid black lines denote the high/low temperature fits of the inverse susceptibility data.

FIGURE 6
www.frontiersin.org

Figure 6. Field dependence of the magnetization for single crystals of TlYbS2 oriented H⊥c over a range of 0 to 7 T at 0.42 K.

The broad deviation from Curie-Weiss behavior at 50 K was investigated further by measuring inverse temperature dependent susceptibility perpendicular to the field at low temperatures (2–50 K) at fields from 0.1 to 5 T (Figure 7A). A clear deviation from linearity can be seen in the 3 T data at 4 K, with a more pronounced deviation at the same temperature for the 5 T data. Supplementary temperature dependent susceptibility measurements (Figure 7B) taken at low temperatures (0.4–2 K) and variable fields (0.2–7 T) reveal a gradual inversion of slope, with the 7 T data appearing as a nearly flat line. The absence of saturation in the field dependence magnetization measurements made in Figure 6 at 7 T suggests the nearly linear susceptibility at 7 T in Figure 7 does not correspond to true saturation.

FIGURE 7
www.frontiersin.org

Figure 7. (A) Temperature dependence (2–8 K) of inverse magnetic susceptibility for single crystals of TlYbS2 oriented H⊥c under applied fields ranging from 0.1 to 5 T. A wider temperature range (2–50 K) of the sample plot is shown in the inset. (B) Temperature dependence (0.4–1.8 K) of susceptibility under applied fields from 0.2 to 7 T.

Similar, albeit more subtle, features were observed in low temperature (0.4–1.4 K) susceptibility measurements (Figure 8) at variable fields (0.5–3 T), along the c axis. The presence of such very subtle transitions and an even weaker magnetic susceptibility perpendicular to c are unsurprising, as magnetic interactions are expected to reside primarily within the triangular lattice ab plane. The large size of non-magnetic Tl atoms that separate the triangular lattices reasonably limit inter-layer magnetic interaction pathways. Although a complete description of this complex behavior is not possible with the current data, similar behavior has been recently reported in Os0.55Cl2 (McGuire et al., 2019) and MErSe2 (M = Na, K) (Xing et al., 2019a). In both reports, no long-range ordering and a spin-liquid-like behavior is observed in triangular lattice frameworks similar to TlYbS2, suggesting that it may also be a candidate for this exotic behavior.

FIGURE 8
www.frontiersin.org

Figure 8. Temperature dependence (0.4–1.4 K) of magnetic susceptibility for single crystals of TlYbS2 oriented H||c under applied fields from 0.5 to 3 T.

Conclusions

In summary, for the first time large single crystals of TlYbS2 were prepared from a molten flux growth technique, and a reinvestigation of the nuclear structure was conducted. The single crystal structure was determined to crystallize in the hexagonal P63/mmc β-RbScO2 type structure, in contrast to the previously reported trigonal R3¯m α-NaFeO2 structure for polycrystalline samples. Anisotropic measurements along the ab plane and along c ranging from 0.42–350 K and 0.1–7 T displayed clear short-range magnetism and metamagnetic behavior. The complete absence of long-range order and consistency of TlYbS2 magnetic behavior to the two-dimensional triangular lattice Heisenberg antiferromagnet model indicates its possible candidacy as a quantum spin liquid.

Data Availability Statement

The datasets generated for this study can be found in the Cambridge Crystallographic Data Centre (https://www.ccdc.cam.ac.uk/structures/) under the identifier 1965470. Alternatively, email data_request@ccdc.cam.ac.uk, or by contact The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB21EZ, UK; fax +44 1223 336033.

Author Contributions

TF was responsible for making the materials in polycrystalline and single crystal forms; additionally, he drafted this manuscript. JX was responsible for the magnetic measurements taken, and LS was responsible for the structural characterization. AS oversaw the experiments and helped in finalizing this manuscript.

Funding

This research, conducted at Oak Ridge National Laboratory (ORNL), was supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), Materials Sciences and Engineering Division (MSE).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Radu Custelcean is gratefully acknowledged for permitted use of his single crystal X-ray diffractometer.

References

Anderson, P. W. (1973). Resonating valence bonds: a new kind of insulator. Mater. Res. Bull. 8:153. doi: 10.1016/0025-5408(73)90167-0

CrossRef Full Text | Google Scholar

Applegate, R., Hayre, N. R., Singh, R. R. P., Lin, T., Day, A. G. R., and Gingras, M. J. P. (2012). Vindication of Yb2Ti2O7 as a model exchange quantum spin ice. Phys. Rev. Lett. 109:097205. doi: 10.1103/PhysRevLett.109.097205

PubMed Abstract | CrossRef Full Text | Google Scholar

Baenitz, M., Schlender, P., Sichelschmidt, J., Onykiienko, Y. A., Zangeneh, Z., Ranjith, K. M., et al. (2018). NaYbS2: a planar spin-1/2 triangular-lattice magnet and putative spin liquid. Phys. Rev. B. 98:220409. doi: 10.1103/PhysRevB.98.220409

CrossRef Full Text | Google Scholar

Balents, L. (2010). Spin liquids in frustrated magnets. Nature 464, 199–208. doi: 10.1038/nature08917

PubMed Abstract | CrossRef Full Text | Google Scholar

Bordelon, M. M., Kenney, E., Liu, C., Hogan, T., Posthuma, L., Kavand, M., et al. (2019). Field-tunable quantum disordered ground state in the triangular-lattice antiferromagnet NaYbO2. Nat. Phys. 15, 1058–1064. doi: 10.1038/s41567-019-0594-5

CrossRef Full Text | Google Scholar

Bruker (2015). APEX3, SAINT+, and SADABS. Madison, WI: Bruker AXS Inc.

Google Scholar

Bugaris, D. E., and zur Loye, H.-C. (2012). Materials discovery by flux crystal growth: quaternary and higher order oxides. Angew. Chem. Int. Ed. 51, 3780–3811. doi: 10.1002/anie.201102676

PubMed Abstract | CrossRef Full Text | Google Scholar

Chubokov, A. V., and Golosov, D. I. (1991). Quantum theory of an antiferromagnet on a triangular lattice in a magnetic field. J. Phys. Condens. Matter. 3:69. doi: 10.1088/0953-8984/3/1/005

CrossRef Full Text | Google Scholar

Curnoe, S. H. (2008). Structural disortion and the spin liquid state in Tb2Ti2O7. Phys. Rev. B. 78:094418. doi: 10.1103/PhysRevB.78.094418

CrossRef Full Text | Google Scholar

Duczmal, M., and Pawlak, L. (1994). Magnetic properties of TlLnS2 compounds (Ln = Nd, Gd, Dy, Er, and Yb). J. Alloys Compd. 209, 271–274. doi: 10.1016/0925-8388(94)91112-6

CrossRef Full Text | Google Scholar

Dun, Z. L., Trinh, J., Li, K., Lee, M., Chen, K. W., Baumbach, R., et al. (2016). Magnetic ground states of the rare-earth tripod kagome lattice Mg2RE2Sb3O14 (RE = Gd, Dy, Er). Phys. Rev. Lett. 116:157201. doi: 10.1103/PhysRevLett.116.157201

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardner, J. S., Gingras, M. J. P., and Greedan, J. E. (2010). Magnetic pyrochlore oxides. Rev. Mod. Phys. 82:53. doi: 10.1103/RevModPhys.82.53

CrossRef Full Text | Google Scholar

Hashimoto, Y., Wakeshima, M., and Hinatsu, Y. (2003). Magnetic properties of ternary sodium oxides NaLnO2 (Ln - rare earths). J. Solid State Chem. 176, 266–272. doi: 10.1016/j.jssc.2003.08.001

CrossRef Full Text | Google Scholar

Helton, J. S., Matan, K., Shores, M. P., Nytko, E. A., Bartlett, B. M., Yoshida, Y., et al. (2007). Spin dynamics of the spin−1/2 Kagome Lattice antiferromagnet ZnCu3(OH)6Cl2. Phys. Rev. Lett. 98:107204. doi: 10.1103/PhysRevLett.98.107204

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W.-J., Gong, S.-S., Zhu, W., and Sheng, D. N. (2015). Competing spin-liquid states in the spin-1/2 Heisenberg model on the triangular lattice. Phys. Rev. B. 92:140403. doi: 10.1103/PhysRevB.92.140403

CrossRef Full Text | Google Scholar

Huse, D. A., and Elser, V. (1988). Simple variational wave functions for two-dimensional heisenberg spin-1/2 antiferromagnets. Phys. Rev. Lett. 60:2531. doi: 10.1103/PhysRevLett.60.2531

PubMed Abstract | CrossRef Full Text | Google Scholar

Klepov, V. V., and zur Loye, H.-C. (2018). Complex topologies from simple building blocks: uranium (IV) thiophosphates. Inorg. Chem. 57, 11175–11183. doi: 10.1021/acs.inorgchem.8b01733

PubMed Abstract | CrossRef Full Text | Google Scholar

Laconis, J., Liu, C., Halaz, G. B., and Balents, L. (2018). Spin liquid versus spin orbit coupling on the triangular lattice. Sci. Post Phys. 4:003. doi: 10.21468/SciPostPhys.4.1.003

CrossRef Full Text | Google Scholar

Lee, P. A., Nagaosa, N., and Wen, X.-G. (2006). Doping a Mott insulator: physics of high-temperature superconductivity. Rev. Mod. Phys. 78:17. doi: 10.1103/RevModPhys.78.17

CrossRef Full Text | Google Scholar

Lee, S.-S., and Lee, P. A. (2005). U(1) gauge theory of the hubbard model: spin liquid states and possible application to k-(BEDT-TTF)2Cu2(CN)3. Phys. Rev. Lett. 95:036403. doi: 10.1103/PhysRevLett.95.036403

CrossRef Full Text | Google Scholar

Li, Y., Adroja, D., Bewley, R. I., Voneshen, D., Tsirlin, A. A., Gegenwart, P., et al. (2017a). Crystalline electric-field randomness in the triangular lattice spin-liquid YbMgGaO4. Phys. Rev. Lett. 118:107202. doi: 10.1103/PhysRevLett.118.107202

CrossRef Full Text | Google Scholar

Li, Y., Adroja, D., Voneshen, D., Bewley, R. I., Zhang, Q., Tsirlin, A. A., et al. (2017b). Nearest-neighbour resonating valence bonds in YbMgGaO4 Nat. Commun. 8:15814. doi: 10.1038/ncomms15814

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Chen, G., Tong, W., Pi, L., Liu, J., Yang, Z., et al. (2015). Rare-earth triangular lattice spin liquid: a single-crystal study of YbMgGaO4. Phys. Rev. Lett. 115:167203. doi: 10.1103/PhysRevLett.115.167203

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y.-D., Wang, X., and Chen, G. (2016). Anisotropic spin model of strong spin-orbit-coupled triangular antiferromagnets. Phys. Rev. B. 94:035107. doi: 10.1103/PhysRevB.94.035107

CrossRef Full Text | Google Scholar

Liu, W., Zhang, Z., Ji, J., Liu, Y., Li, J., Wang, X., et al. (2018). Rare-earth calcogenides: a large family of triangular lattice spin liquid candidates. Chin. Phys. Lett. 35:117501. doi: 10.1088/0256-307X/35/11/117501

CrossRef Full Text | Google Scholar

Lu, L., Song, M., Liu, W., Reyes, A. P., Kuhns, P., Lee, H. O., et al. (2017). Magnetism and local symmetry breaking in a Mott insulator with strong spin orbit interactions. Nat. Commun. 8:14407. doi: 10.1038/ncomms14407

PubMed Abstract | CrossRef Full Text | Google Scholar

McGuire, M. A., Zheng, Q., Yan, J., and Sales, B. C. (2019). Chemical disorder and spin-liquid-like magnetism in the van der Waals layered 5d transition metal halide Os0.55Cl2. Phys. Rev. Lett. 99:214402. doi: 10.1103/PhysRevB.99.214402

CrossRef Full Text | Google Scholar

Okamoto, Y., Nohara, M., Aruga-Katori, H., and Takagi, H. (2007). Spin-liquid state in the S = 1/2 hyperkagome antiferromagnet Na4Ir3O8. Phys. Rev. Lett. 99:137207. doi: 10.1103/PhysRevLett.99.137207

PubMed Abstract | CrossRef Full Text | Google Scholar

Onoda, S., and Tanaka, Y. (2010). Quantum melting of spin ice: emergent cooperative quadrupole and chirality. Phys. Rev. Lett. 105:047201. doi: 10.1103/PhysRevLett.105.047201

PubMed Abstract | CrossRef Full Text | Google Scholar

Paddison, J. A. M., Daum, M., Dun, Z., Ehlers, G., Liu, Y., Stone, M. B., et al. (2017). Continuous excitations of the triangular-lattice quantum spin liquid YbMgGaO4. Nat. Phys. 13, 117–122. doi: 10.1038/nphys3971

CrossRef Full Text | Google Scholar

Ranjith, K. M., Dmytriieva, D. D., Khim, S., Sichelschmidt, J., Luther, S., Ehlers, D., et al. (2019). Field-induced instability of the quantum spin liquid ground state in the Jeff = 1/2 triangular-lattice compound NaYbO2. Phys. Rev. B. 99:180401. doi: 10.1103/PhysRevB.99.180401

CrossRef Full Text | Google Scholar

Savary, L., and Balents, L. (2016). Quantum spin liquids: a review. Rep. Prog. Phys. 80:016502. doi: 10.1088/0034-4885/80/1/016502

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheldrick, G. M. (2008). A short history of SHELX. Acta Crystallogr. A Found. Crystallogr. 64, 112–122. doi: 10.1107/S0108767307043930

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Y., Dong Li, Y., Wo, H., Li, Y., Shen, S., Pan, B., et al. (2016). Evidence for a spinon Fermi surface in a triangular-lattice quantum-spin-liquid candidate. Nature 540, 559–562. doi: 10.1038/nature20614

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirata, Y., Tanaka, H., Matsuo, A., and Kindo, K. (2012). Experimental realization of a spin-1/2 triangular-lattice Heisenberg antiferromagnet. Phys. Rev. Lett. 108:057205. doi: 10.1103/PhysRevLett.108.057205

PubMed Abstract | CrossRef Full Text | Google Scholar

Starykh, O. A. (2015). Unusual ordered phases of highly frustrated magnetism: a review. Rep. Prog. Phys. 78:052502. doi: 10.1088/0034-4885/78/5/052502

CrossRef Full Text | Google Scholar

Stowe, K. (1997). The crystal strctures of ErSeI and NaErSe2 Z. Anog. All. Chem. 623, 1639–1643. doi: 10.1002/zaac.19976231026

CrossRef Full Text | Google Scholar

Tsujimoto, Y., Juillerat, C. A., Zhang, W., Fujii, K., Yashima, M., Halasyamani, P. S., et al. (2018). Function of tetrahedral ZnS3O building blocks in the formation of SrZn2S2O: a phase matchable polar oxysulfide with a large second harmonic generation response. Chem. Mater. 30, 6486–6493. doi: 10.1021/acs.chemmater.8b02967

CrossRef Full Text | Google Scholar

Usman, M., Morrison, G., and zur Loye, H.-C. (2019b). La2USe3S2: a serendipiously grown lanthanide/actinide chalcogenide from a eutectic halide flux. J. Chem. Cryst. 49, 169–173. doi: 10.1007/s10870-019-00780-x

CrossRef Full Text | Google Scholar

Usman, M., Smith, M. D., Morrison, G., Klepov, V. V., Zhang, W., Halasyamani, P. S., et al. (2019a). Molten alkali halide flux growth of an extensive family of noncentrosymmetric rare earth sulfides: structure and magnetic optical (SHG) properties. Inorg. Chem. 58, 8541–8550. doi: 10.1021/acs.inorgchem.9b00849

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, J., Sanjeewa, K. D., Kim, J., Stewart, G. R., Du, M.-H., Reboredo, F. A., et al. (2019b). Crystal synthesis and frustrated magnetism in triangular lattice CsRESe2 (RE = La - Lu): quantum spin liquid candidates CsCeSe2 and CsYbSe2. ACS Materials Lett. 2, 71–75. doi: 10.1021/acsmaterialslett.9b00464

CrossRef Full Text | Google Scholar

Xing, J., Sanjeewa, L. D., Kim, J., Meier, W. R., May, A. F., Zheng, Q., et al. (2019a). Synthesis, magnetization and heat capacity of triangular lattice materials NaErSe2 and KErSe2. Phys. Rev. Mater. 3:114413. doi: 10.1103/PhysRevMaterials.3.114413

CrossRef Full Text | Google Scholar

Xu, Y., Zhang, J., Li, Y. S., Yu, Y. J., Hong, X. C., Zhang, Q.M., et al. (2016). Absence of magnetic thermal conductivity in the quantum spin-liquid candidate YbMgGaO4. Phys. Rev. Lett. 117:267202. doi: 10.1103/PhysRevLett.117.267202

CrossRef Full Text | Google Scholar

Yamamoto, D., Marmorini, G., and Danshita, I. (2015). Microscopic model calculatons for the magnetization process of layered triangular-lattice quantum antiferromagnets. Phys. Rev. Lett. 114:027201. doi: 10.1103/PhysRevLett.114.027201

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamashita, M., Nakkata, N., Senshu, Y., Nagata, M., Yamamoto, H. M., Kato, R., et al. (2010). Highly mobile gapless excitations in a two-dimensional candidate quantum spin liquid. Science. 328, 1246–1248. doi: 10.1126/science.1188200

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Takigawa, M., Yoshida, H., Okamoto, Y., and Hiroi, Z. (2009). Phase diagram and spin dynamics in volborthite with a distorted kagome lattice. Phys. Rev. Lett. 103:077207. doi: 10.1103/PhysRevLett.103.077207

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H. D., Choi, E. S., Li, G., Balicas, L., Wiebe, C. R., Qiu, Y., et al. (2011). Spin liquid state in the S = 1/2 triangular lattice Ba3CuSb2O9. Phys. Rev. Lett. 106:147204. doi: 10.1103/PhysRevLett.106.147204

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Maksimov, P. A., White, S. R., and Chernyshev, A. L. (2018). Topography of spin liquids on a triangular lattice. Phys. Rev. Lett. 120:207203. doi: 10.1103/PhysRevLett.120.207203

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Maksimov, P. A., White, S. R., and Cheryshev, A. L. (2017). Disorder-induced mimicry of a spin liquid in YbMgGaO4. Phys. Rev. Lett. 119:157201. doi: 10.1103/PhysRevLett.119.157201

CrossRef Full Text | Google Scholar

Keywords: flux growth, single crystal, frustration, quantum materials, magnetism, triangular lattice

Citation: Ferreira T, Xing J, Sanjeewa LD and Sefat AS (2020) Frustrated Magnetism in Triangular Lattice TlYbS2 Crystals Grown via Molten Flux. Front. Chem. 8:127. doi: 10.3389/fchem.2020.00127

Received: 15 November 2019; Accepted: 12 February 2020;
Published: 26 February 2020.

Edited by:

Katsuya Teshima, Shinshu University, Japan

Reviewed by:

Akitaka Ito, Kochi University of Technology, Japan
Hiroaki Iguchi, Tohoku University, Japan

Copyright © 2020 Ferreira, Xing, Sanjeewa and Sefat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Timothy Ferreira, ferreirate@ornl.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.