Emergent honeycomb lattice in LiZn2Mo3O8

Chemistry of Quantum Spin Liquids

Quantum spin liquids are an exciting playground for exotic physical phenomena and emergent many-body quantum states. The realization and discovery of quantum spin liquid candidate materials and associated phenomena lie at the intersection of solid-state chemistry, condensed matter physics, and materials science and engineering. In this review, we provide the current status of the crystal chemistry, synthetic techniques, physical properties, and research methods in the field of quantum spin liquids. We highlight a number of specific quantum spin liquid candidate materials and their structure-property relationships, elucidating their fascinating behavior and connecting it to the intricacies of their structures. Furthermore, we share our thoughts on defects and their inevitable presence in materials, of which quantum spin liquids are no exception, which can complicate the interpretation of characterization of these materials, and urge the community to extend their attention to materials preparation and data analysis, cognizant of the impact of defects. This review was written with the intention of providing guidance on improving the materials design and growth of quantum spin liquids, and to paint a picture of the beauty of the underlying chemistry of this exciting class of materials.

Exotic topological point and line nodes in the plaquette excitations of a frustrated Heisenberg antiferromagnet on the honeycomb lattice

A number of topological nodes including Dirac, quadratic and three-band touching points as well as a pair of degenerate Dirac line nodes are found to emerge in the triplet plaquette excitations of the frustrated spin-1/2J1-J2antiferromagnetic Heisenberg honeycomb model when the ground state of the system lies in a spin-disordered plaquette-valence-bond-solid phase. A six-spin plaquette operator theory of this honeycomb model has been developed for this purpose by using the eigenstates of an isolated Heisenberg hexagonal plaquette. Spin-1/2 operators are thus expressed in the Fock space spanned by the plaquette operators those are obtained in terms of exact analytic form of eigenstates for a single frustrated Heisenberg hexagon. Ultimately, an effective interacting boson model of this system is obtained on the basis of low energy singlets and triplets plaquette operators by employing a mean-field approximation. The values of ground state energy and spin gap of this system have been estimated and the validity of this formalism has been tested upon comparison with the known results. Emergence of topological point and line nodes on the basis of spin-disordered ground state noted in this investigation is very rare on any frustrated system as well as the presence of triplet flat band. Evolution of those topological nodes is studied throughout the full frustrated regime. Finally, emergence of topological phases has been reported upon adding a time-reversal-symmetry breaking term to the Hamiltonian. Coexistence of spin gap with either topological nodes or phases turns this honeycomb model an interesting one.

Randomly Hopping Majorana Fermions in the Diluted Kitaev System α-Ru_{0.8}Ir_{0.2}Cl_{3}

dc and ac magnetic susceptibility, magnetization, specific heat, and Raman scattering measurements are combined to probe low-lying spin excitations in α-Ru_{1-x}Ir_{x}Cl_{3} (x≈0.2), which realizes a disordered spin liquid. At intermediate energies (ℏω>3  meV), Raman spectroscopy evidences linearly ω-dependent Majorana-like excitations, obeying Fermi statistics. This points to robustness of a Kitaev paramagnetic state under spin vacancies. At low energies below 3 meV, we observe power-law dependences and quantum-critical-like scalings of the thermodynamic quantities, implying the presence of a weakly divergent low-energy density of states. This scaling phenomenology is interpreted in terms of the random hoppings of Majorana fermions. Our results demonstrate an emergent hierarchy of spin excitations in a diluted Kitaev honeycomb system subject to spin vacancies and bond randomness.

Nontrivial topological phases on the stuffed honeycomb lattice

We report the appearance of nontrivial topological phases in a tight-binding model on the stuffed honeycomb lattice. The model contains nearest neighbor and next nearest neighbor hopping terms coupled with an additional phase depending on the direction of hopping. Chern insulating and semi-metallic phases emerge with the change of hopping parameters. Nonzero Chern numbers characterizing the bands and the existence of topologically protected edge states in the gap between the relevant bands confirm the presence of those phases. We show that adding an extra basis to Haldane's honeycomb model can lead to an additional topological phase characterized by Chern number [Formula: see text]. Transition between different topological phases driven by the hopping parameters has been illustrated in the topological phase diagram of the system. Zero temperature Hall conductivity along with density of states is evaluated. Topological properties of another tight-binding model on the stuffed square lattice are also reported in this article.

Quantum magnetisms in uniform triangular lattices Li2AMo3O8 (A = In, Sc)

