Source author record

T. Papenbrock

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

48works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

48 published item(s)

preprint2022arXiv

Angular-momentum projection in coupled-cluster theory: structure of $^{34}$Mg

Single-reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for $^8$Be and $^{20}$Ne using a two-nucleon interaction from chiral effective field theory and for $pf$-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus $^{34}$Mg and find agreement with data.

preprint2022arXiv

Effective field theory of pairing rotations

Pairing rotations are the low-energy excitations of finite superfluid systems, connecting systems that differ in their number of Cooper pairs. This paper presents a model-independent derivation of pairing rotations within an effective theory that exploits the emergent breaking of U(1) phase symmetries. The symmetries are realized nonlinearly and the Nambu-Goldstone modes depend only on time because the system is finite. Semi-magic nuclei exhibit pairing rotational bands while the pairing spectrum becomes an elliptical paraboloid for open-shell nuclei. Model-independent relations between double charge-exchange reactions and alpha particle capture or knockout in open-shell nuclei are in analogy to the pair transfer reactions in a single superfluid. Odd semi-magic nuclei are described by coupling a fermion to the superfluid. The leading-order theories reproduce data for pairing rotational bands within uncertainty estimates.

preprint2021arXiv

Nuclear charge radii of Na isotopes: A tale of two theories

The accuracy of atomic theory calculations limits the extraction of nuclear charge radii from isotope shift measurements of odd-proton nuclei. For Na isotopes, though precise spectroscopic measurements have existed since more than half a century, calculations by different methods offer a wide range of values. Here, we present accurate atomic calculations to reliably extract the Na charge radii. By combining experimental matter radii with nuclear coupled-cluster calculations based on nucleon-nucleon and three-nucleon forces, we constrain the parameters obtained from the atomic calculations. Therefore, this study guides atomic theory and highlights the importance of using accurate atomic and nuclear computations in our understanding of the size of light nuclei.

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

preprint2021arXiv

Two-pion exchange as a leading-order contribution in chiral effective field theory

Pion exchange is the central ingredient to nucleon-nucleon interactions used in nuclear structure calculations, and one pion exchange (OPE) enters at leading order in chiral effective field theory. In the $^{2S+1}L_J={^1S_0}$ partial wave, however, OPE and a contact term needed for proper renormalization fail to produce the qualitative, and quantitative, features of the scattering phase shifts. Cutoff variation also revealed a surprisingly low breakdown momentum of about 330 MeV in this partial wave. Here we show that potentials consisting of OPE, two pion exchange (TPE), and a single contact address these problems and yield accurate and renormalization group (RG) invariant phase shifts in the $^1S_0$ partial wave. We demonstrate that a leading-order potential with TPE can be systematically improved by adding a contact quadratic in momenta. For momentum cutoffs below about 500 MeV, the removal of relevant physics from TPE loops needs to be compensated by additional contacts to keep RG invariance. Inclusion of the $Δ$ isobar degree of freedom in the potential does not change the strong contributions of TPE.

preprint2020arXiv

A two-neutron halo is unveiled in $^{29}$F

We report the measurement of reaction cross sections ($σ_R^{\rm ex}$) of $^{27,29}$F with a carbon target at RIKEN. The unexpectedly large $σ_R^{\rm ex}$ and derived matter radius identify $^{29}$F as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the $2p_{3/2}$ orbital, thereby vanishing the shell closure associated with the neutron number $N = 20$. The results are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of $^{27}$F but are challenged for $^{29}$F.

preprint2019arXiv

Coherent elastic neutrino-nucleus scattering on 40Ar from first principles

Coherent elastic neutrino scattering on the 40Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin, and estimate systematic uncertainties. The neutron-skin thickness of 40Ar40 is consistent with results from density functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.

preprint2016arXiv

Effective field theory for vibrations in odd-mass nuclei

Heavy even-even nuclei exhibit low-energy collective excitations that are separated in scale from the microscopic (fermion) degrees of freedom. This separation of scale allows us to approach nuclear vibrations within an effective field theory (EFT). In odd-mass nuclei collective and single-particle properties compete at low energies, and this makes their description more challenging. In this article we describe odd-mass nuclei with ground-state spin $I=\sfrac{1}{2}$ by means of an EFT that couples a fermion to the collective degrees of freedom of an even-even core. The EFT relates observables such as energy levels, electric quadrupole ($E2$) transition strengths, and magnetic dipole ($M1$) moments of the odd-mass nucleus to those of its even-even neighbor, and allows us to quantify theoretical uncertainties. For isotopes of rhodium and silver the theoretical description is consistent with data within experimental and theoretical uncertainties. Several testable predictions are made.

