Search

Quantum spin models with a large number of interaction partners per spin are frequently used to describe modern many-body quantum optical systems like arrays of Rydberg atoms, atom-cavity systems or trapped ion crystals. For theoretical analysis the mean-field (MF) ansatz is routinely applied. However, besides special cases of all-to-all or strong long range interactions, the MF ansatz provides only approximate results. Here we present a systematic correction to MF theory based on diagrammatic perturbation theory for quantum spin correlators in thermal equilibrium. Our analytic results are universally applicable for any lattice geometry and spin-length S. We provide pre-computed and easy-to-use building blocks for Ising, Heisenberg and transverse field Ising models in the symmetry-unbroken regime. We showcase the quality and simplicity of the method by computing magnetic phase boundaries and excitations gaps. We also treat the Dicke-Ising model of ground-state superradiance where we show that corrections to the MF phase boundary vanish., Comment: minor changes, references added

In frustrated magnetism, the empirically found quantum-to-classical correspondence (QCC) matches the real-space static susceptibility pattern of a quantum spin-$1/2$ model with its classical counterpart computed at a certain elevated temperature. This puzzling relation was observed via bold line diagrammatic Monte Carlo simulations in dimensions two and three. The matching was within error bars and seemed valid down to the lowest accessible temperatures $T$ about an order of magnitude smaller than the exchange coupling $J$. Here we employ resummed spin diagrammatic perturbation theory to show analytically that the QCC breaks at fourth order in $J/T$ and provide the approximate mapping between classical and quantum temperatures. Our treatment further reveals that QCC is an indication of the surprising accuracy with which static correlators can be approximated by a simple renormalized mean-field form. We illustrate this for all models discussed in the context of QCC so far, including a recent example of the $S=1$ material $\mathrm{K}_2\mathrm{Ni}_2(\mathrm{SO}_4)_3$. The success of the mean-field form is traced back to partial diagrammatic cancellations., Comment: 6 pages, 5 figures, replaced figure 4, removed supplement, added end matter

We implement the temperature flow scheme first proposed by Honerkamp and Salmhofer in Phys.~Rev.~B 64, 184516 (2001) into the pseudo-Majorana functional renormalization group method for quantum spin systems. Since the renormalization group parameter in this approach is a physical quantity -- the temperature $T$ -- the numerical efficiency increases significantly compared to more conventional renormalization group parameters, especially when computing finite temperature phase diagrams. We first apply this method to determine the finite temperature phase diagram of the $J_1$-$J_2$ Heisenberg model on the simple cubic lattice where our findings support claims of a vanishingly small nonmagnetic phase around the high frustration point $J_2=0.25J_1$. Perhaps most importantly, we find the temperature flow scheme to be advantageous in detecting finite temperature phase transitions as, by construction, a phase transition is never encountered at an artificial, unphysical cutoff parameter. Finally, we apply the temperature flow scheme to the dipolar XXZ model on the square lattice where we find a rich phase diagram with a large non-magnetic regime down to the lowest accessible temperatures. Wherever a comparison with error-controlled (quantum) Monte Carlo methods is applicable, we find excellent quantitative agreement with less than $5\%$ deviation from the numerically exact results., Comment: 12 pages, 7 figures

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 of $m$-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., Comment: Review article for Reports on Progress in Physics. Comments are welcome! (43+7) pages, (16+4) figures, 3 tables

Motivated by a recent experiment on a square-lattice Rydberg atom array realizing a long-range dipolar XY model [Chen et al., Nature (2023)], we numerically study the model's equilibrium properties. We obtain the phase diagram, critical properties, entropies, variance of the magnetization, and site-resolved correlation functions. We consider both ferromagnetic and antiferromagnetic interactions and apply quantum Monte Carlo and pseudo-Majorana functional renormalization group techniques, generalizing the latter to a U(1) symmetric setting. Our simulations perform extensive thermometry for the first time in dipolar Rydberg atom arrays and establish conditions for adiabaticity and thermodynamic equilibrium. On the ferromagnetic side of the experiment, we determine the entropy per particle S/N~0.5, close to the one at the critical temperature, S_c/N = 0.585(15). The simulations suggest the presence of an out-of-equilibrium plateau at large distances in the correlation function, thus motivating future studies on the non-equilibrium dynamics of the system., Comment: 17 pages

