Researcher profile

A. Schwenk

A. Schwenk contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
8topics
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

17 published item(s)

preprint2022arXiv

Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions

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

preprint2022arXiv

Importance truncation for the in-medium similarity renormalization group

Ab initio nuclear many-body frameworks require extensive computational resources, especially when targeting heavier nuclei. Importance-truncation (IT) techniques allow to significantly reduce the dimensionality of the problem by neglecting unimportant contributions to the solution of the many-body problem. In this work, we apply IT methods to the nonperturbative in-medium similarity renormalization group (IMSRG) approach and investigate the induced errors for ground-state energies in different mass regimes based on different nuclear Hamiltonians. We study various importance measures, which define the IT selection, and identify two measures that perform best, resulting in only small errors to the full IMSRG(2) calculations even for sizable compression ratios. The neglected contributions are accounted for in a perturbative way and serve as an estimate of the IT-induced error. Overall we find that the IT-IMSRG(2) performs well across all systems considered, while the largest compression ratios for a given error can be achieved when using soft Hamiltonians and for large single-particle bases.

preprint2022arXiv

Least-square approach for singular value decompositions of scattering problems

It was recently observed that chiral two-body interactions can be efficiently represented using matrix factorization techniques such as the singular value decomposition. However, the exploitation of these low-rank structures in a few- or many-body framework is nontrivial and requires reformulations that explicitly utilize the decomposition format. In this work, we present a general least-square approach that is applicable to different few- and many-body frameworks and allows for an efficient reduction to a low number of singular values in the least-square iteration. We verify the feasibility of the least-square approach by solving the Lippmann-Schwinger equation in factorized form. The resulting low-rank approximations of the $T$ matrix are found to fully capture scattering observables. Potential applications of the least-square approach to other frameworks with the goal of employing tensor factorization techniques are discussed.

preprint2021arXiv

Comparing different density-matrix expansions for long-range pion exchange

Empirical energy density functionals (EDFs) are generally successful in describing nuclear properties across the table of nuclides. But their limitations motivate using the density-matrix expansion (DME) to embed long-range pion interactions into a Skyrme functional. Recent results on the impact of the pion were both encouraging and puzzling, necessitating a careful re-examination of the DME implementation. Here we take the first steps, focusing on two-body scalar terms in the DME. Exchange energies with long-range one-pion contributions are well approximated by all DME implementations considered, with preference for variants that do not truncate at two derivatives in every EDF term. The use of the DME for chiral pion contributions is therefore supported by this investigation. For scalar-isovector energies it is important to treat neutrons and protons separately. The results are found to apply under broad conditions, although self-consistency is not yet tested.

preprint2021arXiv

Natural orbitals for many-body expansion methods

The nuclear many-body problem for medium-mass systems is commonly addressed using wave-function expansion methods that build upon a second-quantized representation of many-body operators with respect to a chosen computational basis. While various options for the computational basis are available, perturbatively constructed natural orbitals recently have been shown to lead to significant improvement in many-body applications yielding faster model-space convergence and lower sensitivity to basis set parameters in large-scale no-core shell model diagonalizations. This work provides a detailed comparison of single-particle basis sets and a systematic benchmark of natural orbitals in nonperturbative many-body calculations using the in-medium similarity renormalization group approach. As a key outcome we find that the construction of natural orbitals in a large single-particle basis enables for performing the many-body calculation in a reduced space of much lower dimension, thus offering significant computational savings in practice that help extend the reach of ab initio methods towards heavier masses and higher accuracy.

preprint2021arXiv

New equations of state constrained by nuclear physics, observations, and QCD calculations of high-density nuclear matter

