The classification and construction of symmetry-protected topological (SPT) phases in interacting boson and fermion systems have become a fascinating theoretical direction in recent years. It has been shown that (generalized) group cohomology theory or cobordism theory gives rise to a complete classification of SPT phases in interacting boson or spin systems. The construction and classification of SPT phases in interacting fermion systems are much more complicated, especially in three dimensions. In this work, we revisit this problem based on an equivalence class of fermionic symmetric local unitary transformations. We construct very general fixed-point SPT wave functions for interacting fermion systems. We naturally reproduce the partial classifications given by special group supercohomology theory, and we show that with an additional $\tilde{B}{H}^2(G_b,{\mathbb{Z}}_2)$ structure [the so-called obstruction-free subgroup of $H^2(G_b,{\mathbb{Z}}_2)$], a complete classification of SPT phases for three-dimensional interacting fermion systems with a total symmetry group $G_f=G_b\times {\mathbb{Z}}_2^f$ can be obtained for unitary symmetry group $G_b$. We also discuss the procedure for deriving a general group supercohomology theory in arbitrary dimensions.

We study Abelian braiding statistics of loop excitations in three-dimensional gauge theories with fermionic particles and the closely related problem of classifying 3D fermionic symmetry-protected topological (FSPT) phases with unitary symmetries. It is known that the two problems are related by turning FSPT phases into gauge theories through gauging the global symmetry of the former. We show that there exist certain types of Abelian loop braiding statistics that are allowed only in the presence of fermionic particles, which correspond to 3D “intrinsic” FSPT phases, i.e., those that do not stem from bosonic SPT phases. While such intrinsic FSPT phases are ubiquitous in 2D systems and in 3D systems with antiunitary symmetries, their existence in 3D systems with unitary symmetries was not confirmed previously due to the fact that strong interaction is necessary to realize them. We show that the simplest unitary symmetry to support 3D intrinsic FSPT phases is ${\mathbb{Z}}_2\times {\mathbb{Z}}_4$. To establish the results, we first derive a complete set of physical constraints on Abelian loop braiding statistics. Solving the constraints, we obtain all possible Abelian loop braiding statistics in 3D gauge theories, including those that correspond to intrinsic FSPT phases. Then, we construct exactly soluble state-sum models to realize the loop braiding statistics. These state-sum models generalize the well-known Crane-Yetter and Dijkgraaf-Witten models.

We compute the effects of generic short-range interactions on gapless electrons residing at the quantum critical point separating a two-dimensional Dirac semimetal and a symmetry-preserving band insulator. The electronic dispersion at this critical point is anisotropic ($E_{\mathbf{k}}=\pm \sqrt{v^2{k}_x^2+b^2{k}_y^{2n}}$ with $n=2$), which results in unconventional scaling of thermodynamic and transport quantities. Because of the vanishing density of states [$\varrho (E)\sim |E{|}^{1/n}$], this anisotropic semimetal (ASM) is stable against weak short-range interactions. However, for stronger interactions, the direct Dirac-semimetal to band-insulator transition can either (i) become a fluctuation-driven first-order transition (although unlikely in a particular microscopic model considered here, the anisotropic honeycomb lattice extended Hubbard model) or (ii) get avoided by an intervening broken-symmetry phase. We perform a controlled renormalization group analysis with the small parameter $\epsilon =1/n$, augmented with a $1/n$ expansion (parametrically suppressing quantum fluctuations in the higher dimension) by perturbing away from the one-dimensional limit, realized by setting $\epsilon =0$ and $n\to \infty$. We identify charge density wave (CDW), antiferromagnet (AFM), and singlet $s$-wave superconductivity as the three dominant candidates for broken symmetry. The onset of any such order at strong coupling $(\sim \epsilon )$ takes place through a continuous quantum phase transition across an interacting multicritical point, where the ordered phase, band insulator, Dirac, and anisotropic semimetals meet. We also present the phase diagram of an extended Hubbard model for the ASM, obtained via the controlled deformation of its counterpart in one dimension. The latter displays spin-charge separation and instabilities to CDW, spin density wave, and Luther-Emery liquid phases at arbitrarily weak coupling. The spin density wave and Luther-Emery liquid phases deform into pseudospin SU(2)-symmetric quantum critical points separating the ASM from the AFM and superconducting orders, respectively. Our phase diagram shows an intriguing interplay among CDW, AFM, and $s$-wave paired states that can be germane for a uniaxially strained optical honeycomb lattice for ultracold fermion atoms, or the organic compound $\alpha \text{−}(\mathrm{BEDT}\text{−}\mathrm{TTF}{)}_2{\mathrm{I}}_3$.

