Pseudo-fermion functional renormalization group for spin models.

For decades, frustrated quantum magnets have been a seed for scientific progress and innovation in condensed matter. As much as the numerical tools for low-dimensional quantum magnetism have thrived and improved in recent years due to breakthroughs inspired by quantum information and quantum computation, higher-dimensional quantum magnetism can be considered as the final frontier, where strong quantum entanglement, multiple ordering channels, and manifold ways of paramagnetism culminate. At the same time, efforts in crystal synthesis have induced a significant increase in the number of tangible frustrated magnets which are generically three-dimensional in nature, creating an urgent need for quantitative theoretical modeling. We review the pseudo-fermion (PF) and pseudo-Majorana (PM) functional renormalization group (FRG) and their specific ability to address higher-dimensional frustrated quantum magnetism. First developed more than a decade ago, the PFFRG interprets a Heisenberg model Hamiltonian in terms of Abrikosov pseudofermions, which is then treated in a diagrammatic resummation scheme formulated as a renormalization group flow ofm-particle pseudofermion vertices. The article reviews the state of the art of PFFRG and PMFRG and discusses their application to exemplary domains of frustrated magnetism, but most importantly, it makes the algorithmic and implementation details of these methods accessible to everyone. By thus lowering the entry barrier to their application, we hope that this review will contribute towards establishing PFFRG and PMFRG as the numerical methods for addressing frustrated quantum magnetism in higher spatial dimensions.

Topological Fracton Quantum Phase Transitions by Tuning Exact Tensor Network States.

Gapped fracton phases of matter generalize the concept of topological order and broaden our fundamental understanding of entanglement in quantum many-body systems. However, their analytical or numerical description beyond exactly solvable models remains a formidable challenge. Here we employ an exact 3D quantum tensor-network approach that allows us to study a Z_{N} generalization of the prototypical X cube fracton model and its quantum phase transitions between distinct topological states via fully tractable wave function deformations. We map the (deformed) quantum states exactly to a combination of a classical lattice gauge theory and a plaquette clock model, and employ numerical techniques to calculate various entanglement order parameters. For the Z_{N} model we find a family of (weakly) first-order fracton confinement transitions that in the limit of N→∞ converge to a continuous phase transition beyond the Landau-Ginzburg-Wilson paradigm. We also discover a line of 3D conformal quantum critical points (with critical magnetic flux loop fluctuations) which, in the N→∞ limit, appears to coexist with a gapless deconfined fracton state.

Nishimori's Cat: Stable Long-Range Entanglement from Finite-Depth Unitaries and Weak Measurements.

In the field of monitored quantum circuits, it has remained an open question whether finite-time protocols for preparing long-range entangled states lead to phases of matter that are stable to gate imperfections, that can convert projective into weak measurements. Here, we show that in certain cases, long-range entanglement persists in the presence of weak measurements, and gives rise to novel forms of quantum criticality. We demonstrate this explicitly for preparing the two-dimensional Greenberger-Horne-Zeilinger cat state and the three-dimensional toric code as minimal instances. In contrast to random monitored circuits, our circuit of gates and measurements is deterministic; the only randomness is in the measurement outcomes. We show how the randomness in these weak measurements allows us to track the solvable Nishimori line of the random-bond Ising model, rigorously establishing the stability of the glassy long-range entangled states in two and three spatial dimensions. Away from this exactly solvable construction, we use hybrid tensor network and Monte Carlo simulations to obtain a nonzero Edwards-Anderson order parameter as an indicator of long-range entanglement in the two-dimensional scenario. We argue that our protocol admits a natural implementation in existing quantum computing architectures, requiring only a depth-3 circuit on IBM's heavy-hexagon transmon chips.

Scalable Neural Decoder for Topological Surface Codes.

With the advent of noisy intermediate-scale quantum (NISQ) devices, practical quantum computing has seemingly come into reach. However, to go beyond proof-of-principle calculations, the current processing architectures will need to scale up to larger quantum circuits which will require fast and scalable algorithms for quantum error correction. Here, we present a neural network based decoder that, for a family of stabilizer codes subject to depolarizing noise and syndrome measurement errors, is scalable to tens of thousands of qubits (in contrast to other recent machine learning inspired decoders) and exhibits faster decoding times than the state-of-the-art union find decoder for a wide range of error rates (down to 1%). The key innovation is to autodecode error syndromes on small scales by shifting a preprocessing window over the underlying code, akin to a convolutional neural network in pattern recognition approaches. We show that such a preprocessing step allows to effectively reduce the error rate by up to 2 orders of magnitude in practical applications and, by detecting correlation effects, shifts the actual error threshold up to fifteen percent higher than the threshold of conventional error correction algorithms such as union find or minimum weight perfect matching, even in the presence of measurement errors. An in situ implementation of such a machine learning-assisted quantum error correction will be a decisive step to push the entanglement frontier beyond the NISQ horizon.

