Search

Spontaneous symmetry breaking can persist at all temperatures in certain biconical $\mathrm{O}(N)\times \mathbb{Z}_2$ vector models when the underlying field theories are ultraviolet complete. So far, the existence of such theories has been established in fractional dimensions for local but nonunitary models or in 2+1 dimensions but for nonlocal models. Here, we study local models at zero and finite temperature directly in 2+1 dimensions employing functional methods. At zero temperature, we establish that our approach describes the quantum critical behaviour with high accuracy for all $N\geq 2$. We then exhibit the mechanism of discrete symmetry breaking from $\mathrm{O}(N)\times \mathbb{Z}_2\to \mathrm{O}(N)$ for increasing temperature near the biconical critical point when $N$ is finite but large. We calculate the corresponding finite-temperature phase diagram and further show that the Hohenberg-Mermin-Wagner theorem is fully respected within this approach, i.e., symmetry breaking only occurs in the $\mathbb{Z}_2$ sector. Finally, we determine the critical $N$ above which this phenomenon can be observed to be $N_c \approx 15$., Comment: 6 + 5 pages, 3 + 2 figures, comments welcome

Interaction-induced charge orders with electronic origin occur as states of spontaneously broken symmetry in several materials platforms. An electronic mechanism for charge order requires an attractive component in the effective charge vertex. We put forward such a mechanism for the formation of unconventional charge density waves in a metal. These states result from the condensation of particle-hole pairs with finite wave vector and non-zero angular momentum and correspond to bond or loop current order on a lattice. The mechanism we describe can be viewed as Kohn Luttinger analysis in the particle-hole channel with finite transferred momentum. It incorporates one-loop spin and pairing correctionsn, which are then used as an input for a summation in the charge channel triggering an instability. We extend our analysis to a spin-fluctuation approach, where the effective charge interaction is dressed by the particle-hole ladder with exchanged momentum. We argue that this mechanism works for weakly-interacting metals with nested Fermi surface and a large number of fermion flavors. We apply the Kohn-Luttinger-like approach to square- and triangular-lattice Hubbard models with SU($N_f$) flavour symmetry and show that it leads to different types of $p$-wave charge density waves. We also study effects beyond weak coupling at and away from Van Hove filling in terms of a phenomenological model with additional exchange interaction. In the vicinity of Van Hove filling, we obtain $d$-wave charge density waves with wave vectors determined by nesting as leading instabilities. In addition, we find another charge density wave with wave vector $K/4$ on the triangular lattice on both sides of Van Hove filling. We demonstrate that this $K/4$ instability can win the competition against pairing for $N_f=4$ via an unbiased functional renormalisation group calculation., Comment: 12+5 pages, 8+1 figures

The understanding of phenomena falling outside the Ginzburg-Landau paradigm of phase transitions represents a key challenge in condensed matter physics. A famous class of examples is constituted by the putative deconfined quantum critical points between two symmetry-broken phases in layered quantum magnets, such as pressurised SrCu$_2$(BO$_3$)$_2$. Experiments find a weak first-order transition, which simulations of relevant microscopic models can reproduce. The origin of this behaviour has been a matter of considerable debate for several years. In this work, we demonstrate that the nature of the deconfined quantum critical point can be best understood in terms of a novel dynamical mechanism, termed Nordic walking. Nordic walking denotes a renormalisation group flow arising from a beta function that is flat over a range of couplings. This gives rise to a logarithmic flow that is faster than the well-known walking behaviour, associated with the annihilation and complexification of fixed points, but still significantly slower than the generic running of couplings. The Nordic-walking mechanism can thus explain weak first-order transitions, but may also play a role in high-energy physics, where it could solve hierarchy problems. We analyse the Wess-Zumino-Witten field theory pertinent to deconfined quantum critical points with a topological term in 2+1 dimensions. To this end, we construct an advanced functional renormalisation group approach based on higher-order regulators. We thereby calculate the beta function directly in 2+1 dimensions and provide evidence for Nordic walking., Comment: 12 + 11 pages (typeset with sn-jnl.cls), 4 + 1 figures; comments welcome

Evidence for relativistic quantum criticality of antiferromagnetism and superconductivity in two-dimensional Dirac fermion systems has been found in large-scale quantum Monte Carlo simulations. However, the corresponding ($2+1$)-dimensional Gross--Neveu--Yukawa field theory with $N_f=2$ four-component Dirac fermions coupled to two triplets of order parameters does not exhibit a renormalization group fixed point at one-loop order. Instead, the theory only features a critical point for a large or very small fractional number of fermion flavors $N_f$, which disappears for a broad range of flavor numbers around the physical case, $N_f=2$, due to fixed-point annihilation. This raises the question on how to explain the observed scaling collapse in the quantum Monte Carlo data. Here, we extend previous renormalization-group analyses by studying a generalized model at two-loop order in $4-\epsilon$ spacetime dimensions. We determine the $\epsilon$ correction to the upper and lower critical flavor numbers for the fixed-point annihilation and find that they both go towards the physical case $N_f=2$. However, this only happens very slowly, such that an extrapolation to $\epsilon=1$ still suggests the absence of criticality in $2+1$ dimensions. Thereby, we consolidate the finding that the continuum field theory does not feature a stable renormalization-group fixed point and no true quantum criticality would be expected for the considered system. We briefly discuss a possible reconciliation in terms of a complex conformal field theory. Further, we also explore the fixed-point structure in an enlarged theory space and identify a candidate stable fixed-point solution., Comment: 13 pages, 6 figures, 2 tables

Two-dimensional van der Waals materials have recently been established experimentally as a highly-tunable condensed matter platform, facilitating the controlled manipulation of band structures and interactions. In several of these materials, Dirac cones are present in the low-energy regime near the Fermi level. Thus, fermionic excitations emerging in these materials close to the Dirac cones have a linear dispersion relation near the Fermi surface as massless relativistic Dirac fermions. Here, we study low-energy fermionic excitations of such Dirac materials in the presence of a mass gap that may be generated by symmetry breaking. Introducing a dynamical Fermi velocity and/or time-dependent mass gap for the Dirac quasiparticles, we exhibit the emergence of an analog of cosmological fermion pair production in terms of observables such as the expected occupation number or two-point correlation functions. We find that it is necessary and sufficient for quasiparticle production that only the ratio between the mass gap and the Fermi velocity is time-dependent. In this way, we establish that highly-tunable Dirac materials can serve as analog models for cosmological spacetime geometries, in particular, for Friedmann-Lema\^itre-Robertson-Walker expanding cosmologies. We briefly discuss possibilities for experimental realization, Comment: 15 pages, 4 figures