The strongly correlated insulator ${\mathrm{Ca}}_2{\mathrm{RuO}}_4$ is considered as a paradigmatic realization of both spin-orbital physics and a band-Mott insulating phase, characterized by orbitally selective coexistence of a band and a Mott gap. We present a high resolution oxygen $K$-edge resonant inelastic x-ray scattering study of the antiferromagnetic Mott insulating state of ${\mathrm{Ca}}_2{\mathrm{RuO}}_4$. A set of low-energy (about 80 and 400 meV) and high-energy (about 1.3 and 2.2 eV) excitations are reported, which show strong incident light polarization dependence. Our results strongly support a spin-orbit coupled band-Mott scenario and explore in detail the nature of its exotic excitations. Guided by theoretical modeling, we interpret the low-energy excitations as a result of composite spin-orbital excitations. Their nature unveils the intricate interplay of crystal-field splitting and spin-orbit coupling in the band-Mott scenario. The high-energy excitations correspond to intra-atomic singlet-triplet transitions at an energy scale set by Hund’s coupling. Our findings give a unifying picture of the spin and orbital excitations in the band-Mott insulator ${\mathrm{Ca}}_2{\mathrm{RuO}}_4$.

When a mobile hole is moving in an antiferromagnet it distorts the surrounding Néel order and forms a magnetic polaron. Such interplay between hole motion and antiferromagnetism is believed to be at the heart of high-temperature superconductivity in cuprates. In this article, we study a single hole described by the $t\text{−}{J}_z$ model with Ising interactions between the spins in two dimensions. This situation can be experimentally realized in quantum gas microscopes with Mott insulators of Rydberg-dressed bosons or fermions, or using polar molecules. We work at strong couplings, where hole hopping is much larger than couplings between the spins. In this regime we find strong theoretical evidence that magnetic polarons can be understood as bound states of two partons, a spinon and a holon carrying spin and charge quantum numbers, respectively. Starting from first principles, we introduce a microscopic parton description which is benchmarked by comparison with results from advanced numerical simulations. Using this parton theory, we predict a series of excited states that are invisible in the spectral function and correspond to rotational excitations of the spinon-holon pair. This is reminiscent of mesonic resonances observed in high-energy physics, which can be understood as rotating quark-antiquark pairs carrying orbital angular momentum. Moreover, we apply the strong-coupling parton theory to study far-from-equilibrium dynamics of magnetic polarons observable in current experiments with ultracold atoms. Our work supports earlier ideas that partons in a confining phase of matter represent a useful paradigm in condensed-matter physics and in the context of high-temperature superconductivity in particular. While direct observations of spinons and holons in real space are impossible in traditional solid-state experiments, quantum gas microscopes provide a new experimental toolbox. We show that, using this platform, direct observations of partons in and out of equilibrium are now possible. Extensions of our approach to the $t\text{−}J$ model are also discussed. Our predictions in this case are relevant to current experiments with quantum gas microscopes for ultracold atoms.