Identification of Non-Fermi Liquid Physics in a Quantum Critical Metal via Quantum Loop Topography.

Non-Fermi liquid physics is ubiquitous in strongly correlated metals, manifesting itself in anomalous transport properties, such as a T-linear resistivity in experiments. However, its theoretical understanding in terms of microscopic models is lacking, despite decades of conceptual work and attempted numerical simulations. Here we demonstrate that a combination of sign-problem-free quantum Monte Carlo sampling and quantum loop topography, a physics-inspired machine-learning approach, can map out the emergence of non-Fermi liquid physics in the vicinity of a quantum critical point (QCP) with little prior knowledge. Using only three parameter points for training the underlying neural network, we are able to robustly identify a stable non-Fermi liquid regime tracing the fans of metallic QCPs at the onset of both spin-density wave and nematic order. In particular, we establish for the first time that a spin-density wave QCP commands a wide fan of non-Fermi liquid region that funnels into the quantum critical point. Our study thereby provides an important proof-of-principle example that new physics can be detected via unbiased machine-learning approaches.

In memory of Achim Trebst (1929-2017): a pioneer of photosynthesis research.

The life and work of Achim Trebst (1929-2017) was dedicated to photosynthesis, involving a wide span of seminal contributions which cumulated in more than five decades of active research: Major topics include the separation of light and dark phases in photosynthesis, the elucidation of photosynthesis by the use of inhibitors, the identification of the three-dimensional structure of photosystem II and its degradation, and an explanation of singlet oxygen formation. For this tribute, which has been initiated by Govindjee, twenty-two personal tributes by former coworkers, scientific friends, and his family have been compiled and combined with an introduction tracing the different stages of Achim Trebst's scientific life.

Quantum Spin Liquids in Frustrated Spin-1 Diamond Antiferromagnets.

Motivated by the recent synthesis of the spin-1 A-site spinel NiRh_{2}O_{4}, we investigate the classical to quantum crossover of a frustrated J_{1}-J_{2} Heisenberg model on the diamond lattice upon varying the spin length S. Applying a recently developed pseudospin functional renormalization group approach for arbitrary spin-S magnets, we find that systems with S≥3/2 reside in the classical regime, where the low-temperature physics is dominated by the formation of coplanar spirals and a thermal (order-by-disorder) transition. For smaller local moments S=1 or S=1/2, we find that the system evades a thermal ordering transition and forms a quantum spiral spin liquid where the fluctuations are restricted to characteristic momentum-space surfaces. For the tetragonal phase of NiRh_{2}O_{4}, a modified J_{1}-J_{2}^{-}-J_{2}^{⊥} exchange model is found to favor a conventionally ordered Néel state (for arbitrary spin S), even in the presence of a strong local single-ion spin anisotropy, and it requires additional sources of frustration to explain the experimentally observed absence of a thermal ordering transition.

Competing Orders in a Nearly Antiferromagnetic Metal.

We study the onset of spin-density wave order in itinerant electron systems via a two-dimensional lattice model amenable to numerically exact, sign-problem-free determinantal quantum Monte Carlo simulations. The finite-temperature phase diagram of the model reveals a dome-shaped d-wave superconducting phase near the magnetic quantum phase transition. Above the critical superconducting temperature, an extended fluctuation regime manifests itself in the opening of a gap in the electronic density of states and an enhanced diamagnetic response. While charge density wave fluctuations are moderately enhanced in the proximity of the magnetic quantum phase transition, they remain short ranged. The striking similarity of our results to the phenomenology of many unconventional superconductors points a way to a microscopic understanding of such strongly coupled systems in a controlled manner.

Spin-Peierls Instability of Three-Dimensional Spin Liquids with Majorana Fermi Surfaces.

Three-dimensional (3D) variants of the Kitaev model can harbor gapless spin liquids with a Majorana Fermi surface on certain tricoordinated lattice structures such as the recently introduced hyperoctagon lattice. Here, we investigate Fermi surface instabilities arising from additional spin exchange terms (such as a Heisenberg coupling) which introduce interactions between the emergent Majorana fermion degrees of freedom. We show that independent of the sign and structure of the interactions, the Majorana surface is always unstable. Generically, the system spontaneously doubles its unit cell at exponentially small temperatures and forms a spin liquid with line nodes. Depending on the microscopics, further symmetries of the system can be broken at this transition. These spin-Peierls instabilities of a 3D spin liquid are closely related to BCS instabilities of fermions.

