Search

Identifying experimental probes capable of diagnosing extreme quantum behavior is widely regarded as one of the foremost challenges in modern condensed matter physics. Here, we propose a novel approach for detecting chiral Kitaev spin liquid states through measurements of the local dynamical spin structure factor on the boundary using scanning tunneling microscopy (STM). We specifically focus on unpaired ("dangling") Majorana fermions, which naturally emerge along boundaries of Kitaev spin liquids, and can serve as indicators of chiral boundary modes under broad conditions, thereby offering a clear signature of these exotic quantum states.

Recent experimental work indicates that the Kitaev spin liquid may be realizable close to zero magnetic field in exfoliated $\alpha$-RuCl$_3$ flakes, thus providing a more versatile setting for studying non-Abelian anyons. Here, we propose a robust nanoscale signature of the Kitaev spin liquid that results from its edge states and manifests in the low-energy spin dynamics. In particular, we highlight a singular peak in the dynamical spin structure factor of a zigzag edge whose energy scales linearly with only one component of the magnetic field. This sharp feature in the local spin dynamics directly reflects the bond-directional Kitaev spin interactions and, more generally, the projective symmetries of the Kitaev spin liquid. We demonstrate that our proposed signature survives in the presence of edge disorder as well as non-Kitaev interactions, and provide detailed guidelines for experimentally observing it using inelastic electron tunneling spectroscopy and color-center relaxometry., Comment: 6+2 pages, 4 figures

The excitonic insulator, a quantum mechanical state arising from exciton condensation, was proposed theoretically many years ago but has yet to be experimentally confirmed. The discovery of correlated transition metal oxides based on $4d$ and $5d$ elements, where the on-site Coulomb repulsion is comparable to the dominant hopping amplitude, presents a unique opportunity to study exciton condensation. By constructing an effective mean field Raman operator for the Hubbard model, we derive the low-energy electronic Raman scattering cross section, demonstrating Raman spectroscopy as a powerful tool for detecting exciton condensation. Here, we demonstrate that Raman scattering directly reveals exciton condensation in the bilayer iridate Sr$_3$Ir$_2$O$_7$ under pressure., Comment: The main text is 6 pages long, including 3 figures, and the supplemental material is 8 pages long, including 1 figure

Multi-Q magnetic structures on triangular lattices, with their two-dimensional topological spin texture, have attracted significant interest. However, unambiguously confirming their formation by excluding the presence of three equally-populated single-Q domains remains challenging. In the metallic triangular lattice antiferromagnet Co1/3TaS2, two magnetic ground states have been suggested at different temperature ranges, with the low-temperature phase being a triple-Q structure corresponding to the highest-density Skyrmion lattice. Using inelastic neutron scattering (INS) and advanced spin dynamics simulations, we demonstrate a clear distinction in the excitation spectra between the single-Q and triple-Q phases of Co1/3TaS2 and, more generally, a triangular lattice. First, we refined the spin Hamiltonian by fitting the excitation spectra measured in its paramagnetic phase, allowing us to develop an unbiased model independent of magnetic ordering. Second, we observed that the two magnetically ordered phases in Co1/3TaS2 exhibit markedly different behaviors in their long-wavelength Goldstone modes. Our spin model, derived from the paramagnetic phase, confirms that these behaviors originate from the single-Q and triple-Q nature of the respective ordered phases, providing unequivocal evidence of the single-Q to triple-Q phase transition in Co1/3TaS2. Importantly, we propose that the observed contrast in the long-wavelength spin dynamics between the single-Q and triple-Q orderings is universal, offering a potentially unique way to distinguish a generic triple-Q ordering on a triangular lattice from its multi-domain single-Q counterparts. We describe its applicability with examples of similar hexagonal systems forming potential triple-Q orderings. (For the full abstract, please refer to the manuscript), Comment: 14 pages, 6 figures excluding Appendices

