Chiral Interactions up to Next-to-Next-to-Next-to-Leading Order and Nuclear Saturation

Wavefunction matching for solving quantum many-body problems

Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions1-3 to quantum chemistry4-6 and from atomic and molecular systems7-9 to nuclear physics10-14. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations15,16 of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions17,18 and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations19,20.

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M⊙ neutron star to be [Formula: see text] km.

Strongly interacting matter exhibits deconfined behavior in massive neutron stars

Neutron-star cores contain matter at the highest densities in our Universe. This highly compressed matter may undergo a phase transition where nuclear matter melts into deconfined quark matter, liberating its constituent quarks and gluons. Quark matter exhibits an approximate conformal symmetry, predicting a specific form for its equation of state (EoS), but it is currently unknown whether the transition takes place inside at least some physical neutron stars. Here, we quantify this likelihood by combining information from astrophysical observations and theoretical calculations. Using Bayesian inference, we demonstrate that in the cores of maximally massive stars, the EoS is consistent with quark matter. We do this by establishing approximate conformal symmetry restoration with high credence at the highest densities probed and demonstrating that the number of active degrees of freedom is consistent with deconfined matter. The remaining likelihood is observed to correspond to EoSs exhibiting phase-transition-like behavior, treated as arbitrarily rapid crossovers in our framework.

QCD Equation of State of Dense Nuclear Matter from a Bayesian Analysis of Heavy-Ion Collision Data

Bayesian methods are used to constrain the density dependence of the QCD equation of state (EOS) for dense nuclear matter using the data of mean transverse kinetic energy and elliptic flow of protons from heavy ion collisions (HICs), in the beam energy range sqrt[s_{NN}]=2-10  GeV. The analysis yields tight constraints on the density dependent EOS up to 4 times the nuclear saturation density. The extracted EOS yields good agreement with other observables measured in HIC experiments and constraints from astrophysical observations both of which were not used in the inference. The sensitivity of inference to the choice of observables is also discussed.

Degeneracy in the Inference of Phase Transitions in the Neutron Star Equation of State from Gravitational Wave Data

Gravitational wave (GW) detections of binary neutron star inspirals will be crucial for constraining the dense matter equation of state (EOS). We demonstrate a new degeneracy in the mapping from tidal deformability data to the EOS, which occurs for models with strong phase transitions. We find that there exists a new family of EOS with phase transitions that set in at different densities and that predict neutron star radii that differ by up to ∼500  m but that produce nearly identical tidal deformabilities for all neutron star masses. Next-generation GW detectors and advances in nuclear theory may be needed to resolve this degeneracy.

Plausible presence of new state in neutron stars with masses above 0.98MTOV

