Search

The reduced basis method is used to construct a "universal" basis of Dirac orbitals that may be applicable throughout the nuclear chart to calibrate covariant energy density functionals. Relative to our earlier work using the non-relativistic Schr\"odinger equation, the Dirac equation adds an extra layer of complexity due to the existence of negative energy states. However, once this problem is mitigated, the resulting reduced basis is able to accurately and efficiently reproduce the high-fidelity model at a fraction of the computational cost. We are confident that the resulting reduced basis will serve as a foundational element in developing rapid and accurate emulators. In turn, these emulators will play a critical role in the Bayesian optimization of covariant energy density functionals., Comment: 9 pages, 7 figures, to be submitted to Physical Review C

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the symmetry energy and its slope parameter at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Our Bayesian framework is publicly available, so practitioners can readily use and extend our results., Comment: close to the published version; extended analysis and minor changes; 31 pages, 14 figures, 5 tables

The historical detection of gravitational waves from the binary neutron star merger GW170817 is providing fundamental new insights into the astrophysical site for the creation of the heaviest elements in the cosmos and on the equation of state of neutron-rich matter. Shortly after this historical detection, electromagnetic observations of neutron stars together with measurements of the properties of neutron-rich nuclei at terrestrial facilities have placed additional constraints on the dynamics of neutron-rich matter. It is this unique synergy between heaven and earth that is the focus of this article., Comment: 9 pages and 6 figures. Manuscript submitted to the proceedings of the XLV Symposium on Nuclear Physics in Cocoyoc, Mexico

The question of how and which physical observables or thermodynamic parameters can best predict the onset of a possible phase transition in the inner crust of neutron stars remains largely unresolved. Using semiclassical Monte Carlo simulations, we investigate the isothermal compressibility and density fluctuations in a region of relevance to the dynamics of the inner crust. We show that the isothermal compressibility serves as a robust observable to characterize the transition from the non-uniform crust to the uniform core for proton fractions over 0.2. Moreover, we show explicitly how the two-component isothermal compressibility, computed using the Kirkwood-Buff theory, is directly connected to the fluctuations in the number density, recorded in the grand canonical ensemble by monitoring the number of particles in a small volume located at the center of the simulation box. That is, we compute mean-square particle fluctuations and compare them against the isothermal compressibility for different proton fractions. Although our results show that the mean-square particle fluctuations are proportional to the isothermal compressibility, the lack of a perfect correlation is attributed to the relatively small number of particles included in the simulations. The non-unity slope observed in the dimensionless isothermal compressibility-total nucleon fluctuation variance relationship suggests that the inner crust of neutron stars is composed of anisotropic and inhomogeneous matter., Comment: 10 pages, 5 figures

The recently published CREX results suggest a rather peculiar picture for the density dependence of the symmetry energy. Whereas PREX favors a large neutron skin thickness in $^{208}$Pb, thereby suggesting a stiff equation of state, CREX suggests instead a much softer equation of state. This discrepancy has caused a large spur in the theoretical community since no model has been able to simultaneously reproduce within $1\sigma$ the PREX and CREX results. Motivated by a novel correlation between a CREX observable and a combination of bulk symmetry energy parameters, we calibrate three new covariant energy density functionals that reproduce binding energies and charge radii of spherical nuclei - and also accommodate the constraints imposed by PREX and CREX. Given that these models suggest a stiff equation of state at high densities, predictions for neutron star properties are also discussed., Comment: 10 pages, resubmitted to PRC

The last five years have seen remarkable progress in our quest to determine the equation of state of neutron rich matter. Recent advances across the theoretical, experimental, and observational landscape have been incorporated in a Bayesian framework to refine existing covariant energy density functionals previously calibrated by the properties of finite nuclei. In particular, constraints on the maximum neutron star mass from pulsar timing, on stellar radii from the NICER mission, on tidal deformabilities from the LIGO-Virgo collaboration, and on the dynamics of pure neutron matter as predicted from chiral effective field theories, have resulted in significant refinements to the models, particularly to those predicting a stiff symmetry energy. Still, even after these improvements, we find challenging to reproduce simultaneously the neutron skin thickness of both ${}^{208}$Pb and ${}^{48}$Ca recently reported by the PREX/CREX collaboration., Comment: 16 pages, 7 Figures, 3 Tables

We study the detectability of gravitational-wave signals from sub-solar mass binary neutron star systems by the current generation of ground-based gravitational-wave detectors. We find that finite size effects from large tidal deformabilities of the neutron stars and lower merger frequencies can significantly impact the sensitivity of the detectors to these sources. By simulating a matched-filter based search using injected binary neutron star signals with tidal deformabilities derived from physically motivated equations of state, we calculate the reduction in sensitivity of the detectors. We conclude that the loss in sensitive volume can be as high as $78.4 \%$ for an equal mass binary system of chirp mass $0.17 \, \textrm{M}_{\odot}$, in a search conducted using binary black hole template banks. We use this loss in sensitive volume, in combination with the results from the search for sub-solar mass binaries conducted on data collected by the LIGO-Virgo observatories during their first three observing runs, to obtain a conservative upper limit on the merger rate of sub-solar mass binary neutron stars. Since the discovery of a low-mass neutron star would provide new insight into formation mechanisms of neutron stars and further constrain the equation of state of dense nuclear matter, our result merits a dedicated search for sub-solar mass binary neutron star signals., Comment: 12 pages, 7 figures, supplemental materials at https://github.com/sugwg/sub-solar-ns-detectability

Neutron stars -- compact objects with masses similar to that of our Sun but radii comparable to the size of a city -- contain the densest form of matter in the universe that can be probed in terrestrial laboratories as well as in earth- and space-based observatories. The historical detection of gravitational waves from a binary neutron star merger has opened the brand new era of multimessenger astronomy and has propelled neutron stars to the center of a variety of disciplines, such as astrophysics, general relativity, nuclear physics, and particle physics. The main input required to study the structure of neutron stars is the pressure support generated by its constituents against gravitational collapse. These include neutrons, protons, electrons, and perhaps even more exotic constituents. As such, nuclear physics plays a prominent role in elucidating the fascinating structure, dynamics, and composition of neutron stars., Comment: 15 pages, 6 figures, submitted to Oxford Research Encyclopedia of Physics