Over the past two decades, the enigma of the deconfined quantum critical point (DQCP) has attracted broad attention across the condensed matter, quantum field theory, and high-energy physics communities, as it is expected to offer a new paradigm in theory, experiment, and numerical simulations that goes beyond the Landau-Ginzburg-Wilson framework of symmetry breaking and phase transitions. However, the nature of DQCP has been controversial. For instance, in the square-lattice spin-1/2 $J$-$Q$ model, believed to realize the DQCP between N\'eel and valence bond solid states, conflicting results, such as first-order versus continuous transition, and critical exponents incompatible with conformal bootstrap bounds, have been reported. The enigma of DQCP is exemplified in its anomalous logarithmic subleading contribution in its entanglement entropy (EE), which was discussed in recent studies. In the current work, we demonstrate that similar anomalous logarithmic behavior persists in a class of models analogous to the DQCP. We systematically study the quantum EE of square-lattice SU($N$) DQCP spin models. Based on large-scale quantum Monte Carlo computation of the EE, we show that for a series of $N$ smaller than a critical value, the anomalous logarithmic behavior always exists in the EE, which implies that the previously determined DQCPs in these models do not belong to conformal fixed points. In contrast, when $N\ge N_c$ with a finite $N_c$ that we evaluate to lie between $7$ and $8$, the DQCPs are consistent with conformal fixed points that can be understood within the Abelian Higgs field theory with $N$ complex components., Comment: The revised version focuses on the anomalous logarithmic correction to the EE arising from the smooth boundary, rather than corners. And the critical $N_c$ is determined based on the anomalous EE signal from smooth boundary

The discovery of correlated insulating states in moir\'e heterostructures has renewed the interest in strongly-coupled electron systems where spin and valley (or layer) degrees of freedom are intertwined. In the strong-coupling limit, such systems can be effectively described by SU(4) spin-valley models akin to Kugel-Khomskii models long studied in the context of spin-orbit coupled materials. However, typical moir\'e heterostructures also exhibit interactions that break the SU(4) symmetry down to SU(2)${}_{\mathrm{spin}}\otimes$U(1)${}_{\mathrm{valley}}$. Here we investigate the impact of such symmetry-breaking couplings on the magnetic phase diagram for triangular superlattices considering a filling of two electrons (or holes) per moir\'e unit cell. We explore a broad regime of couplings -- including XXZ anisotropies, Dzyaloshinskii-Moriya exchange and on-site Hund's couplings -- using semi-classical Monte Carlo simulations. We find a multitude of classically ordered phases, including (anti-)ferromagnetic, incommensurate, and stripe order, manifesting in different sectors of the spin-valley model's parameter space. Zooming in on the regimes where quantum fluctuations are likely to have an effect, we employ pseudo-fermion functional renormalization group (pf-FRG) calculations to resolve quantum disordered ground states such as spin-valley liquids, which we indeed find for certain parameter regimes. As a concrete example, we discuss the case of trilayer graphene aligned with hexagonal boron nitride using material-specific parameters., Comment: 20 pages, 16 figures

We study quantum multicritical behavior in a (2+1)-dimensional Gross-Neveu-Yukawa field theory with eight-component Dirac fermions coupled to two triplets of order parameters that act as Dirac masses, and transform as $(1,0) + (0,1)$ representation under the SO(4)$\simeq$SO(3)$\times$SO(3) symmetry group. This field theory is relevant to spin-1/2 fermions on honeycomb or $\pi$-flux lattices, for example, near the transition points between an $s$-wave superconductor and a charge-density wave, on one side, and N\'eel order, on the other. Two triplets of such order parameters always allow for a common pair of two other order parameters that would complete them to the maximal set of compatible (anticommuting) orders of five. We first derive a unitary transformation in the Nambu (particle-hole) space which maps any two such triplets, possibly containing some superconducting orders, onto purely insulating order parameters. This allows one to consider a universal SO(4) Gross-Neveu-Yukawa description of the multicriticality without any Nambu doubling. We then proceed to derive the renormalization-group flow of the coupling constants at one-loop order in $4-\epsilon$ space-time dimensions, allowing also a more general set of order parameters transforming under SO($n_a$)$\times$SO($n_b$). While for $n_a=n_b > 2 $ in the bosonic sector and with fermions decoupled there is a stable fixed point of the flow, the Yukawa coupling to fermions quickly leads to its elimination by a generic fixed-point collision in the relevant range of fermion flavor numbers $N_f$. This suggests the replacement of the critical behavior by a runaway flow in the physical case $n_a=n_b=3$. The structure of the RG flow at $n_a\neq n_b$ is also discussed, and some non-perturbative arguments in favor of the stability of the decoupled critical point when $n_a=3$ and $n_b=1$ in $D=2+1$ are provided., Comment: 13 pages, 4 figures

We study competing orders of spinless fermions in the triangular-lattice Hubbard model with nearest-neighbor interaction. We calculate the effective, momentum-resolved two-particle vertex in an unbiased way in terms of the functional renormalization group method and compare two different schemes for the momentum discretization, one based on dividing the Fermi surface into patches and one based on a channel decomposition. We study attractive and repulsive nearest-neighbor interaction and find a competition of pairing and charge instabilities. In the attractive case, a Pomeranchuk instability occurs at Van Hove filling and $f$-wave and $p$-wave pairing emerge when the filling is reduced. In the repulsive case, we obtain a charge density wave at Van Hove filling and extended $p$-wave pairing with reduced filling. The $p$-wave pairing solution is doubly degenerate and can realize chiral $p+ip$ superconductivity with different Chern numbers in the ground state. We discuss implications for strongly correlated spin-orbit coupled hexagonal electron systems such as moir\'e heterostructures.

