tware architecture is developed that separates the physical model, numerical algorithms, and post-processing analysis into independent computational layers. The framework incorporates multi-start fixed-point searches, numerical Jacobian construction, eigenvalue and eigenvector analysis, residual-based convergence diagnostics, and automated extraction of universal critical exponents. The architecture is designed to accommodate higher-order truncations, alternative regulator choices, and continuation methods without requiring substantial modifications to the numerical infrastructure. The present implementation investigates scalar field models using quartic and higher-order polynomial truncations with a self-consistent treatment of the anomalous dimension within an LPA′-type framework. Validation procedures include Gaussian fixed-point recovery, truncation reduction tests, numerical convergence studies, and stability analysis. The resulting fixed-point structure and critical exponents are benchmarked against established Functional Renormalization Group calculations, providing a quantitative assessment of the computational implementation and its current level of agreement. Beyond the immediate numerical results, this work establishes a reproducible computational platform for systematic investigations of regulator dependence, truncation convergence, and dimensional continuation. Within the broader SDEL program, the framework serves as a foundation for exploring whether scale-dependent emergent dynamics can be formulated consistently alongside established renormalization-group methods. The paper therefore contributes both a validated computational methodology and a structured pathway for investigating future theoretical extensions while maintaining clear separation between established FRG techniques and proposed SDEL-specific developments.
tware architecture is developed that separates the physical model, numerical algorithms, and post-processing analysis into independent computational layers. The framework incorporates multi-start fixed-point searches, numerical Jacobian construction, eigenvalue and eigenvector analysis, residual-based convergence diagnostics, and automated extraction of universal critical exponents. The architecture is designed to accommodate higher-order truncations, alternative regulator choices, and continuation methods without requiring substantial modifications to the numerical infrastructure. The present implementation investigates scalar field models using quartic and higher-order polynomial truncations with a self-consistent treatment of the anomalous dimension within an LPA′-type framework. Validation procedures include Gaussian fixed-point recovery, truncation reduction tests, numerical convergence studies, and stability analysis. The resulting fixed-point structure and critical exponents are benchmarked against established Functional Renormalization Group calculations, providing a quantitative assessment of the computational implementation and its current level of agreement. Beyond the immediate numerical results, this work establishes a reproducible computational platform for systematic investigations of regulator dependence, truncation convergence, and dimensional continuation. Within the broader SDEL program, the framework serves as a foundation for exploring whether scale-dependent emergent dynamics can be formulated consistently alongside established renormalization-group methods. The paper therefore contributes both a validated computational methodology and a structured pathway for investigating future theoretical extensions while maintaining clear separation between established FRG techniques and proposed SDEL-specific developments.
Continuous phase transitions provide a natural setting for understanding universality. Near criticality, the correlation length diverges, microscopic details become increasingly irrelevant, and physical observables obey scaling laws characterized by universal critical exponents. Within the Functional Renormalization Group (FRG), this scale dependence is described through the Effective Average Action, which progressively incorporates fluctuations as the momentum scale is lowered. A first part of this thesis is developed within the Local Potential Approximation (LPA), where we study multicritical fixed points of scalar field theories. Their universal properties are compared with results obtained at higher orders in the derivative expansion. Particular attention is devoted to the relation between global fixed-point potentials and their local representations through finite polynomial truncations. We show that the convergence of such expansions is constrained by the analytic structure of the global solution, and that expanding around a suitable nonzero stationary point significantly improves the reconstruction of multicritical branches. Beyond the structure of individual fixed points, a central part of this thesis concerns RG irreversibility. In two dimensions, Zamolodchikov’s c-theorem describes the flow from the ultraviolet to the infrared through a monotonically decreasing c-function. Within the FRG framework, we develop a variational procedure to reconstruct RG trajectories between fixed points and integrate the corresponding c-function along them. The resulting quantity displays monotonic behaviour and becomes increasingly stable as the truncation order is raised, providing a manifestation of RG irreversibility within the approximation considered. Finally, close to the upper critical dimensions, the linearized FRG flow around the Gaussian fixed point is naturally organized in terms of Hermite polynomials. Using the epsilon-expansion, we derive a curvature-based candidate c-function and obtain an analytic explanation for the alternating curvature of successive multicritical fixed points and for their progressive accumulation around the Gaussian fixed point.
We introduce and develop Functional Dimensional Regularization (FDR), a novel functional renormalization group (RG) framework that extends the principles of dimensional regularization beyond the perturbative ε-expansion. The key insight is to define the FDR beta functions, valid in any continuous dimension d and characterized by threshold functions, as a sum over all critical dimensions of the corresponding scheme-independent beta functions computed in DR. We demonstrate that this formalism exhibits all the hallmarks of a fully-fledged functional RG. Within the local potential approximation, we perform a general fixed point analysis, recovering the ε-expansion for all multi-critical models and establishing a direct connection to functional scaling solutions, spike plots, and eigen-perturbation spectra. Extending the framework to the second order of the derivative expansion, we compute the critical exponents for the Ising universality class in d=3, providing estimates for the anomalous dimension η and RG spectrum that converge rapidly and show favorable agreement with state-of-the-art results. We present a detailed comparison with non-perturbative and proper-time RG approaches, highlighting structural similarities and key differences. In particular, we introduce a ``distillation'' procedure that allows one to systematically derive the FDR flow from any other one-loop-exact RG scheme. Our results establish FDR as a self-contained and competitive RG framework that combines the technical simplicity of dimensional regularization with the versatility of the functional approach.
We present a functional renormalization group (fRG) study of the two dimensional Hubbard model, performed with an algorithmic implementation which lifts some of the common approximations made in fRG calculations. In particular, in our fRG flow; (i) we take explicitly into account the momentum and the frequency dependence of the vertex functions; (ii) we include the feedback effect of the self-energy; (iii) we implement the recently introduced multiloop extension which allows us to sum up all the diagrams of the parquet approximation with their exact weight. Due to its iterative structure based on successive one-loop computations, the loop convergence of the fRG results can be obtained with an affordable numerical effort. In particular, focusing on the analysis of the physical response functions, we show that the results become independent from the chosen cutoff scheme and from the way the fRG susceptibilities are computed, i.e., either through flowing couplings to external fields, or through a “post-processing” contraction of the interaction vertex at the end of the flow. The presented substantial refinement of fRG-based computation schemes paves a promising route towards future quantitative fRG analyses of more challenging systems and/or parameter regimes.
Renormalization group methods are well-established tools for the (numerical) investigation of the low-energy properties of correlated quantum many-body systems, allowing us to capture their scale-dependent nature. The functional renormalization group (FRG) allows us to continuously evolve a microscopic model action to an effective low-energy action as a function of decreasing energy scales via an exact functional flow equation, which is then approximated by some truncation scheme to facilitate computation. Here, we report on our implementation of multiloop FRG, an extended truncation scheme recently developed for electronic FRG calculations, within the pseudofermion functional renormalization group (pf-FRG) framework for interacting quantum spin systems. We discuss in detail the conceptual intricacies of the flow equations generated by the multiloop truncation, as well as essential refinements to the integration scheme for the resulting integrodifferential equations. To benchmark our approach, we analyze antiferromagnetic Heisenberg models on the pyrochlore, simple cubic, and face-centered cubic lattice, discussing the convergence of physical observables for higher-loop calculations and comparing with existing results where available. Combined, these methodological refinements systematically improve the pf-FRG approach to one of the numerical tools of choice when exploring frustrated quantum magnetism in higher spatial dimensions.
The functional renormalization group (FRG) is a powerful tool that has facilitated the investigation of various strongly interacting theories, encompassing systems characterized by both bosonic and fermionic variables. Specifically, this project focuses on the analysis of theories with O(N) symmetry with bosonic scalar fields. These theories are of particular interest due to their broad applicability in real physical systems. For example, Z2 symmetry (the O(1) model) describes the well-known Ising universality class, which in turn models the liquid-gas transition. On the other hand, the O(2) model belongs to the universality class of the XY model, used to describe the transition from fluid to superfluid in 4 He, while the O(3) model, known as the Heisenberg model, describes the ferromagnetic transition in isotropic materials. Finally the O(0) model is related to the self-avoiding walk (SAW). To characterize these theories, various methods have been implemented within the framework of the renormalization group. Perturbation theory has been applied since the early days of the subject, specifically in the form of the ε-expansion, which has reached high loop order in recent years. Alternatively, a non-perturbative approach can be chosen, which is known as the non-perturbative renormalization group (NPRG). Within the NPRG, there are several approximation schemes, with the derivative expansion being particularly noteworthy. The ultimate goal of these efforts is to calculate the critical exponents, which define the given universality class. The standard non-perturbative method involves solving an equation that depends on a regulator or cut-off, which represents one of the main challenges of this approach. Although the regulators are designed so that the theory does not depend on them at the scales of interest, the results for the critical exponents can vary depending on the regulator used once approximations are made. In this project, we propose an approach that eliminates the need for an explicit mass cut-off significantly simplifying the calculation of critical exponents while mitigating the regulator dependence of the results. In particular, we show that it is possible to use dimensional regularization (DR) beyond the ε-expansion in the context of RG calculations of critical properties. Based on this we propose a new functional RG scheme called Functional Dimensional Regularization (FDR) and apply it to the O(N) model in three dimension, finding excellent agreement with state-of-the-art computations
The functional renormalization group (fRG) is an established tool in the treatment of correlated electron systems, notably for the description of competing instabilities. In recent years, methodological advancements led to the multiloop extension of the fRG, which systematically includes loop corrections beyond the conventional one-loop truncation and yields a quantitatively accurate description of two-dimensional lattice systems. At the same time, the single-boson exchange (SBE) decomposition of the two-particle vertex has been shown to offer both computational and interpretative advantages paving the way to more affordable approximation schemes. We here apply their combination coined as multiloop SBE fRG to the two-dimensional Hubbard model at weak coupling. After providing a detailed account of the underlying formalism in physical channels, we analyze the results for the frequency- and momentum-dependent vertex functions. We find that the SBE approximation, i.e., neglecting the flow of the multi-boson exchange contributions, accurately reproduces the parquet approximation at loop convergence. The presented algorithmic improvement opens the route for the treatment of more challenging parameter regimes and more realistic models.