preprint2016arXiv

Effective field theory in the harmonic oscillator basis

We develop interactions from chiral effective field theory (EFT) that are tailored to the harmonic oscillator basis. As a consequence, ultraviolet convergence with respect to the model space is implemented by construction and infrared convergence can be achieved by enlarging the model space for the kinetic energy. In oscillator EFT, matrix elements of EFTs formulated for continuous momenta are evaluated at the discrete momenta that stem from the diagonalization of the kinetic energy in the finite oscillator space. By fitting to realistic phase shifts and deuteron data we construct an effective interaction from chiral EFT at next-to-leading order. Many-body coupled-cluster calculations of nuclei up to 132Sn exhibit a fast convergence of ground-state energies and radii in feasible model spaces.

preprint2016arXiv

Electric dipole polarizability from first principles calculations

The electric dipole polarizability quantifies the low-energy behaviour of the dipole strength and is related to critical observables such as the radii of the proton and neutron distributions. Its computation is challenging because most of the dipole strength lies in the scattering continuum. In this paper we combine integral transforms with the coupled-cluster method and compute the dipole polarizability using bound-state techniques. Employing different interactions from chiral effective field theory, we confirm the strong correlation between the dipole polarizability and the charge radius, and study its dependence on three-nucleon forces. We find good agreement with data for the 4He, 40Ca, and 16O nuclei, and predict the dipole polarizability for the rare nucleus 22O.

preprint2016arXiv

Emergent properties of nuclei from ab initio coupled-cluster calculations

Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. This endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. This paper reviews some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO$_{\rm sat}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to $^{56}$Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon-nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. The coupling to the continuum impacts the energies of the $J^π= {1/2}^-,{3/2}^-,{7/2}^-,{3/2}^+$ states in $^{17,23,25}$O, and - contrary to naive shell-model expectations - the level ordering of the $J^π= {3/2}^+,{5/2}^+,{9/2}^+$ states in $^{53,55,61}$Ca.

preprint2016arXiv

Structure of $^{78}$Ni from first principles computations

Doubly magic nuclei have a simple structure and are the cornerstones for entire regions of the nuclear chart. Theoretical insights into the supposedly doubly magic $^{78}$Ni and its neighbors are challenging because of the extreme neutron-to-proton ratio and the proximity of the continuum. We predict the $J^π=2_1^+$ state in $^{78}$Ni from a correlation with the $J^π=2_1^+$ state in $^{48}$Ca using chiral nucleon-nucleon and three-nucleon interactions. Our results confirm that $^{78}$Ni is doubly magic, and the predicted low-lying states of $^{79,80}$Ni open the way for shell-model studies of many more rare isotopes.

preprint2016arXiv

Unexpectedly large charge radii of neutron-rich calcium isotopes

Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei.

preprint2015arXiv

Accurate nuclear radii and binding energies from a chiral interaction

With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO$_{\rm sat}$, yield accurate binding energies and radii of nuclei up to $^{40}$Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective $J^π=3^-$ states in $^{16}$O and $^{40}$Ca are described accurately, while spectra for selected $p$- and $sd$-shell nuclei are in reasonable agreement with experiment.

preprint2015arXiv

Charge, neutron, and weak size of the atomic nucleus

What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist and the size of neutron stars. We present an ab initio calculation of the neutron distribution of the neutron-rich nucleus $^{48}$Ca. We show that the neutron skin (difference between radii of neutron and proton distributions) is significantly smaller than previously thought. We also make predictions for the electric dipole polarizability and the weak form factor; both quantities are currently targeted by precision measurements. Based on ab initio results for $^{48}$Ca, we provide a constraint on the size of a neutron star.

preprint2015arXiv

Effective field theory for deformed atomic nuclei

We present an effective field theory (EFT) for a model-independent description of deformed atomic nuclei. In leading order this approach recovers the well-known results from the collective model by Bohr and Mottelson. When higher-order corrections are computed, the EFT accounts for finer details such as the variation of the moment of inertia with the band head and the small magnitudes of interband $E2$ transitions. For rotational bands with a finite spin of the band head, the EFT is equivalent to the theory of a charged particle on the sphere subject to a magnetic monopole field.