We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030 + 0451, PSR J0740 + 6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory (χEFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forward neural network model embedded in the Bayesian nonparametric inference, the structure of NS matter's sound speed c_s is explored in a model-agnostic way. It is found that a peak structure is common in the c_s² posterior, locating at 2.4-4.8ρ_sat (nuclear saturation density) and c_s² exceeds c²/3 at 90% credibility. The non-monotonic behavior suggests evidence of the state deviating from the hadronic matter inside the very massive NSs. Assuming the new/exotic state is featured as it is softer than typical hadronic models or even with hyperons, we find that a sizable (⩾10^-3M_⊙) exotic core, likely made of quark matter, is plausible for the NS with a gravitational mass above about 0.98M_TOV, where M_TOV represents the maximum gravitational mass of a non-rotating cold NS. The inferred M_TOV=2.18_-0.13^+0.27M_⊙ (90% credibility) is well consistent with the value of 2.17_-0.12^+0.15M_⊙ estimated independently with GW170817/GRB 170817A/AT2017gfo assuming a temporary supramassive NS remnant formed after the merger. PSR J0740 + 6620, the most massive NS detected so far, may host a sizable exotic core with a probability of ≈0.36.

Gravitational Wave Signal for Quark Matter with Realistic Phase Transition

The cores of neutron stars (NSs) near the maximum mass can realize a transitional change to quark matter (QM). Gravitational waves from binary NS mergers are expected to convey information about the equation of state (EOS) sensitive to the QM transition. Here, we present the first results of gravitational wave simulation with the realistic EOS that is consistent with ab initio approaches of χEFT and pQCD and is assumed to go through smooth crossover. We compare them to results obtained with another EOS with a first-order hadron-quark phase transition. Our results suggest that early collapse to a black hole in the post-merger stage after NS merger robustly signifies softening of the EOS associated with the QM onset in the crossover scenario. The nature of the hadron-quark phase transition can be further constrained by the condition that electromagnetic counterparts should be energized by the material left outside the remnant black hole.

Nuclear Equation of State for Arbitrary Proton Fraction and Temperature Based on Chiral Effective Field Theory and a Gaussian Process Emulator

We calculate the equation of state of asymmetric nuclear matter at finite temperature based on chiral effective field theory interactions to next-to-next-to-next-to-leading order. Our results assess the theoretical uncertainties from the many-body calculation and the chiral expansion. Using a Gaussian process emulator for the free energy, we derive the thermodynamic properties of matter through consistent derivatives and use the Gaussian process to access arbitrary proton fraction and temperature. This enables a first nonparametric calculation of the equation of state in beta equilibrium, and of the speed of sound and the symmetry energy at finite temperature. Moreover, our results show that the thermal part of the pressure decreases with increasing densities.

The Symmetry Energy: Current Status of Ab Initio Predictions vs. Empirical Constraints

Infinite nuclear matter is a suitable laboratory to learn about nuclear forces in many-body systems. In particular, modern theoretical predictions of neutron-rich matter are timely because of recent and planned experiments aimed at constraining the equation of state of isospin-asymmetric matter. For these reasons, we have taken a broad look at the equation of state of neutron-rich matter and the closely related symmetry energy, which is the focal point of this article. Its density dependence is of paramount importance for a number of nuclear and astrophysical systems, ranging from neutron skins to the structure of neutron stars. We review and discuss ab initio predictions in relation to recent empirical constraints. We emphasize and demonstrate that free-space nucleon–nucleon data pose stringent constraints on the density dependence of the neutron matter equation of state, which essentially determines the slope of the symmetry energy at saturation.

Ab initio predictions link the neutron skin of 208Pb to nuclear forces

Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 109 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape.

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars1-9 and from heavy-ion collisions of gold nuclei at relativistic energies10,11 with microscopic nuclear theory calculations12-17 to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission5-8,18. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

How Perturbative QCD Constrains the Equation of State at Neutron-Star Densities

We demonstrate in a general and analytic way how high-density information about the equation of state (EOS) of strongly interacting matter obtained using perturbative quantum chromodynamics constrains the same EOS at densities reachable in physical neutron stars. Our approach is based on utilizing the full information of the thermodynamic potentials at the high-density limit together with thermodynamic stability and causality. This requires considering the pressure as a function of chemical potential p(μ) instead of the commonly used pressure as a function of energy density p(ε). The results can be used to propagate the perturbative quantum chromodynamics calculations reliable around 40n_{s} to lower densities in the most conservative way possible. We constrain the EOS starting from only a few times the nuclear saturation density n≳2.2n_{s}, and at n=5n_{s} we exclude at least 65% of otherwise allowed area in the ε-p plane. This provides information complementary to astrophysical observations that should be taken into account in any complete statistical inference study of the EOS. These purely theoretical results are independent of astrophysical neutron-star input, and hence, they can also be used to test theories of modified gravity and beyond the standard model physics in neutron stars.

The Equation of State of Neutron-Rich Matter at Fourth Order of Chiral Effective Field Theory and the Radius of a Medium-Mass Neutron Star

We report neutron star predictions based on our most recent equations of state. These are derived from chiral effective field theory, which allows for a systematic development of nuclear forces, order by order. We utilize high-quality two-nucleon interactions and include all three-nucleon forces up to fourth order in the chiral expansion. Our ab initio predictions are restricted to the domain of applicability of chiral effective field theory. However, stellar matter in the interior of neutron stars can be up to several times denser than normal nuclear matter at saturation, and its composition is essentially unknown. Following established practices, we extend our microscopic predictions to higher densities matching piecewise polytropes. The radius of the average-size neutron star, about 1.4 solar masses, is sensitive to the pressure at normal densities, and thus it is suitable to constrain ab initio theories of the equation of state. For this reason, we focus on the radius of medium-mass stars. We compare our results with other theoretical predictions and recent constraints.

Holographic dense QCD in the Veneziano limit

Solving the properties of dense QCD matter is an extremely challenging problem because standard theoretical tools do not work at intermediate densities. The gauge/gravity duality may help to provide answers in this region. I give a brief review of recent progress in this field, focusing on the V-QCD model, which is one of the most sophisticated holographic models of QCD. I discuss predictions for the phase diagram, the equation of state, and properties of baryons. I apply these results to analyze the properties of neutron stars and to quark matter production in neutron star mergers.

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ^{208}Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283±0.071  fm for the neutron-skin thickness of ^{208}Pb, implying a slope parameter L=106±37  MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S_{0}, L, and R_{skin}^{^{208}Pb} directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields S_{0}=33.0_{-1.8}^{+2.0}  MeV, L=53_{-15}^{+14}  MeV, and R_{skin}^{^{208}Pb}=0.17_{-0.04}^{+0.04}  fm.

Global Microscopic Description of Nucleon-Nucleus Scattering with Quantified Uncertainties

We develop for the first time a microscopic global nucleon-nucleus optical potential with quantified uncertainties suitable for analyzing nuclear reaction experiments at next-generation rare-isotope beam facilities. Within the improved local density approximation and without any adjustable parameters, we begin by computing proton-nucleus and neutron-nucleus optical potentials from a set of five nuclear forces from chiral effective field theory for 1800 target nuclei in the mass range 12≤A≤242 for energies between 0  MeV<E≲150  MeV. We then parameterize a global optical potential for each chiral force that depends smoothly on the projectile energy as well as the target nucleus mass number and isospin asymmetry. Uncertainty bands for elastic scattering observables are generated from a full covariance analysis of the parameters entering in the description of our global optical potential and benchmarked against existing experimental data for stable target nuclei. Since our approach is purely microscopic, we anticipate a similar quality of the model for nucleon scattering on unstable isotopes.

Hyperons in Neutron Stars

I review the issues related to the appearance of hyperons in neutron star matter, focusing in particular on the problem of the maximum mass supported by hyperonic equations of state. I discuss the general mechanism that leads to the formation of hyperons in the core of neutron stars and I review the main techniques and many-body methods used to construct an appropriate equation of state to describe the strongly interacting system of hadrons hosted in the core of neutron stars. I outline the consequences on the structure and internal composition of neutron stars and also discuss the possible signatures of the presence of hyperons in astrophysical dynamical systems like supernova explosions and binary neutron star mergers. Finally, I briefly report about the possible important role played by hyperons in the transport properties of neutron star matter and on the consequences of neutron star cooling and gravitational wave instabilities induced by the presence of hyperons.

Unified Equation of State for Neutron Stars Based on the Gogny Interaction

The effective Gogny interactions of the D1 family were established by D. Gogny more than forty years ago with the aim to describe simultaneously the mean field and the pairing field corresponding to the nuclear interaction. The most popular Gogny parametrizations, namely D1S, D1N and D1M, describe accurately the ground-state properties of spherical and deformed finite nuclei all across the mass table obtained with Hartree–Fock–Bogoliubov (HFB) calculations. However, these forces produce a rather soft equation of state (EoS) in neutron matter, which leads to predict maximum masses of neutron stars well below the observed value of two solar masses. To remove this limitation, we built new Gogny parametrizations by modifying the density dependence of the symmetry energy predicted by the force in such a way that they can be applied to the neutron star domain and can also reproduce the properties of finite nuclei as good as their predecessors. These new parametrizations allow us to obtain stiffer EoS’s based on the Gogny interactions, which predict maximum masses of neutron stars around two solar masses. Moreover, other global properties of the star, such as the moment of inertia and the tidal deformability, are in harmony with those obtained with other well tested EoSs based on the SLy4 Skyrme force or the Barcelona–Catania–Paris–Madrid (BCPM) energy density functional. Properties of the core-crust transition predicted by these Gogny EoSs are also analyzed. Using these new Gogny forces, the EoS in the inner crust is obtained with the Wigner–Seitz approximation in the Variational Wigner–Kirkwood approach along with the Strutinsky integral method, which allows one to estimate in a perturbative way the proton shell and pairing corrections. For the outer crust, the EoS is determined basically by the nuclear masses, which are taken from the experiments, wherever they are available, or by HFB calculations performed with these new forces if the experimental masses are not known.

Normalizing Flows for Microscopic Many-Body Calculations: An Application to the Nuclear Equation of State

Normalizing flows are a class of machine learning models used to construct a complex distribution through a bijective mapping of a simple base distribution. We demonstrate that normalizing flows are particularly well suited as a Monte Carlo integration framework for quantum many-body calculations that require the repeated evaluation of high-dimensional integrals across smoothly varying integrands and integration regions. As an example, we consider the finite-temperature nuclear equation of state. An important advantage of normalizing flows is the ability to build highly expressive models of the target integrand, which we demonstrate enables precise evaluations of the nuclear free energy and its derivatives. Furthermore, we show that a normalizing flow model trained on one target integrand can be used to efficiently calculate related integrals when the temperature, density, or nuclear force is varied. This work will support future efforts to build microscopic equations of state for numerical simulations of supernovae and neutron star mergers that employ state-of-the-art nuclear forces and many-body methods.

Nuclear Physics and Astrophysics Constraints on the High Density Matter Equation of State

(1) This review has been written in memory of Steven Moszkowski who unexpectedly passed away in December 2020. It has been inspired by our many years of discussions. Steven’s enthusiasm, drive and determination to understand atomic nuclei in simple terms of basic laws of physics was infectious. He sought the fundamental origin of nuclear forces in free space, and their saturation and modification in nuclear medium. His untimely departure left our job unfinished but his legacy lives on. (2) Focusing on the nuclear force acting in nuclear matter of astrophysical interest and its equation of state (EoS), we take several typical snapshots of evolution of the theory of nuclear forces. We start from original ideas in the 1930s moving through to its overwhelming diversity today. The development is supported by modern observational and terrestrial data and their inference in the multimessenger era, as well as by novel mathematical techniques and computer power. (3) We find that, despite the admirable effort both in theory and measurement, we are facing multiple models dependent on a large number of variable correlated parameters which cannot be constrained by data, which are not yet accurate, nor sensitive enough, to identify the theory closest to reality. The role of microphysics in the theories is severely limited or neglected, mostly deemed to be too difficult to tackle. (4) Taking the EoS of high-density matter as an example, we propose to develop models, based, as much as currently possible, on the microphysics of the nuclear force, with a minimal set of parameters, chosen under clear physical guidance. Still somewhat phenomenological, such models could pave the way to realistic predictions, not tracing the measurement, but leading it.

Ab Initio Limits of Atomic Nuclei

We predict the limits of existence of atomic nuclei, the proton and neutron drip lines, from the light through medium-mass regions. Starting from a chiral two- and three-nucleon interaction with good saturation properties, we use the valence-space in-medium similarity renormalization group to calculate ground-state and separation energies from helium to iron, nearly 700 isotopes in total. We use the available experimental data to quantify the theoretical uncertainties for our ab initio calculations towards the drip lines. Where the drip lines are known experimentally, our predictions are consistent within the estimated uncertainty. For the neutron-rich sodium to chromium isotopes, we provide predictions to be tested at rare-isotope beam facilities.

How Well Do We Know the Neutron-Matter Equation of State at the Densities Inside Neutron Stars? A Bayesian Approach with Correlated Uncertainties

We introduce a new framework for quantifying correlated uncertainties of the infinite-matter equation of state derived from chiral effective field theory (χEFT). Bayesian machine learning via Gaussian processes with physics-based hyperparameters allows us to efficiently quantify and propagate theoretical uncertainties of the equation of state, such as χEFT truncation errors, to derived quantities. We apply this framework to state-of-the-art many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χEFT expansion. This produces the first statistically robust uncertainty estimates for key quantities of neutron stars. We give results up to twice nuclear saturation density for the energy per particle, pressure, and speed of sound of neutron matter, as well as for the nuclear symmetry energy and its derivative. At nuclear saturation density, the predicted symmetry energy and its slope are consistent with experimental constraints.

Symmetric Nuclear Matter from the Strong Interaction

We study the equation of state of symmetric nuclear matter at zero temperature over a wide range of densities using two complementary theoretical approaches. At low densities, up to twice nuclear saturation density, we compute the energy per particle based on modern nucleon-nucleon and three-nucleon interactions derived within chiral effective field theory. For higher densities, we derive for the first time constraints in a Fierz-complete setting directly based on quantum chromodynamics using functional renormalization group techniques. We find remarkable consistency of the results obtained from both approaches as they come together in density and the natural emergence of a maximum in the speed of sound c_{S} at supranuclear densities. The presence of this maximum appears tightly connected to the formation of a diquark gap. Notably, this maximum is observed to exceed the asymptotic value c_{S}^{2}=1/3 while its exact position in terms of the density cannot yet be determined conclusively.

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of ^{48}Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. We use our scheme to compute the matrix element that governs the neutrinoless double beta decay of ^{48}Ca to ^{48}Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as ^{76}Ge, ^{130}Te, and ^{136}Xe.

Equation of State Effects in Core-Collapse Supernovae

We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the different effects in changing these properties from the Lattimer and Swesty to the Shen et al. equations of state, the two most commonly used equations of state in simulations. In particular, we trace the contraction behavior to the effective mass, which determines the thermal nucleonic contributions to the equation of state. Larger effective masses lead to lower pressures at nuclear densities and a lower thermal index. This results in a more rapid contraction of the proto-neutron star and consequently higher neutrino energies, which aids the shock evolution to a faster explosion.

Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus

We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus ^{16}O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.032501]. With this method, we compute the binding energy and charge radius of ^{16}O at more than 10^{6} different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the ^{3}S_{1} wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.

Harmonic Potential Theorem: Extension to Spin-, Velocity-, and Density-Dependent Interactions

One of the few exact results for the description of the time evolution of an inhomogeneous, interacting many-particle system is given by the harmonic potential theorem (HPT). The relevance of this theorem is that it sets a tight constraint on time-dependent many-body approximations. In this contribution, we show that the original formulation of the HPT is valid also for the case of spin-, velocity-, and density-dependent interactions. This result is completely general and relevant, among the rest, for nuclear structure theory both in the case of ab initio and of more phenomenological approaches. As an example, we report on a numerical implementation by testing the small-amplitude limit of the time-dependent Hartree-Fock-also known as the random phase approximation-for the translational frequencies of a neutron system trapped in a harmonic potential.