BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. Piekarewicz

J. Piekarewicz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thnucl-exastro-ph.SRastro-ph.HEastro-phcond-matcond-mat.stat-mechcond-mat.str-elhep-latphysics.atom-ph

Catalog footprint

What is connected

47works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2026arXiv

Experimental study of $^{53}$Cr via the $(d,pγ)$ reaction

Excited states in $^{53}$Cr were studied via the $^{52}$Cr$(d,pγ)$ reaction up to the neutron-separation threshold. Proton-$γ$ angular correlations and $γ$ decay branching ratios were measured in particle-$γ$ coincidences between the Super-Enge Split-Pole Spectrograph (SE-SPS) and CeBr$_3$ Array (CeBrA) demonstrator of the John D. Fox Accelerator Laboratory at Florida State University. Previous spin-parity assignments from a $(d,p)$ singles experiment at the SE-SPS are supported and $γ$-ray transitions in $^{53}$Cr reported. We firmly assign higher-lying excited states to $^{53}$Cr because overlapping excited states and contaminants could be identified better due to the complementary $γ$-decay information. We also correct some of the previously reported excitation energies and present a reanalysis of previously measured $^{52}$Cr$(d,p){}^{53}${Cr} angular distributions guided by the complementary $γ$-ray information. Based on this reanalysis, the fragmentation of the neutron $2p_{3/2}$, $2p_{1/2}$, $1f_{5/2}$, $1g_{9/2}$, and $2d_{5/2}$ single-particle strengths is reassessed for $^{53}$Cr. A comparison to the corresponding strengths in $^{55}$Fe is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the possible existence of a soft dipole mode in 8He

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Compressional-mode resonances in the molybdenum isotopes: Emergence of softness in open-shell nuclei near A=90

"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($α,α^\prime$) reaction at $E_α= 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$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

GW190814: Impact of a 2.6 solar mass neutron star on nucleonic equations of state

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sensitivity of the fusion cross section to the density dependence of the symmetry energy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The power of two: Assessing the impact of a second measurement of the weak-charge form factor of 208Pb