Molecular based spin-1/2 triangular lattice systems such as LiZn2Mo3O8 have attracted research interest. Distortions, defects, and intersite disorder are suppressed in such molecular-based magnets, and intrinsic geometrical frustration gives rise to unconventional and unexpected ground states. Li2AMo3O8 (A = In or Sc) is such a compound where spin-1/2 Mo3O13 clusters in place of Mo ions form the uniform triangular lattice. Their ground states are different according to the A site. Li2InMo3O8 undergoes conventional 120° long-range magnetic order below TN = 12 K whereas isomorphic Li2ScMo3O8 exhibits no long-range magnetic order down to 0.5 K. Here, we report exotic magnetisms in Li2InMo3O8 and Li2ScMo3O8 investigated by muon spin rotation (μSR) and inelastic neutron scattering (INS) spectroscopies using polycrystalline samples. Li2InMo3O8 and Li2ScMo3O8 show completely different behaviors observed in both μSR and INS measurements, representing their different ground states. Li2InMo3O8 exhibits spin wave excitation which is quantitatively described by the nearest neighbor anisotropic Heisenberg model based on the 120° spin structure. In contrast, Li2ScMo3O8 undergoes short-range magnetic order below 4 K with quantum-spin-liquid-like magnetic fluctuations down to the base temperature. Origin of the different ground states is discussed in terms of anisotropies of crystal structures and magnetic interactions.

Correlated Partial Disorder in a Weakly Frustrated Quantum Antiferromagnet

Partial disorder-the microscopic coexistence of long-range magnetic order and disorder-is a rare phenomenon that has been experimentally and theoretically reported in some Ising- or easy plane-spin systems, driven by entropic effects at finite temperatures. Here, we present an analytical and numerical analysis of the S=1/2 Heisenberg antiferromagnet on the sqrt[3]×sqrt[3]-distorted triangular lattice, which shows that its quantum ground state has partial disorder in the weakly frustrated regime. This state has a 180° Néel ordered honeycomb subsystem coexisting with disordered spins at the hexagon center sites. These central spins are ferromagnetically aligned at short distances, as a consequence of a Casimir-like effect originated by the zero-point quantum fluctuations of the honeycomb lattice.

Scaling and data collapse from local moments in frustrated disordered quantum spin systems

Recently measurements on various spin–1/2 quantum magnets such as H₃LiIr₂O₆, LiZn₂Mo₃O₈, ZnCu₃(OH)₆Cl₂ and 1T-TaS₂—all described by magnetic frustration and quenched disorder but with no other common relation—nevertheless showed apparently universal scaling features at low temperature. In particular the heat capacity C[H, T] in temperature T and magnetic field H exhibits T/H data collapse reminiscent of scaling near a critical point. Here we propose a theory for this scaling collapse based on an emergent random-singlet regime extended to include spin-orbit coupling and antisymmetric Dzyaloshinskii-Moriya (DM) interactions. We derive the scaling C[H, T]/T ~ H^−γF_q[T/H] with F_q[x] = x^q at small x, with q ∈ {0, 1, 2} an integer exponent whose value depends on spatial symmetries. The agreement with experiments indicates that a fraction of spins form random valence bonds and that these are surrounded by a quantum paramagnetic phase. We also discuss distinct scaling for magnetization with a q-dependent subdominant term enforced by Maxwell’s relations.

There are many proposals for new forms of quantum matter in frustrated magnets but in practice disorder prevents the realisation of theoretically-tractable idealised models. Kimchi et al. show that recently observed scaling behavior common to several disordered quantum magnets can be understood as the emergence of a universal random-singlet regime.

Physical realization of 2D spin liquid state by ab initio design and strain engineering in FeX3

So far, no real physical two dimensional (2D) monolayer materials with spin liquid (SL) state have been identified, although the SL state has been analytically predicted to be present in 2D triangular, honeycomb and kagome model lattices. Identifying a realistic monolayer, 2D SL material thus enables one to more clearly observe and understand the physics of the fractional quantum Hall effect, high-temperature superconductivity and magnetic monopole. In this work, we have performed first principles calculations within density functional theory to investigate the magnetic phase diagram in monolayer FeCl₃, which reveals that the SL state may exist near the quantum phase transition between different antiferromagnetic (AFM) phases. Fundamentally, under a biaxial in-plane strain, the quantum phase transition appear and the ratio between the second neighboring exchange interaction (J ₂) and the first neighboring interaction (J ₁) falls into the SL range (0.21  ~  0.28). Through comparison, the exchange energy ratio of monolayer FeBr₃ (0.327  ~  0.565) is beyond the SL range and thus FeBr₃ could not exhibit SL characters. During the quantum phase transition induced by strain, the magnetic ground state transforms from AFM-Néel phase to AFM-stripy. Near the critical point, the semiconductor FeCl₃ transforms from indirect to direct band gap. Our findings provide insights that the monolayer FeCl₃ is a realistic 2D SL prototype material and the SL state could be observed by strain engineering.