Reduced basis methods provide a powerful framework for building efficient and accurate emulators. Although widely applied in many fields to simplify complex models, reduced basis methods have only been recently introduced into nuclear physics. In this letter we build an emulator to study the single-particle structure of atomic nuclei. By scaling a suitable mean-field Hamiltonian, a "universal" reduced basis is constructed capable of accurately and efficiently reproduce the entire single-particle spectrum of a variety of nuclei. Indeed, the reduced basis model reproduces both ground- and excited-state energies as well as the associated wave-functions with remarkable accuracy. Our results bode well for more demanding applications that use Bayesian optimization to calibrate nuclear energy density functionals., Comment: 5 pages, 5 figures, to be submitted to Physical Review C

Background: Coherent {\pi}0 photoproduction on heavy nuclei has been suggested as a reliable tool to infer neutron skin thicknesses. To this aim, various experiments have been performed, especially on 208Pb. Purpose: We analyze the sensitivity of that reaction process to the nucleonic density, and especially to the neutron skin thickness, for 12C, 40Ca and 208Pb, for which reliable data exist, and on 116,124Sn, for which measurements have been performed in Mainz. We study also the role played by the {\pi}0-nucleus final-state interaction. Method: A model of the reaction is developed at the impulse approximation considering either plane waves or distorted waves to describe the {\pi}0-nucleus scattering in the outgoing channel. Results: Our calculations are in good agreement with existing data, especially for 208Pb. The sensitivity of the theoretical cross sections to the choice of the nucleonic density is small, and below the experimental resolution. Conclusions: Coherent {\pi}0 photoproduction is mostly an isoscalar observable that bares no practical sensitivity to the neutron skin thickness. To infer that structure observable it should be coupled to other reaction measurements, such as electron scattering, or by comparing experiments performed on isotopes of the same chemical element., Comment: Accepted for publication in Phys. Rev. C. One new Ref.[21]

The historical detection of gravitational waves emitted from the binary neutron star merger GW170817 has opened the new era of multi-messenger astronomy. Since then, many other significant discoveries -- both on heaven and earth -- are providing new clues into the behavior of neutron-rich matter. It is the goal of this article to illustrate how the remarkable progress made during the last few years is spearheading the field into the golden age of neutron-star physics., Comment: 9 pages, 5 figures, submitted to the proceedings of the XLIV Brazilian Workshop on Nuclear Physics

With an extreme neutron-to-proton ratio of N/Z=3, 8He provides an ideal laboratory for the study of a variety of exotic phenomena, such as the emergence of a soft dipole mode that is dominated by transitions into the continuum. In this contribution, a covariant density functional theory (DFT) framework is used to compute ground-state properties and the dipole response of 8He. Although 8He is admittedly too light for DFT to be applicable, the great merit of the approach is that the spurious contamination associated with the center-of-mass motion is guaranteed to decouple from the physical response. Given that a strong mixing between the isoscalar and isovector dipole modes is expected for a system with such a large neutron-proton asymmetry as 8He, the narrow structures that emerged at low energies in the isovector dipole response are attributed to the shift of the spurious strength to zero (or near zero) excitation energy. Thus, the theoretical framework implemented here disfavors the emergence of a soft dipole mode in 8He., Comment: 8 pages, 7 figures, and two tables. Submitted to the Physical Review C

The poorly known crustal equation of state plays a critical role in many observational phenomena associated with a neutron star. Using semi-classical Monte Carlo simulations, we explore the possible configurations of the inner edge of the neutron-star crust for a variety of baryon densities and proton fractions. Applying the Kirkwood--Buff theory to these two-component systems, we observe how the isothermal compressibility reaches a maximum when isolated non-symmetric clusters are formed in an extremely dilute neutron gas. To determine the neutron fraction, we suggest a geometrical model based on the behavior of the proton-neutron pair correlation function. Accordingly, the equation of state of the inner crust is calculated, illustrating that the nuclear energy in beta-equilibrium follows a power-law behavior with baryon density. As a possible astrophysical outcome of this study, our results could help refine the mass-radius relation. Finally, our results pave the way towards further investigations of the impact of the proton-neutron pair correlation function on transport properties within the neutron-star crust., Comment: 11 pages, 6 figures

The recent announcement by the PREX collaboration of an unanticipated thick neutron skin in 208Pb (Rskin208) has challenged our understanding of neutron rich matter in the vicinity of nuclear saturation density. Whereas earlier constraints indicate that the symmetry energy is relatively soft, the PREX-2 result seems to suggest the opposite. We confront constraints on the symmetry energy obtained from measurements of the electric dipole polarizability against those informed by the PREX-2 measurement of Rskin208 and by the correlations that it entails. Covariant energy density functionals informed by the properties of finite nuclei are used to compute the electric dipole response of 48Ca, 68Ni, 132Sn, and 208Pb. The set of functionals used in this work are consistent with experimental data, yet are flexible enough in that they span a wide range of values of Rskin208. It is found that theoretical predictions of the electric dipole polarizability that are consistent with the PREX-2 measurement systematically overestimate the corresponding values extracted from the direct measurements of the distribution of electric dipole strength. The neutron skin thickness of 208Pb extracted from parity violating electron scattering and the electric dipole polarizability measured in photoabsorption experiments are two of the cleanest experimental tools used to constrain the symmetry energy around nuclear saturation density. However, the recent value of Rskin208 that suggests a fairly stiff symmetry energy stands in stark contrast to the conclusions derived from the electric dipole polarizability. At present, we offer no solution to this dilemma., Comment: 11 pages and 9 figures

