Motivated by experimental studies that have found signatures of a quantum spin liquid phase in organic crystals whose structure is well described by the two-dimensional triangular lattice, we study the Hubbard model on this lattice at half filling using the infinite-system density matrix renormalization group (iDMRG) method. On infinite cylinders with finite circumference, we identify an intermediate phase between observed metallic behavior at low interaction strength and Mott insulating spin-ordered behavior at strong interactions. Chiral ordering from spontaneous breaking of time-reversal symmetry, a fractionally quantized spin Hall response, and characteristic level statistics in the entanglement spectrum in the intermediate phase provide strong evidence for the existence of a chiral spin liquid in the full two-dimensional limit of the model.

Gauge theories are the cornerstone of our understanding of fundamental interactions among elementary particles. Their properties are often probed in dynamical experiments, such as those performed at ion colliders and high-intensity laser facilities. Describing the evolution of these strongly coupled systems is a formidable challenge for classical computers and represents one of the key open quests for quantum simulation approaches to particle physics phenomena. In this work, we show how recent experiments done on Rydberg atom chains naturally realize the real-time dynamics of a lattice gauge theory at system sizes at the boundary of classical computational methods. We prove that the constrained Hamiltonian dynamics induced by strong Rydberg interactions maps exactly onto the one of a U(1) lattice gauge theory. Building on this correspondence, we show that the recently observed anomalously slow dynamics corresponds to a string-inversion mechanism, reminiscent of the string breaking typically observed in gauge theories. This underlies the generality of this slow dynamics, which we illustrate in the context of one-dimensional quantum electrodynamics on the lattice. Within the same platform, we propose a set of experiments that generically show long-lived oscillations, including the evolution of particle-antiparticle pairs, and discuss how a tunable topological angle can be realized, further affecting the dynamics following a quench. Our work shows that the state of the art for quantum simulation of lattice gauge theories is at 51 qubits and connects the recently observed slow dynamics in atomic systems to archetypal phenomena in particle physics.

A Kitaev quantum spin liquid is a prime example of novel quantum magnetism of spin-orbit entangled pseudospin-$1/2$ moments in a honeycomb lattice. Most candidate materials such as ${\mathrm{Na}}_2{\mathrm{IrO}}_3$ have many competing exchange interactions beyond the minimal Kitaev-Heisenberg model whose small variations in the strength of the interactions produce huge differences in low-energy dynamics. Our incomplete knowledge of dynamic spin correlations hampers identification of a minimal model and quantification of the proximity to the Kitaev quantum spin-liquid phase. Here, we report momentum- and energy-resolved magnetic excitation spectra in a honeycomb lattice ${\mathrm{Na}}_2{\mathrm{IrO}}_3$ measured using a resonant inelastic x-ray scattering spectrometer capable of 12 meV resolution. Measured spectra at a low temperature show that the dynamic response lacks resolution-limited coherent spin waves in most parts of the Brillouin zone but has a discernible dispersion and spectral weight distribution within the energy window of 60 meV. A systematic investigation using the exact diagonalization method and direct comparison of high-resolution experimental spectra and theoretical simulations allow us to confine a parameter regime in which the extended Kitaev-Heisenberg model reasonably reproduces the main feature of the observed magnetic excitations. Hidden Kitaev quantum spin-liquid and Heisenberg phases found in the complex parameter space are used as references to propose the picture of renormalized magnons as explaining the incoherent nature of magnetic excitations. Magnetic excitation spectra are taken at elevated temperatures to follow the temperature evolution of the resonant inelastic x-ray scattering dynamic response in the paramagnetic state. Whereas the low-energy excitation progressively diminishes as the zigzag order disappears, the broad high-energy excitation maintains its spectral weight up to a much higher temperature of 160 K. We suggest that the dominant nearest-neighbor interactions keep short-range correlations up to quite high temperatures with a specific short-range dynamics which has a possible connection to a proximate spin-liquid phase.

We describe the phase diagram of electrons on a fully connected lattice with random hopping, subject to a random Heisenberg spin exchange interaction between any pair of sites and a constraint of no double occupancy. A perturbative renormalization group analysis yields a critical point with fractionalized excitations at a nonzero critical value $p_c$ of the hole doping $p$ away from the half-filled insulator. We compute the renormalization group to two loops, but some exponents are obtained to all loop order. We argue that the critical point $p_c$ is flanked by confining phases: a disordered Fermi liquid with carrier density $1+p$ for $p>p_c$ and a metallic spin glass with carrier density $p$ for $p<p_c$. Additional evidence for the critical behavior is obtained from a large-$M$ analysis of a model which extends the SU(2) spin symmetry to $\mathrm{SU}(M)$. We discuss the relationship of the vicinity of this deconfined quantum critical point to key aspects of cuprate phenomenology.