In the field of quantum many-body physics, the spectral (or Lehmann) representation simplifies the calculation of Matsubara n-point correlation functions if the eigensystem of a Hamiltonian is known. It is expressed via a universal kernel function and a system- and correlator-specific product of matrix elements. Here we provide the kernel functions in full generality, for arbitrary n, arbitrary combinations of bosonic or fermionic operators and an arbitrary number of anomalous terms. As an application, we consider bosonic 3- and 4-point correlation functions for the fermionic Hubbard atom and a free spin of length S, respectively., Comment: close to published version

The pseudo-fermion representation for $S=1/2$ quantum spins introduces unphysical states in the Hilbert space which can be projected out using the Popov-Fedotov trick. However, state-of-the-art implementation of the functional renormalization group method for pseudo-fermions have so far omitted the Popov-Fedotov projection. Instead, restrictions to zero temperature were made and absence of unphysical contributions to the ground-state was assumed. We question this belief by exact diagonalization of several small-system counterexamples where unphysical states do contribute to the ground state. We then introduce Popov-Fedotov projection to pseudo-fermion functional renormalization, enabling finite temperature computations with only minor technical modifications to the method. At large and intermediate temperatures, our results are perturbatively controlled and we confirm their accuracy in benchmark calculations. At lower temperatures, the accuracy degrades due to truncation errors in the hierarchy of flow equations. Interestingly, these problems cannot be alleviated by switching to the parquet approximation. We introduce the spin projection as a method-intrinsic quality check. We also show that finite temperature magnetic ordering transitions can be studied via finite-size scaling., Comment: 14 pages, 8 figures; added acknowledgement

We employ a recently developed variant of the functional renormalization group method for spin systems, the so-called pseudo Majorana functional renormalization group, to investigate three-dimensional spin-1/2 Heisenberg models at finite temperatures. We study unfrustrated and frustrated Heisenberg systems on the simple cubic and pyrochlore lattices. Comparing our results with other quantum many-body techniques, we demonstrate a high quantitative accuracy of our method. Particularly, for the unfrustrated simple cubic lattice antiferromagnet ordering temperatures obtained from finite-size scaling of one-loop data deviate from error controlled quantum Monte Carlo results by $\sim5\%$ and we further confirm the established values for the critical exponent $\nu$ and the anomalous dimension $\eta$. As the PMFRG yields results in good agreement with QMC, but remains applicable when the system is frustrated, we next treat the pyrochlore Heisenberg antiferromagnet as a paradigmatic magnetically disordered system and find nearly perfect agreement of our two-loop static homogeneous susceptibility with other methods. We further investigate the broadening of pinch points in the spin structure factor as a result of quantum and thermal fluctuations and confirm a finite width in the extrapolated limit $T\rightarrow0$. While extensions towards higher loop orders $\ell$ seem to systematically improve our approach for magnetically disordered systems we also discuss subtleties when increasing $\ell$ in the presence of magnetic order. Overall, the pseudo Majorana functional renormalization group is established as a powerful many-body technique in quantum magnetism with a wealth of possible future applications., Comment: 15 (+9 appendix) pages, 11 figures

The combination of two-dimensional Dirac surface states with s-wave superconductivity is expected to generate localized topological Majorana zero modes in vortex cores. Putative experimental signatures of these modes have been reported for heterostructures of proximitized topological insulators, iron-based superconductors or certain transition metal dichalcogenides. Despite these efforts, the Majorana nature of the observed excitation is still under debate. We propose to identify the presence of Majorana vortex modes using a nonlocal transport measurement protocol originally employed for one-dimensional settings. In the case of an isolated subgap state, the protocol provides a spatial map of the ratio of local charge- and probability-density which offers a clear distinction between Majorana and ordinary fermionic modes. We show that these distinctive features survive in the experimentally relevant case of hybridizing vortex core modes., Comment: 12 pages, 6 Figures

The pseudofermion functional renormalization group (PFFRG) method has proven to be a powerful numerical approach to treat frustrated quantum spin systems. In its usual implementation, however, the complex fermionic representation of spin operators introduces unphysical Hilbert space sectors which renders an application at finite temperatures inaccurate. In this work, we formulate a general functional renormalization group approach based on Majorana fermions to overcome these difficulties. We, particularly, implement spin operators via an $SO(3)$ symmetric Majorana representation which does not introduce any unphysical states and, hence, remains applicable to quantum spin models at finite temperatures. We apply this scheme, dubbed pseudo Majorana functional renormalization group (PMFRG) method, to frustrated Heisenberg models on small spin clusters as well as square and triangular lattices. Computing the finite temperature behavior of spin correlations and thermodynamic quantities such as free energy and heat capacity, we find good agreement with exact diagonalization and the high-temperature series expansion down to moderate temperatures. We, particularly, observe a significantly enhanced accuracy of the PMFRG compared to the PFFRG at finite temperatures. More generally, we conclude that the development of functional renormalization group approaches with Majorana fermions considerably extends the scope of applicability of such methods., Comment: 16 pages, 6 figures. v2: Typos corrected. v3: substantially simplified flow equations. v4: Small quantitative corrections to figs. 3 and 4