Transmon platform for quantum computing challenged by chaotic fluctuations.

From the perspective of many-body physics, the transmon qubit architectures currently developed for quantum computing are systems of coupled nonlinear quantum resonators. A certain amount of intentional frequency detuning ('disorder') is crucially required to protect individual qubit states against the destabilizing effects of nonlinear resonator coupling. Here we investigate the stability of this variant of a many-body localized phase for system parameters relevant to current quantum processors developed by the IBM, Delft, and Google consortia, considering the cases of natural or engineered disorder. Applying three independent diagnostics of localization theory - a Kullback-Leibler analysis of spectral statistics, statistics of many-body wave functions (inverse participation ratios), and a Walsh transform of the many-body spectrum - we find that some of these computing platforms are dangerously close to a phase of uncontrollable chaotic fluctuations.

Realization of nearly dispersionless bands with strong orbital anisotropy from destructive interference in twisted bilayer MoS2.

Recently, the twist angle between adjacent sheets of stacked van der Waals materials emerged as a new knob to engineer correlated states of matter in two-dimensional heterostructures in a controlled manner, giving rise to emergent phenomena such as superconductivity or correlated insulating states. Here, we use an ab initio based approach to characterize the electronic properties of twisted bilayer MoS2. We report that, in marked contrast to twisted bilayer graphene, slightly hole-doped MoS2 realizes a strongly asymmetric px-py Hubbard model on the honeycomb lattice, with two almost entirely dispersionless bands emerging due to destructive interference. The origin of these dispersionless bands, is similar to that of the flat bands in the prototypical Lieb or Kagome lattices and co-exists with the general band flattening at small twist angle due to the moiré interference. We study the collective behavior of twisted bilayer MoS2 in the presence of interactions, and characterize an array of different magnetic and orbitally-ordered correlated phases, which may be susceptible to quantum fluctuations giving rise to exotic, purely quantum, states of matter.

Orbital ordering in e(g) orbital systems: ground states and thermodynamics of the 120° model.

Orbital degrees of freedom shape many of the properties of a wide class of Mott insulating, transition metal oxides with partially filled 3d shells. Here we study orbital ordering transitions in systems where a single electron occupies the e(g) orbital doublet and the spatially highly anisotropic orbital interactions can be captured by an orbital-only model, often called the 120° model. Our analysis of both the classical and quantum limits of this model in an extended parameter space shows that the 120° model is in close proximity to several T=0 phase transitions and various competing ordered phases. We characterize the orbital order of these nearby phases and their associated thermal phase transitions by extensive numerical simulations and perturbative arguments.

Emergence of a field-driven U(1) spin liquid in the Kitaev honeycomb model.

In the field of quantum magnetism, the exactly solvable Kitaev honeycomb model serves as a paradigm for the fractionalization of spin degrees of freedom and the formation of [Formula: see text] quantum spin liquids. An intense experimental search has led to the discovery of a number of spin-orbit entangled Mott insulators that realize its characteristic bond-directional interactions and, in the presence of magnetic fields, exhibit no indications of long-range order. Here, we map out the complete phase diagram of the Kitaev model in tilted magnetic fields and report the emergence of a distinct gapless quantum spin liquid at intermediate field strengths. Analyzing a number of static, dynamical, and finite temperature quantities using numerical exact diagonalization techniques, we find strong evidence that this phase exhibits gapless fermions coupled to a massless U(1) gauge field. We discuss its stability in the presence of perturbations that naturally arise in spin-orbit entangled candidate materials.

Machine learning quantum phases of matter beyond the fermion sign problem.

State-of-the-art machine learning techniques promise to become a powerful tool in statistical mechanics via their capacity to distinguish different phases of matter in an automated way. Here we demonstrate that convolutional neural networks (CNN) can be optimized for quantum many-fermion systems such that they correctly identify and locate quantum phase transitions in such systems. Using auxiliary-field quantum Monte Carlo (QMC) simulations to sample the many-fermion system, we show that the Green's function holds sufficient information to allow for the distinction of different fermionic phases via a CNN. We demonstrate that this QMC + machine learning approach works even for systems exhibiting a severe fermion sign problem where conventional approaches to extract information from the Green's function, e.g. in the form of equal-time correlation functions, fail.

Numerical stabilization of entanglement computation in auxiliary-field quantum Monte Carlo simulations of interacting many-fermion systems.