We show how a large family of interacting nonequilibrium phases of matter can arise from the presence of multiple time-translation symmetries, which occur by quasiperiodically driving an isolated, quantum many-body system with two or more incommensurate frequencies. These phases are fundamentally different from those realizable in time-independent or periodically driven (Floquet) settings. Focusing on high-frequency drives with smooth time dependence, we rigorously establish general conditions for which these phases are stable in a parametrically long-lived “preheating” regime. We develop a formalism to analyze the effect of the multiple time-translation symmetries on the dynamics of the system, which we use to classify and construct explicit examples of the emergent phases. In particular, we discuss time quasicrystals which spontaneously break the time-translation symmetries, as well as time-translation symmetry-protected topological phases.

Epitaxial strain provides important pathways to control the magnetic and electronic states in transition-metal oxides. However, the large strain is usually accompanied by a strong reduction of the oxygen-vacancy formation energy, which hinders the direct manipulation of their intrinsic properties. Here, using a postdeposition ozone annealing method, we obtain a series of oxygen stoichiometric ${\mathrm{SrCoO}}_3$ thin films with the tensile strain up to 3.0%. We observe a robust ferromagnetic ground state in all strained thin films, while interestingly the tensile strain triggers a distinct metal-to-insulator transition along with the increase of the tensile strain. The persistent ferromagnetic state across the electrical transition therefore suggests that the magnetic state is directly correlated with the localized electrons, rather than the itinerant ones, which then calls for further investigation of the intrinsic mechanism of this magnetic compound beyond the double-exchange mechanism.

Nonequilibrium phenomena are ubiquitous in nature and in a wide range of systems, including cold atomic gases and solid-state materials. While these phenomena are challenging to describe both theoretically and experimentally, they are essential for the fundamental understanding of many-body systems and practical devices. In the context of spintronics, when a magnetic insulator (MI) is subjected to a thermal gradient, a pure spin current is generated in the form of magnons without the presence and dissipation of a charge current—attractive for reducing energy consumption and central to the emerging field of spin caloritronics. However, the experimental methods for directly quantifying a spin current in insulators and for probing local phonon-magnon nonequilibrium and the associated magnon chemical potential are largely missing. Here, we apply a heating laser to generate a thermal gradient in the MI yttrium iron garnet (YIG), ${\mathrm{Y}}_3{\mathrm{Fe}}_5{\mathrm{O}}_{12}$, and evaluate two components of the spin current, driven by temperature and chemical potential gradients, respectively. The experimental method and theory approach for evaluating quasiparticle chemical potential can be applied for analogous phenomena in other many-body systems.

The quench dynamics of a system involving two competing orders is investigated using a Ginzburg-Landau theory with relaxational dynamics. We consider the scenario where a pump rapidly heats the system to a high temperature, after which the system cools down to its equilibrium temperature. We study the evolution of the order parameter amplitude and fluctuations in the resulting time-dependent free-energy landscape. Exponentially growing thermal fluctuations dominate the dynamics. The system typically evolves into the phase associated with the faster-relaxing order parameter, even if it is not the global free-energy minimum. This theory offers a natural explanation for the widespread experimental observation that metastable states may be induced by laser-induced collapse of a dominant equilibrium order parameter.

Phonon polaritons localized in polar nanoresonators and superlattices are being actively investigated as promising platforms for midinfrared nanophotonics. Here we show that the nonlocal nature of the phonon response can strongly modify their nanoscale physics. Using a nonlocal dielectric approach, we study dielectric nanospheres and thin dielectric films taking into account optical phonons dispersion. We discover a rich nonlocal phenomenology, qualitatively different from the one of plasmonic systems. Our theory allows us to explain the recently reported discrepancy between theory and experiments in atomic-scale superlattices, and it provides a practical tool for the design of phonon-polariton nanodevices.

