DRalgo: Dimensional-Reduction EFT Toolkit
- DRalgo is an algorithmic toolkit that performs dimensional reduction by integrating out nonzero Matsubara modes to produce effective 3D field theories for thermal phase transitions.
- It automatically computes next-to-leading order corrections—including two-loop thermal masses, beta functions, and Debye screening—ensuring refined and controlled matching.
- The package exports matching formulas and effective potentials in Mathematica or Python formats, facilitating seamless integration with phase-transition and gravitational-wave analyses.
DRalgo is an algorithmic implementation and Mathematica package for constructing an effective, dimensionally reduced, high-temperature field theory for generic models. In the thermal phase-transition literature, it denotes a workflow that integrates out nonzero Matsubara modes and, when required, temporal gauge scalars to obtain three-dimensional effective field theories at the hard, soft, and ultrasoft scales; performs matching to next-to-leading order; computes beta functions, anomalous dimensions, Debye screening masses, and perturbative thermal potentials; and exports these results for numerical phase-transition, bubble-nucleation, and gravitational-wave calculations (Ekstedt et al., 2022, Bertenstam et al., 2 Jan 2025).
1. Definition, scope, and development
The original DRalgo release presented the package as “a package for effective field theory approach for thermal phase transitions” and described it as an algorithmic implementation that constructs an effective, dimensionally reduced, high-temperature field theory for generic models. Its stated core capabilities are automatic matching to next-to-leading order, including two-loop thermal corrections to scalar and Debye masses as well as one-loop thermal corrections to couplings; the possibility of integrating out additional heavy scalars; leading-order beta functions for general gauge-charges and fermion-families in both the fundamental and effective theory; and computation of the finite-temperature effective potential within the effective theory (Ekstedt et al., 2022).
Subsequent work expanded the package from super-renormalizable three-dimensional operators to the one-loop matching of generic three-dimensional dimension-five and dimension-six operators for arbitrary models containing scalars, fermions, and gauge fields. In that extension, DRalgo is described as a Mathematica-based toolkit that provides the complete three-dimensional operator bases up to mass dimension six, generates the tensor structure of each Wilson coefficient in terms of the original four-dimensional couplings, computes the one-loop thermal sum-integrals in a general background-field gauge, and solves the matching relations for any subset of operators in a fully symbolic framework (Bernardo et al., 14 May 2026).
A possible source of confusion is purely terminological. In the supplied material, unrelated descriptions also attach the string “DRalgo” to a distributed augmented Lagrangian method and to a distributed retrieval-augmented generation framework. In the arXiv phase-transition literature considered here, however, DRalgo denotes the dimensional-reduction package and its extensions.
2. Dimensional-reduction framework
The package is built around high-temperature dimensional reduction. At temperatures any zero-temperature mass, a four-dimensional gauge theory can be reduced to a three-dimensional Euclidean effective field theory by integrating out nonzero Matsubara modes and, optionally, additional heavy fields whose masses satisfy . The result is a super-renormalizable three-dimensional theory for the static modes (Ekstedt et al., 2022).
The scale hierarchy is organized explicitly:
| Scale | Parametric size | DR step |
|---|---|---|
| hard | full Matsubara tower | |
| soft | after integrating out nonzero modes | |
| ultrasoft | after integrating out temporal scalars |
This separation is central to the package’s treatment of infrared sensitivity. In the standard four-dimensional high- expansion, bosonic zero modes produce infrared divergences and motivate daisy resummation. The dimensional-reduction formulation replaces that ad hoc reorganization with a systematic EFT construction in which ultraviolet and infrared physics are separated cleanly. In the high- regime, the dimensional-reduction procedure is stated to be equivalent—order by order in a combined coupling/high-0 expansion—to thermal-mass daisy resummation, but with the advantage of a controlled EFT framework (Brdar et al., 7 May 2025).
The reduced theory can be written schematically as
1
with tree-level potential
2
The one-loop correction is UV-finite in three dimensions,
3
where the sum runs over the light three-dimensional bosonic degrees of freedom, including scalars and longitudinal gauge fields (Athron et al., 2024).
3. Matching, running, and effective potentials
DRalgo implements hard-to-soft and soft-to-ultrasoft dimensional reduction for a general renormalizable four-dimensional Lagrangian. In the original package description, the matching order is next-to-leading order: all two-loop thermal corrections to scalar masses and to Debye masses, and all one-loop thermal corrections to scalar quartics, gauge couplings, Yukawas, and 4 self-couplings (Ekstedt et al., 2022).
The package also computes all required four-dimensional and three-dimensional beta functions, thermal masses, anomalous dimensions, and two-loop thermal corrections to scalar masses. It automatically generates analytical expressions for the three-dimensional parameters, commonly denoted by hats, such as 5, 6, and 7, together with the NLO contribution to the three-dimensional potential,
8
and supports export of the matching formulae and three-dimensional potentials in Mathematica or Python form (Bertenstam et al., 2 Jan 2025).
A central mapping relates the three-dimensional EFT back to the thermal four-dimensional effective potential:
9
In the leptoquark application, the same relation is written as
0
This mapping is what allows DRalgo-generated EFTs to be passed directly into bounce-action and gravitational-wave solvers (Bertenstam et al., 2 Jan 2025).
The supercooled-phase-transition analysis of the Abelian Higgs model makes explicit how DRalgo combines matching with renormalisation-group improvement. There, consistent RGE evolution of the couplings is described as essential for a meaningful interpretation of the results, and the four-dimensional parameters are matched onto three-dimensional EFT parameters order by order in a loop expansion. Representative one-loop matchings are
1
2
while the scalar mass and Debye mass are matched at two loops (Christiansen et al., 4 Nov 2025).
The same study distinguishes several three-dimensional loop-ordering schemes: “3D 1-L (LO),” “3D 1-L (NLO),” “3D 2-L (NLO),” and “3D 2-L (Mixed).” This suggests that DRalgo is not merely a code generator for one fixed perturbative truncation, but a framework for controlled comparisons between different EFT expansion schemes when residual scale dependence and convergence are under scrutiny.
4. Software architecture and user workflow
The package is implemented in Mathematica and exposes a model-building and matching workflow in terms of named routines. A typical usage pattern begins by loading or defining the four-dimensional model, for example through ImportModelDRalgo[...] after tensor allocation and construction of invariant structures. The hard-to-soft step is then executed with PerformDRhard[], after which the package can print soft EFT couplings and masses through calls such as PrintCouplings[], PrintScalarMass["LO"], PrintScalarMass["NLO"], PrintDebyeMass["LO"], PrintDebyeMass["NLO"], and BetaFunctions4D[]. If desired, PerformDRsoft[{heavyIndices}] integrates out soft heavy scalars or temporal gauge fields to produce an ultrasoft EFT, with corresponding accessors such as PrintCouplingsUS[], PrintScalarMassUS[...], and BetaFunctions3DUS[]. The perturbative three-dimensional effective potential is then assembled via CalculatePotentialUS[] and retrieved with PrintEffectivePotential["LO"/"NLO"/"NNLO"] (Ekstedt et al., 2022).
The internal logic of later pipelines makes the same structure explicit in algorithmic terms. In the leptoquark study, the DRalgo-based thermal potential is constructed by first RG-evolving four-dimensional parameters to 3, typically with 4; then performing hard-to-soft matching; then choosing an ultrasoft matching scale 5, with 6 by default; then integrating out temporal fields; and finally building the three-dimensional effective potential before վերադարձing the four-dimensional form 7 (Bertenstam et al., 2 Jan 2025).
DRalgo’s output format is intentionally designed for downstream use. The package can export matching relations and potentials in Mathematica or Python form, and the higher-dimensional-operator extension continues this design by returning Mathematica rule lists for Wilson coefficients. In that extension, Dimension5Matching[...] and Dimension6Matching[...] solve the hard-mode matching problem symbolically, while HardThermal1LoopInt["B",s,\alpha,d] and HardThermal1LoopInt["F",s,\alpha,d] expose the underlying bosonic and fermionic thermal master integrals (Bernardo et al., 14 May 2026).
5. Integration with phase-transition and gravitational-wave toolchains
DRalgo has become part of several end-to-end software stacks for cosmological phase transitions. The most explicit interfaces described in the literature are Dratopi, PhaseTracer2, and PT2GWFinder.
| Interface | How DRalgo is used | Notable detail |
|---|---|---|
| Dratopi | reads DRalgo-generated Python code and builds ultrasoft EFT parameters | wraps a modified CosmoTransitions |
| PhaseTracer2 | hand-translates DRalgo Mathematica output into C++ model methods | uses solveBetas() and get3d_parameters() |
| PT2GWFinder / DRTools | imports DRalgo output and reconstructs 8 | ComputeDRPotential[...] returns 9 |
Dratopi is described as “a thin Python layer” that reads in the DRalgo-generated Python code, numerically solves the four-dimensional RGEs from an input scale 0 to 1 while enforcing user-specified constraints such as perturbativity and bounded-from-below conditions, uses the DRalgo matching routines to build the ultrasoft three-dimensional EFT parameters 2, provides an auto-diagonalization of field-dependent mass matrices for the NLO thermal potential, and wraps CosmoTransitions with minor patches so that the potential
3
is called at nonzero 4 and the condition 5 defines the temperature range. It also offers “sanity checks,” notably the ratio 6, to monitor high-7 perturbativity (Bertenstam et al., 2 Jan 2025).
PhaseTracer2 incorporates dimensionally reduced effective potentials for models obtained from DRalgo. Rather than file I/O, it ingests the output through two user-written member functions in a model class deriving from EffectivePotential::Potential: solveBetas(x,t), which solves the one-loop 8 beta functions from the input scale 9 to 0, and get3d_parameters(T), which evaluates the four-dimensional couplings at 1{m_32,\lambda_3,g_32,\dots}A_0$5
imports raw DRalgo expressions, solves the four-dimensional RG equations numerically over the chosen range, applies the matching formulae at each temperature via DRStep, reconstructs the analytic three-dimensional potential via DRPotentialN, and optionally rescales it to the four-dimensional potential. The primary returned object is the closed-form function
$\mu_{4d}\approx \pi T$2
which can then be passed directly to the package’s transition finder (Brdar et al., 7 May 2025).
6. Applications, validation, and later extensions
The package has been applied to models with qualitatively different thermal histories. In a minimal leptoquark extension of the Standard Model that explains active neutrino oscillation data while satisfying current flavor physics constraints, the DRalgo–Dratopi chain automated the construction of a three-dimensional high-$\mu_{4d}\approx \pi T$3 EFT from the original four-dimensional model, computed the resulting three-dimensional thermal potential at NLO, and fed that potential into bubble-nucleation and gravitational-wave routines. The study reports diverse phase-transition patterns, including color symmetry-breaking scenarios in the early Universe, and states that a detectable signal in the $\mu_{4d}\approx \pi T$4–$\mu_{4d}\approx \pi T$5 frequency range features color-restoration and leptoquark masses near $\mu_{4d}\approx \pi T$6 (Bertenstam et al., 2 Jan 2025).
That application also provides explicit performance and robustness observations. The DRalgo+Dratopi chain was able to scan $\mu_{4d}\approx \pi T$7 model points, producing $\mu_{4d}\approx \pi T$8 viable phase transitions with $\mu_{4d}\approx \pi T$9 in a matter of hours on a modern multicore workstation. Varying the hard-matching prefactor $\mu_{3d}^{s}\approx gT$0 changed the predicted GW peak amplitude by $\mu_{3d}^{s}\approx gT$1, which is presented as evidence that the explicit removal of large logarithms and inclusion of NLO terms shrinks theoretical uncertainties compared to standard four-dimensional high-$\mu_{3d}^{s}\approx gT$2 potentials. The high-$\mu_{3d}^{s}\approx gT$3 perturbativity diagnostic $\mu_{3d}^{s}\approx gT$4 clustered around $\mu_{3d}^{s}\approx gT$5–$\mu_{3d}^{s}\approx gT$6 for all transitions, and imposing $\mu_{3d}^{s}\approx gT$7 eliminated only very few points on the edge (Bertenstam et al., 2 Jan 2025).
Independent validation comes from software integrations. In PhaseTracer2, DR-EFT potentials are described as simpler than full four-dimensional one-loop potentials because they avoid the special functions $\mu_{3d}^{s}\approx gT$8 and $\mu_{3d}^{s}\approx gT$9 and instead enter through algebraic $\mu_{3d}^{Us}\approx g^2T$0 terms. In a typical two-field example, the paper reports a $\mu_{3d}^{Us}\approx g^2T$1–$\mu_{3d}^{Us}\approx g^2T$2 speed-up in the phase-structure scan relative to the full four-dimensional Coleman–Weinberg plus daisy potential, negligible overhead in RG running at $\mu_{3d}^{Us}\approx g^2T$3 per temperature point, agreement of critical temperatures and latent heats at $\mu_{3d}^{Us}\approx g^2T$4 with the full four-dimensional $\mu_{3d}^{Us}\approx g^2T$5 plus Parwani result for benchmark singlet-model points, and improved stability in strong supercooling because infrared modes are integrated out analytically (Athron et al., 2024). In PT2GWFinder, the same $\mu_{3d}^{Us}\approx g^2T$6 grid scan in $\mu_{3d}^{Us}\approx g^2T$7 that took about $\mu_{3d}^{Us}\approx g^2T$8 with the Coleman–Weinberg plus daisy implementation ran in roughly $\mu_{3d}^{Us}\approx g^2T$9 with the DR approach, while yielding qualitatively similar and in many regions quantitatively similar phase-transition observables (Brdar et al., 7 May 2025).
The supercooled dark-sector study sharpened the package’s role in precision comparisons between approximation schemes. There, the three-dimensional EFT with consistent expansion in the four-dimensional parameters is reported to have significantly reduced scale dependence. For the benchmark “dark photon” model, varying scales in the four-dimensional high-temperature one-loop scheme produces a fractional variation $A_0$0, while the three-dimensional two-loop scheme yields $A_0$1. The paper further states that the four-dimensional high-temperature one-loop and three-dimensional two-loop schemes give nearly identical $A_0$2 curves and agree on the nucleation temperature to within a few percent over the range $A_0$3, while the propagated nanohertz gravitational-wave peak amplitude and frequency differ by less than $A_0$4; by contrast, the simple analytic parametrisation can deviate by tens of percent in the limit of large supercooling (Christiansen et al., 4 Nov 2025).
The 2026 higher-dimensional-operator extension broadens DRalgo’s remit from super-renormalizable three-dimensional EFTs to complete operator bases through dimension six. This extension is motivated by the statement that strong first-order cosmological transitions often require subleading hard-mode effects to be under analytic control, and that these higher-dimensional operators are crucial for quantifying the convergence of the high-temperature expansion and stabilizing nonperturbative lattice studies of three-dimensional EFTs. The package now automates matching for arbitrary models of scalars, fermions, and gauge fields, but still reports coefficients in a redundant off-shell basis with explicit gauge-parameter dependence; elimination of redundancies and associated field redefinitions are described as being left to the user or to future DRalgo updates (Bernardo et al., 14 May 2026).
Taken together, these developments place DRalgo at the intersection of thermal EFT construction, perturbative control of infrared physics, renormalisation-group improvement, and end-to-end phenomenology of cosmological first-order phase transitions. A plausible implication is that its main scientific role is not only to produce reduced theories, but to standardize the passage from an arbitrary renormalizable four-dimensional Lagrangian to numerically tractable three-dimensional actions whose uncertainties can be diagnosed through matching-scale variation, loop-order comparisons, and higher-dimensional-operator estimates.