Two-dimensional (2D) Dirac fermions are a central paradigm of modern condensed matter physics, describing low-energy excitations in graphene, in certain classes of superconductors, and on surfaces of 3D topological insulators. At zero energy E=0, Dirac fermions with mass m are band insulators, with the Chern number jumping by unity at m=0. This observation lead Ludwig et al [Phys. Rev. B 50, 7526 (1994)] to conjecture that the transition in 2D disordered Dirac fermions (DDF) and the integer quantum Hall transition (IQHT) are controlled by the same fixed point and possess the same universal critical properties. Given the far-reaching implications for the emerging field of the quantum anomalous Hall effect, modern condensed matter physics and our general understanding of disordered critical points, it is surprising that this conjecture has never been tested numerically. Here, we report the results of extensive numerics on the phase diagram and criticality of 2D-DDF in the unitary class. We find a critical line at m=0, with energy-dependent localization length exponent. At large energies, our results for the DDF are consistent with state-of-the-art numerical results \nu_IQH = 2.56 - 2.62 from models of the IQHT. At E=0 however, we obtain \nu_0 =2.30-2.36 incompatible with \nu_IQH. This result challenges conjectured relations between different models of the IQHT, and several interpretations are discussed., Comment: close to published version

We propose and investigate numerically a one-dimensional model which exhibits a non-Anderson disorder-driven transition. Such transitions have recently been attracting a great deal of attention in the context of Weyl semimetals, one-dimensional systems with long-range hopping and high-dimensional semiconductors. Our model hosts quasiparticles with the dispersion $\pm |k|^\alpha \mathrm{sign} k$ with $\alpha<1/2$ near two points (nodes) in momentum space and includes short-range-correlated random potential which allows for scattering between the nodes and near each node. In contrast with the previously studied models in dimensions $d<3$, the model considered here exhibits a critical scaling of the Thouless conductance which allows for {an accurate} determination of the critical properties of the non-Anderson transition, with a precision significantly exceeding the results obtained from the critical scaling of the density of states, usually simulated at such transitions. We find that in the limit of the vanishing parameter $\varepsilon=2\alpha-1$ the correlation-length exponent $\nu=2/(3|\varepsilon|)$ at the transition is inconsistent with the prediction $\nu_{RG}=1/|\varepsilon|$ of the perturbative renormalisation-group analysis. Our results allow for a numerical verification of the convergence of $\varepsilon$-expansions for non-Anderson disorder-driven transitions and, in general, interacting field theories near critical dimensions., Comment: 12 pages

We consider the disorder-induced correction to the minimal conductance of an anisotropic two-dimensional Dirac node or a three-dimensional Weyl node. An analytical expression is derived for the correction $\delta G$ to the conductance of a finite-size sample by an arbitrary potential, without taking the disorder average, in second-order perturbation theory. Considering a generic model of a short-range disorder potential, this result is used to compute the probability distribution $P(\delta G)$, which is compared to the numerically exact distribution obtained using the scattering matrix approach. We show that $P(\delta G)$ is Gaussian when the sample has a large width-to-length ratio, and study how the expectation value, the standard deviation, and the probability of finding $\delta G < 0$ depend on the anisotropy of the dispersion., Comment: 10+ pages, 6 figures

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., Comment: v2: Added conductance distributions, improved readability

The self-consistent Born approximation quantitatively fails to capture disorder effects in semimetals. We present an alternative, simple-to-use non-perturbative approach to calculate the disorder induced self-energy. It requires a sufficient broadening of the quasiparticle pole and the solution of a differential equation on the imaginary frequency axis. We demonstrate the performance of our method for various paradigmatic semimetal Hamiltonians and compare our results to exact numerical reference data. For intermediate and strong disorder, our approach yields quantitatively correct momentum resolved results. It is thus complementary to existing RG treatments of weak disorder in semimetals., Comment: 5+6 pages

We present a functional renormalization group approach to interacting topological Green function invariants with a focus on the nature of transitions. The method is applied to chiral symmetric fermion chains in the Mott limit that can be driven into a Haldane phase. We explicitly show that the transition to this phase is accompanied by a zero of the fermion Green function. Our results for the phase boundary are quantitatively benchmarked against DMRG data.