Tunable Quantum Spin Liquidity in the 1/6th-Filled Breathing Kagome Lattice

We present measurements on a series of materials, Li_{2}In_{1-x}Sc_{x}Mo_{3}O_{8}, that can be described as a 1/6th-filled breathing kagome lattice. Substituting Sc for In generates chemical pressure which alters the breathing parameter nonmonotonically. Muon spin rotation experiments show that this chemical pressure tunes the system from antiferromagnetic long range order to a quantum spin liquid phase. A strong correlation with the breathing parameter implies that it is the dominant parameter controlling the level of magnetic frustration, with increased kagome symmetry generating the quantum spin liquid phase. Magnetic susceptibility measurements suggest that this is related to distinct types of charge order induced by changes in lattice symmetry, in line with the theory of Chen et al. [Phys. Rev. B 93, 245134 (2016)PRBMDO2469-995010.1103/PhysRevB.93.245134]. The specific heat for samples at intermediate Sc concentration, which have the minimum breathing parameter, show consistency with the predicted U(1) quantum spin liquid.

Emergent Chiral Spin State in the Mott Phase of a Bosonic Kane-Mele-Hubbard Model

Recently, the frustrated XY model for spins 1/2 on the honeycomb lattice has attracted a lot of attention in relation with the possibility to realize a chiral spin liquid state. This model is relevant to the physics of some quantum magnets. Using the flexibility of ultracold atom setups, we propose an alternative way to realize this model through the Mott regime of the bosonic Kane-Mele-Hubbard model. The phase diagram of this model is derived using bosonic dynamical mean-field theory. Focusing on the Mott phase, we investigate its magnetic properties as a function of frustration. We do find an emergent chiral spin state in the intermediate frustration regime. Using exact diagonalization we study more closely the physics of the effective frustrated XY model and the properties of the chiral spin state. This gapped phase displays a chiral order, breaking time-reversal and parity symmetry, but is not topologically ordered (the Chern number is zero).

Orphan Spins in the S=5/2 Antiferromagnet CaFe_{2}O_{4}

CaFe_{2}O_{4} is an anisotropic S=5/2 antiferromagnet with two competing A (↑↑↓↓) and B (↑↓↑↓) magnetic order parameters separated by static antiphase boundaries at low temperatures. Neutron diffraction and bulk susceptibility measurements, show that the spins near these boundaries are weakly correlated and a carry an uncompensated ferromagnetic moment that can be tuned with a magnetic field. Spectroscopic measurements find these spins are bound with excitation energies less than the bulk magnetic spin waves and resemble the spectra from isolated spin clusters. Localized bound orphaned spins separate the two competing magnetic order parameters in CaFe_{2}O_{4}.

Rearrangement of van der Waals stacking and formation of a singlet state at T = 90 K in a cluster magnet

A change of van der Waals stacking occurs spontaneously at 90 K in a cluster magnet.

Direct assignment of molecular vibrations via normal mode analysis of the neutron dynamic pair distribution function technique

For over a century, vibrational spectroscopy has enhanced the study of materials. Yet, assignment of particular molecular motions to vibrational excitations has relied on indirect methods. Here, we demonstrate that applying group theoretical methods to the dynamic pair distribution function analysis of neutron scattering data provides direct access to the individual atomic displacements responsible for these excitations. Applied to the molecule-based frustrated magnet with a potential magnetic valence-bond state, LiZn2Mo3O8, this approach allows direct assignment of the constrained rotational mode of Mo3O13 clusters and internal modes of MoO6 polyhedra. We anticipate that coupling this well known data analysis technique with dynamic pair distribution function analysis will have broad application in connecting structural dynamics to physical properties in a wide range of molecular and solid state systems.

Molecular quantum magnetism in LiZn2Mo3O8

Inelastic neutron scattering at low temperatures T≤30  K from a powder of LiZn2Mo3O8 demonstrates this triangular-lattice antiferromagnet hosts collective magnetic excitations from spin-1/2 Mo3O13 molecules. Apparently gapless (Δ<0.2  meV) and extending at least up to 2.5 meV, the low-energy magnetic scattering cross section is surprisingly broad in momentum space and involves one-third of the spins present above 100 K. The data are compatible with the presence of valence bonds involving nearest-neighbor and next-nearest-neighbor spins forming a disordered or dynamic state.