Insights into the equation of state of neutron rich matter obtained from the neutron skin thickness of 208Pb are in sharp conflict with earlier measurements of the electric dipole polarizability. We use a set of accurately calibrated energy density functionals to highlight the tension and to articulate how plans for a highly-intense Gamma Factory may alleviate the tension., Comment: 11 pages and 3 figures. Submitted to Annals of Physics as a contribution to the Special Issue on Physics Opportunities with the Gamma Factory

Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a "density ladder" that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of 208Pb (Rskin) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In turn, an improved value of Rskin has been reported recently by the PREX collaboration. Exploiting the strong correlation between Rskin and the slope of the symmetry energy L within a specific class of relativistic energy density functionals, we report a value of L=(106 +/- 37)MeV -- that systematically overestimates current limits based on both theoretical approaches and experimental measurements. The impact of such a stiff symmetry energy on some critical neutron-star observables is also examined., Comment: 5 pages, 4 figures. Manuscript accepted for publication in the Physical Review Letters

The saturation density of nuclear matter $\rho_0$ is a fundamental nuclear physics property that is difficult to predict from fundamental principles. The saturation density is closely related to the interior density of a heavy nucleus, such as $^{208}$Pb. We use parity violating electron scattering to determine the average interior weak charge and baryon densities in $^{208}$Pb. This requires not only measuring the weak radius $R_{\rm wk}$ but also determining the surface thickness of the weak charge density $a$. We obtain $\rho_0=0.150\pm0.010$ fm$^{-3}$, where the 7\% error has contributions form the PREX error on the weak radius, an assumed 10\% uncertainty in the surface thickness $a$, and from the extrapolation to infinite nuclear matter. These errors can be improved with the upcoming PREX II results and with a new parity violating electron scattering experiment, at a somewhat higher momentum transfer, to determine $a$., Comment: 6 pages, 3 figures, minor changes, Phys. Rev. C in press

Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system [R. Abbott et al., ApJ Lett., 896, L44 (2020)]? This is the central question animating this letter. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (a) account for a 2.6 Msun neutron star, (b) satisfy the original constraint on the tidal deformability of a 1.4 Msun neutron star, and (c) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support super-massive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Thus, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 Msun compact object is likely to be the lightest black hole ever discovered., Comment: 6 pages and 3 figures

Uncertainty quantification has emerged as a rapidly growing field in nuclear science. Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of limited value. Also, one often deals with theoretical constructs that involve fundamental quantities that are not accessible to experiment or observation. Particularly relevant in this context is the pressure of pure neutron matter. In this contribution we develop an analytic framework to answer the question of "How can new data reduce uncertainties of current theoretical models?" [P.-G. Reinhard and W. Nazarewicz, Phys. Rev. C81, 051303(R) (2010)]. Simple and insightful expressions are obtained to quantify the impact of one or two new observables on theoretical uncertainties in two critical quantities: the slope of the symmetry energy at saturation density and the pressure of pure neutron matter at twice nuclear matter saturation density., Comment: 6 pages and 2 figures

"Why are the tin isotopes soft?" has remained, for the past decade, an open problem in nuclear structure physics: models which reproduce the isoscalar giant monopole resonance (ISGMR) in the "doubly-closed shell" nuclei, $^{90}$Zr and $^{208}$Pb, overestimate the ISGMR energies of the open-shell tin and cadmium nuclei, by as much as 1 MeV. In an effort to shed some light onto this problem, we present results of detailed studies of the ISGMR in the molybdenum nuclei, with the goal of elucidating where--and how--the softness manifests itself between $^{90}$Zr and the cadmium and tin isotopes. The experiment was conducted using the $^{94,96,98,100}$Mo($\alpha,\alpha^\prime$) reaction at $E_\alpha = 386$ MeV. A comparison of the results with relativistic, self-consistent Random-Phase Approximation calculations indicates that the ISGMR response begins to show softness in the molybdenum isotopes beginning with $A=92$., Comment: Accepted for publication to Physics Letters B

Background: The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce---or rather re-introduced---the "Dirac Oscillator", a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. Purpose: To assess---to our knowledge for the first time---the power and flexibility of the Dirac Oscillator basis in the solution of nuclear structure problems within the framework of covariant density functional theory. Methods: Self-consistent calculations of binding energies and ground-state densities for a selected set of doubly-magic magic are performed using the Dirac-oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results: Results obtained using the Dirac oscillator basis reproduced with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Conclusions: Although the three-dimensional harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry as required in constrained calculations of nuclear excitations., Comment: 10 pages, 2 figures

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. Yet, a coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review., Comment: 23 pages, 5 figures. When citing this paper, please use the following: Yang J, Piekarewicz J. Covariant Density Functional Theory in Nuclear Physics and Astrophysics. Annual Review of Nuclear and Particle Science Volume 70: Submitted. DOI: 10.1146/annurev-nucl-101918-023608

Elastic electron scattering has been used to paint the most accurate picture of the proton distribution in atomic nuclei. Spurred by new experimental developments, it is now possible to gain valuable insights into the neutron distribution using exclusively electroweak probes. Our goal is to assess the information content and complementarity of the following three electroweak experiments in constraining the neutron distribution of atomic nuclei: (a) parity violating elastic electron scattering, (b) coherent elastic neutrino nucleus scattering, and (c) elastic electron scattering of unstable nuclei. Relativistic mean-field models informed by the properties of finite nuclei and neutron stars are used to compute ground state densities and form factors of a variety of nuclei. All the models follow the same fitting protocol, except for the assumed and presently unknown value of the neutron skin thickness of 208Pb. This enables one to tune the density dependence of the symmetry energy without compromising the success in reproducing well known physical observables. We found that the ongoing PREX-II and upcoming CREX campaigns at Jefferson Lab will play a vital role in constraining the weak form factor of xenon and argon, liquid noble gases that are used for the detection of both neutrinos and dark matter particles. We concluded that remarkable new advances in experimental physics have opened a new window into ground state densities of atomic nuclei using solely electroweak probes. The diversity and versatility of these experiments reveal powerful correlations that impose important nuclear structure constraints. In turn, these constraints provide quantitative theoretical uncertainties that are instrumental in searches for new physics and insights into the behavior of dense matter., Comment: 18 pages, 11 figures