[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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nuclear Mass Predictions for the Crustal Composition of Neutron Stars: A Bayesian Neural Network Approach

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the compactness of neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic density functional theory for finite nuclei and neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Searching for isovector signatures in the neutron-rich oxygen and calcium isotopes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The neutron skin thickness from the measured electric dipole polarizability in $^{68}$Ni, $^{120}$Sn, and $^{208}$Pb

The information on the symmetry energy and its density dependence is deduced by comparing the available data on the electric dipole polarizability $α_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 $α_D$ are used to validate different correlations involving $α_D$, the symmetry energy at the saturation density $J$, the corresponding slope parameter $L$ and the neutron skin thickness $Δ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 $Δr_{\!np}$. The resulting intervals are: $J\!=\!30 \text{-}35$ MeV, $L\!=\!20 \text{-} 66$ MeV; and the values for $Δ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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A way forward in the study of the symmetry energy: experiment, theory, and observation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Building relativistic mean field models for finite nuclei and neutron stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Information and statistics: a new paradigm in theoretical nuclear physics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear breathing mode in neutron-rich Nickel isotopes: sensitivity to the symmetry energy and the role of the continuum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pulsar Glitches: The Crust may be Enough

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electric Dipole Polarizability in ${}^{208}$Pb: insights from the Droplet Model

We study the electric dipole polarizability $α_D$ in ${}^{208}$Pb based on the predictions of a large and representative set of relativistic and non-relativistic nuclear mean field models. We adopt the droplet model as a guide to better understand the correlations between $α_D$ and other isovector observables. Insights from the droplet model suggest that the product of $α_D$ and the nuclear symmetry energy at saturation density $J$ is much better correlated with the neutron skin thickness $Δr_{np}$ of ${}^{208}$Pb than the polarizability alone. Correlations of $α_D J$ with $Δr_{np}$ and with the symmetry energy slope parameter $L$ suggest that $α_D J$ is a strong isovector indicator. Hence, we explore the possibility of constraining the isovector sector of thenuclear energy density functional by comparing our theoretical predictions against measurements of both $α_D$ and the parity-violating asymmetry in ${}^{208}$Pb. We find that the recent experimental determination of $α_D$ in ${}^{208}$Pb in combination with the range for the symmetry energy at saturation density $J=[31\pm (2)_{\rm est.}]$\,MeV suggests $Δr_{np}({}^{208}{\rm Pb}) = 0.165 \pm (0.009)_{\rm exp.} \pm (0.013)_{\rm theo.} \pm (0.021)_{\rm est.} {\rm fm}$ and $L= 43 \pm(6)_{\rm exp.} \pm (8)_{\rm theo.}\pm(12)_{\rm est.}$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Giant monopole energies from a constrained relativistic mean-field approach

Background: Average energies of nuclear collective modes may be efficiently and accurately computed using a non-relativistic constrained approach without reliance on a random phase approximation (RPA). Purpose: To extend the constrained approach to the relativistic domain and to establish its impact on the calibration of energy density functionals. Methods: Relativistic RPA calculations are compared against the predictions of the corresponding constrained approach using two accurately calibrated energy density functionals. Results: We find excellent agreement-at the 2% level or better-between the predictions of the relativistic RPA and the corresponding constrained approach for magic (or semi-magic) nuclei ranging from 16O to 208Pb. Conclusions: An efficient and accurate method is proposed for incorporating nuclear collective excitations into the calibration of future energy density functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Has a thick neutron skin in ${}^{208}$Pb been ruled out?

The Lead Radius Experiment (PREX) has provided the first model-independent evidence in favor of a neutron-rich skin in ${}^{208}$Pb. Although the error bars are large, the reported large central value of 0.33\,fm is particularly intriguing. To test whether such a thick neutron-skin in ${}^{208}$Pb is already incompatible with laboratory experiments or astrophysical observations, we employ relativistic models with neutron-skin thickness in ${}^{208}$Pb ranging from 0.16 to 0.33 fm to compute ground state properties of finite nuclei, their collective monopole and dipole response, and mass-{\sl vs}-radius relations for neutron stars. No compelling reason was found to rule out models with large neutron skins in ${}^{208}$Pb from the set of observables considered in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Information content of the weak-charge form factor

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron Skins and Neutron Stars

The neutron-skin thickness of heavy nuclei provides a fundamental link to the equation of state of neutron-rich matter, and hence to the properties of neutron stars. The Lead Radius Experiment ("PREX") at Jefferson Laboratory has recently provided the first model-independence evidence on the existence of a neutron-rich skin in 208Pb. In this contribution we examine how the increased accuracy in the determination of neutron skins expected from the commissioning of intense polarized electron beams may impact the physics of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pygmies, giants, and skins as laboratory constraints on the equation of state of neutron-rich matter

Laboratory experiments sensitive to the density dependence of the symmetry energy may place stringent constraints on the equation of state of neutron-rich matter and, thus, on the structure, dynamics, and composition of neutron stars. Understanding the equation of state of neutron-rich matter is a central goal of nuclear physics that cuts across a variety of disciplines. In this contribution I focus on how laboratory experiments on neutron skins and on both Pygmy and Giant resonances can help us elucidate the structure of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relativistic mean field plus exact pairing approach to open shell nuclei

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?''

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Symmetry Energy Constraints from Giant Resonances: A Theoretical Overview

Giant resonances encapsulate the dynamic response of the nuclear ground state to external perturbations. As such, they offer a unique view of the nucleus that is often not accessible otherwise. Although interesting in their own right, giant resonances are also enormously valuable in providing stringent constraints on the equation of state of asymmetric matter. We this view in mind, we focus on two modes of excitation that are essential in reaching this goal: the isoscalar giant monopole resonance (GMR) and the isovector giant dipole resonance (GDR). GMR energies in heavy nuclei are sensitive to the symmetry energy because they probe the incompressibility of neutron-rich matter. Unfortunately, access to the symmetry energy is hindered by the relatively low neutron-proton asymmetry of stable nuclei. Thus, the measurement of GMR energies in exotic nuclei is strongly encouraged. In the case of the GDR, we find the electric dipole polarizability of paramount importance. Indeed, the electric dipole polarizability appears as one of two laboratory observables -- with the neutron-skin thickness being the other -- that are highly sensitive to the density dependence of the symmetry energy. Finally, we identify the softness of skin and the nature of the pygmy resonance as important unsolved problems in nuclear structure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Nuclear Physics of Neutron Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Charged Ising Model of Neutron Star Matter

Background: The inner crust of a neutron star is believed to consist of Coulomb-frustrated complex structures known as "nuclear pasta" that display interesting and unique low-energy dynamics. Purpose: To elucidate the structure and composition of the neutron-star crust as a function of temperature, density, and proton fraction. Methods: A new lattice-gas model, the "Charged-Ising Model" (CIM), is introduced to simulate the behavior of neutron-star matter. Preliminary Monte Carlo simulations on 30^3 lattices are performed for a variety of temperatures, densities, and proton fractions. Results: Results are obtained for the heat capacity, pair-correlation function, and static structure factor for a variety of conditions appropriate to the inner stellar crust. Conclusions: Although relatively simple, the CIM captures the essence of Coulomb frustration that is required to simulate the subtle dynamics of the inner stellar crust. Moreover, the computationally demanding long-range Coulomb interactions have been pre-computed at the appropriate lattice sites prior to the start of the simulation resulting in enormous computational gains. This work demonstrates the feasibility of future CIM simulations involving a large number of particles as a function of density, temperature, and proton fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electric dipole polarizability and the neutron skin

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Giant Monopole Resonance in even-A Cd isotopes, the asymmetry term in nuclear incompressibility, and the "softness" of Sn and Cd nuclei

The isoscalar giant monopole resonance (ISGMR) in even-A Cd isotopes has been studied by inelastic $α$-scattering at 100 MeV/u and at extremely forward angles, including 0deg. The asymmetry term in the nuclear incompressibility extracted from the ISGMR in Cd isotopes is found to be $K_τ = -555 \pm 75$ MeV, confirming the value previously obtained from the Sn isotopes. ISGMR strength has been computed in relativistic RPA using NL3 and FSUGold effective interactions. Both models significantly overestimate the centroids of the ISGMR strength in the Cd isotopes. Combined with other recent theoretical effort, the question of the "softness" of the open-shell nuclei in the tin region remains open still.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of spin-orbit currents on the electroweak skin of neutron-rich nuclei

Background: Measurements of neutron radii provide important constraints on the isovector sector of nuclear density functionals and offer vital guidance in areas as diverse as atomic parity violation, heavy-ion collisions, and neutron-star structure. Purpose: To assess the impact of spin-orbit currents on the electromagnetic- and weak-charge radii of a variety of nuclei. Special emphasis is placed on the experimentally accessible electroweak skin, defined as the difference between weak-charge and electromagnetic-charge radii. Methods: Two accurately calibrated relativistic mean field models are used to compute proton, neutron, charge, and weak-charge radii of a variety of nuclei. Results: We find that spin-orbit contributions to the electroweak skin of light neutron-rich nuclei, such as 22O and 48Ca, are significant and result in a substantial increase in the size of the electroweak skin relative to the neutron skin. Conclusions: Given that spin-orbit contributions to both the charge and weak-charge radii of nuclei are often as large as present or anticipated experimental error bars, future calculations must incorporate spin-orbit currents in the calculation of electroweak form factors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutron skins and neutron stars

Background: The neutron skin of a heavy nucleus as well as many neutron-star properties are highly sensitive to the poorly constrained density dependence of the symmetry energy. Purpose: To provide for the first time meaningful theoretical errors and to assess the degree of correlation between the neutron-skin thickness of ${}^{208}$Pb and several neutron-star properties. Methods: A proper covariance analysis based on the predictions of an accurately-calibrated relativistic functional {\sl "FSUGold"} is used to quantify theoretical errors and correlation coefficients. Results: We find correlation coefficients of nearly one (or minus one) between the neutron-skin thickness of ${}^{208}$Pb and a host of observables of relevance to the structure, dynamics, and composition of neutron stars. Conclusions: We suggest that a follow-up PREX measurement, ideally with a 0.5% accuracy, could significantly constrain the equation of state of neutron-star matter.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pygmies, Giants, and Skins

Understanding the equation of state (EOS) of neutron-rich matter is a central goal of nuclear physics that cuts across a variety of disciplines. Indeed, the limits of nuclear existence, the collision of energetic heavy ions, the structure of neutron stars, and the dynamics of core-collapse supernova all depend critically on the nuclear-matter EOS. In this contribution I focus on the EOS of cold baryonic matter with special emphasis on its impact on the structure, dynamics, and composition of neutron stars. In particular, I discuss how laboratory experiments on neutron skins as well as on Pygmy and Giant resonances can help us elucidate the structure of these fascinating objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sensitivity of the electric dipole polarizability to the neutron skin thickness in ${}^{208}$Pb

The static dipole polarizability, $α_{\rm D}$, in ${}^{208}$Pb has been recently measured with high-resolution via proton inelastic scattering at the Research Center for Nuclear Physics (RCNP). This observable is thought to be intimately connected with the neutron skin thickness, $r_{\rm skin}$, of the same nucleus and, more fundamentally, it is believed to be associated with the density dependence of the nuclear symmetry energy. The impact of $r_{\rm skin}$ on $α_{\rm D}$ in ${}^{208}$Pb is investigated and discussed on the basis of a large and representative set of relativistic and non-relativistic nuclear energy density functionals (EDF).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accurate calibration of relativistic mean-field models: correlating observables and providing meaningful theoretical uncertainties

Theoretical uncertainties in the predictions of relativistic mean-field models are estimated using a chi-square minimization procedure that is implemented by studying the small oscillations around the chi-square minimum. By diagonalizing the matrix of second derivatives, one gains access to a wealth of information---in the form of powerful correlations---that would normally remain hidden. We illustrate the power of the {\sl covariance analysis} by using two relativistic mean-field models: (a) the original linear Walecka model and (b) the accurately calibrated FSUGold parametrization. In addition to providing meaningful theoretical uncertainties for both model parameters and predicted observables, the covariance analysis establishes robust correlations between physical observables. In particular, we show that whereas the correlation coefficient between the slope of the symmetry energy and the neutron-skin thickness of Lead is indeed very large, a 1% measurement of the neutron radius of Lead may only be able to constrain the slope of the symmetry energy to about 30%.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Influence of the nuclear symmetry energy on the structure and composition of the outer crust

We review and extend with nonrelativistic nuclear mean field calculations a previous study of the impact of the nuclear symmetry energy on the structure and composition of the outer crust of nonaccreting neutron stars. First, we develop a simple "toy model" to understand the most relevant quantities determining the structure and composition of the outer crust: the nuclear symmetry energy and the pressure of the electron gas. While the latter is a well determined quantity, the former ---specially its density dependence--- still lacks an accurate characterization. We thus focus on the influence of the nuclear symmetry energy on the crustal composition. For that, we employ different nuclear models that are accurate in the description of terrestrial nuclei. We show that those models with stiffer symmetry energies ---namely, those that generate thicker neutron skins in heavy nuclei and have smaller symmetry energies at subnormal nuclear densities--- generate more exotic isotopes in the stellar crust than their softer counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Proton fraction in the inner neutron-star crust

Monte Carlo simulations of neutron-rich matter of relevance to the inner neutron-star crust are performed for a system of A=5,000 nucleons. To determine the proton fraction in the inner crust, numerical simulations are carried out for a variety of densities and proton fractions. We conclude---as others have before us using different techniques---that the proton fraction in the inner stellar crust is very small. Given that the purported "nuclear pasta" phase in stellar crusts develops as a consequence of the long-range Coulomb interaction among protons, we question whether pasta formation is possible in such proton-poor environments. To answer this question, we search for physical observables sensitive to the transition between spherical nuclei and exotic pasta structures. Of particular relevance is the static structure factor S(k)---an observable sensitive to density fluctuations. However, no dramatic behavior was observed in S(k). We regard the identification of physical observables sensitive to the existence---or lack-thereof---of a pasta phase in proton-poor environments as an open problem of critical importance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pygmy Resonances and Neutron Skins

Motivated by a recent experiment, the distribution of electric dipole strength in the neutron-rich 68Ni isotope was computed using a relativistic random phase approximation with a set of effective interactions that - although well calibrated - predict significantly different values for the neutron-skin thickness in 208Pb. The emergence of low-energy "Pygmy" strength that exhausts about 5-8% of the energy weighted sum rule (EWSR) is clearly identified. In addition to the EWSR, special emphasis is placed on the dipole polarizability. In particular, our results suggest a strong correlation between the dipole polarizability of 68Ni and the neutron-skin thickness of 208Pb. Yet we find a correlation just as strong and an even larger sensitivity between the neutron-skin thickness of 208Pb and the fraction of the dipole polarizability exhausted by the Pygmy resonance. These findings suggest that the dipole polarizability may be used as a proxy for the neutron skin.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic effective interaction for nuclei, giant resonances, and neutron stars

Nuclear effective interactions are useful tools in astrophysical applications especially if one can guide the extrapolations to the extremes regions of isospin and density that are required to simulate dense, neutron-rich systems. Isospin extrapolations may be constrained in the laboratory by measuring the neutron skin thickness of a heavy nucleus, such as 208Pb. Similarly, future observations of massive neutron stars will constrain the extrapolations to the high-density domain. In this contribution we introduce a new relativistic effective interaction that is simultaneously constrained by the properties of finite nuclei, their collective excitations, and neutron-star properties. By adjusting two of the empirical parameters of the theory, one can efficiently tune the neutron skin thickness of 208Pb and the maximum neutron star mass. We illustrate this procedure in response to the recent interpretation of X-ray observations by Steiner, Lattimer, and Brown that suggests that the FSUGold effective interaction predicts neutron star radii that are too large and a maximum stellar mass that is too small. The new effective interaction is fitted to a neutron skin thickness in 208Pb of only R_n - R_p = 0.16 fm and a moderately large maximum neutron star mass of 1.94 Msun.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Relativistic models of the neutron-star matter equation of state

Motivated by a recent astrophysical measurement of the pressure of cold matter above nuclear-matter saturation density, we compute the equation of state of neutron star matter using accurately calibrated relativistic models. The uniform stellar core is assumed to consist of nucleons and leptons in beta equilibrium; no exotic degrees of freedom are included. We found the predictions of these models to be in fairly good agreement with the measured equation of state. Yet the Mass-vs-Radius relations predicted by these same models display radii that are consistently larger than the observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sensitivity of the Moment of Inertia of Neutron Stars to the Equation of State of Neutron-Rich Matter

The sensitivity of the stellar moment of inertia to the neutron-star matter equation of state is examined using accurately-calibrated relativistic mean-field models. We probe this sensitivity by tuning both the density dependence of the symmetry energy and the high density component of the equation of state, properties that are at present poorly constrained by existing laboratory data. Particularly attractive is the study of the fraction of the moment of inertia contained in the solid crust. Analytic treatments of the crustal moment of inertia reveal a high sensitivity to the transition pressure at the core-crust interface. This may suggest the existence of a strong correlation between the density dependence of the symmetry energy and the crustal moment of inertia. However, no correlation was found. We conclude that constraining the density dependence of the symmetry energy - through, for example, the measurement of the neutron skin thickness in 208Pb - will place no significant bound on either the transition pressure or the crustal moment of inertia.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Do we understand the incompressibility of neutron-rich matter?

The ``breathing mode'' of neutron-rich nuclei is our window into the incompressibility of neutron-rich matter. After much confusion on the interpretation of the experimental data, consistency was finally reached between different models that predicted both the distribution of isoscalar monopole strength in finite nuclei and the compression modulus of infinite matter. However, a very recent experiment on the Tin isotopes at the Research Center for Nuclear Physics(RCNP) in Japan has again muddled the waters. Self-consistent models that were successful in reproducing the energy of the giant monopole resonance (GMR) in nuclei with various nucleon asymmetries (such as 90Zr, 144Sm, and 208Pb) overestimate the GMR energies in the Tin isotopes. As important, the discrepancy between theory and experiment appears to grow with neutron excess. This is particularly problematic as models artificially tuned to reproduce the rapid softening of the GMR in the Tin isotopes become inconsistent with the behavior of dilute neutron matter. Thus, we regard the question of ``why is Tin so soft?'' as an important open problem in nuclear structure.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Atomic Parity Non-Conservation, Neutron Radii, and Effective Field Theories of Nuclei

Accurately calibrated effective field theories are used to compute atomic parity non-conserving (APNC) observables. Although accurately calibrated, these effective field theories predict a large spread in the neutron skin of heavy nuclei. While the neutron skin is strongly correlated to a large number of physical observables, in this contribution we focus on its impact on new physics through APNC observables. The addition of an isoscalar-isovector coupling constant to the effective Lagrangian generates a wide range of values for the neutron skin of heavy nuclei without compromising the success of the model in reproducing well constrained nuclear observables. Earlier studies have suggested that the use of isotopic ratios of APNC observables may eliminate their sensitivity to atomic structure. This leaves nuclear structure uncertainties as the main impediment for identifying physics beyond the standard model. We establish that uncertainties in the neutron skin of heavy nuclei are at present too large to measure isotopic ratios to better than the 0.1% accuracy required to test the standard model. However, we argue that such uncertainties will be significantly reduced by the upcoming measurement of the neutron radius in 208Pb at the Jefferson Laboratory.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Dynamical response of the nuclear pasta in neutron star crusts

The nuclear pasta -- a novel state of matter having nucleons arranged in a variety of complex shapes -- is expected to be found in the crust of neutron stars and in core-collapse supernovae at subnuclear densities of about $10^{14}$ g/cm$^3$. Due to frustration, a phenomenon that emerges from the competition between short-range nuclear attraction and long-range Coulomb repulsion, the nuclear pasta displays a preponderance of unique low-energy excitations. These excitations could have a strong impact on many transport properties, such as neutrino propagation through stellar environments. The excitation spectrum of the nuclear pasta is computed via a molecular-dynamics simulation involving up to 100,000 nucleons. The dynamic response of the pasta displays a classical plasma oscillation in the 1-2 MeV region. In addition, substantial strength is found at low energies. Yet this low-energy strength is missing from a simple ion model containing a single-representative heavy nucleus. The low-energy strength observed in the dynamic response of the pasta is likely to be a density wave involving the internal degrees of freedom of the clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Color screening in a constituent quark model of hadronic matter

The effect of color screening on the formation of a heavy quark-antiquark ($Q\bar{Q}$) bound state--such as the $J/ψ$ meson--is studied using a constituent-quark model. The response of the nuclear medium to the addition of two color charges is simulated directly in terms of its quark constituents via a string-flip potential that allows for quark confinement within hadrons yet enables the hadrons to separate without generating unphysical long-range forces. Medium modifications to the properties of the heavy meson, such as its energy and its mean-square radius, are extracted by solving Schrödinger's equation for the $Q\bar{Q}$ pair in the presence of a (screened) density-dependent potential. The density dependence of the heavy-quark potential is in qualitative agreement with earlier studies of its temperature dependence extracted from lattice calculations at finite temperature. In the present model it is confirmed that abrupt changes in the properties of the $J/ψ$-meson in the hadronic medium ({\it plasma}), correlate strongly with the deconfining phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Perturbation Theory for Spin Ladders Using Angular-Momentum Coupled Bases

We compute bulk properties of Heisenberg spin-1/2 ladders using Rayleigh-Schrödinger perturbation theory in the rung and plaquette bases. We formulate a method to extract high-order perturbative coefficients in the bulk limit from solutions for relatively small finite clusters. For example, a perturbative calculation for an isotropic $2\times 12$ ladder yields an eleventh-order estimate of the ground-state energy per site that is within 0.02% of the density-matrix-renormalization-group (DMRG) value. Moreover, the method also enables a reliable estimate of the radius of convergence of the perturbative expansion. We find that for the rung basis the radius of convergence is $λ_c\simeq 0.8$, with $λ$ defining the ratio between the coupling along the chain relative to the coupling across the chain. In contrast, for the plaquette basis we estimate a radius of convergence of $λ_c\simeq 1.25$. Thus, we conclude that the plaquette basis offers the only currently available perturbative approach which can provide a reliable treatment of the physically interesting case of isotropic $(λ=1)$ spin ladders. We illustrate our methods by computing perturbative coefficients for the ground-state energy per site, the gap, and the one-magnon dispersion relation.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

A Plaquette Basis for the Study of Heisenberg Ladders

We employ a plaquette basis-generated by coupling the four spins in a $2\times2$ lattice to a well-defined total angular momentum-for the study of Heisenberg ladders with antiferromagnetic coupling. Matrix elements of the Hamiltonian in this basis are evaluated using standard techniques in angular-momentum (Racah) algebra. We show by exact diagonalization of small ($2\times4$ and $2\times6$) systems that in excess of 90% of the ground-state probability is contained in a very small number of basis states. These few basis states can be used to define a severely truncated basis which we use to approximate low-lying exact eigenstates. We show how, in this low-energy basis, the isotropic spin-1/2 Heisenberg ladder can be mapped onto an anisotropic spin-1 ladder for which the coupling along the rungs is much stronger than the coupling between the rungs. The mapping thereby generates two distinct energy scales which greatly facilitates understanding the dynamics of the original spin-1/2 ladder. Moreover, we use these insights to define an effective low-energy Hamiltonian in accordance to the newly developed COntractor REnormalization group (CORE) method. We show how a simple range-2 CORE approximation to the effective Hamiltonian to be used with our truncated basis reproduces the low-energy spectrum of the exact $2\times6$ theory at the $\alt 1%$ level.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Relativistic treatment of the charged-current reaction ${}^{12}{\rm C}(ν_μ,μ^{-})$ near threshold

We compute the cross section for the charged-current reaction ${}^{12}{\rm C}(ν_μ,μ^{-})$ near threshold using a relativistic mean-field formalism. A reduced value of the nucleon mass in the medium --- coupled to the finite muon mass --- leads to a 30\% reduction in the cross section relative to a free Fermi-gas estimate. Isovector RPA correlations, which are strongly repulsive at these momentum transfers, reduce the cross section by an additional 15-30\%. Hence, relativistic nuclear-structure effects can account for the more than a factor-of-two reduction in the cross section recently reported by the LSND collaboration.

0 citations0 reviewsNo code linkedOpen access