We develop a systematic large-$N$ expansion based on the Schwinger boson representation of SU(4) coherent states of dimers for the paradigmatic spin-$1/2$ bilayer square lattice Heisenberg antiferromagnet. This system exhibits a quantum phase transition between a quantum paramagnetic state and a N\'eel order state, driven by the coupling constant $g = J'/J$, which is defined as the ratio between the inter-dimer $J'$ and intra-dimer $J$ exchange interactions. We demonstrate that this approach accurately describes static and dynamic properties on both sides of the quantum phase transition. The critical coupling constant $g_c \approx 0.42$ and the dynamic spin structure factor reproduce quantum Monte Carlo results with high precision. Notably, the $1/N$ corrections reveal the longitudinal mode of the magnetically ordered phase along with the overdamping caused by its decay into the two-magnon continuum. The present large-$N$ $SU(N)$ Schwinger boson theory can be extended to more general cases of quantum paramagnets that undergo a quantum phase transition into magnetically ordered states.

By utilizing a multi-orbital periodic Anderson model with parameters obtained from \textit{ab initio} band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both non-magnetic $4f^0$ and magnetic $4f^2$ virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn$_3$. The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of $4f$ materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution., Comment: 30 pages, 16 figures

Anderson's groundbreaking ideas of resonant valence bond (RVB) liquid and spin-charge separation initiated a transformative shift in modern physics. Extensive implications influenced a broad spectrum of fields, from high-temperature superconductors to quantum computing, and gave birth to key concepts in physics, such as quantum spin liquids, emergent gauge symmetries, topological order, and fractionalisation. Despite extensive efforts to demonstrate the existence of an RVB phase in the Hubbard model, a definitive realisation has proven elusive. Here, we present a concise, realistic, and elegant solution to this longstanding problem by demonstrating analytically that an RVB spin liquid, exhibiting spin-charge separation, emerges as the ground state of doped Mott insulators on corner-sharing tetrahedral lattices with frustrated hopping near half-filling -- a manifestation of the counter-Nagaoka effect. We confirm numerically that this result holds for finite-size systems, finite dopant density, and small exchange interactions. While much attention has been devoted to the emergence of new states from geometrically frustrated interactions, our work demonstrates that kinetic energy frustration in doped Mott insulators may be pivotal to stabilise robust, topologically ordered states in real materials.

We report neutron scattering, pressure-dependent AC calorimetry, and AC magnetic susceptibility measurements of triangular lattice NaYbSe$_2$. We observe a continuum of scattering, which is reproduced by matrix product simulations, and no phase transition is detected in any bulk measurements. Comparison to heat capacity simulations suggest the material is within the Heisenberg spin liquid phase. AC Susceptibility shows a significant 23~mK downturn, indicating a gap in the magnetic spectrum. The combination of a gap with no detectable magnetic order, comparison to theoretical models, and comparison to other $A$YbSe$_2$ compounds all strongly indicate NaYbSe$_2$ is within the quantum spin liquid phase. The gap also allows us to rule out a gapless Dirac spin liquid, with a gapped $\mathbb{Z}_2$ liquid the most natural explanation., Comment: 5 pages, 4 figures; 7 pages and 13 figures supplemental materials

We employ a classical limit grounded in SU(4) coherent states to investigate the temperature-dependent dynamical spin structure factor of the $S = 1/2$ ladder consisting of weakly coupled dimers. By comparing the outcomes of this classical approximation with density matrix renormalization group and exact diagonalization calculations in finite size ladders, we demonstrate that the classical dynamics offers an accurate approximation across the entire temperature range when the interdimer coupling is weak and a good approximation in the high temperature regime even when the interdimer coupling is strong. This agreement is achieved after appropriately rescaling the temperature axis and renormalizing expectation values to satisfy a quantum sum rule, following D. Dahlbom et al. [Phys. Rev. B 109, 014427 (2024)]. We anticipate the method will be particularly effective when applied to 2D and 3D lattices composed of weakly-coupled dimers, situations that remain challenging for alternative numerical methods., Comment: 15 pages, 8 figures

We investigate the ground state magnetic configurations of a Fibonacci chain of classical spins with nearest-neighbor ferromagnetic and monoaxial Dzyaloshinskii-Moriya exchange interactions. Our analysis reveals a diverse array of magnetic textures induced by an external magnetic field perpendicular to the Dzyaloshinskii-Moriya vector. These textures exhibit a spectrum ranging from a quasi-fully polarized non-collinear state under high magnetic fields, capable of maintaining metastable chiral soliton topological defects, to a variety of quasicrystalline chiral soliton lattices below a critical field $H_c$. For a range of magnetic fields below $H_c$, the ground state spin textures result from the interplay between an effective quasiperiodic potential influencing the solitons and their repulsive interactions. At lower magnetic fields, the system experiences a commensurate-incommensurate transition, signified by the appearance of discommensurations in the quasicrystalline soliton lattice. In the absence of an external magnetic field, the ground state assumes a helical configuration with a quasiperiodic pitch angle., Comment: 14 pages, 12 figures