Despite a length-scale difference of 18 orders of magnitude, the internal structure of neutron stars and the spatial distribution of neutrons in atomic nuclei are profoundly connected., Comment: Manuscript published in final form in Physics Today, July 2019, page 30

The complex nature of the nuclear forces generates a broad range and diversity of observational phenomena. Heavy nuclei, though orders of magnitude less massive than neutron stars, are governed by the same underlying physics, which is enshrined in the nuclear equation of state. Heavy nuclei are expected to develop a neutron-rich skin where many neutrons collect near the surface. Such a skin thickness is strongly sensitive to the poorly-known density dependence of the symmetry energy near saturation density. An accurate and model-independent determination of the neutron-skin thickness of heavy nuclei would provide a significant first constraint on the density dependence of the nuclear symmetry energy. The determination of the neutron-skin thickness of heavy nuclei has far reaching consequences in many areas of physics as diverse as heavy-ion collisions, polarized electron and proton scattering off nuclei, precision tests of the standard model using atomic parity violation, and nuclear astrophysics. While a systematic and concerted experimental effort has been made to measure the neutron-skin thickness of heavy nuclei, a precise and model-independent determination remains elusive. How to move forward at a time when many new facilities are being commissioned and how to strengthen the synergy with other areas of physics are primary goals of this review., Comment: Submitted for publication to Journal of Physics G

The first detection of a binary neutron star merger has opened the brand new era of multimessenger astronomy. This historic detection has been instrumental in providing constraints on the tidal polarizability of neutron stars. In turn, the tidal polarizability has been used to impose limits on stellar radii and ultimately on the equation of state (EOS). The tidal polarizability is also sensitive to the second tidal Love number k2. It is the main purpose of this work to perform a detailed study of k2 which, for a given compactness parameter, encodes the entire sensitivity of the tidal polarizability to the EOS. In particular, we examine the role that the crustal component of the EOS plays in the determination of k2. A set of realistic models of the equation of state that yield an accurate description of the properties of finite nuclei and support neutron stars of two solar masses is used. Given that the tidal polarizability scales as the fifth power of the compactness parameter, a universal relation exists among the tidal polarizability and the compactness parameter that is highly insensitive to the underlying EOS. Thus, besides an extraction of the tidal polarizabilities, a measurement of the individual stellar masses is also required to impact the mass-radius relation. However, we observe a strong sensitivity of k2 to the EOS, particularly to the contribution from the inner crust. Although by itself the tidal polarizability can not contribute to the determination of the mass-radius relation, future detections of binary neutron star mergers are poised to provide significant constraints on both the tidal polarizabilities and masses of the individual stars, and thus ultimately on the mass-radius relation. Yet, subleading corrections to the tidal polarizability are encoded in the second Love number k2 which displays a large sensitivity to the entire (crust-plus-core) EOS., Comment: 15 pages and 5 figures; submitted to Physical Review C

On August 17, 2017 the LIGO-Virgo collaboration detected for the first time gravitational waves from the binary merger of two neutron stars (GW170817). Unlike the merger of two black holes, the associated electromagnetic radiation was also detected by a host of telescopes operating over a wide range of frequencies---opening a brand new era of multimessenger astronomy. This historical detection is providing fundamental new insights into the astrophysical site for the r-process and on the nature of dense matter. In this contribution we examine the impact of GW170817 on the equation of state of neutron rich matter, particularly on the density dependence of the symmetry energy. Limits on the tidal polarizability extracted from GW170817 seem to suggest that the symmetry energy is soft, thereby excluding models that predict overly large stellar radii., Comment: 16 pages, 4 figures, Invited Lecture presented at the Zakopane Conference on Nuclear Physics "Extremes of the Nuclear Landscape", Zakopane, Poland, August 26-September 2, 2018. To be published in ACTA PHYSICA POLONICA B Volume 50

Background: The "proton radius puzzle" refers to an eight-year old problem that highlights major inconsistencies in the extraction of the charge radius of the proton from muonic Lamb-shift experiments as compared against experiments using elastic electron scattering. For the latter, the determination of the charge radius involves an extrapolation of the experimental form factor to zero momentum transfer. Purpose: To estimate the proton radius by introducing a novel non-parametric approach to model the electric form factor of the proton. Methods: Within a Bayesian paradigm, we develop a model flexible enough to fit the data without any parametric assumptions on the form factor. The Bayesian estimation is guided by imposing only two physical constraints on the form factor: (a) its value at zero momentum transfer (normalization) and (b) its overall shape, assumed to be a monotonically decreasing function of the momentum transfer. Variants of these assumptions are explored to assess the impact of these constraints. Results: So far our results are inconclusive in regard to the proton puzzle, as they depend on both, the assumed constrains and the range of experimental data used. For example, if only low momentum-transfer data is used, adopting only the normalization constraint provides a value compatible with the smaller muonic result, while imposing only the shape constraint favors the larger electronic value. Conclusions: We have presented a novel technique to estimate the proton radius from electron scattering data based on a non-parametric Gaussian process. We have shown the impact of the physical constraints imposed on the form factor and of the range of experimental data used. In this regard, we are hopeful that as this technique is refined and with the anticipated new results from the PRad experiment, we will get closer to resolve of the puzzle.

Neutron stars are unique cosmic laboratories for the exploration of matter under extreme conditions of density and neutron-proton asymmetry. Due to their enormous dynamic range, neutron stars display a myriad of exotic states of matter that are impossible to recreate under normal laboratory conditions. In these three lectures I will discuss how the strong synergy that has developed between nuclear physics and astrophysics will uncover some of the deepest secrets behind these fascinating objects. In particular, I will highlight the enormous impact that the very first detection of gravitational waves from the binary neutron-star merger GW170817 is having in constraining the composition, structure, and dynamics of neutron stars., Comment: Lectures presented at the XIV International Workshop on Hadron Physics, Florianopolis, Brazil, March 2018. Few typos corrected in revised version