preprint2015arXiv

Effective field theory for nuclear vibrations with quantified uncertainties

We develop an effective field theory (EFT) for nuclear vibrations. The key ingredients - quadrupole degrees of freedom, rotational invariance, and a breakdown scale around the three-phonon level - are taken from data. The EFT is developed for spectra and electromagnetic moments and transitions. We employ tools from Bayesian statistics for the quantification of theoretical uncertainties. The EFT consistently describes spectra and electromagnetic transitions for $^{62}$Ni, $^{98,100}$Ru, $^{106,108}$Pd, $^{110,112,114}$Cd, and $^{118,120,122}$Te within the theoretical uncertainties. This suggests that these nuclei can be viewed as anharmonic vibrators.

preprint2015arXiv

Effective Field Theory of Emergent Symmetry Breaking in Deformed Atomic Nuclei

Spontaneous symmetry breaking in non-relativistic quantum systems has previously been addressed in the framework of effective field theory. Low-lying excitations are constructed from Nambu-Goldstone modes using symmetry arguments only. We extend that approach to finite systems. The approach is very general. To be specific, however, we consider atomic nuclei with intrinsically deformed ground states. The emergent symmetry breaking in such systems requires the introduction of additional degrees of freedom on top of the Nambu-Goldstone modes. Symmetry arguments suffice to construct the low-lying states of the system. In deformed nuclei these are vibrational modes each of which serves as band head of a rotational band.

preprint2015arXiv

Effective theory for the non-rigid rotor in an electromagnetic field: Toward accurate and precise calculations of E2 transitions in deformed nuclei

We present a model-independent approach to electric quadrupole transitions of deformed nuclei. Based on an effective theory for axially symmetric systems, the leading interactions with electromagnetic fields enter as minimal couplings to gauge potentials, while subleading corrections employ gauge-invariant non-minimal couplings. This approach yields transition operators that are consistent with the Hamiltonian, and the power counting of the effective theory provides us with theoretical uncertainty estimates. We successfully test the effective theory in homonuclear molecules that exhibit a large separation of scales. For ground-state band transitions of rotational nuclei, the effective theory describes data well within theoretical uncertainties at leading order. In order to probe the theory at subleading order, data with higher precision would be valuable. For transitional nuclei, next-to-leading order calculations and the high-precision data are consistent within the theoretical uncertainty estimates. We also study the faint inter-band transitions within the effective theory and focus on the $E2$ transitions from the $0^+_2$ band (the "$β$ band") to the ground-state band. Here, the predictions from the effective theory are consistent with data for several nuclei, thereby proposing a solution to a long-standing challenge.

preprint2015arXiv

Giant and pigmy dipole resonances in 4He, 16,22O, and 40Ca from chiral nucleon-nucleon interactions

We combine the coupled-cluster method and the Lorentz integral transform for the computation of inelastic reactions into the continuum. We show that the bound-state-like equation characterizing the Lorentz integral transform method can be reformulated based on extensions of the coupled-cluster equation-of-motion method, and we discuss strategies for viable numerical solutions. Starting from a chiral nucleon-nucleon interaction at next-to-next-to-next-to-leading order, we compute the giant dipole resonances of 4He, 16,22O and 40Ca, truncating the coupled-cluster equation-of-motion method at the two-particle-two-hole excitation level. Within this scheme, we find a low-lying E1 strength in the neutron-rich 22O nucleus, which compares fairly well with data from [Leistenschneider et al. Phys. Rev. Lett. 86, 5442 (2001)]. We also compute the electric dipole polariziability in 40Ca. Deficiencies of the employed Hamiltonian lead to overbinding, too small charge radii and a too small electric dipole polarizability in 40Ca.

preprint2015arXiv

Infrared extrapolations of quadrupole moments and transitions

We study the convergence of bound-state quadrupole moments in finite harmonic oscillator spaces. We derive an expression for the infrared extrapolation for the quadrupole moment of a nucleus and benchmark our results using different model interactions for the deuteron. We find good agreement between the analytically derived and numerically obtained convergence behavior. We also derive an extrapolation formula for electric quadrupole transitions and find good agreement with the numerical calculation of a simple system.