The S = 1/2 triangular lattice antiferromagnet (TLAF) is a paradigmatic example of frustrated quantum magnetism. An ongoing challenge involves understanding the influence of exchange anisotropy on the collective behavior within such systems. Using inelastic neutron scattering (INS) and advanced calculation techniques, we have studied the low and high-temperature spin dynamics of Ba2La2CoTe2O12 (BLCTO): a Co2+-based Jeff = 1/2 TLAF that exhibits 120{\deg} order below TN = 3.26 K. We determined the spin Hamiltonian by fitting the energy-resolved paramagnetic excitations measured at T > TN, revealing exceptionally strong easy-plane XXZ anisotropy. Below TN, the excitation spectrum exhibits a high energy continuum having a larger spectral weight than the single-magnon modes, suggesting a scenario characterized by a spinon confinement length that markedly exceeds the lattice spacing. We conjecture that this phenomenon arises from the proximity to a quantum melting point, even under strong easy-plane XXZ anisotropy. Finally, we highlight characteristic flat features in the excitation spectrum, which are connected to higher-order van Hove singularities in the magnon dispersion directly induced by easy-plane XXZ anisotropy. Our results provide a rare experimental insight into the nature of highly anisotropic S = 1/2 TLAFs between the Heisenberg and XY limits., Comment: 5 figures, Supplementary Information can be found on the journal's website

Calculating dynamical spin correlations is essential for matching model magnetic exchange Hamiltonians to momentum-resolved spectroscopic measurements. A major numerical bottleneck is the diagonalization of the dynamical matrix, especially in systems with large magnetic unit cells, such as those with incommensurate magnetic structures or quenched disorder. In this paper, we demonstrate an efficient scheme based on the kernel polynomial method for calculating dynamical correlations of relevance to inelastic neutron scattering experiments. This method reduces the scaling of numerical cost from cubic to linear in the magnetic unit cell size., Comment: Submission to SciPost. 27 pages, 6 figures. Typos corrected, added additional minor comments to derivations

The asymmetric distribution of chiral objects with opposite chirality is of great fundamental interests ranging from molecular biology to particle physics. In quantum materials, chiral states can build on inversion-symmetry-breaking lattice structures or emerge from spontaneous magnetic ordering induced by competing interactions. Although the handedness of a chiral state can be changed through external fields, a spontaneous chirality flipping has yet to be discovered. In this letter, we present experimental evidence of chirality flipping via changing temperature in a topological magnet EuAl$_4$, which features orthogonal spin and charge density waves (SDW/CDW). Using circular dichroism of Bragg peaks in the resonant magnetic x-ray scattering, we find that the chirality of the helical SDW flips through a first order phase transition with modified SDW wavelength. Intriguingly, we observe that the CDW couples strongly with the SDW and displays a rare commensurate-to-incommensurate transition at the chirality flipping temperature. Combining with first principles calculations and angle resolved photoemission spectroscopy, we establish the Fermi surface origin of the helical SDW with intertwined spin, charge, and lattice degrees of freedom in EuAl$_4$. Our results reveal an unprecedented spontaneous chirality flipping and lays the groundwork for a new functional manipulation of chirality through momentum dependent spin-charge-lattice interactions., Comment: Supplementary materials are available from the corresponding author upon request

Simulating quantum spin systems at finite temperatures is an open challenge in many-body physics. This work studies the temperature-dependent spin dynamics of a pivotal compound, FeI$_2$, to determine if universal quantum effects can be accounted for by a phenomenological renormalization of the dynamical spin structure factor $S(\mathbf{q}, \omega)$ measured by inelastic neutron scattering. Renormalization schemes based on the quantum-to-classical correspondence principle are commonly applied at low temperatures to the harmonic oscillators describing normal modes. However, it is not clear how to extend this renormalization to arbitrarily high temperatures. Here we introduce a temperature-dependent normalization of the classical moments, whose magnitude is determined by imposing the quantum sum rule, i.e. $\int d\omega d\mathbf{q} S(\mathbf{q}, \omega) = N_S S (S+1)$ for $N_S$ dipolar magnetic moments. We show that this simple renormalization scheme significantly improves the agreement between the calculated and measured $S(\mathbf{q}, \omega)$ for FeI$_{2}$ at all temperatures. Due to the coupled dynamics of dipolar and quadrupolar moments in that material, this renormalization procedure is extended to classical theories based on SU(3) coherent states, and by extension, to any SU(N) coherent state representation of local multipolar moments., Comment: Associated source code for reproducing calculations available at: https://github.com/SunnySuite/SunnyContributed