We investigate the spin of a multielectron GaAs quantum dot in a sequence of nine charge occupancies, by exchange coupling the multielectron dot to a neighboring two-electron double quantum dot. For all nine occupancies, we make use of a leakage spectroscopy technique to reconstruct the spectrum of spin states in the vicinity of the interdot charge transition between a single- and a multielectron quantum dot. In the same regime we also perform time-resolved measurements of coherent exchange oscillations between the single- and multielectron quantum dot. With these measurements, we identify distinct characteristics of the multielectron spin state, depending on whether the dot’s occupancy is even or odd. For three out of four even occupancies, we do not observe any exchange interaction with the single quantum dot, indicating a spin-0 ground state. For the one remaining even occupancy, we observe an exchange interaction that we associate with a spin-1 multielectron quantum dot ground state. For all five of the odd occupancies, we observe an exchange interaction associated with a spin-$1/2$ ground state. For three of these odd occupancies, we clearly demonstrate that the exchange interaction changes sign in the vicinity of the charge transition. For one of these, the exchange interaction is negative (i.e., triplet preferring) beyond the interdot charge transition, consistent with the observed spin-1 for the next (even) occupancy. Our experimental results are interpreted through the use of a Hubbard model involving two orbitals of the multielectron quantum dot. Allowing for the spin correlation energy (i.e., including a term favoring Hund’s rules) and different tunnel coupling to different orbitals, we qualitatively reproduce the measured exchange profiles for all occupancies.

We put the theory of interacting topological crystalline phases on a systematic footing. These are topological phases protected by space-group symmetries. Our central tool is an elucidation of what it means to “gauge” such symmetries. We introduce the notion of a crystalline topological liquid and argue that most (and perhaps all) phases of interest are likely to satisfy this criterion. We prove a crystalline equivalence principle, which states that in Euclidean space, crystalline topological liquids with symmetry group $G$ are in one-to-one correspondence with topological phases protected by the same symmetry $G$, but acting internally, where if an element of $G$ is orientation reversing, it is realized as an antiunitary symmetry in the internal symmetry group. As an example, we explicitly compute, using group cohomology, a partial classification of bosonic symmetry-protected topological phases protected by crystalline symmetries in ($3+1$) dimensions for 227 of the 230 space groups. For the 65 space groups not containing orientation-reversing elements (Sohncke groups), there are no cobordism invariants that may contribute phases beyond group cohomology, so we conjecture that our classification is complete.

Uniform triangular crystals are the ground state for particles that interact isotropically in two dimensions. However, when immersed in an external potential, for example, one arising from an electric field, a flow field, or gravity, the resulting phases are significantly distorted in a way reminiscent of conformal transformations of planar lattices. We study these “conformal crystals” using colloidal experiments and molecular dynamics simulations. By establishing a projection from these self-assembled inhomogeneous crystals to homogeneous crystals on curved surfaces, we are able to both predict the distribution of defects and establish that defects are an almost inevitable part of the ground state. We determine how the inherent geometry emerges from an interplay between the confining potential and the interparticle interactions. Using molecular dynamics simulations, we demonstrate the generic behavior of this emergent geometry and the resulting defect structures throughout a variety of physical systems.

Anyons are exotic quasiparticles with fractional charge that can emerge as fundamental excitations of strongly interacting topological quantum phases of matter. Unlike ordinary fermions and bosons, they may obey non-Abelian statistics—a property that would help realize fault-tolerant quantum computation. Non-Abelian anyons have long been predicted to occur in the fractional quantum Hall (FQH) phases that form in two-dimensional electron gases in the presence of a large magnetic field, such as the $\nu =5/2$ FQH state. However, direct experimental evidence of anyons and tests that can distinguish between Abelian and non-Abelian quantum ground states with such excitations have remained elusive. Here, we propose a new experimental approach to directly visualize the structure of interacting electronic states of FQH states with the STM. Our theoretical calculations show how spectroscopy mapping with the STM near individual impurity defects can be used to image fractional statistics in FQH states, identifying unique signatures in such measurements that can distinguish different proposed ground states. The presence of locally trapped anyons should leave distinct signatures in STM spectroscopic maps, and enables a new approach to directly detect—and perhaps ultimately manipulate—these exotic quasiparticles.

The mixedness of a quantum state is usually seen as an adversary to topological quantization of observables. For example, exact quantization of the charge transported in a so-called Thouless adiabatic pump is lifted at any finite temperature in symmetry-protected topological insulators. Here, we show that certain directly observable many-body correlators preserve the integrity of topological invariants for mixed Gaussian quantum states in one dimension. Our approach relies on the expectation value of the many-body momentum-translation operator and leads to a physical observable—the “ensemble geometric phase” (EGP)—which represents a bona fide geometric phase for mixed quantum states, in the thermodynamic limit. In cyclic protocols, the EGP provides a topologically quantized observable that detects encircled spectral singularities (“purity-gap” closing points) of density matrices. While we identify the many-body nature of the EGP as a key ingredient, we propose a conceptually simple, interferometric setup to directly measure the latter in experiments with mesoscopic ensembles of ultracold atoms.

