Functional Renormalization Research Articles

Relativistic BEC extracted from a complex FRG flow equation

Abstract Based on the functional renormalization group (FRG) under the local potential approximation, we analyze the Bose-Einstein condensation (BEC) in the relativistic complex scalar theory. This framework leads to a complex flow equation of the effective potential, even with the well-known Litim regulator. In order to evaluate the condensate from such a complex effective potential, we impose a condition between chemical potential and mass, analogously to those in the free theory or the mean field theory. We elucidate that for the strongly (weakly) coupled theory, the phase diagrams computed from the FRG are more (less) deviated from that under the mean field approximation. This result implies that quantum fluctuations strongly affect the nonperturbative formation of the BEC.

Magic angle of Sr2RuO4 : Optimizing correlation-driven superconductivity

Understanding of unconventional superconductivity is crucial for engineering materials with specific order parameters or elevated superconducting transition temperatures. However, for many materials, the pairing mechanism and symmetry of the order parameter remain unclear: reliable and efficient methods of predicting the order parameter and its response to tuning parameters are lacking. Here, we investigate the response of superconductivity in Sr2RuO4 to structural distortions via the random phase approximation (RPA) and functional renormalization group (FRG), starting from realistic models of the electronic structure. Our results suggest that RPA misses the interplay of competing fluctuation channels. FRG reproduces key experimental findings. We predict a magic octahedral rotation angle, maximizing the superconducting Tc and a dominant dx2−y2 pairing symmetry. To enable experimental verification, we provide calculations of the phase-referenced Bogoliubov Quasiparticle Interference imaging. Our work demonstrates a designer approach to tuning unconventional superconductivity with relevance and applicability for a wide range of quantum materials. Published by the American Physical Society 2024

Effective action and black hole solutions in asymptotically safe quantum gravity

We derive the quantum effective action and the respective quantum equations of motion from multi-graviton correlation functions in asymptotically safe quantum gravity. The fully momentum-dependent couplings of three- and four-graviton scatterings are computed within the functional renormalization group approach and the effective action is reconstructed from these vertices. The resulting quantum equations of motion are solved numerically for quantum black hole geometries. Importantly, the black hole solutions show signatures of quantum gravity outside the classical horizon, which manifest in the behavior of the temporal and radial components of the metric. Three different types of solutions with distinct causal structures are identified and the phase structure of the solution space is investigated. Published by the American Physical Society 2024

Asymptotic safe nonassociative quantum gravity with star R-flux products, Goroff–Sagnotti counter-terms, and geometric flows

Nonassociative modifications of general relativity, GR, defined by star products with R-flux deformations in string gravity consist an important subclass of modified gravity theories, MGTs. A longstanding criticism for elaborating quantum gravity, QG, argue that the asymptotic safety does not survive once certain perturbative terms (in general, nonassociative and noncommutative) are included in the projection space. The goal of this work is to prove that a generalized asymptotic safety scenario allows us to formulate physically viable nonassociative QG theories using effective models defined by generic off-diagonal solutions and nonlinear symmetries in nonassociative geometric flow and gravity theories. We elaborate on a new nonholonomic functional renormalization techniques with parametric renormalization group, RG, flow equations for effective actions supplemented by certain canonical two-loop counter-terms. The geometric constructions and quantum deformations are performed for nonassociative phase spaces modelled as R-flux deformed cotangent Lorentz bundles. Our results prove that theories involving nonassociative modifications of GR can be well defined both as classical nonassociative MGTs and QG models. Such theories are characterized by generalized G. Perelman thermodynamic variables which are computed for certain examples of nonassociative geometric and RG flows.

New Universality Class Describes Vicsek's Flocking Phase in Physical Dimensions.

The Vicsek simulation model of flocking together with its theoretical treatment by Toner and Tu in 1995 were two foundational cornerstones of active matter physics. However, despite the field's tremendous progress, the actual universality class (UC) governing the scaling behavior of Viscek's "flocking" phase remains elusive. Here, we use nonperturbative, functional renormalization group methods to analyze, numerically and analytically, a simplified version of the Toner-Tu model, and uncover a novel UC with scaling exponents that agree remarkably well with the values obtained in a recent simulation study by Mahault etal. [Phys. Rev. Lett. 123, 218001 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.218001], in both two and three spatial dimensions. We therefore believe that there is strong evidence that the UC uncovered here describes Vicsek's flocking phase.

