Researcher profile

D. R. Phillips

D. R. Phillips contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2022arXiv

Analyzing rotational bands in odd-mass nuclei using Effective Field Theory and Bayesian methods

We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of $^{99}$Tc, ${}^{155,157}$Gd, ${}^{159}$Dy, ${}^{167, 169}$Er, ${}^{167, 169}$Tm, ${}^{183}$W, ${}^{235}$U and ${}^{239}$Pu. The error model in our Bayesian analysis includes both experimental and EFT truncation uncertainties. It also accounts for the fact that low-energy constants (LECs) at even and odd orders are expected to have different sizes. We use Markov Chain Monte Carlo (MCMC) sampling to explore the joint posterior of the EFT and error-model parameters and show both the LECs and the expansion parameter, $Q$, can be reliably determined. We extract the LECs up to fourth order in the EFT and find that, provided we correctly account for EFT truncation errors in our likelihood, results for lower-order LECs are stable as we go to higher orders. LEC results are also stable with respect to the addition of higher-energy data. We extract the expansion parameter for all the nuclei listed above and find a clear correlation between the extracted and the expected $Q$ based on the single-particle and vibrational energy scales. However, the $Q$ that actually determines the convergence of the EFT expansion is markedly smaller than would be naively expected based on those scales.

preprint2022arXiv

Does Bayesian Model Averaging improve polynomial extrapolations? Two toy problems as tests

We assess the accuracy of Bayesian polynomial extrapolations from small parameter values, x, to large values of x. We consider a set of polynomials of fixed order, intended as a proxy for a fixed-order effective field theory (EFT) description of data. We employ Bayesian Model Averaging (BMA) to combine results from different order polynomials (EFT orders). Our study considers two "toy problems" where the underlying function used to generate data sets is known. We use Bayesian parameter estimation to extract the polynomial coefficients that describe these data at low x. A "naturalness" prior is imposed on the coefficients, so that they are O(1). We Bayesian-Model-Average different polynomial degrees by weighting each according to its Bayesian evidence and compare the predictive performance of this Bayesian Model Average with that of the individual polynomials. The credibility intervals on the BMA forecast have the stated coverage properties more consistently than does the highest evidence polynomial, though BMA does not necessarily outperform every polynomial.

preprint2022arXiv

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution for the ${}^3$H binding energy, the ${}^4$He binding energy and radius, and the ${}^3$H $β$-decay rate can only be obtained if $χ$EFT truncation errors are included in the analysis. All of these except the $β$-decay rate give essentially degenerate constraints on the 3NF low-energy constants, so it is crucial for estimating these parameters. We use eigenvector continuation for fast and accurate emulation of No-Core Shell Model calculations of the considered few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the hyperparameters that quantify the truncation error. We find a $χ$EFT expansion parameter of $Q=0.33 \pm 0.06$ for these observables.

preprint2021arXiv

How well do we know the neutron-matter equation of state at the densities inside neutron stars? A Bayesian approach with correlated uncertainties

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

preprint2021arXiv

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory ($χ$EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the $χ$EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. We then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium $χ$EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

preprint2020arXiv

An effective field theory approach to fermionic rotational bands in nuclei

We extend an effective field theory developed to describe rotational bands in even-even nuclei to the odd-mass case. This organizes Bohr and Mottelson's treatment of a particle coupled to a rotor as a model-independent expansion in powers of the angular velocity of the overall system. We carry out this expansion up to fourth order in the angular velocity and present results for $^{99}$Tc, ${}^{159}$Dy, ${}^{167, 169}$Er, ${}^{167, 169}$Tm, ${}^{183}$W, ${}^{235}$U and ${}^{239}$Pu. In each case, the accuracy and breakdown scale of the effective field theory can be understood based on the single-particle and vibrational energy scales in that nucleus.

preprint2020arXiv

Finite-Size Effects in Heavy Halo Nuclei from Effective Field Theory

Halo/Cluster Effective Field Theory describes halo/cluster nuclei in an expansion in the small ratio of the size of the core(s) to the size of the system. Even in the point-particle limit, neutron halo nuclei have a finite charge radius, because their center of mass does not coincide with their center of charge. This point-particle contribution decreases as $1/A_c$, where $A_c$ is the mass number of the core, and diminishes in importance compared to other effects, e.g., the size of the core to which the neutrons are bound. Here we propose that for heavy cores the EFT expansion should account for the small factors of $1/A_c$. As a specific example, we discuss the implications of this organizational scheme for the inclusion of finite-size effects in expressions for the charge radii of halo nuclei. We show in particular that a short-range operator could be the dominant effect in the charge radius of one-neutron halos bound by a P-wave interaction. The point-particle contribution remains the leading piece of the charge radius for one-proton halos, and so Halo EFT has more predictive power in that case.

preprint2019arXiv

$S$-factor and scattering-parameter extractions from ${}^{3}\mathrm{He} +{}^{4}\mathrm{He} \rightarrow {}^{7}\mathrm{Be} + γ$

Previous studies of the reaction ${}^{3}\mathrm{He} +{}^{4}\mathrm{He} \rightarrow {}^{7}\mathrm{Be} + γ$ have focused on providing the best central value and error bar for the $S$ factor at solar energies. Measurements of this capture reaction, the ${}^{3}\mathrm{He}$-${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$, have been used to constrain employed theoretical models. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for the reaction to perform the extrapolation. At this order the EFT describes the reaction using an s-wave scattering length and effective range, the asymptotic properties of the final bound states, and short-distance contributions to the $E1$ capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total $S$ factor $S(0)= 0.578^{+0.015}_{-0.016}$ keV~b, while the branching ratio for excited- to ground-state capture at zero energy, $Br(0)=0.406^{+0.013}_{-0.011}$, both at 68\% degree of belief. This $S(0)$ is broadly consistent with other recent evaluations, including the previously recommended value $S(0)=0.56 \pm 0.03$ eV b, but has a smaller error bar. We also find significant constraints on the scattering parameters, and we obtain constraints on the $S(E)$'s angular dependence. The path forward seems to lie with better measurements of the scattering phase shift and $S(E)$'s angular dependence, together with better understanding of the asymptotic normalization coefficients of the ${}^7$Be bound states' wave functions. Data on these could further reduce $S(0)$'s uncertainty.