The nuclear energy density functional (EDF) method is the only microscopic framework that reaches essentially the entire nuclear chart, yet its functionals are phenomenological: there is no identified expansion parameter, no systematic order-by-order improvement, no quantified breakdown scale, and no controlled treatment of the fact that nuclei are self-bound systems with no external potential. The source paper (Furnstahl, Eur. Phys. J. A 56, 85 (2020)) asks whether the EDF method can be turned into a proper effective field theory (EFT) and proposes a path based on an effective-action formulation with BRST (Becchi–Rouet–Stora–Tyutin) handling of zero modes. The question remains open. Later work resolved pieces of it — hybrid functionals show that long-range pion terms improve the EDF without added refitting flexibility, and projection-based emulators substantially advance the feasibility strand — but the core is untouched: no in-medium power counting with a demonstrated expansion parameter for finite nuclei exists, and the BRST treatment of self-bound nuclei has not been demonstrated beyond toy models.

Background

Nuclear structure theory describes medium-mass and heavy nuclei mainly through the energy density functional (EDF) method: the ground-state energy of a nucleus is modeled as a functional of nucleon densities and currents (Skyrme-type zero-range, Gogny-type finite-range, or covariant functionals), minimized in a Hartree–Fock–Bogoliubov mean-field approximation — a self-consistent one-particle potential supplemented by a pairing field. Fitted to data with a dozen or so parameters, these functionals reproduce masses and radii across the table of nuclides and capture emergent phenomena such as nuclear saturation (the approximately constant density of nuclear matter) and pairing. Their weaknesses are equally well known: there is no principled rule for adding or improving terms, extrapolations to the driplines (the edges of the chart of nuclides beyond which nuclei are no longer bound), superheavy nuclei, or neutron-star matter come with unquantified model dependence, and single-particle energies and currents are extracted without a controlled many-body justification.

An effective field theory (EFT) is a field theory written in the degrees of freedom active at the resolution of interest, with every term allowed by the symmetries organized by a power counting — a rule ranking terms in powers of a small expansion parameter, so that calculations can be improved order by order with quantifiable truncation errors until the breakdown scale of the expansion is reached. For nuclei the reference EFT is chiral EFT, built from nucleons and pions with long-range pion exchange fixed by chiral symmetry; a pionless EFT, valid below the pion mass scale, treats only contact interactions. Hints that empirical EDFs already encode an EFT include the naturalness of their coupling constants (they scale as naive dimensional analysis would suggest) and the fact that only a few parameter combinations carry most of the physics. The density matrix expansion (DME), a technique that maps finite-range pion-exchange physics into a local functional of the density and its gradients, is one established bridge between chiral EFT and Skyrme-like functionals.

What is missing for a proper EFT of the EDF method was crystallized by Furnstahl (Eur. Phys. J. A 56, 85 (2020)). First, the correct degrees of freedom are unresolved: pion-including, pionless, or quasinucleon densities and mean fields built around the nuclear liquid. Second, no in-medium expansion parameter is known: free-space Weinberg power counting (the chiral-EFT rule that ranks diagrams by powers of the ratio of a typical momentum to the breakdown scale) does not transplant to finite density (saturation in chiral EFT is driven by three-nucleon forces that only enter at third order, whereas a proper EDF EFT should saturate at leading order), Pauli blocking reshuffles diagram counting, and loop integrals at finite density show significant regulator artifacts that depend on cutoff choices. Third, nuclei are self-bound: no external potential anchors the density, so the mean-field solution spontaneously breaks translational (and often rotational) symmetry, the center-of-mass (COM) motion contaminates intrinsic observables, and exact projection onto good quantum numbers is costly. Furnstahl's proposed framework is the effective action — the Legendre transform of a generating functional with auxiliary fields, the field-theoretic underpinning of density functional theory and its Kohn-Sham decoupling — in which broken symmetry appears as the zero-mode problem of a loop expansion. There, collective coordinates can be gauge-fixed and the associated spurious degrees of freedom cancelled by BRST symmetry: one enlarges the Hilbert space with ghost fields and a BRST charge whose invariance eliminates zero modes and infrared divergences, yields a unique (projected) ground state, and survives approximation without generating zero modes. A pre-existing operator-based BRST literature for nuclear collective motion demonstrated the idea on algebraic models, but never within an effective-action EDF for realistic nuclei. The practicality question — whether such a construction can be evaluated at a cost comparable to present EDFs, possibly with emulators (fast surrogate models, such as eigenvector-continuation or reduced-basis emulators, that interpolate expensive parameter-dependent calculations) — completes the list of open issues.

Problem Statement

Construct and demonstrate a proper effective field theory of the nuclear energy density functional (EDF) method. The common objective uniting all requested results is a single, systematically improvable EDF grounded in EFT principles and applicable to nuclei across the nuclear chart. Specifically:

(i) Degrees of freedom and in-medium power counting: identify the appropriate degrees of freedom (pion-including chiral EFT fields, pionless/contact constructions, or quasinucleon densities and mean fields) and establish a power counting for the EDF valid at finite density — a demonstrable expansion parameter together with a breakdown scale — that organizes every term of the functional, including its density dependence, particle-hole and pairing channels, and higher-body (three- and more-body) contributions, with regulator and cutoff artifacts brought under control so that order-by-order results are cutoff independent within the quoted errors.