The functional renormalization group (FRG) has been used widely to investigate phase diagrams, in particular the one of the two-dimensional Hubbard model. So far, the study of one-dimensional models has not attracted as much attention. We use the FRG to investigate the phases of a one-dimensional spinless tight-binding chain with nearest and next-nearest neighbor interactions at half filling. The phase diagram of this model has already been established with other methods, and phase transitions from a metallic phase to ordered phases take place at intermediate to strong interactions. The model is thus well suited to analyze the potential and the limitations of the FRG in this regime of interactions. We employ flow equations that are exact up to second order in the interaction, which implies that we take into account the frequency dependence of the two-particle vertex as well as the feedback of the dynamic self-energy. For intermediate nearest neighbor interactions, our scheme captures the phase transition from a metallic phase to a charge density wave with alternating occupation. The critical interaction, at which this transition occurs, is underestimated due to our approximations. Similarly, for intermediate next-nearest neighbor interactions, we observe a transition to a charge density wave with occupation pattern ..00110011... We show that taking into account a feedback of the two-particle vertex in the flow equation is essential for the detection of those phases., Comment: 15 pages, 13 figures

Magnetotransport experiments on Weyl semimetals are essential for investigating the intriguing topological and low-energy properties of Weyl nodes. If the transport direction is perpendicular to the applied magnetic field, experiments have shown a large positive magnetoresistance. In this work, we present a theoretical scattering matrix approach to transversal magnetotransport in a Weyl node. Our numerical method confirms and goes beyond the existing perturbative analytical approach by treating disorder exactly. It is formulated in real space and is applicable to mesoscopic samples as well as in the bulk limit. In particular, we study the case of clean and strongly disordered samples., Comment: 10 pages, 4 figures

We devise a functional renormalization group treatment for a chain of interacting spinless fermions which is correct up to second order in the interaction strength. We treat both inhomogeneous systems in real-space as well as the translational invariant case in a k-space formalism. The strengths and shortcomings of the different schemes as well as technical details of their implementation are discussed. We use the method to study two proof-of-principle problems in the realm of Luttinger liquid physics, namely reflection at interfaces and power laws in the occupation number as a function of crystal momentum., Comment: 9 pages

Disorder effects are especially pronounced around nodal points in linearly dispersing bandstructures as present in graphene or Weyl semimetals. Despite the enormous experimental and numerical progress, even a simple quantity like the average density of states cannot be assessed quantitatively by analytical means. We demonstrate how this important problem can be solved employing the functional renormalization group method and, for the two dimensional case, demonstrate excellent agreement with reference data from numerical simulations based on tight-binding models. In three dimensions our analytic results also improve drastically on existing approaches., Comment: Erratum prepended (1 + 4 + 3 pages)

Although Lorentz invariance forbids the presence of a term that tilts the energy-momentum relation in the Weyl Hamiltonian, a tilted dispersion is not forbidden and, in fact, generic for condensed matter realizations of Weyl semimetals. We here investigate the combined effect of such a tilted Weyl dispersion and the presence of potential disorder. In particular, we address the influence of a tilt on the disorder-induced phase transition between a quasi-ballistic phase at weak disorder, in which the disorder is an irrelevant perturbation, and a diffusive phase at strong disorder. Our main result is that the presence of a tilt leads to a reduction of the critical disorder strength for this transition or, equivalently, that increasing the tilt at fixed disorder strength drives the system through the phase transition to the diffusive strong-disorder phase. Notably this obscures the tilt induced Lifshitz transition to an over tilted type-II Weyl phase at any finite disorder strength. Our results are supported by analytical calculations using the self-consistent Born approximation and numerical calculations of the density of states and of transport properties., Comment: 8 pages, 8 figures

We study Weyl semimetals in the presence of generic disorder, consisting of a random vector potential as well as a random scalar potential. We derive renormalization group flow equations to second order in the disorder strength. These flow equations predict a disorder-induced phase transition between a pseudo-ballistic weak-disorder phase and a diffusive strong-disorder phase for sufficiently strong random scalar potential or for a pure three-component random vector potential. We verify these predictions using a numerical study of the density of states near the Weyl point and of quantum transport properties at the Weyl point. In contrast, for a pure single-component random vector potential the diffusive strong-disorder phase is absent., Comment: published version with minor changes