The triangular lattice antiferromagnet (TLAF) has been the standard paradigm of frustrated magnetism for several decades. The most common magnetic ordering in insulating TLAFs is the 120 structure. However, a new triple-Q chiral ordering can emerge in metallic TLAFs, representing the short wavelength limit of magnetic skyrmion crystals. We report the metallic TLAF Co1/3TaS2 as the first example of tetrahedral triple-Q magnetic ordering with the associated topological Hall effect (non-zero {\sigma}_{xy}(H=0)). We also present a theoretical framework that describes the emergence of this magnetic ground state, which is further supported by the electronic structure measured by angle-resolved photoemission spectroscopy. Additionally, our measurements of the inelastic neutron scattering cross section are consistent with the calculated dynamical structure factor of the tetrahedral triple-Q state., Comment: 23 pages, 4 figures, SI not included. Accepted for publication in Nature Communications

Magnetic materials with noncoplanar magnetic structures can show unusual physical properties driven by nontrivial topology. Topologically-active states are often multi-q structures, which are challenging to stabilize in models and to identify in materials. Here, we use inelastic neutron-scattering experiments to show that the insulating double perovskites Ba2YRuO6 and Ba2LuRuO6 host a noncoplanar 3-q structure on the face-centered cubic lattice. Quantitative analysis of our neutron-scattering data reveals that these 3-q states are stabilized by biquadratic interactions. Our study identifies double perovskites as a highly promising class of materials to realize topological magnetism, elucidates the stabilization mechanism of the 3-q state in these materials, and establishes neutron spectroscopy on powder samples as a valuable technique to distinguish multi-q from single-q states, facilitating the discovery of topologically-nontrivial magnetic materials., Comment: 5 pages, 4 figures

We consider the problem of extracting an effective low-energy spin model from a Kondo Lattice Model (KLM) with classical localized moments. The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange $J$ excludes the possibility of using perturbation theory beyond the second order Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction at zero temperature. Here we introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions by integrating out the conduction electrons of the original KLM. The resulting effective spin model reproduces the phase diagram obtained with the original KLM as a function of magnetic field and easy-axis anisotropy and reveals the effective four-spin interactions that are responsible for the field induced skyrmion crystal phase. Moreover, this minimal spin model enables an efficient computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with the effective model and the original KLM reveals a good agreement for the magnon dispersion despite the fact that this information was not included in the training data set.

We consider a class of Hubbard-Stratonovich transformations suitable for treating Hubbard interactions in the context of quantum Monte Carlo simulations. A tunable parameter $p$ allows us to continuously vary from a discrete Ising auxiliary field ($p=\infty$) to a compact auxiliary field that couples to electrons sinusoidally ($p=0$). In tests on the single-band Hubbard model, we find that the severity of the sign problem decreases systematically with increasing $p$. Selecting $p$ finite, however, enables continuous sampling methods like the Langevin or Hamiltonian Monte Carlo methods. We explore the tradeoffs between various simulation methods through numerical benchmarks.