preprint2015arXiv

Infrared length scale and extrapolations for the no-core shell model

We precisely determine the infrared (IR) length scale of the no-core shell model (NCSM). In the NCSM, the $A$-body Hilbert space is truncated by the total energy, and the IR length can be determined by equating the intrinsic kinetic energy of $A$ nucleons in the NCSM space to that of $A$ nucleons in a $3(A-1)$-dimensional hyper-radial well with a Dirichlet boundary condition for the hyper radius. We demonstrate that this procedure indeed yields a very precise IR length by performing large-scale NCSM calculations for $^{6}$Li. We apply our result and perform accurate IR extrapolations for bound states of $^{4}$He, $^{6}$He, $^{6}$Li, $^{7}$Li. We also attempt to extrapolate NCSM results for $^{10}$B and $^{16}$O with bare interactions from chiral effective field theory over tens of MeV.

preprint2014arXiv

Effective Field Theory for Finite Systems with Spontaneously Broken Symmetry

We extend effective field theory to the case of spontaneous symmetry breaking in genuinely finite quantum systems such as small superfluid systems, molecules or atomic nuclei, and focus on deformed nuclei. In finite superfluids, symmetry arguments alone relate the spectra of systems with different particle numbers. For systems with non-spherical intrinsic ground states such as atomic nuclei or molecules, symmetry arguments alone yield the universal features of the low-lying excitations as vibrations that are the heads of rotational bands. The low-lying excitations in deformed nuclei differ from those in molecules because of symmetry properties caused by pairing.

preprint2014arXiv

Effects of three-nucleon forces and two-body currents on Gamow-Teller strengths

We optimize chiral interactions at next-to-next-to leading order to observables in two- and three-nucleon systems, and compute Gamow-Teller transitions in carbon-14, oxygen-22 and oxygen-24 using consistent two-body currents. We compute spectra of the daughter nuclei nitrogen-14, fluorine-22 and fluorine-24 via an isospin-breaking coupled-cluster technique, with several predictions. The two-body currents reduce the Ikeda sum rule, corresponding to a quenching factor q^2 ~ 0.84-0.92 of the axial-vector coupling. The half life of carbon-14 depends on the energy of the first excited 1+ state, the three-nucleon force, and the two-body current.

preprint2014arXiv

Infrared extrapolations for atomic nuclei

Harmonic oscillator model-space truncations introduce systematic errors to the calculation of binding energies and other observables. We identify the relevant infrared scaling variable and give values for this nucleus-dependent quantity. We consider isotopes of oxygen computed with the coupled-cluster method from chiral nucleon-nucleon interactions at next-to-next-to-leading order and show that the infrared component of the error is sufficiently understood to permit controlled extrapolations. By employing oscillator spaces with relatively large frequencies, well above the energy minimum, the ultraviolet corrections can be suppressed while infrared extrapolations over tens of MeVs are accurate for ground-state energies. However, robust uncertainty quantification for extrapolated quantities that fully accounts for systematic errors is not yet developed.

preprint2014arXiv

Systematic expansion for infrared oscillator basis extrapolations

Recent work has demonstrated that the infrared effects of harmonic oscillator basis truncations are well approximated by imposing a partial-wave Dirichlet boundary condition at a properly identified radius L. This led to formulas for extrapolating the corresponding energy E_L and other observables to infinite L and thus infinite basis size. Here we reconsider the energy for a two-body system with a Dirichlet boundary condition at L to identify and test a consistent and systematic expansion for E_L that depends only on observables. We also generalize the energy extrapolation formula to nonzero angular momentum, and apply it to the deuteron. Formulas given previously for extrapolating the radius are derived in detail.

preprint2014arXiv

Ultraviolet extrapolations in finite oscillator bases

The use of finite harmonic oscillator spaces in many-body calculations introduces both infrared (IR) and ultraviolet (UV) errors. The IR effects are well approximated by imposing a hard-wall boundary condition at a properly identified radius L_eff. We show that duality of the oscillator implies that the UV effects are equally well described by imposing a sharp momentum cutoff at a momentum Lambda_eff complementary to L_eff. By considering two-body systems with separable potentials, we show that the UV energy corrections depend on details of the potential, in contrast to the IR energy corrections, which depend only on the S-matrix. An adaptation of the separable treatment to more general interactions is developed and applied to model potentials as well as to the deuteron with realistic potentials. The previous success with a simple phenomenological form for the UV error is also explained. Possibilities for controlled extrapolations for A > 2 based on scaling arguments are discussed.