Origami structures are characterized by a network of folds and vertices joining unbendable plates. For applications to mechanical design and self-folding structures, it is essential to understand the interplay between the set of folds in the unfolded origami and the possible 3D folded configurations. When deforming a structure that has been folded, one can often linearize the geometric constraints, but the degeneracy of the unfolded state makes a linear approach impossible there. We derive a theory for the second-order infinitesimal rigidity of an initially unfolded triangulated origami structure and use it to study the set of nearly unfolded configurations of origami with four boundary vertices. We find that locally, this set consists of a number of distinct “branches” which intersect at the unfolded state, and that the number of these branches is exponential in the number of vertices. We find numerical and analytical evidence that suggests that the branches are characterized by choosing each internal vertex to either “pop up” or “pop down.” The large number of pathways along which one can fold an initially unfolded origami structure strongly indicates that a generic structure is likely to become trapped in a “misfolded” state. Thus, new techniques for creating self-folding origami are likely necessary; controlling the popping state of the vertices may be one possibility.

We study the topological properties of superconductors with paired $j=\frac{3}{2}$ quasiparticles. Higher spin Fermi surfaces can arise, for instance, in strongly spin-orbit coupled band-inverted semimetals. Examples include the Bi-based half-Heusler materials, which have recently been established as low-temperature and low-carrier density superconductors. Motivated by this experimental observation, we obtain a comprehensive symmetry-based classification of topological pairing states in systems with higher angular momentum Cooper pairing. Our study consists of two main parts. First, we develop the phenomenological theory of multicomponent (i.e., higher angular momentum) pairing by classifying the stationary points of the free energy within a Ginzburg-Landau framework. Based on the symmetry classification of stationary pairing states, we then derive the symmetry-imposed constraints on their gap structures. We find that, depending on the symmetry quantum numbers of the Cooper pairs, different types of topological pairing states can occur: fully gapped topological superconductors in class DIII, Dirac superconductors, and superconductors hosting Majorana fermions. Notably, we find a series of nematic fully gapped topological superconductors, as well as double- and triple-Dirac superconductors, with quadratic and cubic dispersion, respectively. Our approach, applied here to the case of $j=\frac{3}{2}$ Cooper pairing, is rooted in the symmetry properties of pairing states, and can therefore also be applied to other systems with higher angular momentum and high-spin pairing. We conclude by relating our results to experimentally accessible signatures in thermodynamic and dynamic probes.

The modern semiclassical theory of a Bloch electron in a magnetic field now encompasses the orbital magnetic moment and the geometric phase. These two notions are encoded in the Bohr-Sommerfeld quantization condition as a phase ($\lambda$) that is subleading in powers of the field; $\lambda$ is measurable in the phase offset of the de Haas–van Alphen oscillation, as well as of fixed-bias oscillations of the differential conductance in tunneling spectroscopy. In some solids and for certain field orientations, $\lambda /\pi$ are robustly integer valued, owing to the symmetry of the extremal orbit; i.e., they are the topological invariants of magnetotransport. Our comprehensive symmetry analysis identifies solids in any (magnetic) space group for which $\lambda$ is a topological invariant, as well as the symmetry-enforced degeneracy of Landau levels. The analysis is simplified by our formulation of ten (and only ten) symmetry classes for closed, Fermi-surface orbits. Case studies are discussed for graphene, transition metal dichalcogenides, 3D Weyl and Dirac metals, and crystalline and ${\mathbb{Z}}_2$ topological insulators. In particular, we point out that a $\pi$ phase offset in the fundamental oscillation should not be viewed as a smoking gun for a 3D Dirac metal.