We present a high-field study of the strongly anisotropic easy-plane square lattice $S$ = 2 quantum magnet Ba$_{2}$FeSi$_{2}$O$_{7}$. This compound is a rare high-spin antiferromagnetic system with very strong easy-plane anisotropy, such that the interplay between spin level crossings and antiferromagnetic order can be studied. We observe a magnetic field-induced spin level crossing occurring within an ordered state. This spin level crossing appears to preserve the magnetic symmetry while producing a non-monotonic dependence the order parameter magnitude. The resulting temperature-magnetic field phase diagram exhibits two dome-shaped regions of magnetic order overlapping around 30 T. The ground state of the lower-field dome is predominantly a linear combination of $| S^{z} = 0 \rangle$ and $ |S^{z} = 1 \rangle$ states, while the ground state of the higher-field dome can be approximated by a linear combination of $| S^{z} = 1 \rangle $ and $ | S^{z} = 2\rangle$ states. At 30 T, where the spin levels cross, the magnetization exhibits a slanted plateau, {\color {black}the magnetocaloric effect shows a broad hump, and the electric polarization shows a weak slope change}. We determined the detailed magnetic phase boundaries and the spin level crossings using measurements of magnetization, electric polarization, and the magnetocaloric effect in pulsed magnetic fields to 60 T. We calculate these properties using a mean field theory based on direct products of SU(5) coherent states and find good agreement. Finally, we measure and calculate the magnetically-induced electric polarization that reflects magnetic ordering and spin level crossings. This multiferroic behavior provides another avenue for detecting phase boundaries and symmetry changes., Comment: 9 pages, 5 figures

We study the BCS-BEC crossover phenomenon of triplet exciton condensation using a half-filled bilayer Hubbard model. We calculate the dynamical spin structure factor and the exciton wave function on the phase boundary between the antiferromagnetically ordered and disordered phases. In the BCS regime, the formation and condensation of particle-hole bound states, namely excitons, occur simultaneously at the phase transition point, and the exciton wave function is extended in real space. In the BEC regime, on the other hand, bound states are well defined over the whole Brillouin zone also in the disordered phase, and the size of particle-hole pairs is smaller than a lattice space. Quantum criticality of charge-spin-orbital entangled states can emerge from intermediate coupling materials in the crossover regime., Comment: 6 pages, 3 figures, accepted in JPS Conference Proceedings (LT29)

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude $S$ have $N=2S+1$ levels. Alternatively, the spin state is fully characterized by a set of $N^{2}-1$ local physical observables, which we interpret as generalized spin components. For example, a spin with $S=1$ has three dipole components and five quadrupole components. These components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU($N$) coherent states in the Schr\"odinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a non-equilibrium relaxation process that creates CP$^{2}$ Skyrmions, which are topological defects with both dipole and quadrupole character.

Electrons at the border of localization generate exotic states of matter across all classes of strongly correlated electron materials and many other quantum materials with emergent functionality. Heavy electron metals are a model example, in which magnetic interactions arise from the opposing limits of localized and itinerant electrons. This remarkable duality is intimately related to the emergence of a plethora of novel quantum matter states such as unconventional superconductivity, electronic-nematic states, hidden order and most recently topological states of matter such as topological Kondo insulators and Kondo semimetals and putative chiral superconductors. The outstanding challenge is that the archetypal Kondo lattice model that captures the underlying electronic dichotomy is notoriously difficult to solve for real materials. Here we show, using the prototypical strongly-correlated antiferromagnet CeIn$_3$, that a multi-orbital periodic Anderson model embedded with input from ab initio bandstructure calculations can be reduced to a simple Kondo-Heisenberg model, which captures the magnetic interactions quantitatively. We validate this tractable Hamiltonian via high-resolution neutron spectroscopy that reproduces accurately the magnetic soft modes in CeIn$_3$, which are believed to mediate unconventional superconductivity. Our study paves the way for a quantitative understanding of metallic quantum states such as unconventional superconductivity.

Quantum spin liquids (QSL) are theoretical states of matter with long-range entanglement and exotic quasiparticles. However, they generally elude quantitative theory, rendering their underlying phases mysterious and hampering efforts to identify experimental QSL states. Here we study triangular lattice resonating valence bond QSL candidate materials KYbSe$_2$ and NaYbSe$_2$. We measure the magnon modes in their 1/3 plateau phase, where quantitative theory is tractable, using inelastic neutron scattering and fit them using nonlinear spin wave theory. We also fit the KYbSe$_2$ heat capacity using high temperature series expansion. Both KYbSe$_2$ fits yield the same magnetic Hamiltonian to within uncertainty, confirming previous estimates and showing the Heisenberg $J_2/J_1$ to be an accurate model for these materials. Most importantly, comparing KYbSe$_2$ and NaYbSe$_2$ shows that smaller $A$-site Na$^+$ ion has a larger $J_2/J_1$ ratio. However, hydrostatic pressure applied to KYbSe$_2$ increases the ordering temperature (a result consistent with density functional theory calculations), indicating that pressure decreases $J_2/J_1$. These results show how periodic table and hydrostatic pressure can tune the $A$YbSe$_2$ materials in a controlled way., Comment: 7 pages, 7 figures; 4 pages and 7 additional figures of supplemental information