In the absence of spin-orbit coupling, the conventional dogma of Anderson localization asserts that all states localize in two dimensions, with a glaring exception: the quantum Hall plateau transition (QHPT). In that case, the localization length diverges and interference-induced quantum-critical spatial fluctuations appear at all length scales. Normally, QHPT states occur only at isolated energies; accessing them therefore requires fine-tuning of the electron density or magnetic field. In this paper we show that QHPT states can be realized throughout an energy continuum, i.e., as an “energy stack” of critical states wherein each state in the stack exhibits QHPT phenomenology. The stacking occurs without fine-tuning at the surface of a class AIII topological phase, where it is protected by U(1) and (anomalous) chiral or time-reversal symmetries. Spectrum-wide criticality is diagnosed by comparing numerics to universal results for the longitudinal Landauer conductance and wave function multifractality at the QHPT. Results are obtained from an effective 2D surface field theory and from a bulk 3D lattice model. We demonstrate that the stacking of quantum-critical QHPT states is a robust phenomenon that occurs for AIII topological phases with both odd and even winding numbers. The latter conclusion may have important implications for the still poorly understood logarithmic conformal field theory believed to describe the QHPT.

Quantized energy levels are one of the hallmarks of quantum mechanics at the atomic level. The manifestation of quantization in macroscopic physical systems has showcased important quantum phenomena, such as quantized conductance in (fractional) quantum Hall effects and quantized vortices in superconductors. Here we report the first experimental observation of quantized exciton energies in a macroscopic system with strong Coulomb interaction, monolayer ${\mathrm{WSe}}_2$ crystal under a strong magnetic field. Employing helicity-resolved magnetoreflectance spectroscopy, we observe a striking ladder of plateaus as a function of the gate voltage for both exciton resonance in one valley and exciton-polariton branch in the opposite valley, thanks to the inter-Landau levels transitions governed by unique valley-selective selection rules. The observed quantized excitation energy level spacing sensitively depends on the doping level, indicating strong many-body effects. Our work will inspire the study of intriguing quantum phenomena originating from the interplay between Landau levels and many-body interactions in two-dimension monolayer crystals.

A microwave amplifier combining noise performances as close as possible to the quantum limit with large bandwidth and high saturation power is highly desirable for many solid-state quantum technologies. Here, we introduce a new traveling-wave parametric amplifier based on superconducting quantum interference devices. It displays a 3-GHz bandwidth, a $-100$-dBm saturation (1-dB compression) point and added noise near the quantum limit. Compared to the previous state of the art, it is an order of magnitude more compact, its characteristic impedance is in situ tunable, and its fabrication process requires only two lithography steps. The key is the engineering of a gap in the dispersion relation of the transmission line. This is obtained using a periodic modulation of the SQUID size, similarly to what is done with photonic crystals. Moreover, we provide a new theoretical treatment to describe the nontrivial interplay between nonlinearity and such periodicity. Our approach provides a path to cointegration with other quantum devices such as qubits given the low footprint and easy fabrication of our amplifier.

We introduce a complex-plane generalization of the consecutive level-spacing ratio distribution used to distinguish regular from chaotic quantum spectra. Our approach features the distribution of complex-valued ratios between nearest- and next-to-nearest-neighbor spacings. We show that this quantity can successfully detect the chaotic or regular nature of complex-valued spectra, which is done in two steps. First, we show that, if eigenvalues are uncorrelated, the distribution of complex spacing ratios is flat within the unit circle, whereas random matrices show a strong angular dependence in addition to the usual level repulsion. The universal fluctuations of Gaussian unitary and Ginibre unitary universality classes in the large-matrix-size limit are shown to be well described by Wigner-like surmises for small-size matrices with eigenvalues on the circle and on the two-torus, respectively. To study the latter case, we introduce the toric unitary ensemble, characterized by a flat joint eigenvalue distribution on the two-torus. Second, we study different physical situations where non-Hermitian matrices arise: dissipative quantum systems described by a Lindbladian, nonunitary quantum dynamics described by non-Hermitian Hamiltonians, and classical stochastic processes. We show that known integrable models have a flat distribution of complex spacing ratios, whereas generic cases, expected to be chaotic, conform to random matrix theory predictions. Specifically, we are able to clearly distinguish chaotic from integrable dynamics in boundary-driven dissipative spin-chain Liouvillians and in the classical asymmetric simple exclusion process and to differentiate localized from delocalized regimes in a non-Hermitian disordered many-body system.