Massless $(2+1)\mathrm{D}$ Dirac fermions arise in a variety of systems from graphene to the surfaces of topological insulators, where generating a mass is typically associated with breaking a symmetry. However, with strong interactions, a symmetric gapped phase can arise for multiples of eight Dirac fermions. A continuous quantum phase transition from the massless Dirac phase to this massive phase, which we term symmetric mass generation, is necessarily beyond the Landau paradigm and is hard to describe even at the conceptual level. Nevertheless, such transition has been consistently observed in several numerical studies recently. Here, we propose a theory for the symmetric mass generation transition which is reminiscent of deconfined criticality and involves emergent non-Abelian gauge fields coupled both to Dirac fermions and to critical Higgs bosons. We motivate the theory using an explicit parton construction and discuss predictions for numerics. Additionally, we show that the fermion Green’s function is expected to undergo a zero-to-pole transition across the critical point.

We consider the ground-state properties of an impurity particle (“polaron”) resonantly interacting with a Bose-Einstein condensate (BEC). Focusing on the equal-mass system, we use a variational wave function for the polaron that goes beyond previous work and includes up to three Bogoliubov excitations of the BEC, thus allowing us to capture both Efimov trimers and associated tetramers. We find that the length scale associated with Efimov trimers (i.e., the three-body parameter) can strongly affect the polaron’s behavior, even at densities where there are no well-defined Efimov states. However, by comparing our results with recent quantum Monte Carlo calculations, we argue that the polaron energy is a universal function of the Efimov three-body parameter for sufficiently low boson densities. We further support this conclusion by showing that the energies of the deepest bound Efimov trimers and tetramers at unitarity are universally related to one another, regardless of the microscopic model. On the other hand, we find that the quasiparticle residue and effective mass sensitively depend on the coherence length $\xi$ of the BEC, with the residue tending to zero as $\xi$ diverges, in a manner akin to the orthogonality catastrophe.

Harnessing the full power of nascent quantum processors requires the efficient management of a limited number of quantum bits with finite coherent lifetimes. Hybrid algorithms, such as the variational quantum eigensolver (VQE), leverage classical resources to reduce the required number of quantum gates. Experimental demonstrations of VQE have resulted in calculation of Hamiltonian ground states, and a new theoretical approach based on a quantum subspace expansion (QSE) has outlined a procedure for determining excited states that are central to dynamical processes. We use a superconducting-qubit-based processor to apply the QSE approach to the ${\mathrm{H}}_2$ molecule, extracting both ground and excited states without the need for auxiliary qubits or additional minimization. Further, we show that this extended protocol can mitigate the effects of incoherent errors, potentially enabling larger-scale quantum simulations without the need for complex error-correction techniques.

Quantum fluctuations of the Néel state of the square lattice antiferromagnet are usually described by a ${\mathbb{CP}}^1$ theory of bosonic spinons coupled to a U(1) gauge field, and with a global SU(2) spin rotation symmetry. Such a theory also has a confining phase with valence bond solid (VBS) order, and upon including spin-singlet charge-2 Higgs fields, deconfined phases with ${\mathbb{Z}}_2$ topological order possibly intertwined with discrete broken global symmetries. We present dual theories of the same phases starting from a mean-field theory of fermionic spinons moving in $\pi$ flux in each square lattice plaquette. Fluctuations about this $\pi$-flux state are described by ($2+1$)-dimensional quantum chromodynamics (${\mathrm{QCD}}_3$) with a SU(2) gauge group and $N_f=2$ flavors of massless Dirac fermions. It has recently been argued by Wang et al. [Deconfined Quantum Critical Points: Symmetries and Dualities, Phys. Rev. X 7, 031051 (2017).] that this ${\mathrm{QCD}}_3$ theory describes the Néel-VBS quantum phase transition. We introduce adjoint Higgs fields in ${\mathrm{QCD}}_3$ and obtain fermionic dual descriptions of the phases with ${\mathbb{Z}}_2$ topological order obtained earlier using the bosonic ${\mathbb{CP}}^1$ theory. We also present a fermionic spinon derivation of the monopole Berry phases in the U(1) gauge theory of the VBS state. The global phase diagram of these phases contains multicritical points, and our results imply new boson-fermion dualities between critical gauge theories of these points.