preprint2013arXiv

An optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order

We optimize the nucleon-nucleon interaction from chiral effective field theory at next-to-next- to-leading order. The resulting new chiral force NNLOopt yields χ^2 \approx 1 per degree of freedom for laboratory energies below approximately 125 MeV. In the A = 3, 4 nucleon systems, the contributions of three-nucleon forces are smaller than for previous parametrizations of chiral interactions. We use NNLOopt to study properties of key nuclei and neutron matter, and demonstrate that many aspects of nuclear structure can be understood in terms of this nucleon-nucleon interaction, without explicitly invoking three-nucleon forces.

preprint2013arXiv

Coupled-cluster calculations of nucleonic matter

Background: The equation of state (EoS) of nucleonic matter is central for the understanding of bulk nuclear properties, the physics of neutron star crusts, and the energy release in supernova explosions. Purpose: This work presents coupled-cluster calculations of infinite nucleonic matter using modern interactions from chiral effective field theory (EFT). It assesses the role of correlations beyond particle-particle and hole-hole ladders, and the role of three-nucleon-forces (3NFs) in nuclear matter calculations with chiral interactions. Methods: This work employs the optimized nucleon-nucleon NN potential NNLOopt at next-to-next-to leading-order, and presents coupled-cluster computations of the EoS for symmetric nuclear matter and neutron matter. The coupled-cluster method employs up to selected triples clusters and the single-particle space consists of a momentum-space lattice. We compare our results with benchmark calculations and control finite-size effects and shell oscillations via twist-averaged boundary conditions. Results: We provide several benchmarks to validate the formalism and show that our results exhibit a good convergence toward the thermodynamic limit. While neutron matter with interactions from chiral EFT is perturbative, symmetric nuclear matter requires nonperturbative approaches. Correlations beyond the standard particle-particle ladder approximation yield non-negligible contributions. The saturation point of symmetric nuclear matter is sensitive to the employed 3NFs and the employed regularization scheme. 3NFs with nonlocal cutoffs exhibit a considerably improved convergence than their local cousins. We are unable to find values for the parameters of the short-range part of the local 3NF that simultaneously yield acceptable values for the saturation point in symmetric nuclear matter and the binding energies of light nuclei.

preprint2013arXiv

Coupled-cluster computations of atomic nuclei

In the past decade, coupled-cluster theory has seen a renaissance in nuclear physics, with computations of neutron-rich and medium-mass nuclei. The method is efficient for nuclei with product-state references, and it describes many aspects of weakly bound and unbound nuclei. This report reviews the technical and conceptual developments of this method in nuclear physics, and the results of coupled-cluster calculations for nucleonic matter, and for exotic isotopes of helium, oxygen, calcium, and some of their neighbors.

preprint2013arXiv

Rotational constants of multi-phonon bands in an effective theory for deformed nuclei

We consider deformed nuclei within an effective theory that exploits the small ratio between rotational and vibrational excitations. For even-even nuclei, the effective theory predicts small changes in the rotational constants of bands built on multi-phonon excitations that are linear in the number of excited phonons. In 166Er and 168Er, this explains the main variations of the rotational constants of the two-phonon gamma vibrational bands. In 232Th, the effective theory correctly explains the trend that the rotational constants decrease with increasing spin of the band head. We also study the effective theory for deformed odd nuclei. Here, time-odd terms enter the Lagrangian and generate effective magnetic forces that yield the high level densities observed in such nuclei.

preprint2013arXiv

Universal properties of infrared oscillator basis extrapolations

Recent work has shown that a finite harmonic oscillator basis in nuclear many-body calculations effectively imposes a hard-wall boundary condition in coordinate space, motivating infrared extrapolation formulas for the energy and other observables. Here we further refine these formulas by studying two-body models and the deuteron. We accurately determine the box size as a function of the model space parameters, and compute scattering phase shifts in the harmonic oscillator basis. We show that the energy shift can be well approximated in terms of the asymptotic normalization coefficient and the bound-state momentum, discuss higher-order corrections for weakly bound systems, and illustrate this universal property using unitarily equivalent calculations of the deuteron.