We present new equations of state for applications in core-collapse supernova and neutron star merger simulations. We start by introducing an effective mass parametrization that is fit to recent microscopic calculations up to twice saturation density. This is important to capture the predicted thermal effects, which have been shown to determine the proto-neutron star contraction in supernova simulations. The parameter range of the energy-density functional underlying the equation of state is constrained by chiral effective field theory results at nuclear densities as well as by functional renormalization group computations at high densities based on QCD. We further implement observational constraints from measurements of heavy neutron stars, the gravitational wave signal of GW170817, and from the recent NICER results. Finally, we study the resulting allowed ranges for the equation of state and for properties of neutron stars, including the predicted ranges for the neutron star radius and maximum mass.

preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

preprint2020arXiv

$\boldsymbol{N=32}$ shell closure below calcium: Low-lying structure of $^{50}$Ar

Low-lying excited states in the $N=32$ isotope $^{50}$Ar were investigated by in-beam $γ$-ray spectroscopy following proton- and neutron-knockout, multi-nucleon removal, and proton inelastic scattering at the RIKEN Radioactive Isotope Beam Factory. The energies of the two previously reported transitions have been confirmed, and five additional states are presented for the first time, including a candidate for a 3$^-$ state. The level scheme built using $γγ$ coincidences was compared to shell-model calculations in the $sd-pf$ model space, and to ab initio predictions based on chiral two- and three-nucleon interactions. Theoretical proton- and neutron-knockout cross sections suggest that two of the new transitions correspond to $2^+$ states, while the previously proposed $4^+$ state could also correspond to a $2^+$ state.

preprint2020arXiv

Ab initio limits of atomic nuclei

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

preprint2020arXiv

Ab initio short-range-correlation scaling factors from light to medium-mass nuclei

High-energy scattering processes, such as deep inelastic scattering (DIS) and quasielastic (QE) scattering provide a wealth of information about the structure of atomic nuclei. The remarkable discovery of the empirical linear relationship between the slope of the European Muon Collaboration (EMC) effect in DIS and the short-range-correlation (SRC) scaling factors $a_2$ in QE kinematics is naturally explained in terms of scale separation in effective field theory. This explanation has powerful consequences, allowing us to calculate and predict SRC scaling factors from ab initio low-energy nuclear theory. We present ab initio calculations of SRC scaling factors for a nucleus $A$ relative to the deuteron $a_2(A/d)$ and relative to $^3\rm He$ $a_2(A/^3\rm He)$ in light and medium-mass nuclei. Our framework further predicts that the EMC effect and SRC scaling factors have minimal or negligible isovector corrections.

preprint2020arXiv

Charge radius of the short-lived $^{68}$Ni and correlation with the dipole polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus $^{68}$Ni at the \mbox{$N=40$} subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α_{\rm D}$ has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus $^{48}$Ca. Three-particle--three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of $^{68}$Ni.

preprint2020arXiv

Constraining the dense matter equation of state with joint analysis of NICER and LIGO/Virgo measurements

The NICER collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. (2019) explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars. In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of neutron stars from the gravitational wave signal of the compact binary merger GW170817. We compare the constraints on the EOS to those set by the recent measurement of a 2.14 solar mass pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a neutron star-black hole merger, which yields weaker constraints on the EOS.

preprint2020arXiv

Dilute Fermi gas at fourth order in effective field theory

Using effective field theory methods, we calculate for the first time the complete fourth-order term in the Fermi-momentum or $k_{\rm F} a_s$ expansion for the ground-state energy of a dilute Fermi gas. The convergence behavior of the expansion is examined for the case of spin one-half fermions and compared against quantum Monte-Carlo results, showing that the Fermi-momentum expansion is well-converged at this order for $| k_{\rm F} a_s | \lesssim 0.5$.

preprint2020arXiv

First glimpse of the $N=82$ shell closure below $Z=50$ from masses of neutron-rich cadmium isotopes and isomers