Excited states of the $N=32$ nucleus $^{54}$Ti have been studied, via both inverse-kinematics proton scattering and one-neutron knockout from $^{55}$Ti by a liquid hydrogen target, using the GRETINA $\gamma$-ray tracking array. Inelastic proton-scattering cross sections and deformation lengths have been determined. A low-lying octupole state has been tentatively identified in $^{54}$Ti for the first time. A comparison of $(p,p')$ results on low-energy octupole states in the neutron-rich Ca and Ti isotopes with the results of Random Phase Approximation calculations demonstrates that the observed systematic behavior of these states is unexpected., Comment: 7 pages, 8 figures

The historical first detection of a binary neutron star merger by the LIGO-Virgo collaboration [B. P. Abbott et al. Phys. Rev. Lett. 119, 161101 (2017)] is providing fundamental new insights into the astrophysical site for the $r$-process and on the nature of dense matter. A set of realistic models of the equation of state (EOS) that yield an accurate description of the properties of finite nuclei, support neutron stars of two solar masses, and provide a Lorentz covariant extrapolation to dense matter are used to confront its predictions against tidal polarizabilities extracted from the gravitational-wave data. Given the sensitivity of the gravitational-wave signal to the underlying EOS, limits on the tidal polarizability inferred from the observation translate into constraints on the neutron-star radius. Based on these constraints, models that predict a stiff symmetry energy, and thus large stellar radii, can be ruled out. Indeed, we deduce an upper limit on the radius of a $1.4\,M_{\odot}$ neutron star of $R_{\star}^{1.4}\!<\!13.76\,{\rm km}$. Given the sensitivity of the neutron-skin thickness of ${}^{208}$Pb to the symmetry energy, albeit at a lower density, we infer a corresponding upper limit of about $R_{\rm skin}^{208}\!\lesssim\!0.25\,{\rm fm}$. However, if the upcoming PREX-II experiment measures a significantly thicker skin, this may be evidence of a softening of the symmetry energy at high densities---likely indicative of a phase transition in the interior of neutron stars., Comment: 6 pages, 4 figures, edited figures and the text with a more correct interpretation of the LIGO-Virgo results

It has been recently suggested that differences in the charge radii of mirror nuclei are proportional to the neutron-skin thickness of neutron-rich nuclei and to the slope of the symmetry energy $L$ [B.A. Brown Phys. Rev. Lett. 119, 122502 (2017)]. The determination of the neutron skin has important implications for nuclear physics and astrophysics. Although the use of electroweak probes provides a largely model-independent determination of the neutron skin, the experimental challenges are enormous. Thus, the possibility that differences in the charge radii of mirror nuclei may be used as a surrogate for the neutron skin is a welcome alternative. To test the validity of this assumption we perform calculations based on a set of relativistic energy density functionals that span a wide region of values of $L$. Our results confirm that the difference in charge radii of various neutron-deficient nickel isotopes and their corresponding mirror nuclei is indeed strongly correlated to both the neutron-skin thickness and $L$. Moreover, given that various neutron-star properties are also sensitive to $L$, a data-to-data relation emerges between the difference in charge radii of mirror nuclei and the radius of low-mass neutron stars., Comment: 5 pages, 3 figures

Nuclear astrophysics centers on the role of nuclear physics in the cosmos. In particular, nuclear masses at the limits of stability are critical in the development of stellar structure and the origin of the elements. In this contribution we test and validate the predictions of recently refined nuclear mass models against the newly published AME2016 compilation. The basic paradigm underlining the recently refined nuclear mass models is based on existing state-of-the-art models that are subsequently refined through the training of an artificial neural network. We observe a significant improvement in the Bayesian Neural Network (BNN) predictions relative to the corresponding "bare" models when compared to the nearly 50 new masses reported in the AME2016 compilation. Further, AME2016 estimates for the handful of impactful isotopes in the determination of r-process abundances are found to be in fairly good agreement with our theoretical predictions. Indeed, the BNN-improved Duflo-Zuker model predicts a root-mean-square deviation relative to experiment of about 400 keV. Given the excellent performance of the BNN refinement in confronting the recently published AME2016 compilation, we are confident of its critical role in our quest for mass models of the highest quality. Moreover, as uncertainty quantification is at the core of the BNN approach, the improved mass models are in a unique position to identify those nuclei that will have the strongest impact in resolving some of the outstanding questions in nuclear astrophysics., Comment: 7 pages, 3 figures, and 2 tables

A high-resolution study of the electromagnetic response of $^{206}$Pb below the neutron separation energy is performed using a ($\vec{\gamma}$,$\gamma'$) experiment at the HI$\vec{\gamma}$S facility. Nuclear resonance fluorescence with 100% linearly polarized photon beams is used to measure spins, parities, branching ratios, and decay widths of excited states in $^{206}$Pb from $4.9$ to $8.1$MeV. The extracted $\Sigma$$B$(E1)$\uparrow$ and $\Sigma$$B$(M1)$\uparrow$ values for the total electric and magnetic dipole strength below the neutron separation energy are 0.9$\pm$0.2e$^2$fm$^2$ and 8.3$\pm$2.0$\mu_{N}^2$, respectively. These measurements are found to be in very good agreement with the predictions from an energy-density functional (EDF) plus quasiparticle phonon model (QPM). Such a detailed theoretical analysis allows to separate the pygmy dipole resonance from both the tail of the giant dipole resonance and multi-phonon excitations. Combined with earlier photonuclear experiments above the neutron separation energy, one extracts a value for the electric dipole polarizability of $^{206}$Pb of $\alpha_{D}\!=\!122\pm10$mb/MeV. When compared to predictions from both the EDF+QPM and accurately calibrated relativistic EDFs, one deduces a range for the neutron-skin thickness of $R_{\rm skin}^{206}\!=\!0.12$-$0.19$fm and a corresponding range for the slope of the symmetry energy of $L\!=\!48$-$60$MeV. This newly obtained information is also used to estimate the Maxwellian-averaged radiative cross section $^{205}$Pb(n,$\gamma$)$^{206}$Pb at 30keV to be $\sigma\!=\!130\!\pm\!25$mb. The astrophysical impact of this measurement--on both the s-process in stellar nucleosynthesis and on the equation of state of neutron-rich matter--is discussed., Comment: 8 pages and 4 figures. Manuscript accepted for publication in Physics Letters B

The isoscalar monopole response of neutron-rich nuclei is sensitive to both the incompressibility coefficient of symmetric nuclear matter and the density dependence of the symmetry energy. The main goal of this paper is to explore the emergence, evolution, and origin of low energy monopole strength along the even-even calcium isotopes: from 40Ca to 60Ca. The distribution of isoscalar monopole strength is computed in a relativistic random phase approximation (RPA) using three effective interactions that have been calibrated to the properties of finite nuclei and neutron stars. A non-spectral approach is adopted that allows for an exact treatment of the continuum without any reliance on discretization. For the stable calcium isotopes, no evidence of low-energy monopole strength is observed, even as the 1f7/2 neutron orbital is being filled and the neutron-skin thickness progressively grows. Further, in contrast to experimental findings, a mild softening of the monopole response with increasing mass number is predicted. Beyond 48Ca, a significant amount of low-energy monopole strength emerges as soon as the weak-binding neutron orbitals (2p and 1f5/2) become populated. The emergence and evolution of low-energy strength is identified with transitions from these weakly-bound states into the continuum. Moreover, given that models with a soft symmetry energy tend to reach the neutron-drip line earlier than their stiffer counterparts, we identify an inverse correlation between the neutron-skin thickness and the inverse energy weighted sum., Comment: 14 pages, 8 figures, and 5 tables

[Background] Besides its intrinsic value as a fundamental nuclear-structure observable, the weak-charge density of 208Pb - a quantity that is closely related to its neutron distribution - is of fundamental importance in constraining the equation of state of neutron-rich matter. [Purpose] To assess the impact that a second electroweak measurement of the weak-charge form factor of 208Pb may have on the determination of its overall weak-charge density. [Methods] Using the two putative experimental values of the form factor, together with a simple implementation of Bayes' theorem, we calibrate a theoretically sound - yet surprisingly little known - symmetrized Fermi function, that is characterized by a density and form factor that are both known exactly in closed form. [Results] Using the charge form factor of 208Pb as a proxy for its weak-charge form factor, we demonstrate that using only two experimental points to calibrate the symmetrized Fermi function is sufficient to accurately reproduce the experimental charge form factor over a significant range of momentum transfers. [Conclusions] It is demonstrated that a second measurement of the weak-charge form factor of 208Pb supplemented by a robust theoretical input in the form of the symmetrized Fermi function, would place significant constraints on the neutron distribution of 208Pb and, ultimately, on the equation of state of neutron-rich matter., Comment: 14 pages, 3 tables, and 6 figures

It is the aim of this paper to discuss the impact of nuclear fusion on the EOS. This is a timely subject given the expected availability of increasingly exotic beams at rare isotope facilities\,\cite{balantekin2014}. In practice, we focus on $^{48}$Ca+$^{48}$Ca fusion. We employ three different approaches to calculate fusion cross-sections for a set of energy density functionals with systematically varying nuclear matter properties. Fusion calculations are performed using frozen densities, using a dynamic microscopic method based on density-constrained time-dependent Hartree-Fock (DC-TDHF) approach, as well as direct TDHF study of above barrier cross-sections. For these studies, we employ a family of Skyrme parametrizations with systematically varied nuclear matter properties. We find a slight preference for forces which deliver a slope of symmetry energy of $L\approx 50$\,MeV that corresponds to a neutron-skin thickness of $^{48}$Ca of $R_\mathrm{skin}\!=\!(0.180\!-\!0.210)$\,fm., Comment: 15 pages, 10 figures

The information on the symmetry energy and its density dependence is deduced by comparing the available data on the electric dipole polarizability $\alpha_D$ of ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb with the predictions of the Random Phase Approximation, using a representative set of nuclear energy density functionals. The calculated values of $\alpha_D$ are used to validate different correlations involving $\alpha_D$, the symmetry energy at the saturation density $J$, the corresponding slope parameter $L$ and the neutron skin thickness $\Delta r_{\!np}$, as suggested by the Droplet Model. A subset of models that reproduce simultaneously the measured polarizabilities in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are employed to predict the values of the symmetry energy parameters at saturation density and $\Delta r_{\!np}$. The resulting intervals are: $J\!=\!30 \text{-}35$ MeV, $L\!=\!20 \text{-} 66$ MeV; and the values for $\Delta r_{\!np}$ in ${}^{68}$Ni, ${}^{120}$Sn, and ${}^{208}$Pb are in the ranges: 0.15\text{-}0.19 fm, 0.12\text{-}0.16 fm, and 0.13\text{-}0.19 fm, respectively. The strong correlation between the electric dipole polarizabilities of two nuclei is instrumental to predict the values of electric dipole polarizabilities in other nuclei., Comment: 12 pages, 2 tables and 6 figures

Besides their intrinsic nuclear-structure value, nuclear mass models are essential for astrophysical applications, such as r-process nucleosynthesis and neutron-star structure. To overcome the intrinsic limitations of existing "state-of-the-art" mass models, we propose a refinement based on a Bayesian Neural Network (BNN) formalism. A novel BNN approach is implemented with the goal of optimizing mass residuals between theory and experiment. A significant improvement (of about 40%) in the mass predictions of existing models is obtained after BNN refinement. Moreover, these improved results are now accompanied by proper statistical errors. Finally, by constructing a "world average" of these predictions, a mass model is obtained that is used to predict the composition of the outer crust of a neutron star. The power of the Bayesian neural network method has been successfully demonstrated by a systematic improvement in the accuracy of the predictions of nuclear masses. Extension to other nuclear observables is a natural next step that is currently under investigation., Comment: 13 pages, 5 figures, submitted to Physical Review C

Recent progress in the determination of both masses and radii of neutron stars are starting to place stringent constraints on the dense matter equation of state. In particular, new theoretical developments together with improved statistical tools seem to favor stellar radii that are significantly smaller than those predicted by models using purely nucleonic equations of state. Given that the underlying equation of state must also account for the observation of $2M_{\odot}$ neutron stars, theoretical approaches to the study of the dense matter equation of state are facing serious challenges. In response to this challenge, we compute in a model-independent way the underlying equation of state associated with an assumed mass-radius template similar to the "common radius" assumption used in recent studies. Once such a mass-radius template is adopted, the equation of state follows directly from the implementation of Lindblom's algorithm; assumptions on the nature or composition of the dense stellar core are not required. By analyzing mass-radius profiles with a maximum mass consistent with observation and common radii in the 8 to 11 km range, a lower limit on the stellar radius of a $1.4M_{\odot}$ neutron star of $R_{\rm NS}\!\gtrsim\!10.7$ km is obtained in order to prevent the equation of state from violating causality., Comment: 5 pages, 4 figures

The main goal of the present contribution is a pedagogical introduction to the fascinating world of neutron stars by relying on relativistic density functional theory. Density functional theory provides a powerful--and perhaps unique--framework for the calculation of both the properties of finite nuclei and neutron stars. Given the enormous densities that may be reached in the core of neutron stars, it is essential that such theoretical framework incorporates from the outset the basic principles of Lorentz covariance and special relativity. After a brief historical perspective, we present the necessary details required to compute the equation of state of dense, neutron-rich matter. As the equation of state is all that is needed to compute the structure of neutron stars, we discuss how nuclear physics--particularly certain kind of laboratory experiments--can provide significant constrains on the behavior of neutron-rich matter., Comment: Contributing chapter to the book "Relativistic Density Functional for Nuclear Structure"; World Scientific Publishing Company (Singapore); Editor Prof. Jie Meng (Peking University)

We search for potential isovector signatures in the neutron-rich oxygen and calcium isotopes within the framework of a relativistic mean-field theory with an exact treatment of pairing correlations. To probe the isovector sector we calibrate a few relativistic density functionals using the same isoscalar constraints but with one differing isovector assumption. It is found that under certain conditions, the isotopic chain in oxygen can be made to terminate at the experimentally observed ${}^{24}$O isotope and in the case of the calcium isotopes at ${}^{60}$Ca. To produce such behavior, the resulting symmetry energy must be soft, with predicted values for the symmetry energy and its slope at saturation density being $J\!=\!(30.92\pm0.47)$ MeV and $L\!=\!(51.0\pm1.5)$ MeV, respectively. As a consequence, the neutron-skin thickness of ${}^{208}$Pb is rather small: $R_{\rm skin}^{208}\!=\!(0.161\pm0.011)$ fm. This same model - labelled "FSUGarnet" - predicts $R_{1.4}\!=\!(13.1\pm0.1)$ km for the radius of a "canonical" 1.4$M_{\odot}$ neutron star, yet is also able to support a two-solar-mass neutron star., Comment: 6 pages, 5 figures

In this new era of radioactive beam facilities, the discovery of novel modes of excitation in nuclei far away from stability represents an area of intense research activity. In addition, these modes of excitation appear to be sensitive to the uncertain density dependence of the symmetry energy. We study the emergence, evolution, and nature of both the soft and giant isoscalar monopole modes as a function of neutron excess in three unstable Nickel isotopes: 56Ni, 68Ni, and 78Ni. The distribution of isoscalar monopole strength is computed in a relativistic random-phase approximation using several accurately calibrated effective interactions. In particular, a non-spectral Green's function approach is adopted that allows for an exact treatment of the continuum without any reliance on discretization. The discretization of the continuum is neither required nor admitted. In the case of 56Ni, the lack of low-energy strength results in a direct correlation between the centroid energy of the giant monopole resonance and the incompressibility coefficient of symmetric nuclear matter. In contrast, the large neutron excess in both 68Ni and 78Ni generates a significant, yet relatively featureless, amount of low-energy strength that is driven by transitions into the continuum. Moreover, the evolution of monopole strength with neutron excess displays sensitivity to the density dependence of the symmetry energy. Our results suggest that future measurements of the distribution of isoscalar monopole strength at radioactive beam facilities using a very long chain of both stable and unstable isotopes could place important constraints on the equation of state of neutron-rich matter and ultimately on the properties of neutron stars. However, given the nature of the low-energy monopole excitations, a proper treatment of the continuum is essential., Comment: Manuscript contains 12 pages, 3 figures, and 3 tables. Submitted to the Physical Review C

Background: Theoretical approaches based on density functional theory provide the only tractable method to incorporate the wide range of densities and isospin asymmetries required to describe finite nuclei, infinite nuclear matter, and neutron stars. Purpose: A relativistic energy density functional (EDF) is developed to address the complexity of such diverse nuclear systems. Moreover, a statistical perspective is adopted to describe the information content of various physical observables. Methods: We implement the model optimization by minimizing a suitably constructed chi-square objective function using various properties of finite nuclei and neutron stars. The minimization is then supplemented by a covariance analysis that includes both uncertainty estimates and correlation coefficients. Results: A new model, FSUGold2, is created that can well reproduce the ground-state properties of finite nuclei, their monopole response, and that accounts for the maximum neutron star mass observed up to date. In particular, the model predicts both a stiff symmetry energy and a soft equation of state for symmetric nuclear matter--suggesting a fairly large neutron-skin thickness in Pb208 and a moderate value of the nuclear incompressibility. Conclusions: We conclude that without any meaningful constraint on the isovector sector, relativistic EDFs will continue to predict significantly large neutron skins. However, the calibration scheme adopted here is flexible enough to create models with different assumptions on various observables. Such a scheme--properly supplemented by a covariance analysis--provides a powerful tool to identify the critical measurements required to place meaningful constraints on theoretical models., Comment: 23 pages, 7 figures

Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of very limited utility. In this contribution we use maximum-likelihood estimation to compute theoretical errors and assess correlations between physical observables. We illustrate the power and elegance of these methods using examples of both pedagogical and realistic interest. In particular, we implement a gaussian approximation to the likelihood function to develop a new relativistic effective interaction constrained by ground-state properties of finite nuclei, their monopole response, and masses of neutron stars., Comment: 17 pages, 5 figures. Submitted to Journal of Physics G Focus issue "Enhancing the interaction between nuclear experiment and theory through information and statistics" (ISNET)

Pulsar glitches-the sudden spin-up in the rotational frequency of a neutron star-suggest the existence of an angular-momentum reservoir confined to the inner crust of the neutron star. Large and regular glitches observed in the Vela pulsar have originally constrained the fraction of the stellar moment of inertia that must reside in the solid crust to about 1.4%. However, crustal entrainment-which until very recently has been ignored-suggests that in order to account for the Vela glitches, the fraction of the moment of inertia residing in the crust must increase to about 7%. This indicates that the required angular momentum reservoir may exceed that which is available in the crust. We explore the possibility that uncertainties in the equation of state provide enough flexibility for the construction of models that predict a large crustal thickness and consequently a large crustal moment of inertia. Given that analytic results suggest that the crustal moment of inertia is sensitive to the transition pressure at the crust-core interface, we tune the parameters of the model to maximize the transition pressure, while still providing an excellent description of nuclear observables. In this manner we are able to obtain fractional moments of inertia as large as 7% for neutron stars with masses below 1.6 solar masses. In particular, we find that if the neutron-skin thickness of 208Pb falls within the (0.20-0.26) fm range, large enough transition pressures can be generated to explain the large Vela glitches without invoking an additional angular-momentum reservoir beyond that confined to the solid crust. Our results suggest that the crust may be enough., Comment: 14 pages, 5 figures, and 4 tables

The symmetry energy describes how the energy of nuclear matter rises as one goes away from equal numbers of neutrons and protons. This is very important to describe neutron rich matter in astrophysics. This article reviews our knowledge of the symmetry energy from theoretical calculations, nuclear structure measurements, heavy ion collisions, and astronomical observations. We then present a roadmap to make progress in areas of relevance to the symmetry energy that promotes collaboration between the astrophysics and the nuclear physics communities., Comment: 19 pages, results from first International Collaborations in Nuclear Theory (ICNT) program at NSCL/FRIB, submitted to J. Phys G

We explore the unique and fascinating structure of neutron stars. Although neutron stars are of interest in many areas of Physics, our aim is to provide an intellectual bridge between Nuclear Physics and Astrophysics. We argue against the naive perception of a neutron star as a uniform assembly of neutrons packed to enormous densities. Rather, by focusing on the many exotic phases that are speculated to exist in a neutron star, we show how the reality is different and far more interesting., Comment: 8 pages and 4 figures. Excerpts from a lecture presented at the Seventh European Summer School on Experimental Nuclear Astrophysics; Catania, IT (September, 2013)

Background: Pairing correlations play a critical role in determining numerous properties of open-shell nuclei. Traditionally, they are included in a mean-field description of atomic nuclei through the approximate Bardeen-Cooper-Schrieffer or Hartree-Fock-Bogoliubov formalism. Purpose: We propose a new hybrid ''relativistic-mean-field-plus-pairing'' approach in which pairing is treated exactly so the number of particles is conserved. To verify the reliability of the formalism, we apply it to the study of both ground-state properties and isoscalar monopole excitations of the Tin isotopes. Methods: Accurately-calibrated relativistic mean-field models supplemented by an exact treatment of pairing correlations are used to compute ground-state observables along the isotopic chain in Tin. In turn, ground-state densities are used as input to the calculation of giant monopole resonances through a constrained-relativistic approach. Results: We compute a variety of ground-state observables sensitive to pairing correlations as well as the evolution of giant monopole energies along the isotopic chain in Tin. Whereas ground-state properties are consistent with experiment, we find that pairing correlations have a minor effect on the giant monopole energies. Conclusions: A new mean-field-plus-pairing approach is introduced to compute properties of open-shell nuclei. The formalism provides an efficient and powerful alternative to the computation of both ground-state properties and monopole energies of open-shell nuclei. We find ground-state properties to be well reproduced in this approach. However, as many have concluded before us, we find that pairing correlations are unlikely to provide an answer to the question of ''why is Tin so soft?'', Comment: 13 pages, 4 figures

Parity-violating electron scattering provides a model-independent determination of the nuclear weak-charge form factor that has widespread implications across such diverse areas as fundamental symmetries, nuclear structure, heavy-ion collisions, and neutron-star structure. We assess the impact of precise measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb on a variety of nuclear observables, such as the neutron skin and the electric-dipole polarizability. We use the nuclear Density Functional Theory with several accurately calibrated non-relativistic and relativistic energy density functionals. To assess the degree of correlation between nuclear observables and to explore systematic and statistical uncertainties on theoretical predictions, we employ the chi-square statistical covariance technique. We find a strong correlation between the weak-charge form factor and the neutron radius, that allows for an accurate determination of the neutron skin of neutron-rich nuclei. We determine the optimal range of the momentum transfer $q$ that maximizes the information content of the measured weak-charge form factor and quantify the uncertainties associated with the strange quark contribution. Moreover, we confirm the role of the electric-dipole polarizability as a strong isovector indicator. Accurate measurements of the weak-charge form factor of ${}^{48}$Ca and ${}^{208}$Pb will have a profound impact on many aspects of nuclear theory and hadronic measurements of neutron skins of exotic nuclei at radioactive-beam facilities., Comment: 10 pages, 4 figures