The delafossite metals ${\mathrm{PdCoO}}_2$, ${\mathrm{PtCoO}}_2$, and ${\mathrm{PdCrO}}_2$ are among the highest conductivity materials known, with low-temperature mean free paths of tens of microns in the best as-grown single crystals. A key question is whether these very low resistive scattering rates result from strongly suppressed backscattering due to special features of the electronic structure or are a consequence of highly unusual levels of crystalline perfection. We report the results of experiments in which high-energy electron irradiation was used to introduce point disorder to the Pd and Pt layers in which the conduction occurs. We obtain the cross section for formation of Frenkel pairs in absolute units, and cross-check our analysis with first-principles calculations of the relevant atomic displacement energies. We observe an increase of resistivity that is linear in defect density with a slope consistent with scattering in the unitary limit. Our results enable us to deduce that the as-grown crystals contain extremely low levels of in-plane defects of approximately 0.001%. This confirms that crystalline perfection is the most important factor in realizing the long mean free paths and highlights how unusual these delafossite metals are in comparison with the vast majority of other multicomponent oxides and alloys. We discuss the implications of our findings for future materials research.

A grand challenge underlies the entire field of topology-enabled quantum logic and information science: how to establish topological control principles driven by quantum coherence and understand the time dependence of such periodic driving. Here we demonstrate a few-cycle THz-pulse-induced phase transition in a Dirac semimetal ${\mathrm{ZrTe}}_5$ that is periodically driven by vibrational coherence due to excitation of the lowest Raman active mode. Above a critical THz-pump field threshold, there emerges a long-lived metastable phase, approximately 100 ps, with unique Raman phonon-assisted topological switching dynamics absent for optical pumping. The switching also manifests itself by distinct features: nonthermal spectral shape, relaxation slowing near the Lifshitz transition where the critical Dirac point occurs, and diminishing signals at the same temperature that the Berry-curvature-induced anomalous Hall effect magnetoresistance vanishes. These results, together with first-principles modeling, identify a mode-selective Raman coupling that drives the system from strong to weak topological insulators with a Dirac semimetal phase established at a critical atomic displacement controlled by the phonon coherent pumping. Harnessing of vibrational coherence can be extended to steer symmetry-breaking transitions, i.e., Dirac to Weyl ones, with implications for THz topological quantum gate and error correction applications.

Two-dimensional semiconductors provide an ideal platform for exploration of linear exciton and polariton physics, primarily due to large exciton binding energy and strong light-matter coupling. These features, however, generically imply reduced exciton-exciton interactions, hindering the realization of active optical devices such as lasers or parametric oscillators. Here, we show that electrical injection of itinerant electrons into monolayer molybdenum diselenide allows us to overcome this limitation: dynamical screening of exciton-polaritons by electrons leads to the formation of new quasiparticles termed polaron-polaritons that exhibit unexpectedly strong interactions as well as optical amplification by Bose-enhanced polaron-electron scattering. To measure the nonlinear optical response, we carry out time-resolved pump-probe measurements and observe polaron-polariton interaction enhancement by a factor of 50 ($0.5\mu \mathrm{eV}\mu {\mathrm{m}}^2$) as compared to exciton-polaritons. Concurrently, we measure a spectrally integrated transmission gain of the probe field of $\gtrsim 2$ stemming from stimulated scattering of polaron-polaritons. We show theoretically that the nonequilibrium nature of optically excited quasiparticles favors a previously unexplored interaction mechanism stemming from a phase-space filling in the screening cloud, which provides an accurate explanation of the strong repulsive interactions observed experimentally. Our findings show that itinerant electron-exciton interactions provide an invaluable tool for electronic manipulation of optical properties, demonstrate a new mechanism for dramatically enhancing polariton-polariton interactions, and pave the way for realization of nonequilibrium polariton condensates.

In contrast to magnetic order formed by electrons’ dipolar moments, ordering phenomena associated with higher-order multipoles (quadrupoles, octupoles, etc.) are more difficult to characterize because of the limited choice of experimental probes that can distinguish different multipolar moments. The heavy-fermion compound ${\mathrm{CeB}}_6$ and its La-diluted alloys are among the best-studied realizations of the long-range-ordered multipolar phases, often referred to as “hidden order.” Previously, the hidden order in phase II was identified as primary antiferroquadrupolar and field-induced octupolar order. Here, we present a combined experimental and theoretical investigation of collective excitations in phase II of ${\mathrm{CeB}}_6$. Inelastic neutron scattering (INS) in fields up to 16.5 T reveals a new high-energy mode above 14 T in addition to the low-energy magnetic excitations. The experimental dependence of their energy on the magnitude and angle of the applied magnetic field is compared to the results of a multipolar interaction model. The magnetic excitation spectrum in a rotating field is calculated within a localized approach using the pseudospin representation for the ${\mathrm{\Gamma }}_8$ states. We show that the rotating-field technique at fixed momentum can complement conventional INS measurements of the dispersion at a constant field and holds great promise for identifying the symmetry of multipolar order parameters and the details of intermultipolar interactions that stabilize hidden-order phases.