Detwinning of magnetic (nematic) domains in Fe-based superconductors has so far only been obtained through mechanical straining, which considerably perturbs the ground state of these materials. The recently discovered nonmechanical detwinning in ${\mathrm{EuFe}}_2{\mathrm{As}}_2$ by ultralow magnetic fields offers an entirely different, nonperturbing way to achieve the same goal. However, this way seemed risky due to the lack of a microscopic understanding of the magnetically driven detwinning. Specifically, the following issues remained unexplained: (i) ultralow value of the first detwinning field of approximately 0.1 T, two orders of magnitude below that of ${\mathrm{BaFe}}_2{\mathrm{As}}_2$, and (ii) reversal of the preferential domain orientation at approximately 1 T and restoration of the low-field orientation above 10–15 T. In this paper, we present, using published as well as newly measured data, a full theory that quantitatively explains all the observations. The key ingredient of this theory is a biquadratic coupling between Fe and Eu spins, analogous to the Fe-Fe biquadratic coupling that drives the nematic transition in this family of materials.

While charge density wave (CDW) instabilities are ubiquitous to superconducting cuprates, the different ordering wave vectors in various cuprate families have hampered a unified description of the CDW formation mechanism. Here, we investigate the temperature dependence of the low-energy phonons in the canonical CDW-ordered cuprate ${\mathrm{La}}_{1.875}{\mathrm{Ba}}_{0.125}{\mathrm{CuO}}_4$. We discover that the phonon softening wave vector associated with CDW correlations becomes temperature dependent in the high-temperature precursor phase and changes from a wave vector of 0.238 reciprocal lattice units (r.l.u.) below the ordering transition temperature to 0.3 r.l.u. at 300 K. This high-temperature behavior shows that “214”-type cuprates can host CDW correlations at a similar wave vector to previously reported CDW correlations in non-214-type cuprates such as ${\mathrm{YBa}}_2{\mathrm{Cu}}_3{\mathrm{O}}_{6+\delta }$. This indicates that cuprate CDWs may arise from the same underlying instability despite their apparently different low-temperature ordering wave vectors.

Made from stacks of two-dimensional materials, van der Waals heterostructures exhibit unique light-matter interactions and are promising for novel optoelectronic devices. The performance of such devices is governed by near-field coupling through, e.g., interlayer charge and/or energy transfer. New concepts and experimental methodologies are needed to properly describe two-dimensional heterointerfaces. Here, we report an original study of interlayer charge and energy transfer in atomically thin metal-semiconductor [i.e., graphene-transition metal dichalcogenide (TMD, here molybdenum diselenide, ${\mathrm{MoSe}}_2$)] heterostructures using a combination of microphotoluminescence and Raman scattering spectroscopies. The photoluminescence intensity in $\mathrm{graphene}/{\mathrm{MoSe}}_2$ is quenched by more than 2 orders of magnitude and rises linearly with the incident photon flux, demonstrating a drastically shortened (about 1 ps) room-temperature ${\mathrm{MoSe}}_2$ exciton lifetime. Key complementary insights are provided from a comprehensive analysis of the graphene and ${\mathrm{MoSe}}_2$ Raman modes, which reveals net photoinduced electron transfer from ${\mathrm{MoSe}}_2$ to graphene and hole accumulation in ${\mathrm{MoSe}}_2$. Remarkably, the steady-state Fermi energy of graphene saturates at $290\pm 15\mathrm{meV}$ above the Dirac point. This reproducible behavior is observed both in ambient air and in vacuum and is discussed in terms of intrinsic factors (i.e., band offsets) and environmental effects. In this saturation regime, balanced photoinduced flows of electrons and holes may transfer to graphene, a mechanism that effectively leads to energy transfer. Using a broad range of incident photon fluxes and diverse environmental conditions, we find that the presence of net photoinduced charge transfer has no measurable impact on the near-unity photoluminescence quenching efficiency in $\mathrm{graphene}/{\mathrm{MoSe}}_2$. This absence of correlation strongly suggests that energy transfer to graphene (either in the form of electron exchange or dipole-dipole interaction) is the dominant interlayer coupling mechanism between atomically thin TMDs and graphene.