Numerical accuracy of the derivative-expansion-based functional renormalization group

Abstract We investigate the precision of the numerical implementation of the functional renormalization group based on extracting the eigenvalues from the linearized renormalization group transformation. For this purpose, we implement the local potential approximation and O ( ∂ 2 ) orders of the derivative expansion for the three-dimensional O(N) models with N ∈ { 1 , 2 , 3 } . We identify several categories of numerical error and devise simple tests to track their magnitude as functions of numerical parameters. Our numerical schemes converge properly and are characterized by errors of several orders of magnitude smaller than the error bars of the derivative expansion for these models. We highlight situations in which our methods cease to converge, most often due to rounding errors. In particular, we observe an impaired convergence of the discretization scheme when the ρ ~ grid is cut off at the value ρ ~ Max smaller than 3.5 times the local potential minimum. The program performing the numerical calculations for this study is shared as an open-source library accessible for review and reuse.

Momentum-dependent quasiparticle properties of the Fermi polaron from the functional renormalization group

We study theoretically the lifetimes of attractive and repulsive Fermi polarons, as well as the molecule at finite momentum in three dimensions. To this end, we develop a technique that allows for the computation of Green's functions in the whole complex frequency plane using exact analytical continuation within the functional renormalization group. The improved numerical stability and reduced computational cost of this method yield access to previously inaccessible momentum-dependent quasiparticle properties of low-lying excited states. While conventional approaches like the non-self-consistent T-matrix approximation method cannot determine these lifetimes, we are able to find the momentum-dependent lifetime at different interaction strengths of both the attractive and repulsive polaron as well as the molecule. At weak coupling our results confirm predictions made from effective Fermi liquid theory regarding the decay of the attractive polaron, and we demonstrate that Fermi liquidlike behavior extends far into the strong-coupling regime where the attractive polaron and molecule exhibit a p4 momentum scaling in their decay widths. Our results offer an intriguing insight into the momentum-dependent quasiparticle properties of the Fermi polaron problem, which can be measured using techniques such as Raman transfer and Ramsey interferometry. Published by the American Physical Society 2024

Criticality of the O(N) universality via global solutions to nonperturbative fixed-point equations

Fixed-point equations in the functional renormalization group approach are integrated from large to vanishing field, where an asymptotic potential in the limit of large field is implemented as initial conditions. This approach allows us to obtain a global fixed-point potential with high numerical accuracy, that incorporates the correct asymptotic behavior in the limit of large field. Our calculated global potential is in good agreement with the Taylor expansion in the region of small field, and it also coincides with the Laurent expansion in the regime of large field. Laurent expansion of the potential in the limit of large field for general case, that the spatial dimension d is a continuous variable in the range 2≤d≤4\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2\\le d \\le 4$$\\end{document}, is obtained. Eigenfunctions and eigenvalues of perturbations near the Wilson–Fisher fixed point are computed with the method of eigenperturbations. Critical exponents for different values of d and N of the O(N) universality class are calculated.

Scalar spectral functions from the spectral functional renormalization group

We compute spectral functions in a scalar ϕ4-theory in three spacetime dimensions via the spectral functional renormalization group. This approach allows for the direct, manifestly Lorentz covariant computation of correlation functions in Minkowski spacetime, including a physical on-shell renormalization. We present numerical results for the spectral functions of the two- and four-point correlation functions for different values of the coupling parameter. These results agree very well with those obtained from another functional real-time approach, the spectral Dyson-Schwinger equation. Published by the American Physical Society 2024

No-go for the formation of heavy mass Primordial Black Holes in Single Field Inflation

