Search

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., 6 pages, 4 figures, accepted to Physical Review Letters

Born in the aftermath of core collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yield novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons between nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state., Comment: additional references and minor improvements; 31 pages, 8 figures; invited review prepared for Annu. Rev. Nucl. Part. Sci. (2021)

We add an ensemble of nuclei to the equation of state for homogeneous nucleonic matter to generate a new set of models suitable for astrophysical simulations of core-collapse supernovae and neutron star mergers. We implement empirical constraints from (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. Our model is also guided by microscopic many-body theory calculations based on realistic nuclear forces, including the zero-temperature neutron matter equation of state from quantum Monte Carlo simulations and thermal contributions to the free energy from finite-temperature many-body perturbation theory. We ensure that the parameters of our model can be varied while preserving thermodynamic consistency and the connection to experimental or observational data, thus providing a probability distribution of the astrophysical hot and dense matter equation of state. We compare our results with those obtained from other available equations of state. While our probability distributions indeed represent a large number of possible equations of state, we cannot yet claim to have fully explored all of the uncertainties, especially with regard to the structure of nuclei in the hot and dense medium., Comment: 13 pages, 9 figures

We examine the properties of the isospin-asymmetry expansion of the nuclear equation of state from chiral two- and three-body forces. We focus on extracting the high-order symmetry energy coefficients that consist of both normal terms (occurring with even powers of the isospin asymmetry) as well as terms involving the logarithm of the isospin asymmetry that are formally nonanalytic around the expansion point of isospin-symmetric nuclear matter. These coefficients are extracted from numerically precise perturbation theory calculations of the equation of state coupled with a new set of finite difference formulas that achieve stability by explicitly removing the effects of higher-order terms in the expansion. We consider contributions to the symmetry energy coefficients from both two- and three-body interactions. It is found that the coefficients of the logarithmic terms are generically larger in magnitude than those of the normal terms from second-order perturbation theory diagrams, but overall the normal terms give larger contributions to the ground state energy. The high-order isospin-asymmetry terms are especially relevant at large densities where they affect the proton fraction in beta-equilibrium matter., 12 pages, 11 figures

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\ensuremath{\le}A\ensuremath{\le}242$ for energies between $0\text{ }\text{ }\mathrm{MeV}lE\ensuremath{\lesssim}150\text{ }\text{ }\mathrm{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.

We formulate a microscopic optical potential from chiral two- and three-body forces. The real and imaginary central terms of the optical potential are obtained from the nucleon self-energy in infinite matter, while the real spin–orbit term is extracted from a nuclear energy density functional constructed from the density matrix expansion using the same chiral potential. The density-dependent optical potential is then folded with the nuclear density distributions for selected Calcium isotopes resulting in energy-dependent nucleon–nucleus optical potentials from which we study proton–nucleus elastic scattering cross sections calculated using the TALYS reaction code. We compare the results of the microscopic calculations to phenomenological models and experimental data.

Motivated by the unknown nature of the $2.50-2.67\,M_\odot$ compact object in the binary merger event GW190814, we study the maximum neutron star mass based on constraints from low-energy nuclear physics, neutron star tidal deformabilities from GW170817, and simultaneous mass-radius measurements of PSR J0030+045 from NICER. Our prior distribution is based on a combination of nuclear modeling valid in the vicinity of normal nuclear densities together with the assumption of a maximally stiff equation of state at high densities, a choice that enables us to probe the connection between observed heavy neutron stars and the transition density at which conventional nuclear physics models must break down. We demonstrate that a modification of the highly uncertain supra-saturation density equation of state beyond 2.64 times normal nuclear density is required in order for chiral effective field theory models to be consistent with current neutron star observations and the existence of $2.6\,M_\odot$ neutron stars. We also show that the existence of very massive neutron stars strongly impacts the radii of $\sim 2.0\,M_\odot$ neutron stars (but not necessarily the radii of $1.4\,M_\odot$ neutron stars), which further motivates future NICER radius measurements of PSR J1614-2230 and PSR J0740+6620., Comment: 6 pages, 4 figures, revised figures

We construct a new class of phenomenological equations of state for homogeneous matter for use in simulations of hot and dense matter in local thermodynamic equilibrium. We construct a functional form which respects experimental, observational and theoretical constraints on the nature of matter in various density and temperature regimes. Our equation of state matches (i) the virial coefficients expected from nucleon-nucleon scattering phase shifts, (ii) experimental measurements of nuclear masses and charge radii, (iii) observations of neutron star radii, (iv) theory results on the equation of state of neutron matter near the saturation density, and (v) theory results on the evolution of the EOS at finite temperatures near the saturation density. Our analytical model allows one to compute the variation in the thermodynamic quantities based on the uncertainties in the nature of the nucleon-nucleon interaction. Finally, we perform a correction to ensure the equation of state is causal at all densities, temperatures, and electron fractions., Comment: 16 pages, 12 figures

Three-nucleon (3N) forces are an indispensable ingredient for accurate few-body and many-body nuclear structure and reaction theory calculations. While the direct implementation of chiral 3N forces can be technically very challenging, a simpler approach is given by employing instead a medium-dependent NN interaction V_med that reflects the physics of three-body forces at the two-body normal-ordered approximation. We review the derivation and construction of V_med from the chiral 3N interaction at next-to-next-to-leading order (N2LO), consisting of a long-range $2\pi$-exchange term, a mid-range $1\pi$-exchange component and a short-range contact term. Several applications of V_med to the equation of state of cold nuclear and neutron matter, the nucleon single-particle potential in nuclear matter, and the nuclear quasiparticle interaction are discussed. We also explore differences in using local vs. nonlocal regulating functions on 3N forces and make direct comparisons to exact results at low order in perturbation theory expansions for the equation of state and single-particle potential. We end with a discussion and numerical calculation of the in-medium NN potential V_med from the next-to-next-to-next-to-leading order (N3LO) chiral 3N force, which consists of a series of long-range and short-range terms., Comment: 31 pages, 15 figures, Published in Frontiers in Physics special volume "The Long-Lasting Quest for Nuclear Interactions: The Past, the Present and the Future"

We present predictions for neutron star tidal deformabilities obtained from a Bayesian analysis of the nuclear equation of state, assuming a minimal model at high-density that neglects the possibility of phase transitions. The Bayesian posterior probability distribution is constructed from priors obtained from microscopic many-body theory based on realistic two- and three-body nuclear forces, while the likelihood functions incorporate empirical information about the equation of state from nuclear experiments. The neutron star crust equation of state is constructed from the liquid drop model, and the core-crust transition density is found by comparing the energy per baryon in inhomogeneous matter and uniform nuclear matter. From the cold $\beta$-equilibrated neutron star equation of state, we then compute neutron star tidal deformabilities as well as the mass-radius relationship. Finally, we investigate correlations between the neutron star tidal deformability and properties of finite nuclei., Comment: 15 pages, 17 figures. Invited contribution to the EPJA Topical Issue "First joint gravitational wave and electromagnetic observations: Implications for nuclear and particle physics"

We formulate microscopic optical potentials for nucleon-nucleus scattering from chiral two- and three-nucleon forces. The real and imaginary central terms of the optical potentials are obtained from the nucleon self energy in infinite nuclear matter at a given density and isospin asymmetry, calculated self-consistently to second order in many-body perturbation theory. The real spin-orbit term is extracted from the same chiral potential using an improved density matrix expansion. The density-dependent optical potential is then folded with the nuclear density distributions of 40Ca, 42Ca, 44Ca, and 48Ca from which we study proton-nucleus elastic scattering and total reaction cross sections using the reaction code TALYS. We compare the results of the microscopic calculations to those of phenomenological models and experimental data up to projectile energies of E = 180 MeV. While overall satisfactory agreement with the available experimental data is obtained, we find that the elastic scattering and total reaction cross sections can be significantly improved with a weaker imaginary optical potential, particularly for larger projectile energies., Comment: 12 pages, 10 figures

We investigate neutron star moments of inertia from Bayesian posterior probability distributions of the nuclear equation of state that incorporate information from microscopic many-body theory and empirical data of finite nuclei. We focus on PSR J0737-3039A and predict that for this 1.338 M_sun neutron star the moment of inertia lies in the range $1.04 \times 10^{45}$ g cm$^{2} < I < 1.51 \times 10^{45}$ g cm$^{2}$ at the 95% credibility level, while the most probable value for the moment of inertia is $\tilde I = 1.36 \times 10^{45}$ g cm$^{2}$. Assuming a measurement of the PSR J0737-3039A moment of inertia to 10% precision, we study the implications for neutron star radii and tidal deformabilities. We also determine the crustal component of the moment of inertia and find that for typical neutron star masses of 1.3 M_sun < M < 1.5 M_sun the crust contributes 1% - 6% of the total moment of inertia, below what is needed to explain large pulsar glitches in the scenario of strong neutron entrainment., 12 pages, 12 figures, published in Phys. Rev. C

We confront observational data from gravitational wave event GW170817 with microscopic modeling of the cold neutron star equation of state. We develop and employ a Bayesian statistical framework that enables us to implement constraints on the equation of state from laboratory measurements of nuclei and state-of-the-art chiral effective field theory methods. The energy density functionals constructed from the posterior probability distributions are then used to compute consistently the neutron star equation of state from the outer crust to the inner core, assuming a composition consisting of protons, neutrons, electrons, and muons. In contrast to previous studies, we find that the 95% credibility range of predicted neutron star tidal deformabilities ($136 < \Lambda < 519$) for a 1.4 solar-mass neutron star is already consistent with the upper bound deduced from observations of the GW170817 event. However, we find that lower bounds on the neutron star tidal deformability will very strongly constrain microscopic models of the dense matter equation of state. We also demonstrate a strong correlation between the neutron star tidal deformability and the pressure of beta-equilibrated matter at twice saturation density., Comment: 6 pages, 6 figures, to appear in Physical Review Letters

We present calculations of shell-model effective interactions for both degenerate and non-degenerate model spaces using the Krenciglowa–Kuo (KK) and the extended Krenciglowa–Kuo iteration method recently developed by Okamoto, Suzuki et al. (EKKO). The starting point is the low-momentum nucleon–nucleon interaction V low − k obtained from the N3LO chiral two-nucleon interaction. The model spaces spanned by the sd and sdpf shells are both considered. With a solvable model, we show that both the KK and EKKO methods are convenient for deriving the effective interactions for non-degenerate model spaces. The EKKO method is especially desirable in this situation since the vertex function Z ˆ -box employed therein is well behaved while the corresponding vertex function Q ˆ -box employed in the Lee–Suzuki (LS) and KK methods may have singularities. The converged shell-model effective interactions given by the EKKO and KK methods are equivalent, although the former method is considerably more efficient. The degenerate sd-shell effective interactions given by the LS method are practically identical to those from the EKKO and KK methods. Results of the sd one-shell and sdpf two-shell calculations for 18O, 18F, 19O and 19F using the EKKO effective interactions are compared, and the importance of an inter-space consistency between the sd and sdpf calculations is discussed.

Core-polarization (CP) and Brown-Rho (BR) scaling were among Gerry's most favorite topics. In this contribution, we will discuss some of the early history as well as more recent work associated with these two fascinating phenomena. I (TTSK) will begin with some recollections of Princeton, where I met Gerry for the first time in 1964 and worked as his postdoc. Core polarization was in fact the first topic he assigned to me. JWH started working with Gerry at Stony Brook from 2003 and was Gerry's last student in nuclear physics. We three had teamed up well, working closely on both CP and BR scaling, particularly on the latter's connection to the anomalously-long beta-decay lifetime of carbon-14. We shall here briefly review these topics, including a recently developed new Brown-Rho scaling based on a Skyrmion half-Skyrmion two-phase model., Comment: 8 pages, 7 figures, published in Nuclear Physics A special issue in memory of Gerry Brown

We compute from chiral two- and three-nucleon interactions the energy per particle of symmetric nuclear matter and pure neutron matter at third-order in perturbation theory including self-consistent second-order single-particle energies. Particular attention is paid to the third-order particle-hole ring diagram, which is often neglected in microscopic calculations of the equation of state. We provide semianalytic expressions for the direct terms from central and tensor model-type interactions that are useful as theoretical benchmarks. We investigate uncertainties arising from the order-by-order convergence in both many-body perturbation theory and the chiral expansion. Including also variations in the resolution scale at which nuclear forces are resolved, we provide new error bands on the equation of state, the isospin-asymmetry energy, and its slope parameter. We find in particular that the inclusion of third-order diagrams reduces the theoretical uncertainty at low densities, while in general the largest error arises from omitted higher-order terms in the chiral expansion of the nuclear forces.

We study the ${}^{1}S_0$ proton pairing gap in beta-equilibrated neutron star matter within the framework of chiral effective field theory. We focus on the role of three-body forces, which strongly modify the effective proton-proton spin-singlet interaction in dense matter. We find that three-body forces generically reduce both the size of the pairing gap and the maximum density at which proton pairing may occur. The pairing gap is computed within BCS theory, and model uncertainties are estimated by varying the nuclear potential and the choice of single-particle spectrum in the gap equation. We find that a second-order perturbative treatment of the single-particle spectrum suppresses the proton ${}^{1}S_0$ pairing gap relative to the use of a free spectrum. We estimate the critical temperature for the onset of proton superconductivity to be $T_c = (3.7 - 6.0)\times 10^{9} $ K, which is consistent with previous theoretical results in the literature and marginally within the range deduced from a recent Bayesian analysis of neutron star cooling observations., Comment: 8 pages, 9 figures

We compute from chiral two- and three-body forces the complete quasiparticle interaction in symmetric nuclear matter up to twice nuclear matter saturation density. Second-order perturbative contributions that account for Pauli-blocking and medium polarization are included, allowing for an exploration of the full set of central and noncentral operator structures permitted by symmetries and the long-wavelength limit. At the Hartree-Fock level, the next-to-next-to-leading order three-nucleon force contributes to all noncentral interactions, and their strengths grow approximately linearly with the nucleon density up that of saturated nuclear matter. Three-body forces are shown to enhance the already strong proton-neutron effective tensor interaction, while the corresponding like-particle tensor force remains small. We also find a large isovector cross-vector interaction but small center-of-mass tensor interactions in the isoscalar and isovector channels. The convergence of the expansion of the noncentral quasiparticle interaction in Landau parameters and Legendre polynomials is studied in detail., Comment: 15 pages, 14 figures

We construct a new Skyrme interaction Sk$\chi$m$^*$ by fitting the equation of state and nucleon effective masses in asymmetric nuclear matter from chiral two- and three-body forces as well as the binding energies of finite nuclei. Employing this interaction to study the electric dipole polarizabilities of $^{48}$Ca, $^{68}$Ni, $^{120}$Sn, and $^{208}$Pb in the random-phase approximation, we find that the theoretical predictions are in good agreement with experimentally measured values without additional fine tuning of the Skyrme interaction, thus confirming the usefulness of the new Skyrme interaction in studying the properties of nuclei. We further use this interaction to study the neutron skin thicknesses of $^{48}$Ca and $^{208}$Pb, and they are found to be consistent with the experimental data., Comment: Significantly revised, 7 pages, 4 figures. Published version in PLB

The equations of state of nuclear and neutron matter (and, more generally, neutron-rich matter) play a fundamental role towards the understanding of a broad spectrum of systems, ranging from the skins of neutron-rich nuclei to the structure of compact stars. After a brief review of different approaches to predict the properties of nuclear and neutron matter, we will focus on error quantification in effective field theory. The various sources of uncertainty and their impact in infinite matter will be explored. We will then report on recent calculations of the nuclear and neutron matter equations of state at different orders of the chiral expansion as well as changing resolution scale. We will discuss the significance of such predictions as a foundation for future studies of convergence of the chiral perturbation series.

The isospin-asymmetry dependence of the nuclear matter equation of state obtained from microscopic chiral two- and three-body interactions in second-order many-body perturbation theory is examined in detail. The quadratic, quartic and sextic coefficients in the Maclaurin expansion of the free energy per particle of infinite homogeneous nuclear matter with respect to the isospin asymmetry are extracted numerically using finite differences, and the resulting polynomial isospin-asymmetry parametrizations are compared to the full isospin-asymmetry dependence of the free energy. It is found that in the low-temperature and high-density regime where the radius of convergence of the expansion is generically zero, the inclusion of higher-order terms beyond the leading quadratic approximation leads overall to a significantly poorer description of the isospin-asymmetry dependence. In contrast, at high temperatures and densities well below nuclear saturation density, the interaction contributions to the higher-order coefficients are negligible and the deviations from the quadratic approximation are predominantly from the noninteracting term in the many-body perturbation series. Furthermore, we extract the leading logarithmic term in the isospin-asymmetry expansion of the equation of state at zero temperature from the analysis of linear combinations of finite differences. It is shown that the logarithmic term leads to a considerably improved description of the isospin-asymmetry dependence at zero temperature., 14 pages, 9 figures, 2 tables, some minor changes, references updated, matches published version

The notion of "hadronic freedom" is introduced based on the vector manifestation of hidden local symmetry and is used to suggest that the dileptons measured in relativistic heavy-ion collisions do not provide {\em direct} information on the spontaneous breaking of chiral symmetry and hence on the mechanism for mass generation of light-quark hadrons. We give arguments how the dileptons emitted from those vector mesons whose masses are shifted \`a la Brown-Rho (BR) scaling by the vacuum change in temperature -- as in heavy-ion collisions -- and/or density -- as in cold compressed matter -- could be strongly suppressed in the hadronic free region between the chiral restoration point and the "flash point" at which the vector mesons recover $\gsim 90%$ of their free-space on-shell masses and the full strong coupling strength. It may seem ironical that the very mechanism i.e., the vector manifestation in hidden local symmetry, that is to make the mass drop and modify the spectral function is in turn responsible for the dilepton suppression. A possible falsification of this drastic prediction will be indicated. We also briefly discuss the potential role played by the holographic dimension intrinsic in gravity-gauge duality that provides a unified field theory description of hadrons -- both mesons and baryons -- under extreme conditions. Baryons arise as coherent states of pions and vector mesons of holographic QCD which define the ground state or the "vacuum" on which mesonic excitations could undergo BR scaling., Comment: To revamp and replace arXiv:0804.3196 which was withdrawn due to nontrivial errors

Background: An accurate description of nuclear pairing gaps is extremely important for understanding static and dynamic properties of the inner crusts of neutron stars and to explain their cooling process. Purpose: We plan to study the behavior of the pairing gaps $\Delta_F$ as a function of the Fermi momentum $k_F$ for neutron and nuclear matter in all relevant angular momentum channels where superfluidity is believed to naturally emerge. The calculations will employ realistic chiral nucleon-nucleon potentials with the inclusion of three-body forces and self-energy effects. Methods: The superfluid states of neutron and nuclear matter are studied by solving the BCS gap equation for chiral nuclear potentials using the method suggested by Khodel et al., where the original gap equation is replaced by a coupled set of equations for the dimensionless gap function $\chi(p)$ defined by $\Delta(p) = \Delta_F \chi(p)$ and a non-linear algebraic equation for the gap magnitude $\Delta_F = \Delta(p_F)$ at the Fermi surface. This method is numerically stable even for small pairing gaps, such as that encountered in the coupled $^3PF_2$ partial wave. Results: We have successfully applied Khodel's method to singlet ($S$) and coupled channel ($SD$ and $PF$) cases in neutron and nuclear matter. Our calculations agree with other ab-initio approaches, where available, and provide crucial inputs for future applications in superfluid systems., Comment: 18 pages and 9 figures

Chiral symmetry, first entering in nuclear physics in the 1970's for which Gerry Brown played a seminal role, has led to a stunningly successful framework for describing strongly-correlated nuclear dynamics both in finite and infinite systems. We review how the early germinal idea, conceived with the soft-pion theorems in the pre-QCD era, has evolved into a highly predictive theoretical framework for nuclear physics, aptly assessed by Steven Weinberg: "it (chiral effective field theory) allows one to show in a fairly convincing way that what they (nuclear physicists) have been doing all along... is the correct first step in a consistent approximation scheme." Our review recounts both how the theory presently fares in confronting Nature and how one can understand its extremely intricate workings in terms of the multifaceted aspects of chiral symmetry, namely, chiral perturbation theory, skyrmions, Landau Fermi-liquid theory, the Cheshire cat phenomenon, and hidden local and mended symmetries., 108 pages, 50 figures, submitted to Physics Reports special volume dedicated to Gerry Brown

Realistic nuclear potentials, derived within chiral perturbation theory, are a major breakthrough in modern nuclear structure theory, since they provide a direct link between nuclear physics and its underlying theory, namely the QCD. As a matter of fact, chiral potentials are tailored on the low-energy regime of nuclear structure physics, and chiral perturbation theory provides on the same footing two-nucleon forces as well as many-body ones. This feature fits well with modern advances in ab-initio methods and realistic shell-model. Here, we will review recent nuclear structure calculations, based on realistic chiral potentials, for both finite nuclei and infinite nuclear matter., 10 pages, 8 figures, plenary talk presented at "Nucleus-Nucleus 2015" Conference, 21-26 June 2015, Catania, to be published in the "Conference Proceedings" Series of the Italian Physical Society

We present a description of symmetric nuclear matter within the framework of Landau Fermi liquid theory. The low momentum nucleon-nucleon interaction V(low-k) is used to calculate the effective interaction between quasiparticles on the Fermi surface, from which we extract the quasiparticle effective mass, the nuclear compression modulus, the symmetry energy, and the anomalous orbital gyromagnetic ratio. The exchange of density, spin, and isospin collective excitations is included through the Babu-Brown induced interaction, and it is found that in the absence of three-body forces the self-consistent solution to the Babu-Brown equations is in poor agreement with the empirical values for the nuclear observables. This is improved by lowering the nucleon and meson masses according to Brown-Rho scaling, essentially by including a scalar tadpole contribution to the meson and nucleon masses, as well as by scaling g_A. We suggest that modifying the masses of the exchanged mesons is equivalent to introducing a short-range three-body force, and the net result is that the Brown-Rho double decimation is accomplished all at once., 26 pages, 9 figures

Recent developments in our description of RHIC and related heavy-ion phenomena in terms of hidden local symmetry theories are reviewed with a focus on the novel nearly massless states in the vicinity of -- both below and above -- the chiral restoration temperature T_c. We present complementary and intuitive ways to understand both Harada-Yamawaki's vector manifestation structure and Brown-Rho scaling -- which are closely related -- in terms of "melting" of soft glues observed in lattice calculations and join the massless modes that arise in the vector manifestation (in the chiral limit) just below T_c to tightly bound massless states above T_c. This phenomenon may be interpreted in terms of the Beg-Shei theorem. It is suggested that hidden local symmetry theories arise naturally in holographic dual QCD from string theory, and a clear understanding of what really happens near the critical point could come from a deeper understanding of the dual bulk theory. Other matters discussed are the relation between Brown-Rho scaling and Landau Fermi-liquid fixed point parameters at the equilibrium density, its implications for "low-mass dileptons" produced in heavy-ion collisions, the reconstruction of vector mesons in peripheral collisions, the pion velocity in the vicinity of the chiral transition point, kaon condensation viewed from the VM fixed point, nuclear physics with Brown-Rho scaling, and the generic feature of dropping masses at the RGE fixed points in generalized hidden local symmetry theories.

We explore the use of mean field models to approximate microscopic nuclear equations of state derived from chiral effective field theory across the densities and temperatures relevant for simu- lating astrophysical phenomena such as core-collapse supernovae and binary neutron star mergers. We consider both relativistic mean field theory with scalar and vector meson exchange as well as energy density functionals based on Skyrme phenomenology and compare to thermodynamic equa- tions of state derived from chiral two- and three-nucleon forces in many-body perturbation theory. Quantum Monte Carlo simulations of symmetric nuclear matter and pure neutron matter are used to determine the density regimes in which perturbation theory with chiral nuclear forces is valid. Within the theoretical uncertainties associated with the many-body methods, we find that select mean field models describe well microscopic nuclear thermodynamics. As an additional consistency requirement, we study as well the single-particle properties of nucleons in a hot/dense environment, which affect e.g., charged-current weak reactions in neutron-rich matter. The identified mean field models can be used across a larger range of densities and temperatures in astrophysical simulations than more computationally expensive microscopic models.

The density and temperature dependence of the nuclear symmetry free energy is investigated using microscopic two- and three-body nuclear potentials constructed from chiral effective field theory. The nuclear force models and many-body methods are benchmarked to properties of isospin-symmetric nuclear matter in the vicinity of the saturation density as well as the virial expansion of the neutron matter equation of state at low fugacities. The free energy per particle of isospin-asymmetric nuclear matter is calculated assuming a quadratic dependence of the interaction contributions on the isospin asymmetry. The spinodal instability at subnuclear densities is examined in detail., Comment: 20 pages, 15 figures, 2 tables. additional typos corrected, cf. footnote 2 in arXiv:1603.02935

We calculate the nuclear and neutron matter equations of state from microscopic nuclear forces at different orders in chiral effective field theory and with varying momentum-space cutoff scales. We focus attention on how the order-by-order convergence depends on the choice of resolution scale and the implications for theoretical uncertainty estimates on the isospin asymmetry energy. Specifically we study the equations of state using consistent NLO and N2LO (next-to-next-to-leading order) chiral potentials where the low-energy constants cD and cE associated with contact vertices in the N2LO chiral three-nucleon force are fitted to reproduce the binding energies of 3H and 3He as well as the beta-decay lifetime of 3H. At these low orders in the chiral expansion there is little sign of convergence, while an exploratory study employing the N3LO two-nucleon force together with the N2LO three-nucleon force give first indications for (slow) convergence with low-cutoff potentials and poor convergence with higher-cutoff potentials. The consistent NLO and N2LO potentials described in the present work provide the basis for estimating theoretical uncertainties associated with the order-by-order convergence of nuclear many-body calculations in chiral effective field theory., Comment: 9 pages, 6 figures. Accepted for publication in Physical Review C

The Kirson-Babu-Brown (KBB) induced interaction approach is applied to the derivation of the sd-shell nucleon effective interactions. We start from a low-momentum NN interaction Vlow-k obtained by integrating out the high-momentum components of modern realistic NN potentials beyond a decimation momentum λ. Then a large class of core polarization diagrams are summed up to all orders by way of the KBB self-consistent equations. It is found that the solution of these equations is simplified by the use of Vlow-k, which is energy independent, and by treating them in terms of the Green's fuctions in the particle-particle and the particle-hole channels. The effective interactions calculated with the all-order KBB core polarizations and those with the polarization in second-order perturbation theory are found to be rather similar to each other, typical difference between them being less than 10%.

We compute the isospin-asymmetry dependence of microscopic optical model potentials from realistic chiral two- and three-body interactions over a range of resolution scales $\Lambda \simeq 400-500$\,MeV. We show that at moderate projectile energies, $E_{\rm inv} = 110 - 200$\,MeV, the real isovector part of the optical potential changes sign, a phenomenon referred to as isospin inversion. We also extract the strength and energy dependence of the imaginary isovector optical potential and find no evidence for an analogous phenomenon over the range of energies, $E \leq 200$\,MeV, considered in the present work. Finally, we compute for the first time the leading corrections to the Lane parametrization for the isospin-asymmetry dependence of the optical potential and observe an enhanced importance at low scattering energies., Comment: 8 pages, 7 figures

International audience; We study the impact of three-body forces on the response functions of cold neutron matter. These response functions are determined in the random phase approximation (RPA) from a residual interaction expressed in terms of Landau parameters. Special attention is paid to the non-central part, including all terms allowed by the relevant symmetries. Using Landau parameters derived from realistic nuclear two- and three-body forces grounded in chiral effective field theory, we find that the three-body term has a strong impact on the excited states of the system and in the static and long-wavelength limit of the response functions for which a new exact formula is established.

We investigate the thermodynamic equation of state of isospin-symmetric nuclear matter with microscopic nuclear forces derived within the framework of chiral effective field theory. Two- and three-body nuclear interactions constructed at low resolution scales form the basis for a perturbative calculation of the finite-temperature equation of state. The nuclear force models and many-body methods are benchmarked against bulk properties of isospin-symmetric nuclear matter at zero temperature, which are found to be well reproduced when chiral nuclear interactions constructed at the lowest resolution scales are employed. The calculations are then extended to finite temperatures, where we focus on the liquid-gas phase transition and the associated critical point. The Maxwell construction is applied to construct the physical equation of state, and the value of the critical temperature is determined to be T_c =17.2-19.1 MeV, in good agreement with the value extracted from multifragmentation reactions of heavy ions., 22 pages, 12 figures, 2 tables. v3 matches published version

We calculate perturbatively the energy per nucleon in infinite nuclear matter with a chiral N3LO (next-to-next-to-next-to-leading order) two-body potential plus a N2LO three-body force (3BF). The 3BF low-energy constants which cannot be constrained by two-body observables are chosen such as to reproduce the A=3 binding energies and the triton Gamow-Teller matrix element. This enables to study the nuclear matter equation of state in a parameter-free approach., 7 pages, 5 figures, talk presented at "XIV Convegno su Problemi di Fisica Nucleare Teorica", Cortona October 29-31, 2013. Submitted to Journal of Physics: Conferences Series

In a recent paper we studied the behavior of the pairing gaps $\Delta_F$ as a function of the Fermi momentum $k_F$ for neutron and nuclear matter in all relevant angular momentum channels where superfluidity is believed to naturally emerge. The calculations employed realistic chiral nucleon-nucleon potentials with the inclusion of three-body forces and self-energy effects. In this contribution, after a detailed description of the numerical method we employed in the solution of the BCS equations, we will show a preliminary analysis of the Cooper pair wavefunctions., Comment: 9 pages, 4 figures, contribution to the Proceedings of ICNFP2014 and PANIC14

We calculate neutrino absorption rates due to charged-current reactions ${\ensuremath{\nu}}_{e}+n\ensuremath{\rightarrow}{e}^{\ensuremath{-}}+p$ and ${\overline{\ensuremath{\nu}}}_{e}+p\ensuremath{\rightarrow}{e}^{+}+n$ in the outer regions of a newly born neutron star called the neutrino-sphere. To improve on recent work which has shown that nuclear mean fields enhance the ${\ensuremath{\nu}}_{e}$ cross section and suppress the ${\overline{\ensuremath{\nu}}}_{e}$ cross section, we employ realistic nucleon-nucleon interactions that fit measured scattering phase shifts. Using these interactions we calculate the momentum-, density-, and temperature-dependent nucleon self-energies in the Hartree-Fock approximation. A potential derived from chiral effective field theory and a pseudopotential constructed to reproduce nucleon-nucleon phase shifts at the mean-field level are used to study the equilibrium proton fraction and charged-current rates. We compare our results to earlier calculations obtained using phenomenological mean-field models and to those obtained in the virial expansion valid at low density and high temperature. In the virial regime our results are consistent with previous calculations, and at higher densities relevant for the neutrino-sphere, $\ensuremath{\rho}\ensuremath{\gtrsim}{10}^{12}$ ${\mathrm{g}/\mathrm{cm}}^{3}$, we find the difference between the ${\ensuremath{\nu}}_{e}$ and ${\overline{\ensuremath{\nu}}}_{e}$ absorption rates to be larger than predicted earlier. Our results may have implications for heavy-element nucleosynthesis in supernovae, and for supernova neutrino detection.

We present variational Monte Carlo calculations of the neutron matter equation of state using chiral nuclear forces. The ground-state wavefunction of neutron matter, containing non-perturbative many-body correlations, is obtained from auxiliary-field quantum Monte Carlo simulations of up to about 340 neutrons interacting on a 10^3 discretized lattice. The evolution Hamiltonian is chosen to be attractive and spin-independent in order to avoid the fermion sign problem and is constructed to best reproduce broad features of the chiral nuclear force. This is facilitated by choosing a lattice spacing of 1.5 fm, corresponding to a momentum-space cutoff of Lambda = 414 MeV/c, a resolution scale at which strongly repulsive features of nuclear two-body forces are suppressed. Differences between the evolution potential and the full chiral nuclear interaction (Entem and Machleidt Lambda = 414 MeV) are then treated perturbatively. Our results for the equation of state are compared to previous quantum Monte Carlo simulations which employed chiral two-body forces at next-to-next-to-leading order (N2LO). In addition we include the effects of three-body forces at N2LO, which provide important repulsion at densities higher than 0.02 fm^-3, as well as two-body forces at N3LO.

We compute the energy per particle of infinite symmetric nuclear matter from chiral N3LO (next-to-next-to-next-to-leading order) two-body potentials plus N2LO three-body forces. The low-energy constants of the chiral three-nucleon force that cannot be constrained by two-body observables are fitted to reproduce the triton binding energy and the 3H-3He Gamow-Teller transition matrix element. In this way, the saturation properties of nuclear matter are reproduced in a parameter-free approach. The equation of state is computed up to third order in many-body perturbation theory, with special emphasis on the role of the third-order particle-hole diagram. The dependence of these results on the cutoff scale and regulator function is studied. We find that the inclusion of three-nucleon forces consistent with the applied two-nucleon interaction leads to a reduced dependence on the choice of the regulator only for lower values of the cutoff., Comment: 9 pages, 12 figures, 3 tables, to be published in Physical Review C. arXiv admin note: text overlap with arXiv:1209.5537

We review efforts to describe the approach to chiral symmetry restoration in neutron matter from the low-energy realization of QCD, chiral effective field theory., Comment: Contribution to the Gerry Brown 90th birthday memorial book published by World Scientific. v2 updated references

We calculate the energy per particle in infinite neutron matter perturbatively using chiral N3LO two-body potentials plus N2LO three-body forces. The cutoff dependence of the predictions is investigated by employing chiral interactions with different regulators. We find that the inclusion of three-nucleon forces, which are consistent with the applied two-nucleon interaction, leads to a strongly reduced regulator dependence of the results.