Double Weyl nodes are topologically protected band crossing points which carry chiral charge $\pm2$. They are stabilized by $C_{4}$ point group symmetry and are predicted to occur in $\mathrm{SrSi_{2}}$ or $\mathrm{HgCr_{2}Se_{4}}$. We study their stability and physical properties in the presence of a disorder potential. We investigate the density of states and the quantum transport properties at the nodal point. We find that, in contrast to their counterparts with unit chiral charge, double Weyl nodes are unstable to any finite amount of disorder and give rise to a diffusive phase, in agreement with predictions of Goswami and Nevidomskyy [Phys. Rev. B 92, 214504 (2015)] and Bera, Sau, and Roy [Phys. Rev. B 93, 201302(R) (2016)]. However, for finite system sizes a crossover between pseudodiffusive and diffusive quantum transport can be observed., Comment: Supersedes arXiv:1512.07164

Strong topological insulators may have nonzero weak indices. The nonzero weak indices allow for the existence of topologically protected helical states along line defects of the lattice. If the lattice admits line defects that connect opposite surfaces of a slab of such a "weak-and-strong" topological insulator, these states effectively connect the surface states at opposite surfaces. Depending on the phases accumulated along the dislocation lines, this connection results in a suppression of in-plane transport and the opening of a spectral gap or in an enhanced density of states and an increased conductivity.

Double Weyl nodes are topologically protected band crossing points which carry chiral charge $\pm2$. They are stabilized by $C_4$ point group symmetry and are predicted to occur in $\mathrm{SrSi_2}$ or $\mathrm{HgCr_{2}Se_{4}}$. We study their fate in the presence of quenched disorder by numerical transport calculations and scaling arguments. We find that a double Weyl node is unstable in a disordered environment and splits into a pair of emergent simple Weyl nodes which carry equal chiral charge of $\pm1$. This pair of simple Weyl nodes is robust to disorder up to a certain critical disorder strength, in complete analogy to their `elementary' counterparts., Comment: This manuscript has been withdrawn due to an error in the numerics

Three-dimensional Dirac and Weyl semimetals exhibit a disorder-induced quantum phase transition between a semimetallic phase at weak disorder and a diffusive-metallic phase at strong disorder. Despite considerable effort, both numerically and analytically, the critical exponents $\nu$ and $z$ of this phase transition are not known precisely. Here we report a numerical calculation of the critical exponent $\nu=1.47\pm0.03$ using a minimal single-Weyl node model and a finite-size scaling analysis of conductance. Our high-precision numerical value for $\nu$ is incompatible with previous numerical studies on tight-binding models and with one- and two-loop calculations in an $\epsilon$-expansion scheme. We further obtain $z=1.49\pm0.02$ from the scaling of the conductivity with chemical potential.

We calculate conductance and noise for quantum transport at the nodal point for arbitrarily tilted and anisotropic Dirac or Weyl cones. Tilted and anisotropic dispersions are generic in absence of certain discrete symmetries, such as particle-hole and lattice point group symmetries. Whereas anisotropy affects the conductance g, but leaves the Fano factor F (the ratio of shot noise power and current) unchanged, a tilt affects both g and F. Since F is a universal number in many other situations, this finding is remarkable. We apply our general considerations to specific lattice models of strained graphene and a pyrochlore Weyl semi-metal., Comment: Published. 6 pages, 5 figures

Weyl semimetals are paradigmatic topological gapless phases in three dimensions. We here address the effect of disorder on charge transport in Weyl semimetals. For a single Weyl node with energy at the degeneracy point and without interactions, theory predicts the existence of a critical disorder strength beyond which the density of states takes on a nonzero value. Predictions for the conductivity are divergent, however. In this work, we present a numerical study of transport properties for a disordered Weyl cone at zero energy. For weak disorder our results are consistent with a renormalization group flow towards an attractive pseudoballistic fixed point with zero conductivity and a scale-independent conductance; for stronger disorder diffusive behavior is reached. We identify the Fano factor as a signature that discriminates between these two regimes.

The role of disorder in the field of three-dimensional time reversal invariant topological insulators has become an active field of research recently. However, the computation of Z2 invariants for large, disordered systems still poses a considerable challenge. In this paper we apply and extend a recently proposed method based on the scattering matrix approach, which allows the study of large systems at reasonable computational effort with few-channel leads. By computing the Z2 invariant directly for the disordered topological Anderson insulator, we unambiguously identify the topological nature of this phase without resorting to its connection with the clean case. We are able to efficiently compute the Z2 phase diagram in the mass-disorder plane. The topological phase boundaries are found to be well described by the self consistent Born approximation, both for vanishing and finite chemical potential., Comment: 10 pages