Quantum magnets admit more than one classical limit and $N$-level systems with strong single-ion anisotropy are expected to be described by a classical approximation based on SU($N$) coherent states. Here we test this hypothesis by modeling finite temperature inelastic neutron scattering (INS) data of the effective spin-one antiferromagnet \bfso{}. The measured dynamic structure factor is calculated with a generalized Landau-Lifshitz dynamics for SU($3$) spins. Unlike the traditional classical limit based on SU($2$) coherent states, the results obtained with classical SU($3$) spins are in good agreement with the measured temperature-dependent spectrum. The SU($3$) approach developed here provides a general framework to understand the broad class of materials comprising weakly coupled antiferromagnetic dimers, trimers, or tetramers, and magnets with strong single-ion anisotropy., Comment: 8 pages 4 figures

The Landau-Lifshitz equation describes the time-evolution of magnetic dipoles, and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here, we reformulate the dynamics of these $N^2-1$ generalized spin components as a mean-field Schr\"odinger equation on the $N$-dimensional coherent state. This viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics.

Magnetic skyrmions are nanoscale topological textures that have been recently observed in different families of quantum magnets. These textures are known as CP$^1$ skyrmions because the target manifold of the magnetization field is the 2D sphere $S^2 \cong$ CP$^{1}$. Here we report the emergence of magnetic CP$^2$ skyrmions in realistic spin-$1$ models. Unlike CP$^1$ skyrmions, CP$^2$ skyrmions can also arise as metastable textures of quantum paramagnets, opening a new road to discover emergent topological solitons in non-magnetic materials. The quantum phase diagram of the spin-$1$ models also includes magnetic field-induced CP$^2$ skyrmion crystals that can be detected with regular momentum- (diffraction) and real-space (Lorentz transmission electron microscopy) experimental techniques., Comment: main text: 8 pages, 4 figures. Supplemental Info: 6 pages, 3 figures

Recent inelastic neutron scattering (INS) experiments of the triangular antiferromagnet Ba$_3$CoSb$_2$O$_9$ revealed strong deviations from semiclassical theories. We demonstrate that key features of the INS data are well reproduced by a parton Schwinger boson theory beyond the saddle-point approximation. The measured magnon dispersion is well reproduced by the dispersion of two-spinon bound states (poles of the emergent gauge fields propagator), while the low-energy continuum scattering is reproduced by a quasifree two-spinon continuum, suggesting that a free-spinon gas is a good initial framework to study magnetically ordered states near a quantum melting point.

Excitonic insulators are usually considered to form via the condensation of a soft charge mode of bound electron-hole pairs. This, however, presumes that the soft exciton is of spin-singlet character. Early theoretical considerations have also predicted a very distinct scenario, in which the condensation of magnetic excitons results in an antiferromagnetic excitonic insulator state. Here we report resonant inelastic x-ray scattering (RIXS) measurements of Sr3Ir2O7. By isolating the longitudinal component of the spectra, we identify a magnetic mode that is well-defined at the magnetic and structural Brillouin zone centers, but which merges with the electronic continuum in between these high-symmetry points and which decays upon heating concurrent with a decrease in the material's resistivity. We show that a bilayer Hubbard model, in which electron-hole pairs are bound by exchange interactions, consistently explains all the electronic and magnetic properties of Sr3Ir2O7 indicating that this material is a realization of the long-predicted antiferromagnetic excitonic insulators phase., Comment: accepted in Nature Communications; 6 pages main materials

We present a systematic study of the formation of skyrmion crystals in a triangular Kondo Lattice model for generic electron filling fractions. Our results indicate that the four-sublattice chiral antiferromagnetic ordering that was reported more than one decade ago can be understood as the dense limit of a sequence of skyrmion crystals whose lattice parameter is dictated by the Fermi wave-vector. This observation has important implications for the ongoing search of skyrmion crystals in metallic materials with localized magnetic moments., Comment: 17 pages, 4 figures