Both helium and ammonia are main components of icy giant planets. While ammonia is very reactive, helium is the most inert element in the universe. It is of great interest whether ammonia and helium can react with each other under planetary conditions, and if so, what kinds of structures and states of matter can form. Here, using crystal structure prediction methods and first-principles calculations, we report three new stable stoichiometries and eight new stable phases of $\mathrm{He}\text{−}{\mathrm{NH}}_3$ compounds under pressures up to 500 GPa. These structures may exhibit perovskitelike structures for ${\mathrm{HeNH}}_3$ and ${\mathrm{He}}_2{\mathrm{NH}}_3$, and a host-guest crystal structure for $\mathrm{He}{({\mathrm{NH}}_3)}_2$. Superionic states are found in all these $\mathrm{He}\text{−}{\mathrm{NH}}_3$ compounds under high pressures and temperatures in which the hydrogen atoms are diffusive while the nitrogen and helium atoms remain fixed. Such dynamical behavior in helium ammonia compounds is quite different from that in helium water compounds, where weakly interacting helium is more diffusive than stronger bound hydrogen. The low-density host-guest phase of space group $I4cm$ is found to be stable at very low pressures (about 3 GPa) and it enters into a plastic state, characterized by freely rotating ammonia molecules. The present results suggest that plastic or superionic helium ammonia compounds may exist under planetary conditions, and helium contributes crucially to the exotic physics and chemistry observed under extreme conditions.

Inspired by recent experiments on graphene, we examine the nondissipative viscoelastic response of anisotropic two-dimensional quantum systems. We pay particular attention to electron fluids with point group symmetries and those with discrete translational symmetry. We start by extending the Kubo formalism for viscosity to systems with internal degrees of freedom and discrete translational symmetry, highlighting the importance of properly considering the role of internal angular momentum. We analyze the Hall components of the viscoelastic response tensor in systems with discrete point group symmetry, focusing on the hydrodynamic implications of the resulting forces. We show that though there are generally six Hall viscosities, there are only three independent contributions to the viscous force density in the bulk. To compute these coefficients, we develop a framework to consistently write down the long-wavelength stress tensor and viscosity for multicomponent lattice systems. We apply our formalism to lattice and continuum models, including a lattice Chern insulator and anisotropic superfluid.

Exploring new parameter regimes to realize and control novel phases of matter has been a main theme in modern condensed matter physics research. The recent discovery of two-dimensional (2D) magnetism in nearly freestanding monolayer atomic crystals has already led to observations of a number of novel magnetic phenomena absent in bulk counterparts. Such intricate interplays between magnetism and crystalline structures provide ample opportunities for exploring quantum phase transitions in this new 2D parameter regime. Here, using magnetic field- and temperature-dependent circularly polarized Raman spectroscopy of phonons and magnons, we map out the phase diagram of chromium triiodide (${\mathrm{CrI}}_3$) that has been known to be a layered antiferromagnet (AFM) in its 2D films and a ferromagnet (FM) in its three-dimensional (3D) bulk. However, we reveal a novel mixed state of layered AFM and FM in 3D ${\mathrm{CrI}}_3$ bulk crystals where the layered AFM survives in the surface layers, and the FM appears in deeper bulk layers. We then show that the surface-layered AFM transits into the FM at a critical magnetic field of 2 T, similar to what was found in the few-layer case. Interestingly, concurrent with this magnetic phase transition, we discover a first-order structural phase transition that alters the crystallographic point group from $C_{3i}$ (rhombohedral) to $C_{2h}$ (monoclinic). Our result not only unveils the complex single-magnon behavior in 3D ${\mathrm{CrI}}_3$, but it also settles the puzzle of how ${\mathrm{CrI}}_3$ transits from a bulk FM to a thin-layered AFM semiconductor, despite recent efforts in understanding the origin of layered AFM in ${\mathrm{CrI}}_3$ thin layers, and reveals the intimate relationship between the layered AFM-to-FM and the crystalline rhombohedral-to-monoclinic phase transitions. These findings further open opportunities for future 2D magnet-based magnetomechanical devices.