(ii) Self-bound nuclei and zero modes: within an effective-action formulation (or an equivalent construction achieving the same requirements), treat the fact that nuclei are self-bound: translational and rotational symmetry breaking, the associated zero modes, and spurious center-of-mass contributions must be consistently eliminated, with symmetry restoration emerging automatically rather than by explicit numerical projection. The BRST-based construction highlighted by the source as the most critical step of its proposed program is the suggested route, but any construction meeting these requirements is admissible.

(iii) Demonstration and systematics: show order-by-order improvement with quantified truncation errors for nuclear observables — ground-state energies, radii, and related bulk properties — on a progression from toy models through simple nuclei to realistic nuclei, at a computational cost compatible with present EDF practice (emulators may be used to achieve it).

The three components are jointly required: a power counting without a correct treatment of self-boundness (or vice versa) does not constitute a proper EFT of the EDF method. The source's own staged path — proofs of principle on toy models, then simple nuclear models, then realistic systems — defines the escalation; the question is answered when such a construction exists, its expansion parameter and breakdown scale are identified and validated, and its systematic improvement is demonstrated at a computational cost compatible with present EDF practice.

References

  1. R. J. Furnstahl, 'Turning the nuclear energy density functional method into a proper effective field theory: reflections', Eur. Phys. J. A 56, 85 (2020), DOI: 10.1140/epja/s10050-020-00095-y, arXiv:1906.00833, https://arxiv.org/abs/1906.00833
  2. R. J. Furnstahl, H.-W. Hammer, A. Schwenk, 'Nuclear Structure at the Crossroads', Few Body Syst. 62, 72 (2021), DOI: 10.1007/s00601-021-01658-5, arXiv:2107.00413, https://arxiv.org/abs/2107.00413
  3. S. Burrello, M. Grasso, C.-J. Yang, 'Towards a power counting in nuclear energy-density-functional theories through a perturbative analysis', Phys. Lett. B 811, 135938 (2020), DOI: 10.1016/j.physletb.2020.135938, arXiv:2010.12339, https://arxiv.org/abs/2010.12339
  4. J. A. Oller, 'An in-medium chiral power-counting scheme for nuclear matter and some applications', J. Phys. G 46, 073001 (2019), DOI: 10.1088/1361-6471/ab0997, arXiv:1902.06065, https://arxiv.org/abs/1902.06065
  5. L. Zurek, E. A. Coello Pérez, S. K. Bogner, R. J. Furnstahl, A. Schwenk, 'Comparing different density-matrix expansions for long-range pion exchange', Phys. Rev. C 103, 014325 (2021), DOI: 10.1103/PhysRevC.103.014325, arXiv:2010.12518, https://arxiv.org/abs/2010.12518
  6. L. Zurek, S. K. Bogner, R. J. Furnstahl, R. Navarro Pérez, N. Schunck, A. Schwenk, 'Optimized nuclear energy density functionals including long-range pion contributions', Phys. Rev. C 109, 014319 (2024), DOI: 10.1103/PhysRevC.109.014319, arXiv:2307.13568, https://arxiv.org/abs/2307.13568
  7. T. Duguet, A. Ekström, R. J. Furnstahl, S. König, D. Lee, 'Colloquium: Eigenvector continuation and projection-based emulators', Rev. Mod. Phys. 96, 031002 (2024), DOI: 10.1103/RevModPhys.96.031002, arXiv:2310.19419, https://arxiv.org/abs/2310.19419
  8. K. Fraboulet, J.-P. Ebran, 'Addressing energy density functionals in the language of path-integrals I: comparative study of diagrammatic techniques applied to the (0+0)-D O(N)-symmetric φ⁴-theory', Eur. Phys. J. A 59, 91 (2023), DOI: 10.1140/epja/s10050-023-00933-9, arXiv:2208.13044, https://arxiv.org/abs/2208.13044
  9. K. Fraboulet, J.-P. Ebran, 'Addressing energy density functionals in the language of path-integrals II: comparative study of functional renormalization group techniques applied to the (0+0)-D O(N)-symmetric φ⁴-theory', Eur. Phys. J. A 60, 45 (2024), DOI: 10.1140/epja/s10050-023-01069-6, arXiv:2210.07748, https://arxiv.org/abs/2210.07748
  10. P. Sharma, R. J. Furnstahl, 'Testing variational perturbation theory for effective actions using the Gaudin-Yang model', Phys. Rev. C 112, 044309 (2025), DOI: 10.1103/tg2k-qg9v, arXiv:2506.02919, https://arxiv.org/abs/2506.02919
  11. T. Duguet, J.-P. Ebran, M. Frosini, H. Hergert, V. Somà, 'Rooting the EDF method into the ab initio framework. PGCM-PT formalism based on MR-IMSRG pre-processed Hamiltonians', Eur. Phys. J. A 59, 13 (2023), DOI: 10.1140/epja/s10050-023-00914-y, arXiv:2209.03424, https://arxiv.org/abs/2209.03424
  12. A. Kievsky, G. Orlandini, M. Gattobigio, 'Many-body energy density functional', Phys. Rev. A 104, L030801 (2021), DOI: 10.1103/PhysRevA.104.L030801, arXiv:2106.13582, https://arxiv.org/abs/2106.13582
  13. M. Pavon Valderrama, 'Renormalization group evolution and power counting in nuclear matter', arXiv:2603.19892 (2026), https://arxiv.org/abs/2603.19892
  14. U. Vernik, K. Hebeler, A. Schwenk, 'Exploring quark-mass-dependent three-nucleon forces in medium-mass nuclei', Phys. Rev. C 113, 064305 (2026), DOI: 10.1103/gcpw-cy14, arXiv:2512.20454, https://arxiv.org/abs/2512.20454