preprint2012arXiv

Condensates of p-wave pairs are exact solutions for rotating two-component Bose gases

We derive exact analytical results for the wave functions and energies of harmonically trapped two-component Bose-Einstein condensates with weakly repulsive interactions under rotation. The isospin symmetric wave functions are universal and do not depend on the matrix elements of the two-body interaction. The comparison with the results from numerical diagonalization shows that the ground state and low-lying excitations consists of condensates of p-wave pairs for repulsive contact interactions, Coulomb interactions, and the repulsive interactions between aligned dipoles.

preprint2012arXiv

Continuum effects and three-nucleon forces in neutron-rich oxygen isotopes

We employ interactions from chiral effective field theory and compute binding energies, excited states, and radii for isotopes of oxygen with the coupled-cluster method. Our calculation includes the effects of three-nucleon forces and of the particle continuum, both of which are important for the description of neutron-rich isotopes in the vicinity of the nucleus O-24. Our main results are the placement of the neutron drip-line at O-24, the assignment of spins, parities and resonance widths for several low-lying states of the drip-line nucleus, and an efficient approximation that incorporates the effects of three-body interactions.

preprint2012arXiv

Evolution of shell structure in neutron-rich calcium isotopes

We employ interactions from chiral effective field theory and compute the binding energies and low-lying excitations of calcium isotopes with the coupled-cluster method. Effects of three-nucleon forces are included phenomenologically as in-medium two-nucleon interactions, and the coupling to the particle continuum is taken into account using a Berggren basis. The computed ground-state energies and the low-lying 2+ states for the isotopes 42,48,50,52Ca are in good agreement with data, and we predict the excitation energy of the first 2+ state in 54Ca at 1.9 MeV, displaying only a weak sub-shell closure. In the odd-mass nuclei 53,55,61Ca we find that the positive parity states deviate strongly from the naive shell model.

preprint2012arXiv

Time-dependent coupled-cluster method for atomic nuclei

We study time-dependent coupled-cluster theory in the framework of nuclear physics. Based on Kvaal's bi-variational formulation of this method [S. Kvaal, arXiv:1201.5548], we explicitly demonstrate that observables that commute with the Hamiltonian are conserved under time evolution. We explore the role of the energy and of the similarity-transformed Hamiltonian under real and imaginary time evolution and relate the latter to similarity renormalization group transformations. Proof-of-principle computations of He-4 and O-16 in small model spaces, and computations of the Lipkin model illustrate the capabilities of the method.

preprint2011arXiv

Ab initio coupled-cluster theory for open-shell nuclei

We develop a new method to describe properties of truly open-shell nuclei. This method is based on single-reference coupled-cluster theory and the equation-of-motion method with extensions to nuclei with $A\pm 2$ nucleons outside a closed shell. We perform proof-of-principle calculations for the ground states of the helium isotopes $^{3-6}$He and the first excited $2^+$ state in $^6$He. The comparison with exact results from matrix diagonalization in small model spaces demonstrates the accuracy of the coupled-cluster methods. Three-particle--one-hole excitations of $^4$He play an important role for the accurate description of $^6$He. For the open-shell nucleus $^6$He, the computational cost of the method is comparable with the coupled-cluster singles-and-doubles approximation while its accuracy is similar to coupled-cluster with singles, doubles and triples excitations.

preprint2011arXiv

Chaos in fermionic many-body systems and the metal-insulator transition

We show that finite Fermi systems governed by a mean field and a few-body interaction generically possess spectral fluctuations of the Wigner-Dyson type and are, thus, chaotic. Our argument is based on an analogy to the metal-insulator transition. We construct a sparse random-matrix ensemble ScE that mimics that transition. Our claim then follows from the fact that the generic random-matrix ensemble modeling a fermionic interacting many-body system is much less sparse than ScE.

preprint2011arXiv

Exploring the anomaly in the interaction cross section and matter radius of 23O

New measurements of the interaction cross sections of 22,23O at 900A MeV performed at the GSI, Darmstadt are reported that address the unsolved puzzle of the large cross section previously observed for 23O. The matter radii for these oxygen isotopes extracted through a Glauber model analysis are in good agreement with the new predictions of the ab initio coupled-cluster theory reported here. They are consistent with a 22O+neutron description of 23O as well.