Evidence for correlated insulating and superconducting phases around regions of high density of states was reported in the strongly spin-orbit coupled van-der Waals material twisted tungsten diselenide (tWSe$_2$). We investigate their origin and interplay by using a functional renormalization group approach that allows to describe superconducting and spin/charge instabilities in an unbiased way. We map out the phase diagram as function of filling and perpendicular electric field, and find that the moir\'e Hubbard model for tWSe$_2$ features mixed-parity superconducting order parameters with $s/f$-wave and topological $d/p$-wave symmetry next to (incommensurate) density wave states. Our work systematically characterizes competing interaction-driven phases in tWSe$_2$ beyond mean-field approximations and provides guidance for experimental measurements by outlining the fingerprint of correlated states in interacting susceptibilities., Comment: 7 pages, 3 figures, supplemental material

A simple yet paradigmatic model for the interplay of strong electronic correlations and geometric frustration is the triangular lattice Hubbard model. Recently it was proposed that moir\'e structures of transition metal dichalcogenides can be used to simulate extended versions that include non-local density-density interactions. We study competing instabilities of interacting electrons in such an extended Hubbard model on the triangular lattice near a filling where the density of states has a Van Hove singularity. We employ a truncated-unity functional renormalization group approach to investigate two cases: a paradigmatic minimally extended Hubbard model and a specific model with parameters that are applicable to hetero-bilayers of transition metal dichalcogenides. We unravel rich phase diagrams, including tendencies to spin-density-wave order and unconventional pairing, which can give rise to topological superconductivity. We classify the symmetry of the superconducting instabilities according to their irreducible representations and show that higher lattice harmonics are dominant when the nearest-neighbor interaction is sizable indicating pair formation between second-nearest neighbors. The phenomenological consequences can be enhanced spin and thermal quantum Hall responses in a topological superconductor., Comment: 11+7 pages, 11+2 figures

Experimental demonstrations of tunable correlation effects in magic-angle twisted bilayer graphene have put two-dimensional moir\'e quantum materials at the forefront of condensed-matter research. Other twisted few-layer graphitic structures, boron-nitride, and homo- or hetero-stacks of transition metal dichalcogenides (TMDs) have further enriched the opportunities for analysis and utilization of correlations in these systems. Recently, within the latter material class, strong spin-orbit coupling or excitonic physics were experimentally explored. The observation of a Mott insulating state and other fascinating collective phenomena such as generalized Wigner crystals, stripe phases and quantum anomalous Hall insulators confirmed the relevance of many-body interactions, and demonstrated the importance of their extended range. Since the interaction, its range, and the filling can be tuned experimentally by twist angle, substrate engineering and gating, we here explore Fermi surface instabilities and resulting phases of matter of hetero-bilayer TMDs. Using an unbiased renormalization group approach, we establish in particular that hetero-bilayer TMDs are unique platforms to realize topological superconductivity with winding number $|\mathcal{N}|=4$. We show that this state reflects in pronounced experimental signatures, such as distinct quantum Hall features., Comment: 8 pages, 3 figures

Moir\'e materials provide a highly tunable environment for the realization of band structures with engineered physical properties. Specifically, moir\'e structures with Fermi surface flat bands - a synthetic environment for the realization of correlated phases - have moir\'e unit cells containing thousands of atoms and tantalizingly complex bands structures. In this paper we show that statistical principles go a long way in explaining universal physical properties of these systems. Our approach builds on three conceptual elements: the presence of quantum chaos caused by the effective irregularity of the atomic configurations on short length scales, Anderson localization in momentum space, and the presence of approximate crystalline symmetries. Which of these principles dominates depends on material parameters such as the extension of the Fermi surface or the strength of the moir\'e lattice potential. The phenomenological consequences of this competition are predictions for the characteristic group velocity of moir\'e bands, a primary indicator for their average flatness. In addition to these generic features, we identify structures outside the statistical context, notably almost flat bands close to the extrema of the unperturbed spectra, and the celebrated zero energy `magic angle' flat bands, where the latter require exceptionally fine tuned material parameters., Comment: 15 pages, 10 figures

Two-dimensional spin-orbital magnets with strong exchange frustration have recently been predicted to facilitate the realization of a quantum critical point in the Gross-Neveu-SO(3) universality class. In contrast to previously known Gross-Neveu-type universality classes, this quantum critical point separates a Dirac semimetal and a long-range-ordered phase, in which the fermion spectrum is only partially gapped out. Here, we characterize the quantum critical behavior of the Gross-Neveu-SO(3) universality class by employing three complementary field-theoretical techniques beyond their leading orders. We compute the correlation-length exponent $\nu$, the order-parameter anomalous dimension $\eta_\phi$, and the fermion anomalous dimension $\eta_\psi$ using a three-loop $\epsilon$ expansion around the upper critical space-time dimension of four, a second-order large-$N$ expansion (with the fermion anomalous dimension obtained even at the third order), as well as a functional renormalization group approach in the improved local potential approximation. For the physically relevant case of $N=3$ flavors of two-component Dirac fermions in 2+1 space-time dimensions, we obtain the estimates $1/\nu = 1.03(15)$, $\eta_\phi = 0.42(7)$, and $\eta_\psi = 0.180(10)$ from averaging over the results of the different techniques, with the displayed uncertainty representing the degree of consistency among the three methods., Comment: 18 pages, 7 figures, 3 tables, 1 ancillary file "GNSO3-exponents.m". v2: additional explanations, minor corrections

The question whether global symmetries can be realized in quantum-gravity-matter-systems has far-reaching phenomenological consequences. Here, we collect evidence that within an asymptotically safe context, discrete global symmetries of the form $\mathbb{Z}_n$, $n>4$, cannot be realized in a near-perturbative regime. In contrast, an effective-field-theory approach to quantum gravity might feature such symmetries, providing a mechanism to generate mass hierarchies in the infrared without the need for additional fine-tuning., Comment: 27 pages, 10 figures. v2: matches published version

Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. Here we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive exchange interaction, we find a new type of spin Pomeranchuk order, in which the spin order parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case., Comment: 11+10 pages, 9+4 figures

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. 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.