We probe the $N=82$ nuclear shell closure by mass measurements of neutron-rich cadmium isotopes with the ISOLTRAP spectrometer at ISOLDE-CERN. The new mass of $^{132}$Cd offers the first value of the $N=82$, two-neutron shell gap below $Z=50$ and confirms the phenomenon of mutually enhanced magicity at $^{132}$Sn. Using the recently implemented phase-imaging ion-cyclotron-resonance method, the ordering of the low-lying isomers in $^{129}$Cd and their energies are determined. The new experimental findings are used to test large-scale shell-model, mean-field and beyond-mean-field calculations, as well as the ab initio valence-space in-medium similarity renormalization group.

preprint2020arXiv

Spin-polarized neutron matter, the maximum mass of neutron stars, and GW170817

We investigate how a phase transition from neutron-star matter to spin-polarized neutron matter affects the equation of state and mass-radius relation of neutron stars. While general extension schemes for the equation of state allow for high pressures inside neutron stars, we find that a phase transition to spin-polarized neutron matter excludes extreme regimes. Hence, such a transition limits the maximum mass of neutron stars to lie below 2.6$-$$2.9 \, M_{\odot}$, depending on the microscopic nuclear forces used, while significantly larger masses could be reached without these constraints. These limits are in good agreement with recent constraints extracted from the neutron-star merger GW170817 and its electromagnetic counterpart. Assuming the description in terms of spin-polarized neutron matter to be valid in the center of neutron stars, we find that stars with a large spin-polarized domain in their core are ruled out by GW170817.

preprint2020arXiv

Testing ab initio nuclear structure in neutron-rich nuclei: lifetime measurements of second 2+ states in 16C and 20O

To test the predictive power of ab initio nuclear structure theory, the lifetime of the second 2+ state in neutron-rich 20O, tau(2+_2 ) = 150(+80-30) fs, and an estimate for the lifetime of the second 2+ state in 16C have been obtained, for the first time. The results were achieved via a novel Monte Carlo technique that allowed us to measure nuclear state lifetimes in the tens-to-hundreds femtoseconds range, by analyzing the Doppler-shifted gamma-transition line shapes of products of low-energy transfer and deep-inelastic processes in the reaction 18O (7.0 MeV/u) + 181Ta. The requested sensitivity could only be reached owing to the excellent performances of the AGATA gamma-tracking array, coupled to the PARIS scintillator array and to the VAMOS++ magnetic spectrometer. The experimental lifetimes agree with predictions of ab initio calculations using two- and three-nucleon interactions, obtained with the valence-space in-medium similarity renormalization group for 20O, and with the no-core shell model for 16C. The present measurement shows the power of electromagnetic observables, determined with high-precision gamma spectroscopy, to assess the quality of first-principles nuclear structure calculations, complementing common benchmarks based on nuclear energies. The proposed experimental approach will be essential for short lifetimes measurements in unexplored regions of the nuclear chart, including r-process nuclei, when intense ISOL-type beams become available.

preprint2019arXiv

Electromagnetic properties of $^{21}$O for benchmarking nuclear Hamiltonians

The structure of exotic nuclei provides valuable tests for state-of-the-art nuclear theory. In particular electromagnetic transition rates are more sensitive to aspects of nuclear forces and many-body physics than excitation energies alone. We report the first lifetime measurement of excited states in $^{21}$O, finding $τ_{1/2^+}=420^{+35}_{-32}\text{(stat)}^{+34}_{-12}\text{(sys)}$\,ps. This result together with the deduced level scheme and branching ratio of several $γ$-ray decays are compared to both phenomenological shell-model and ab initio calculations based on two- and three-nucleon forces derived from chiral effective field theory. We find that the electric quadrupole reduced transition probability of $\rm B(E2;1/2^+ \rightarrow 5/2^+_{g.s.}) = 0.71^{+0.07\ +0.02}_{-0.06\ -0.06}$~e$^2$fm$^4$, derived from the lifetime of the $1/2^+$ state, is smaller than the phenomenological result where standard effective charges are employed, suggesting the need for modifications of the latter in neutron-rich oxygen isotopes. We compare this result to both large-space and valence-space ab initio calculations, and by using multiple input interactions we explore the sensitivity of this observable to underlying details of nuclear forces.