preprint2011arXiv

Occupation number-based energy functional for nuclear masses

We develop an energy functional with shell-model occupations as the relevant degrees of freedom and compute nuclear masses across the nuclear chart. The functional is based on Hohenberg-Kohn theory with phenomenologically motivated terms. A global fit of the 17-parameter functional to nuclear masses yields a root-mean-square deviation of χ= 1.31 MeV. Nuclear radii are computed within a model that employs the resulting occupation numbers.

preprint2010arXiv

Ab initio computation of the 17F proton-halo state and resonances in A = 17 nuclei

We perform coupled-cluster calculations of the energies and lifetimes of single-particle states around the doubly magic nucleus 16O based on chiral nucleon-nucleon interactions at next-to-next- to-next-to-leading order. To incorporate effects from the scattering continuum, we solve the coupled- cluster equations with a Gamow-Hartree-Fock basis. Our calculations for the 1/2+ proton-halo state in 17F and the 1/2+ state in 17O agree well with experiment, while the calculated spin-orbit splitting between d5/2 and d3/2 states is too small due to the lack of three-nucleon forces. We find that continuum effects yield a significant amount of additional binding energy for the 1/2+ and 3/2+ states in 17O and 17F.

preprint2010arXiv

Ab initio coupled-cluster approach to nuclear structure with modern nucleon-nucleon interactions

We perform coupled-cluster calculations for the doubly magic nuclei 4He, 16O, 40Ca and 48Ca, for neutron-rich isotopes of oxygen and fluorine, and employ "bare" and secondary renormalized nucleon-nucleon interactions. For the nucleon-nucleon interaction from chiral effective field theory at order next-to-next-to-next-to leading order, we find that the coupled-cluster approximation including triples corrections binds nuclei within 0.4 MeV per nucleon compared to data. We employ interactions from a resolution-scale dependent similarity renormalization group transformations and assess the validity of power counting estimates in medium-mass nuclei. We find that the missing contributions due to three-nucleon forces are consistent with these estimates. For the unitary correlator model potential, we find a slow convergence with respect to increasing the size of the model space. For the G-matrix approach, we find a weak dependence of ground-state energies on the starting energy combined with a rather slow convergence with respect to increasing model spaces. We also analyze the center-of-mass problem and present a practical and efficient solution.

preprint2010arXiv

Computation of spectroscopic factors with the coupled-cluster method

We present a calculation of spectroscopic factors within coupled-cluster theory. Our derivation of algebraic equations for the one-body overlap functions are based on coupled-cluster equation-of-motion solutions for the ground and excited states of the doubly magic nucleus with mass number $A$ and the odd-mass neighbor with mass $A-1$. As a proof-of-principle calculation, we consider $^{16}$O and the odd neighbors $^{15}$O and $^{15}$N, and compute the spectroscopic factor for nucleon removal from $^{16}$O. We employ a renormalized low-momentum interaction of the $V_{\mathrm{low-}k}$ type derived from a chiral interaction at next-to-next-to-next-to-leading order. We study the sensitivity of our results by variation of the momentum cutoff, and then discuss the treatment of the center of mass.

preprint2006arXiv

Level Repulsion in Constrained Gaussian Random-Matrix Ensembles

Introducing sets of constraints, we define new classes of random-matrix ensembles, the constrained Gaussian unitary (CGUE) and the deformed Gaussian unitary (DGUE) ensembles. The latter interpolate between the GUE and the CGUE. We derive a sufficient condition for GUE-type level repulsion to persist in the presence of constraints. For special classes of constraints, we extend this approach to the orthogonal and to the symplectic ensembles. A generalized Fourier theorem relates the spectral properties of the constraining ensembles with those of the constrained ones. We find that in the DGUEs, level repulsion always prevails at a sufficiently short distance and may be lifted only in the limit of strictly enforced constraints.

preprint1998arXiv

On the Special Role of Symmetric Periodic Orbits in a Chaotic System

After early work of Henon it has become folk knowledge that symmetric periodic orbits are of particular importance. We reinforce this belief by additional studies and we further find that invariant closed symplectic submanifolds caused by discrete symmetries prove to be an important complement to the long known role of the orbits. The latter have particular importance in semi-classics. Based on the structural stability of hyperbolic horseshoes we give an argument that opens an avenue to the understanding of these facts.