The Heisenberg triangular lattice quantum spin liquid and the phase transitions to nearby magnetic orders have received much theoretical attention, but clear experimental manifestations of these states are rare. This work investigates a new spin-half Yb$^{3+}$ delafossite material, KYbSe$_2$, whose inelastic neutron scattering spectra reveal a diffuse continuum with a sharp lower bound. Applying entanglement witnesses to the data reveals significant multipartite entanglement spread between its neighbors, and analysis of its magnetic exchange couplings shows close proximity to the triangular lattice Heisenberg quantum spin liquid. Key features of the data are reproduced by Schwinger-boson theory and tensor network calculations with a significant second-neighbor coupling $J_2$. The strength of the dynamical structure factor at the $K$ point shows a scaling collapse in $\hbar\omega/k_\mathrm{B}T$ down to 0.3 K, indicating a second-order quantum phase transition. Comparing this to previous theoretical work suggests that the proximate phase at larger $J_2$ is a gapped $\mathbb{Z}_2$ spin liquid, resolving a long-debated issue. We thus show that KYbSe$_2$ is close to a spin liquid phase, which in turn sheds light on the theoretical phase diagram itself., Comment: 6 pages main text, 6 pages methods, 7 pages supplemental information

The Schwinger boson theory provides a natural path for treating quantum spin systems with large quantum fluctuations. In contrast to semi-classical treatments, this theory allows us to describe a continuous transition between magnetically ordered and spin liquid states, as well as the continuous evolution of the corresponding excitation spectrum. The square lattice Heisenberg antiferromagnet is one of the first models that was approached with the Schwinger boson theory. Here we revisit this problem to reveal several subtle points that were omitted in previous treatments and that are crucial to further develop this formalism. These points include the freedom for the choice of the saddle point (Hubbard-Stratonovich decoupling and choice of the condensate) and the $1/N$ expansion in the presence of a condensate. A key observation is that the spinon condensate leads to Feynman diagrams that include contributions of different order in $1/N$, which must be accounted to get a qualitatively correct excitation spectrum. We demonstrate that a proper treatment of these contributions leads to an exact cancellation of the single-spinon poles of the dynamical spin structure factor, as expected for a magnetically ordered state. The only surviving poles are the ones arising from the magnons (two-spinon bound states), which are the true collective modes of an ordered magnet., Comment: 25 pages, 12 figures

The search for new elementary particles is one of the most basic pursuits in physics, spanning from subatomic physics to quantum materials. Magnons are the ubiquitous elementary quasiparticle to describe the excitations of fully-ordered magnetic systems. But other possibilities exist, including fractional and multipolar excitations. Here, we demonstrate that strong quantum interactions exist between three flavors of elementary quasiparticles in the uniaxial spin-one magnet FeI2. Using neutron scattering in an applied magnetic field, we observe spontaneous decay between conventional and heavy magnons and the recombination of these quasiparticles into a super-heavy bound-state. Akin to other contemporary problems in quantum materials, the microscopic origin for new physics in FeI2 is the quasi-flat nature of excitation bands and the presence of Kitaev anisotropic magnetic exchange interactions., Comment: 31 pages, 3 main figures

We introduce a classical limit of the dynamics of quantum spin systems based on coherent states of SU($N$), where $N$ is the dimension of the local Hilbert space. This approach, that generalizes the well-known Landau-Lifshitz dynamics from SU(2) to SU($N$), provides a better approximation to the exact quantum dynamics for a large class of realistic spin Hamiltonians, including $S \geq 1$ systems with large single-ion anisotropy and weakly-coupled multi-spin units, such as dimers or trimers. We illustrate this idea by comparing the spin structure factors of a single-ion $S=1$ model that are obtained with the SU(2) and SU(3) classical spin dynamics against the exact solution., Comment: 12 pages, 2 figures

Recent numerical studies indicate that the antiferromagnetic Kitaev honeycomb lattice model undergoes a magnetic-field-induced quantum phase transition into a new spin-liquid phase. This intermediate-field phase has been previously characterized as a gapless spin liquid. By implementing a recently developed variational approach based on the exact fractionalized excitations of the zero-field model, we demonstrate that the field-induced spin liquid is gapped and belongs to Kitaev's 16-fold way. Specifically, the low-field non-Abelian liquid with Chern number $C=\pm 1$ transitions into an Abelian liquid with $C=\pm 4$. The critical field and the field-dependent behaviors of key physical quantities are in good quantitative agreement with published numerical results. Furthermore, we derive an effective field theory for the field-induced critical point which readily explains the ostensibly gapless nature of the intermediate-field spin liquid.