The coupling between fermionic matter and gauge fields plays a fundamental role in our understanding of nature, while at the same time posing a challenging problem for theoretical modeling. In this situation, controlled information can be gained by combining different complementary approaches. Here, we study a confinement transition in a system of $N_f$ flavors of interacting Dirac fermions charged under a U(1) gauge field in 2+1 dimensions. Using Quantum Monte Carlo simulations, we investigate a lattice model that exhibits a continuous transition at zero temperature between a gapless deconfined phase, described by three-dimensional quantum electrodynamics, and a gapped confined phase, in which the system develops valence-bond-solid order. We argue that the quantum critical point is in the universality class of the QED$_3$-Gross-Neveu-XY model. We study this field theory within a $1/N_f$ expansion in fixed dimension as well as a renormalization group analysis in $4-\epsilon$ space-time dimensions. The consistency between numerical and analytical results is revealed from large to intermediate flavor number., Comment: 10 pages, 9 figures; v2: additional data and explanations, corrected erroneous interpretation of dimer scaling dimension

The quantum phase diagram and critical behavior of two-dimensional Dirac fermions coupled to two compatible order-parameter fields with $O(N_1)\oplus O(N_2)$ symmetry is investigated. Recent numerical studies of such systems have reported evidence for non-Landau-Ginzburg-Wilson transitions and emergent $O(N_1+N_2)$ symmetry between the two ordered states, which has been interpreted within a scenario of deconfined quantum criticality in (2+1)-dimensional Dirac materials. Here, we provide two theoretical approaches to refine the phase diagrams of such systems. In the immediate vicinity of the multicritical point between the ordered phases and the semimetallic phase, we employ a non-perturbative field-theoretical analysis based on the functional renormalization group. For the particular case of $N_1=3$, $N_2=1$, we perform a large-scale quantum Monte Carlo analysis of the strong-coupling region, where both orders meet. Our findings support the robust emergence of enhanced symmetry at the multicritical point and suggest the transition between the two ordered phases to take place via a sequence of continuous transitions. In particular, we find that intermediate regimes of coexistence are present in the phase diagram for all values of $N_1$ and $N_2$., Comment: 4+9 pages, 3+4 figures, 1+0 tables