In the absence of a fermion sign problem, auxiliary-field (or determinantal) quantum Monte Carlo (DQMC) approaches have long been the numerical method of choice for unbiased, large-scale simulations of interacting many-fermion systems. More recently, the conceptual scope of this approach has been expanded by introducing ingenious schemes to compute entanglement entropies within its framework. On a practical level, these approaches, however, suffer from a variety of numerical instabilities that have largely impeded their applicability. Here we report on a number of algorithmic advances to overcome many of these numerical instabilities and significantly improve the calculation of entanglement measures in the zero-temperature projective DQMC approach, ultimately allowing us to reach similar system sizes as for the computation of conventional observables. We demonstrate the applicability of this improved DQMC approach by providing an entanglement perspective on the quantum phase transition from a magnetically ordered Mott insulator to a band insulator in the bilayer square lattice Hubbard model at half filling.

Overcoming the slowing down of flat-histogram Monte Carlo simulations: cluster updates and optimized broad-histogram ensembles.

We study the performance of Monte Carlo simulations that sample a broad histogram in energy by determining the mean first-passage time to span the entire energy space of d-dimensional ferromagnetic Ising/Potts models. We first show that flat-histogram Monte Carlo methods with single-spin flip updates such as the Wang-Landau algorithm or the multicanonical method perform suboptimally in comparison to an unbiased Markovian random walk in energy space. For the d = 1, 2, 3 Ising model, the mean first-passage time tau scales with the number of spins N = L(d) as tau proportional N2L(z). The exponent z is found to decrease as the dimensionality d is increased. In the mean-field limit of infinite dimensions we find that z vanishes up to logarithmic corrections. We then demonstrate how the slowdown characterized by z > 0 for finite d can be overcome by two complementary approaches--cluster dynamics in connection with Wang-Landau sampling and the recently developed ensemble optimization technique. Both approaches are found to improve the random walk in energy space so that tau proportional N2 up to logarithmic corrections for the d = 1, 2 Ising model.

Optimizing the ensemble for equilibration in broad-histogram Monte Carlo simulations.

We present an adaptive algorithm which optimizes the statistical-mechanical ensemble in a generalized broad-histogram Monte Carlo simulation to maximize the system's rate of round trips in total energy. The scaling of the mean round-trip time from the ground state to the maximum entropy state for this local-update method is found to be O ( [N ln N](2) ) for both the ferromagnetic and the fully frustrated two-dimensional Ising model with N spins. Our algorithm thereby substantially outperforms flat-histogram methods such as the Wang-Landau algorithm.

Collective states of interacting anyons, edge states, and the nucleation of topological liquids.

Quantum mechanical systems, whose degrees of freedom are so-called su(2)k anyons, form a bridge between ordinary SU(2) quantum magnets (of arbitrary spin-S) and systems of interacting non-Abelian anyons. Anyonic spin-1/2 chains exhibit a topological protection mechanism that stabilizes their gapless ground states and which vanishes only in the limit (k-->infinity) of the ordinary spin-1/2 Heisenberg chain. For anyonic spin-1 chains the phase diagram closely mirrors the one of the biquadratic SU(2) spin-1 chain. Our results describe, at the same time, nucleation of different 2D topological quantum fluids within a "parent" non-Abelian quantum Hall state, arising from a macroscopic occupation with localized, interacting anyons. The edge states between the "nucleated" and the parent liquids are neutral, and correspond precisely to the gapless modes of the anyonic chains.

Local interactions and non-abelian quantum loop gases.

Two-dimensional quantum loop gases are elementary examples of topological ground states with Abelian or non-Abelian anyonic excitations. While Abelian loop gases appear as ground states of local, gapped Hamiltonians such as the toric code, we show that gapped non-Abelian loop gases require nonlocal interactions (or nontrivial inner products). Perturbing a local, gapless Hamiltonian with an anticipated "non-Abelian" ground-state wave function immediately drives the system into the Abelian phase, as can be seen by measuring the Hausdorff dimension of loops. Local quantum critical behavior is found in a loop gas in which all equal-time correlations of local operators decay exponentially.

Collective states of interacting Fibonacci anyons.

We show that chains of interacting Fibonacci anyons can support a wide variety of collective ground states ranging from extended critical, gapless phases to gapped phases with ground-state degeneracy and quasiparticle excitations. In particular, we generalize the Majumdar-Ghosh Hamiltonian to anyonic degrees of freedom by extending recently studied pairwise anyonic interactions to three-anyon exchanges. The energetic competition between two- and three-anyon interactions leads to a rich phase diagram that harbors multiple critical and gapped phases. For the critical phases and their higher symmetry end points we numerically establish descriptions in terms of two-dimensional conformal field theories. A topological symmetry protects the critical phases and determines the nature of gapped phases.