We examine the possibility of Primordial Black Holes (PBHs) formation in single field models of inflation. Using the adiabatic or wave function renormalization scheme in the short range modes, we show that one-loop correction to the power spectrum is free from quadratic UV divergence. We consider a framework in which PBHs are produced during the transition from Slow Roll (SR) to Ultra Slow Roll (USR) followed by the end of inflation. We demonstrate that the renormalized power spectrum soften the contribution of the logarithmic IR divergence and severely restricts the possible mass range of produced PBHs in the said transition, namely, MPBH∼102g\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$M_\ extrm{PBH}\\sim 10^{2}\\,\ extrm{g}$$\\end{document}a^\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hat{\ extrm{a}}$$\\end{document}la a no-go theorem. In particular, we find that the produced PBHs are short lived (tPBHevap∼10-20s\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$t^\ extrm{evap}_\ extrm{PBH}\\sim 10^{-20}\\,\ extrm{s}$$\\end{document}) and the corresponding number of e-folds in the USR region is restricted to ΔNUSR≈2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Delta {\\mathcal {N}}_\ extrm{USR}\\approx 2$$\\end{document}.

Functional renormalisation group approach to the finite-temperature Bose polaron

The functional renormalisation group (FRG) approach is employed to study Bose polarons at finite temperatures in the regime of strong attractive bath-impurity interactions. Both two- and three-dimensional configurations are considered. The appearance of two polaron quasiparticle branches at finite temperatures is revealed, consistent with recent findings by other analytical techniques. Ground-state polaron energies are also reported for selected interactions and temperatures within the gas superfluid phase. The findings of this work present the FRG as a useful tool for studying finite-temperature polarons in quantum gases.

Improved memory truncation scheme for quasi-adiabatic propagator path integral via influence functional renormalization.

Accurately simulating non-Markovian quantum dynamics in system-bath coupled problems remains challenging. In this work, we present a novel memory truncation scheme for the iterative quasi-adiabatic propagator path integral (iQuAPI) method to improve accuracy. Conventional memory truncation in iQuAPI discards all influence functional beyond a certain time interval, which is not effective for problems with a long memory time. Our proposed scheme selectively retains the most significant parts of the influence functional using the density matrix renormalization group algorithm. We validate the effectiveness of our scheme through simulations of the spin-boson model across various parameter sets, demonstrating faster convergence and improved accuracy compared to the conventional scheme. Our findings suggest that the new memory truncation scheme significantly advances the capabilities of iQuAPI for problems with a long memory time.

Dynamics of K2Ni2(SO4)3 governed by proximity to a 3D spin liquid model

Quantum spin liquids (QSLs) have become a key area of research in magnetism due to their remarkable properties, such as long-range entanglement, fractional excitations, and topologically protected phenomena. Recently, the search for QSLs has expanded into the three-dimensional world, despite the suppression of quantum fluctuations due to high dimensionality. A new candidate material, K2Ni2(SO4)3, belongs to the langbeinite family and consists of two interconnected trillium lattices. Although magnetically ordered, it exhibits a highly dynamical and correlated state. In this work, we combine inelastic neutron scattering measurements with density functional theory (DFT), pseudo-fermion functional renormalization group (PFFRG), and classical Monte Carlo (cMC) calculations to study the magnetic properties of K2Ni2(SO4)3, revealing a high level of agreement between experiment and theory. We further reveal the origin of the dynamical state in K2Ni2(SO4)3 to be centred around a magnetic network composed of tetrahedra on a trillium lattice.

Self-organized time crystal in driven-dissipative quantum system

Continuous time crystals (CTCs) are characterized by sustained oscillations that break the time-translation symmetry. Since the ruling out of equilibrium CTCs by no-go theorems, the emergence of such dynamical phases has been observed in various driven-dissipative quantum platforms. The current understanding of CTCs is mainly based on mean-field theories, which fail to address the problem of whether the continuous time-translation symmetry can be broken in noisy, spatially extended systems absent in all-to-all couplings. Here, we propose a CTC realized in a quantum contact model through self-organized bistability. The CTCs stem from the interplay between collective dissipation induced by the first-order absorbing phase transitions and slow constant driving provided by an incoherent pump. The stability of such oscillatory phases in finite dimensions under the action of intrinsic quantum fluctuations is scrutinized by the functional renormalization group method and numerical simulations. Occurring at the edge of many-body synchronization, the CTC phase exhibits an inherent period and amplitude with a coherence time linearly diverging with system size, thus also constituting a boundary time crystal. Our results serve as a solid route towards self-protected CTCs in strongly interacting open systems. Published by the American Physical Society 2024