Twisting moir\'e heterostructures to the flatband regime allows for the formation of strongly correlated quantum states, since the dramatic reduction of the bandwidth can cause the residual electronic interactions to set the principal energy scale. An effective description for such correlated moir\'e heterostructures, derived in the strong-coupling limit at integer filling, generically leads to spin-valley Heisenberg models. Here we explore the emergence and stability of spin liquid behavior in an SU(2)$^{\mathrm{spin}}\otimes$SU(2)$^{\mathrm{valley}}$ Heisenberg model upon inclusion of Hund's-induced and longer-ranged exchange couplings, employing a pseudofermion functional renormalization group approach. We consider two lattice geometries, triangular and honeycomb (relevant to different moir\'e heterostructures), and find, for both cases, an extended parameter regime surrounding the SU(4) symmetric point where no long-range order occurs, indicating a stable realm of quantum spin liquid behavior. For large Hund's coupling, we identify the adjacent magnetic orders, with both antiferromagnetic and ferromagnetic ground states emerging in the separate spin and valley degrees of freedom. For both lattice geometries the inclusion of longer-ranged exchange couplings is found to have both stabilizing and destabilizing effects on the spin liquid regime depending on the sign of the additional couplings., Comment: 5 pages, 4 figures plus appendix (4 pages, 5 figures)

The abelian Higgs model is the textbook example for the superconducting transition and the Anderson-Higgs mechanism, and has become pivotal in the description of deconfined quantum criticality. We study the abelian Higgs model with $n$ complex scalar fields at unprecedented four-loop order in the $4-\epsilon$ expansion and find that the annihilation of the critical and bicritical points occurs at a critical number of $n_c \approx 182.95\left(1 - 1.752\epsilon + 0.798 \epsilon^2 + 0.362\epsilon^3\right) + \mathcal{O}\left(\epsilon^4\right)\nonumber$. Consequently, below $n_c$, the transition turns from second to first order. Resummation of the series to extract the result in three-dimensions provides strong evidence for a critical $n_c(d=3)$ which is significantly below the leading-order value, but the estimates for $n_c$ are widely spread. Conjecturing the topology of the renormalization group flow between two and four dimensions, we obtain a smooth interpolation function for $n_c(d)$ and find $n_c(3)\approx 12.2\pm 3.9$ as our best estimate in three dimensions. Finally, we discuss Miransky scaling occurring below $n_c$ and comment on implications for weakly first-order behavior of deconfined quantum transitions. We predict an emergent hierarchy of length scales between deconfined quantum transitions corresponding to different $n$., Comment: 15 pages, 6 figures, 3 tables

The pseudofermion functional renormalization group (pf-FRG) has been put forward as a semi-analytical scheme that, for a given microscopic spin model, allows to discriminate whether the low-temperature states exhibit magnetic ordering or a tendency towards the formation of quantum spin liquids. However, the precise nature of the putative spin liquid ground state has remained hard to infer from the original (single-site) pf-FRG scheme. Here we introduce a cluster pf-FRG approach, which allows for a more stringent connection between a microscopic spin model and its low-temperature spin liquid ground states. In particular, it allows to calculate spatially structured fermion bilinear expectation values on spatial clusters, which are formed by splitting the original lattice into several sublattices, thereby allowing for the positive identification of a family of bilinear spin liquid states. As an application of this cluster pf-FRG approach, we consider the $J_1$-$J_2$ SU($N$)-Heisenberg model on a square lattice, which is a paradigmatic example for a frustrated quantum magnet exhibiting quantum spin liquid behavior for intermediate coupling strengths. In the well-established large-$N$ limit of this model, we show that our approach correctly captures the emergence of the $\pi$-flux spin liquid state at low temperatures. For small $N$, where the precise nature of the ground state remains controversial, we focus on the widely studied case of $N=2$, for which we determine the low-temperature phase diagram near the strongly-frustrated regime after implementing the fermion number constraint by the flowing Popov-Fedotov method. Our results suggest that the $J_1$-$J_2$-Heisenberg model does not support the formation of a fermion bilinear spin liquid state., Comment: 13 pages, 9 figures

We study the quantum many-body instabilities of interacting electrons with SU(2)$\times$SU(2) symmetry in spin and orbital degrees of freedom on the triangular lattice near van-Hove filling. Our work is motivated by effective models for the flat bands in hexagonal moir\'e heterostructures like twisted bilayer boron nitride and trilayer graphene-boron nitride systems. We consider an extended Hubbard model including onsite Hubbard and Hund's couplings, as well as nearest-neighbor exchange interactions and analyze the different ordering tendencies with the help of an unbiased functional renormalization group approach. We find three classes of instabilities controlled by the filling and bare interactions. For a nested Fermi surface at van-Hove filling, Hund-like couplings induce a weak instability towards spin or orbital density wave phases. An SU(4) exchange interaction moves the system towards a Chern insulator, which is robust with respect to perturbations from Hund-like interactions or deviations from perfect nesting. Further, in an extended range of fillings and interactions, we find topological $d\pm id$ and (spin-singlet)-(orbital-singlet) $f$-wave superconductivity., Comment: 5+4 pages, 4+1 figures

We study collective modes near the quantum critical point of a pair-density-wave (PDW) superconductor in 2+1 dimensions. The fate of gaps of various collective modes is investigated by functional renormalization. For incommensurate PDW superconductors, we show that the gapless Leggett mode, protected by the emergent $U(1)$ symmetry, can induce an exponentially small Higgs mass compared to the superconducting gap. Further, for commensurate PDW superconductors, we find an emergent mass hierarchy in the collective modes, i.e. the masses of Leggett boson, Higgs boson, and the superconducting gap can differ by several magnitudes in the infrared. This may shed light to a mechanism underlying the hierarchy problem in the Standard Model of particle physics., Comment: 8 pages, 5 figures

The QED$_3$-Gross-Neveu model is a (2+1)-dimensional U(1) gauge theory involving Dirac fermions and a critical real scalar field. This theory has recently been argued to represent a dual description of the deconfined quantum critical point between Neel and valence bond solid orders in frustrated quantum magnets. We study the critical behavior of the QED$_3$-Gross-Neveu model by means of an epsilon expansion around the upper critical space-time dimension of $D_c^+=4$ up to the three-loop order. Estimates for critical exponents in 2+1 dimensions are obtained by evaluating the different Pade approximants of their series expansion in epsilon. We find that these estimates, within the spread of the Pade approximants, satisfy a nontrivial scaling relation which follows from the emergent SO(5) symmetry implied by the duality conjecture. We also construct explicit evidence for the equivalence between the QED$_3$-Gross-Neveu model and a corresponding critical four-fermion gauge theory that was previously studied within the 1/N expansion in space-time dimensions 2

Gapless Dirac fermions appear as quasiparticle excitations in various condensed-matter systems. They feature quantum critical points with critical behavior in the 2+1 dimensional Gross-Neveu universality class. The precise determination of their critical exponents defines a prime benchmark for complementary theoretical approaches, such as lattice simulations, the renormalization group and the conformal bootstrap. Despite promising recent developments in each of these methods, however, no satisfactory consensus on the fermionic critical exponents has been achieved, so far. Here, we perform a comprehensive analysis of the Ising Gross-Neveu universality classes based on the recently achieved four-loop perturbative calculations. We combine the perturbative series in $4-\epsilon$ spacetime dimensions with the one for the purely fermionic Gross-Neveu model in $2+\epsilon$ dimensions by employing polynomial interpolation as well as two-sided Pad\'e approximants. Further, we provide predictions for the critical exponents exploring various resummation techniques following the strategies developed for the three-dimensional scalar $O(n)$ universality classes. We give an exhaustive appraisal of the current situation of Gross-Neveu universality by comparison to other methods. For large enough number of spinor components $N\geq 8$ as well as for the case of emergent supersymmetry $N=1$, we find our renormalization group estimates to be in excellent agreement with the conformal bootstrap, building a strong case for the validity of these values. For intermediate $N$ as well as in comparison with recent Monte Carlo results, deviations are found and critically discussed., Comment: 21 pages, 7 figures, 6 tables

Spinless fermions on the honeycomb lattice with repulsive nearest-neighbor interactions are known to harbour a quantum critical point at half-filling, with critical behaviour in the Gross-Neveu (chiral Ising) universality class. The critical interaction strength separates a weak-coupling semimetallic regime from a commensurate charge-density-wave phase. The phase diagram of this basic model of correlated fermions on the honeycomb lattice beyond half-filling is, however, less well established. Here, we perform an analysis of its many-body instabilities using the functional renormalization group method with a basic Fermi surface patching scheme, which allows us to treat instabilities in competing channels on equal footing also away from half-filling. Between half-filling and the van-Hove filling, the free Fermi surface is hole-like and we again find a charge-density wave instability to be dominant at large interactions. Moreover, its characteristics are those of the half-filled case. Directly at the van-Hove filling the nesting property of the free Fermi surface stabilizes a dimerized bond-order phase. At lower filling the free Fermi surface becomes electron-like and a superconducting instability with $f$-wave symmetry is found to emerge from the interplay of intra-unitcell repulsion and collective fluctuations in the proximity to the charge-density wave instability. We estimate the extent of the various phases and extract the corresponding order parameters from the effective low-energy Hamiltonians., Comment: 11 pages, 11 figures

The quantum phase transition to a $\mathbb{Z}_3$-ordered Kekul\'e valence bond solid in two-dimensional Dirac semimetals is governed by a fermion-induced quantum critical point, which renders the putatively discontinuous transition continuous. We study the resulting universal critical behavior in terms of a functional RG approach, which gives access to the scaling behavior on the symmetry-broken side of the phase transition, for general dimension and number of Dirac fermions. In particular, we investigate the emergence of the fermion-induced quantum critical point for space-time dimensions $2

We study the quantum multicritical point in a (2+1)-dimensional Dirac system between the semimetallic phase and two ordered phases that are characterized by anticommuting mass terms with $O(N_1)$ and $O(N_2)$ symmetry, respectively. Using $\epsilon$ expansion around the upper critical space-time dimension of four, we demonstrate the existence of a stable renormalization-group fixed point, enabling a direct and continuous transition between the two ordered phases directly at the multicritical point. This point is found to be characterized by an emergent $O(N_1+N_2)$ symmetry for arbitrary values of $N_1$ and $N_2$ and fermion flavor numbers $N_f$, as long as the corresponding representation of the Clifford algebra exists. Small $O(N)$-breaking perturbations near the chiral $O(N)$ fixed point are therefore irrelevant. This result can be traced back to the presence of gapless Dirac degrees of freedom at criticality, and it is in clear contrast to the purely bosonic $O(N)$ fixed point, which is stable only when $N < 3$. As a by-product, we obtain predictions for the critical behavior of the chiral $O(N)$ universality classes for arbitrary $N$ and fermion flavor number $N_f$. Implications for critical Weyl and Dirac systems in 3+1 dimensions are also briefly discussed., Comment: 5+2 pages, 1 figure, 1 table

In frustrated magnetism, making a stringent connection between microscopic spin models and macroscopic properties of spin liquids remains an important challenge. A recent step towards this goal has been the development of the pseudofermion functional renormalization group approach (pf-FRG) which, building on a fermionic parton construction, enables the numerical detection of the onset of spin liquid states as temperature is lowered. In this work, focusing on the SU(N) Heisenberg model at large N, we extend this approach in a way that allows us to directly enter the low-temperature spin liquid phase, and to probe its character. Our approach proceeds in momentum space, making it possible to keep the truncation minimalistic, while also avoiding the bias introduced by an explicit decoupling of the fermionic parton interactions into a given channel. We benchmark our findings against exact mean-field results in the large-N limit, and show that even without prior knowledge the pf-FRG approach identifies the correct mean-field decoupling channel. On a technical level, we introduce an alternative finite temperature regularization scheme that is necessitated to access the spin liquid ordered phase. In a companion paper arXiv:1711.02182 we present a different set of modifications of the pf-FRG scheme that allow us to study SU(N) Heisenberg models (using a real-space RG approach) for arbitrary values of N, albeit only up to the phase transition towards spin liquid physics., Comment: 12 pages, 5 figures

The link between a modified Higgs self-coupling and the strong first-order phase transition necessary for baryogenesis is well explored for polynomial extensions of the Higgs potential. We broaden this argument beyond leading polynomial expansions of the Higgs potential to higher polynomial terms and to non-polynomial Higgs potentials. For our quantitative analysis we resort to the functional renormalization group, which allows us to evolve the full Higgs potential to higher scales and finite temperature. In all cases we find that a strong first-order phase transition manifests itself in an enhancement of the Higgs self-coupling by at least 50%, implying that such modified Higgs potentials should be accessible at the LHC., Comment: 24 pages, 11 figures; v2: minor changes, refs updated, matches journal version

We study the chiral Ising, the chiral XY and the chiral Heisenberg models at four-loop order with the perturbative renormalization group in $4-\epsilon$ dimensions and compute critical exponents for the Gross-Neveu-Yukawa fixed points to order $\mathcal{O}(\epsilon^4)$. Further, we provide Pad\'e estimates for the correlation length exponent, the boson and fermion anomalous dimension as well as the leading correction to scaling exponent in 2+1 dimensions. We also confirm the emergence of supersymmetric field theories at four loops for the chiral Ising and the chiral XY models with $N=1/4$ and $N=1/2$ fermions, respectively. Furthermore, applications of our results relevant to various quantum transitions in the context of Dirac and Weyl semimetals are discussed, including interaction-induced transitions in graphene and surface states of topological insulators., Comment: 18 pages, 1 figure, 3 supplemental files attached containing the critical exponents for general number of fermion flavors for each of the models

We study the quantum many-body ground states of electrons on the half-filled honeycomb lattice with short- and long-ranged density-density interactions as a model for graphene. To this end, we employ the recently developed truncated-unity functional renormalization group (TU-fRG) approach which allows for a high resolution of the interaction vertex' wavevector dependence. We connect to previous lattice quantum Monte Carlo (QMC) results which predict a stabilization of the semimetallic phase for realistic \emph{ab initio} interaction parameters and confirm that the application of a finite biaxial strain can induce a quantum phase transition towards an ordered ground state. In contrast to lattice QMC simulations, the TU-fRG is not limited in the choice of tight-binding and interaction parameters to avoid the occurrence of a sign problem. Therefore, we also investigate a range of parameters relevant to the realistic graphene material which are not accessible by numerically exact methods. Although a plethora of charge density waves arise under medium-range interactions, we find the antiferromagnetic spin-density wave to be the prevailing instability for long-range interactions. We further explore the impact of an extended tight-binding Hamiltonian with second-nearest neighbor hopping and a finite chemical potential for a more accurate description of the band structure of graphene's $p_z$ electrons., Comment: 17 pages, 12 figures

We establish a scenario where fluctuations of new degrees of freedom at a quantum phase transition change the nature of a transition beyond the standard Landau-Ginzburg paradigm. To this end we study the quantum phase transition of gapless Dirac fermions coupled to a $\mathbb{Z}_3$ symmetric order parameter within a Gross-Neveu-Yukawa model in 2+1 dimensions, appropriate for the Kekul\'e transition in honeycomb lattice materials. For this model the standard Landau-Ginzburg approach suggests a first order transition due to the symmetry-allowed cubic terms in the action. At zero temperature, however, quantum fluctuations of the massless Dirac fermions have to be included. We show that they reduce the putative first-order character of the transition and can even render it continuous, depending on the number of Dirac fermions $N_f$. A non-perturbative functional renormalization group approach is employed to investigate the phase transition for a wide range of fermion numbers. For the first time we obtain the critical $N_f$, where the nature of the transition changes. Furthermore, it is shown that for large $N_f$ the change from the first to second order of the transition as a function of dimension occurs exactly in the physical 2+1 dimensions. We compute the critical exponents and predict sizable corrections to scaling for $N_f =2$., Comment: 12+5 pages, 5 figures

Dirac and Weyl fermions appear as quasi-particle excitations in many different condensed-matter systems. They display various quantum transitions which represent unconventional universality classes related to the variants of the Gross-Neveu model. In this work we study the bosonized version of the standard Gross-Neveu model -- the Gross-Neveu-Yukawa theory -- at three-loop order, and compute critical exponents in $4-\epsilon$ dimensions for general number of fermion flavors. Our results fully encompass the previously known two-loop calculations, and agree with the known three-loop results in the purely bosonic limit of the theory. We also find the exponents to satisfy the emergent super-scaling relations in the limit of a single-component fermion, order by order up to three loops. Finally, we apply the computed series for the exponents and their Pad\'e approximants to several phase transitions of current interest: metal-insulator transitions of spin-1/2 and spinless fermions on the honeycomb lattice, emergent supersymmetric surface field theory in topological phases, as well as the disorder-induced quantum transition in Weyl semimetals. Comparison with the results of other analytical and numerical methods is discussed., Comment: 5+ pages

Over the last years, huge resources of biological and medical data have become available for research. This data offers great chances for machine learning applications in health care, e.g. for precision medicine, but is also challenging to analyze. Typical challenges include a large number of possibly correlated features and heterogeneity in the data. One flourishing field of biological research in which this is relevant is epigenetics. Here, especially large amounts of DNA methylation data have emerged. This epigenetic mark has been used to predict a donor's 'epigenetic age' and increased epigenetic aging has been linked to lifestyle and disease history. In this paper we propose an adaptive model which performs feature selection for each test sample individually based on the distribution of the input data. The method can be seen as partially blind domain adaptation. We apply the model to the problem of age prediction based on DNA methylation data from a variety of tissues, and compare it to a standard model, which does not take heterogeneity into account. The standard approach has particularly bad performance on one tissue type on which we show substantial improvement with our new adaptive approach even though no samples of that tissue were part of the training data., Comment: NIPS 2016 Workshop on Machine Learning for Health, Barcelona, Spain

We study the quantum phase transition of $U(1)$ - charged Dirac fermions Yukawa-coupled to the Kekul\'e valence bond solid order parameter with $Z_3$ symmetry of the honeycomb lattice. The symmetry allows for the presence of the term in the action which is cubic in the Kekul\'e order parameter, and which is expected to prevent the quantum phase transition in question from being continuous. The Gross-Neveu-Yukawa theory for the transition is investigated using a perturbative renormalization group and within the $\epsilon$ expansion close to four space-time dimensions. For a vanishing $U(1)$ charge we show that quantum fluctuations may render the phase transition continuous only sufficiently far away from 3+1 dimensions, where the validity of the conclusions based on the leading order $\epsilon$ expansion appear questionable. In the presence of a fluctuating gauge field, on the other hand, we find quantum critical behavior even at weak coupling to appear close to 3+1 dimensions, that is, within the domain of validity of the perturbation theory. We also determine the renormalization group scaling of the cubic coupling at higher-loop orders and for a large number of Dirac fermions for vanishing charge., Comment: 8 pages, 4 figures

Interacting fermions on the half-filled honeycomb lattice with short-range repulsions have been suggested to host a variety of interesting many-body ground states, e.g., a topological Mott insulator. A number of recent studies of the spinless case in terms of exact diagonalization, the infinite density matrix renormalization group and the functional renormalization group, however, indicate a suppression of the topological Mott insulating phase in the whole range of interaction parameters. Here, we complement the previous studies by investigating the quantum many-body instabilities of the physically relevant case of spin-1/2 fermions with onsite, nearest-neighbor and second-nearest-neighbor repulsion. To this end, we employ the multi-patch functional renormalization group for correlated fermions with refined momentum resolution observing the emergence of an antiferromagnetic spin-density wave and a charge-density wave for dominating onsite and nearest-neighbor repulsions, respectively. For dominating second-nearest neighbor interaction our results favor an ordering tendency towards a charge-modulated ground state over the topological Mott insulating state. The latter evades a stabilization as the leading instability by the additional onsite interaction., Comment: 9 pages, 6 figures

We explore O(N) models in dimensions $4

We study the phase diagram of two-color Quantum Chromodynamics at finite temperature and chemical potential. This is done within an effective low-energy description in terms of quarks, mesons and diquarks. Quantum, thermal and density fluctuations are taken into account with the functional renormalisation group approach. In particular, we establish the phenomenon of pre-condensation, affecting the location of the phase boundary to Bose-Einstein condensation. We also discuss the Silver Blaze property in the context of the functional renormalisation group., Comment: 11+9 pages

The Renormalisation Group is a versatile tool for the study of many systems where scale-dependent behaviour is important. Its functional formulation can be cast into the form of an exact flow equation for the scale-dependent effective action in the presence of an infrared regularisation. The functional RG flow for the scale-dependent effective action depends explicitly on the choice of regulator, while the physics does not. In this work, we systematically investigate three key aspects of how the regulator choice affects RG flows: (i) We study flow trajectories along closed loops in the space of action functionals varying both, the regulator scale and shape function. Such a flow does not vanish in the presence of truncations. Based on a definition of the length of an RG trajectory, we suggest a practical procedure for devising optimised regularisation schemes within a truncation. (ii) In systems with various field variables, a choice of relative cutoff scales is required. At the example of relativistic bosonic two-field models, we study the impact of this choice as well as its truncation dependence. We show that a crossover between different universality classes can be induced and conclude that the relative cutoff scale has to be chosen carefully for a reliable description of a physical system. (iii) Non-relativistic continuum models of coupled fermionic and bosonic fields exhibit also dependencies on relative cutoff scales and regulator shapes. At the example of the Fermi polaron problem in three spatial dimensions, we illustrate such dependencies and show how they can be interpreted in physical terms., Comment: 22 pages, 16 figures

We study the competition of spin- and charge-density waves and their quantum multicritical behavior for the semimetal-insulator transitions of low-dimensional Dirac fermions. Employing the effective Gross-Neveu-Yukawa theory with two order parameters as a model for graphene and a growing number of other two-dimensional Dirac materials allows us to describe the physics near the multicritical point at which the semimetallic and the spin- and charge-density-wave phases meet. With the help of a functional renormalization group approach, we are able to reveal a complex structure of fixed points, the stability properties of which decisively depend on the number of Dirac fermions $N_f$. We give estimates for the critical exponents and observe crucial quantitative corrections as compared to the previous first-order $\epsilon$ expansion. For small $N_f$, the universal behavior near the multicritical point is determined by the chiral Heisenberg universality class supplemented by a decoupled, purely bosonic, Ising sector. At large $N_f$, a novel fixed point with nontrivial couplings between all sectors becomes stable. At intermediate $N_f$, including the graphene case ($N_f = 2$) no stable and physically admissible fixed point exists. Graphene's phase diagram in the vicinity of the intersection between the semimetal, antiferromagnetic and staggered density phases should consequently be governed by a triple point exhibiting first-order transitions., Comment: 13+7 pages, 6 figures

We use the functional renormalization group and the $\epsilon$-expansion concertedly to explore multicritical universality classes for coupled $\bigoplus_i O(N_i)$ vector-field models in three Euclidean dimensions. Exploiting the complementary strengths of these two methods we show how to make progress in theories with large numbers of interactions, and a large number of possible symmetry-breaking patterns. For the three- and four-field models we find a new fixed point that arises from the mutual interaction between different field sectors, and we establish the absence of infrared-stable fixed point solutions for the regime of small $N_i$. Moreover, we explore these systems as toy models for theories that are both asymptotically safe and infrared complete. In particular, we show that these models exhibit complete renormalization group trajectories that begin and end at nontrivial fixed points., Comment: 10 pages, 6 figures; minor changes, as published in EPJ C

We investigate the quantum many-body instabilities of the extended Hubbard model for spinless fermions on the honeycomb lattice with repulsive nearest-neighbor and 2nd nearest-neighbor density-density interactions. Recent exact diagonalization and infinite density matrix renormalization group results suggest that a putative topological Mott insulator phase driven by the 2nd nearest-neighbor repulsion is suppressed, while other numerically exact approaches support the topological Mott insulator scenario. In the present work, we employ the functional renormalization group (fRG) for correlated fermionic systems. Our fRG results hint at a strong suppression of the scattering processes stabilizing the topological Mott insulator. From analyzing the effects of fermionic fluctuations, we obtain a phase diagram which is the result of the competition of various charge ordering instabilities., Comment: 9 pages, 8 figures

We study the multicritical behavior for the semimetal-insulator transitions on graphene's honeycomb lattice using the Gross-Neveu-Yukawa effective theory with two order parameters: the SO(3) (Heisenberg) order parameter describes the antiferromagnetic transition, and the $\mathbb{Z}_2$ (Ising) order parameter describes the transition to a staggered density state. Their coupling induces multicritical behavior which determines the structure of the phase diagram close to the multicritical point. Depending on the number of fermion flavors $N_f$ and working in the perturbative regime in vicinity of three (spatial) dimensions, we observe first order or continuous phase transitions at the multicritical point. For the graphene case of $N_f=2$ and within our low order approximation, the phase diagram displays a tetracritical structure., Comment: 5 pages, 3 figures

In view of the measured Higgs mass of 125 GeV, the perturbative renormalization group evolution of the Standard Model suggests that our Higgs vacuum might not be stable. We connect the usual perturbative approach and the functional renormalization group which allows for a straightforward inclusion of higher-dimensional operators in the presence of an ultraviolet cutoff. In the latter framework we study vacuum stability in the presence of higher-dimensional operators. We find that their presence can have a sizable influence on the maximum ultraviolet scale of the Standard Model and the existence of instabilities. Finally, we discuss how such operators can be generated in specific models and study the relation between the instability scale of the potential and the scale of new physics required to avoid instabilities., Comment: 27 pages, 12 figures

Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account on the ground state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite onsite repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonous finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure., Comment: 12 pages, 15 figures

We study models with three coupled vector fields characterized by $O(N_1)\oplus O(N_2) \oplus O(N_3)$ symmetry. Using the nonperturbative functional renormalization group, we derive $\beta$ functions for the couplings and anomalous dimensions in $d$ dimensions. Specializing to the case of three dimensions, we explore interacting fixed points that generalize the $O(N)$ Wilson-Fisher fixed point. We find a symmetry-enhanced isotropic fixed point, a large class of fixed points with partial symmetry enhancement, as well as partially and fully decoupled fixed point solutions. We discuss their stability properties for all values of $N_1, N_2$, and $N_3$, emphasizing important differences to the related two-field models. For small numbers of field components we find no stable fixed point solutions, and we argue that this can be attributed to the presence of a large class of possible (mixed) couplings in the three-field and multifield models. Furthermore, we contrast different mechanisms for stability interchange between fixed points in the case of the two- and three-field models, which generically proceed through fixed-point collisions., Comment: 14 pages, 15 figures; PRE version, minor changes

We use a functional renormalization group approach to study the instabilities due to electron-electron interactions in a bilayer honeycomb lattice model with AA stacking, as it might be relevant for layered graphene with this structure. Starting with a tight- binding description for the four $\pi$-bands, we integrate out the modes of the dispersion by successively lowering an infrared cutoff and determine the leading tendencies in the effective interactions. The antiferromagnetic spin-density wave is an expected instability for dominant local repulsion among the electrons, but for nonlocal interaction terms also other instabilities occur. We discuss the phase diagrams depending on the model parameters. We compare our results to single-layer graphene and the more common AB-stacked bilayer, both qualitatively and quantitatively., Comment: 9 pages, 3 figures, Annalen der Physik, online available (2014)

We study the impact of electron-phonon interactions on the many-body instabilities of electrons on the honeycomb lattice and their interplay with repulsive local and non-local Coulomb interactions at charge neutrality. To that end, we consider in-plane optical phonon modes with wavevectors close to the $\Gamma$ point as well as to the $K, -K$ points and calculate the effective phonon-mediated electron-electron interaction by integrating out the phonon modes. Ordering tendencies are studied by means of a momentum-resolved functional renormalization group approach allowing for an unbiased investigation of the appearing instabilities. In the case of an exclusive and supercritical phonon-mediated interaction, we find a Kekul\'e and a nematic bond ordering tendency being favored over the $s$-wave superconducting state. The competition between the different phonon-induced orderings clearly shows a repulsive interaction between phonons at small and large wavevector transfers. We further discuss the influence of phonon-mediated interactions on electronically-driven instabilities induced by onsite, nearest neighbor and next-to-nearest neighbor density-density interactions. We find an extension of the parameter regime of the spin density wave order going along with an increase of the critical scales where ordering occurs, and a suppression of competing orders., Comment: 9 pages, 5 figures