We compute the low-energy excitation spectrum and the dynamical spin structure factor of the Kitaev-Heisenberg-Gamma model through a variational approach based on the exact fractionalized excitations of the pure Kitaev honeycomb model. This novel approach reveals the physical reason for the asymmetric stability of the Kitaev spin liquid phases around the ferromagnetic and antiferromagnetic Kitaev limits. Moreover, we demonstrate that the fractionalized excitations form bound states in specific regions of each Kitaev spin liquid phase and that certain phase transitions induced by Heisenberg and Gamma interactions are driven by the condensation of such a bound state. Remarkably, this bound state appears as a sharp mode in the dynamical spin structure factor, while its condensation patterns at the appropriate phase transitions provide a simple explanation for the magnetically ordered phases surrounding each Kitaev spin liquid phase., Comment: 20 pages, 11 figures

We investigate the nature of the magnetic excitations of a bilayer single-orbital Hubbard model in the intermediate-coupling regime. This model exhibits a quantum phase transition (QPT) between a paramagnetic (PM) and an insulating antiferromagnetic (AFM) phase at a critical value of the coupling strength. By using the random phase approximation, we show that the QPT is continuous when the PM state is a band insulator and that the corresponding quantum critical point (QCP) arises from the condensation of preformed excitons. These low-energy excitons reemerge on the other side of the QCP as the transverse and longitudinal modes of the AFM state. In particular, the longitudinal mode remains sharp for the model parameters relevant to Sr$_{3}$Ir$_{2}$O$_{7}$ because of the strong easy-axis anisotropy of this material., Comment: 11 pages, 5 figures

Employing a combination of symmetry analysis, low-energy modeling, and ab initio simulations, we predict the presence of magnetic-field-induced Weyl points close to the Fermi level in CaKFe$_4$As$_4$. Depending on the relative strengths of the magnetic field and of the spin-orbit coupling, the Weyl fermions can carry a topological charge of $\pm1$ or $\pm2$, making CaKFe$_4$As$_4$ a rare realization of a double-Weyl semimetal. We further predict experimental manifestations of these Weyl points, both in bulk properties, such as the anomalous Hall effect, and in surface properties, such as the emergence of prominent Fermi arcs. Because CaKFe$_4$As$_4$ displays unconventional fully-gapped superconductivity below 30 K, our findings open a novel route to investigate the interplay between superconductivity and Weyl fermions.

Geometrically frustrated spin-chain compounds such as Ca3Co2O6 exhibit extremely slow relaxation under a changing magnetic field. Consequently, both low-temperature laboratory experiments and Monte Carlo simulations have shown peculiar out-of-equilibrium magnetization curves, which arise from trapping in metastable configurations. In this work we simulate this phenomenon in a superconducting quantum annealing processor, allowing us to probe the impact of quantum fluctuations on both equilibrium and dynamics of the system. Increasing the quantum fluctuations with a transverse field reduces the impact of metastable traps in out-of-equilibrium samples, and aids the development of three-sublattice ferrimagnetic (up-up-down) long-range order. At equilibrium we identify a finite-temperature shoulder in the 1/3-to-saturated phase transition, promoted by quantum fluctuations but with entropic origin. This work demonstrates the viability of dynamical as well as equilibrium studies of frustrated magnetism using large-scale programmable quantum systems, and is therefore an important step toward programmable simulation of dynamics in materials using quantum hardware.

An ongoing challenge in the study of quantum materials, is to reveal and explain collective quantum effects in spin systems where interactions between different modes types are important. Here we approach this problem through a combined experimental and theoretical study of interacting transverse and longitudinal modes in an easy-plane quantum magnet near a continuous quantum phase transition. Our inelastic neutron scattering measurements of Ba$_{2}$FeSi$_{2}O$_{7}$ reveal the emergence, decay, and renormalization of a longitudinal mode throughout the Brillouin zone. The decay of the longitudinal mode is particularly pronounced at the zone center. To account for the many-body effects of the interacting low-energy modes in anisotropic magnets, we generalize the standard spin-wave theory. The measured mode decay and renormalization is reproduced by including all one-loop corrections. The theoretical framework developed here is broadly applicable to quantum magnets with more than one type of low energy mode., Comment: 29 pages, 8 figures in main text