Nonlocal Effects in Asymptotically Safe Gravity

The asymptotically safe gravity is investigated in the framework of the functional renormalization group method. The low energy region of the model can account for the cosmological behavior, where it is assumed that the nonlocal effects play a crucial role. Using the Wegner–Houghton equation it is shown that the dynamically induced bilocal term modifies the infrared scaling of the model.

Functional renormalization group for fermions on a one dimensional lattice at arbitrary filling

Functional renormalization group for fermions on a one dimensional lattice at arbitrary filling

Local Solutions of RG Flow Equations from the Nash–Moser Theorem

We prove local existence of solutions of a functional renormalisation group equation for the effective action of an interacting quantum field theory, when a suitable local potential approximation is considered. To obtain this equation in a Lorentzian setting a quantum state for the theory is selected and a regulator consisting in a mass is added to the action. The flow equation for mass rescalings is then studied using the renown Nash–Moser theorem.

First order phase transition with the functional renormalization group method

The renormalization group method, more specifically the Wegner-Houghton equation, is used to find first order phase transitions in a simple scalar field theory with a polynomial potential. An improved definition of the running parameters allows us to explore the renormalization group flow down to the IR endpoint and to locate phase transitions. Beyond the expected first order transition further radiative correction generated first and second order transitions are found. The phase diagram is reviewed by a Monte-Carlo simulation of the lattice regulated version of the theory but the serious slow down of the convergence prevents us from obtaining conclusive results from the simulation. Published by the American Physical Society 2024

KeldyshQFT: A C++ codebase for real-frequency multiloop functional renormalization group and parquet computations of the single-impurity Anderson model.

We provide a detailed exposition of our computational framework designed for the accurate calculation of real-frequency dynamical correlation functions of the single-impurity Anderson model in the regime of weak to intermediate coupling. Using quantum field theory within the Keldysh formalism to directly access the self-energy and dynamical susceptibilities in real frequencies, as detailed in our recent publication [Ge et al., Phys. Rev. B 109, 115128 (2024)], the primary computational challenge is the full three-dimensional real-frequency dependence of the four-point vertex. Our codebase provides a fully MPI+OpenMP parallelized implementation of the functional renormalization group (fRG) and the self-consistent parquet equations within the parquet approximation. It leverages vectorization to handle the additional complexity imposed by the Keldysh formalism, using optimized data structures and highly performant integration routines. Going beyond the results shown in the previous publication, the code includes functionality to perform fRG calculations in the multiloop framework, up to arbitrary loop order, including self-consistent self-energy iterations. Moreover, implementations of various regulators, such as hybridization, interaction, frequency, and temperature, are supplied.

On the interplay between boundary conditions and the Lorentzian Wetterich equation

In the framework of the functional renormalization group and of the perturbative, algebraic approach to quantum field theory (pAQFT), in D’Angelo et al. [Ann. Henri Poinc. 25 (2024) 2295–2352] it has been derived a Lorentzian version of a flow equation à la Wetterich, which can be used to study nonlinear, quantum scalar field theories on a globally hyperbolic spacetime. In this work, we show that the realm of validity of this result can be extended to study interacting scalar field theories on globally hyperbolic manifolds with a timelike boundary. By considering the specific examples of half-Minkowski spacetime and of the Poincaré patch of Anti-de Sitter, we show that the form of the Lorentzian Wetterich equation is strongly dependent on the boundary conditions assigned to the underlying field theory. In addition, using a numerical approach, we are able to provide strong evidences that there is a qualitative and not only a quantitative difference in the associated flow and we highlight this feature by considering Dirichlet and Neumann boundary conditions on half-Minkowski spacetime.