Search

Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding is approached near the neutron- and proton driplines, traditional shell-structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe this multiscale physics starting from nuclear forces with roots in quantum chromodynamics. Here we achieve this within a unified and non-perturbative framework that captures both short- and long-range correlations starting from modern nucleon-nucleon and three-nucleon forces from chiral effective field theory. The short-range correlations which accounts for the bulk of the binding energy is included within a symmetry-breaking framework, while long-range correlations (and fine details about the collective structure) are included via symmetry projection. Our calculations accurately reproduce available experimental data for low-lying collective states and the electromagnetic quadrupole transitions in $^{20-30}$Ne. We also reveal coexisting spherical and deformed shapes in $^{30}$Ne, which indicates the breakdown of the magic neutron number $N=20$ as the key nucleus $^{28}$O is approached, and we predict that the dripline nuclei $^{32,34}$Ne are strongly deformed. By developing reduced-order-models for symmetry-projected states, we perform a global sensitivity analysis and find that the subleading singlet S-wave contact and a pion-nucleon coupling strongly impact nuclear deformation in chiral effective-field-theory. The techniques developed in this work clarify how microscopic nuclear forces generate the multiscale physics of nuclei spanning collective phenomena as well as short-range correlations and allow to capture emergent and dynamical phenomena in finite fermion systems., Comment: This work significantly expands upon ideas from arXiv:2305.06955. While it includes some overlap in text to maintain context, substantial new content, methodology, and analysis are presented. 26 pages, 16 figures

We construct efficient emulators for the \emph{ab initio} computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for $> 10^6$ different parametrizations of a chiral interaction model with explicit $\Delta$-isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to \nuc{4}{He}, we use history matching to identify non-implausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling/importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties., Comment: 16 pages, 15 figures, updated Supplemental Material can be found in the PRC version

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation-of-state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of $\Delta$-full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including \nuc{16}{O} observables gives saturation predictions that are more precise than those that only use few-body data., Comment: 5 pages, 3 figures, updated Supplemental Material can be found in PRC version

Ab initio has been used as a label in nuclear theory for over two decades. Its meaning has evolved and broadened over the years. We present our interpretation, briefly review its historical use, and discuss its present-day relation to theoretical uncertainty quantification., Comment: Contribution to Research Topic "Uncertainty Quantification in Nuclear Physics" of Frontiers in Physics

The energy levels of light hypernuclei are experimentally accessible observables that contain valuable information about the interaction between hyperons and nucleons. In this work we study strangeness $S = -1$ systems $^{3,4}_\Lambda$H and $^{4,5}_\Lambda$He using the ab initio no-core shell model (NCSM) with realistic interactions obtained from chiral effective field theory ($\chi$EFT). In particular, we quantify the finite precision of theoretical predictions that can be attributed to nuclear physics uncertainties. We study both the convergence of the solution of the many-body problem (method uncertainty) and the regulator- and calibration data-dependence of the nuclear $\chi$EFT Hamiltonian (model uncertainty). For the former, we implement infrared correction formulas and extrapolate finite-space NCSM results to infinite model space. We then use Bayesian parameter estimation to quantify the resulting method uncertainties. For the latter, we employ a family of 42 realistic Hamiltonians and measure the standard deviation of predictions while keeping the leading-order hyperon-nucleon interaction fixed. Following this procedure we find that model uncertainties of ground-state $\Lambda$ separation energies amount to $\sim 20(100)$ keV in $^3_\Lambda$H($^4_\Lambda$H,He) and $\sim 400$ keV in $^5_\Lambda$He. Method uncertainties are comparable in magnitude for the $^4_\Lambda$H,He $1^+$ excited states and $^5_\Lambda$He, which are computed in limited model spaces, but otherwise much smaller. This knowledge of expected theoretical precision is crucial for the use of binding energies of light hypernuclei to infer the elusive hyperon-nucleon interaction., Comment: 16 pages with 8 figures

We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces -- which are conventionally regarded as next-to-next-to-leading order -- should accompany the two-nucleon force already at leading order (LO) starting with mass number $A\approx 10-20$. We find that this promotion enables the first realistic description of the $^{16}$O ground state based on a renormalization-group-invariant LO interaction. We also performed coupled-cluster calculations of the equation of state for symmetric nuclear matter and our results indicate that LO four-nucleon forces could play a crucial role for describing heavy-mass nuclei. The enhancement mechanism we found is very general and could be important also in other many-body problems., Comment: 6 pages, 4 figures

We make ab initio predictions for the A = 6 nuclear level scheme based on two- and three-nucleon interactions up to next-to-next-to-leading order in chiral effective field theory ($\chi$EFT). We utilize eigenvector continuation and Bayesian methods to quantify uncertainties stemming from the many-body method, the $\chi$EFT truncation, and the low-energy constants of the nuclear interaction. The construction and validation of emulators is made possible via the development of JupiterNCSM -- a new M-scheme no-core shell model code that uses on-the-fly Hamiltonian matrix construction for efficient, single-node computations up to $N_\mathrm{max} = 10$ for ${}^{6}\mathrm{Li}$. We find a slight underbinding of ${}^{6}\mathrm{He}$ and ${}^{6}\mathrm{Li}$, although consistent with experimental data given our theoretical error bars. As a result of incorporating a correlated $\chi$EFT-truncation errors we find more precise predictions (smaller error bars) for separation energies: $S_d({}^{6}\mathrm{Li}) = 0.89 \pm 0.44$ MeV, $S_{2n}({}^{6}\mathrm{He}) = 0.20 \pm 0.60$ MeV, and for the beta decay Q-value: $Q_{\beta^-}({}^{6}\mathrm{He}) = 3.71 \pm 0.65$ MeV. We conclude that our error bars can potentially be reduced further by extending the model space used by JupiterNCSM., Comment: 13 pages, 9 figures

Normal-ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been systematically quantified. In this work we study in particular the normal-ordering two-body approximation applied to a single harmonic-oscillator reference state. We explicate the breaking of translational invariance and demonstrate the magnitude of the approximation error as a function of model space parameters for $^4\rm{He}$ and $^{16}\rm{O}$ by performing full no-core shell-model calculations with and without three-nucleon forces. We combine two different diagnostics to better monitor the breaking of translational invariance. While the center-of-mass effect is shown to become potentially very large for $^4\rm{He}$, it is also shown to be much smaller for $^{16}\rm{O}$ although full convergence is not reached. These tools can be easily implemented in studies using other many-body frameworks and bases., Comment: 13 pages, 10 figures

We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($\chi$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 $\chi$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 $\beta$-decay rate can only be obtained if $\chi$EFT truncation errors are included in the analysis. All of these except the $\beta$-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 $\chi$EFT expansion parameter of $Q=0.33 \pm 0.06$ for these observables., Comment: 15 pages, 7 figures; same as published version

Nuclear charge radii are sensitive probes of different aspects of the nucleon-nucleon interaction and the bulk properties of nuclear matter; thus, they provide a stringent test and challenge for nuclear theory. The calcium region has been of particular interest, as experimental evidence has suggested a new magic number at $N = 32$ [1-3], while the unexpectedly large increases in the charge radii [4,5] open new questions about the evolution of nuclear size in neutron-rich systems. By combining the collinear resonance ionization spectroscopy method with $\beta$-decay detection, we were able to extend the charge radii measurement of potassium ($Z =19$) isotopes up to the exotic $^{52}$K ($t_{1/2}$ = 110 ms), produced in minute quantities. Our work provides the first charge radii measurement beyond $N = 32$ in the region, revealing no signature of the magic character at this neutron number. The results are interpreted with two state-of-the-art nuclear theories. For the first time, a long sequence of isotopes could be calculated with coupled-cluster calculations based on newly developed nuclear interactions. The strong increase in the charge radii beyond $N = 28$ is not well captured by these calculations, but is well reproduced by Fayans nuclear density functional theory, which, however, overestimates the odd-even staggering effect. These findings highlight our limited understanding on the nuclear size of neutron-rich systems, and expose pressing problems that are present in some of the best current models of nuclear theory., Comment: submitted version; revision accepted in Nature Physics

Chiral effective field theory ($\chi$EFT), as originally proposed by Weinberg, promises a theoretical connection between low-energy nuclear interactions and quantum chromodynamics (QCD). However, the important property of renormalization-group (RG) invariance is not fulfilled in current implementations and its consequences for predicting atomic nuclei beyond two- and three-nucleon systems has remained unknown. In this work we present a first and systematic study of recent RG-invariant formulations of $\chi$EFT and their predictions for the binding energies and other observables of selected nuclear systems with mass-numbers up to $A =16$. Specifically, we have carried out ab initio no-core shell-model and coupled cluster calculations of the ground-state energy of $^3$H, $^{3,4}$He, $^{6}$Li, and $^{16}$O using several recent power-counting (PC) schemes at leading order (LO) and next-to-leading order (NLO), where the subleading interactions are treated in perturbation theory. Our calculations indicate that RG-invariant and realistic predictions can be obtained for nuclei with mass number $A \leq 4$. We find, however, that $^{16}$O is either unbound with respect to the four $\alpha$-particle threshold, or deformed, or both. Similarly, we find that the $^{6}$Li ground-state resides above the $\alpha$-deuteron separation threshold. These results are in stark contrast with experimental data and point to either necessary fine-tuning of all relevant counterterms, or that current state-of-the-art RG-invariant PC schemes at LO in $\chi$EFT lack necessary diagrams -- such as three-nucleon forces -- to realistically describe nuclei with mass number $A>4$., Comment: 18 pages, 12 figures, published version

We optimize $\Delta$-full nuclear interactions from chiral effective field theory. The low-energy constants of the contact potentials are constrained by two-body scattering phase shifts, and by properties of bound-state of $A=2$ to $4$ nucleon systems and nuclear matter. The pion-nucleon couplings are taken from a Roy-Steiner analysis. The resulting interactions yield accurate binding energies and radii for a range of nuclei from $A=16$ to $A=132$, and provide accurate equations of state for nuclear matter and realistic symmetry energies. Selected excited states are also in agreement with data., Comment: 7 pages, 6 figures

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., Comment: Published version

Theoretical models of the strong nuclear interaction contain unknown coupling constants (parameters) that must be determined using a pool of calibration data. In cases where the models are complex, leading to time consuming calculations, it is particularly challenging to systematically search the corresponding parameter domain for the best fit to the data. In this paper, we explore the prospect of applying Bayesian optimization to constrain the coupling constants in chiral effective field theory descriptions of the nuclear interaction. We find that Bayesian optimization performs rather well with low-dimensional parameter domains and foresee that it can be particularly useful for optimization of a smaller set of coupling constants. A specific example could be the determination of leading three-nucleon forces using data from finite nuclei or three-nucleon scattering experiments., Comment: 33 pages, 14 figures

Ab initio methods aim to solve the nuclear many-body problem with controlled approximations. Virtually exact numerical solutions for realistic interactions can only be obtained for certain special cases such as few-nucleon systems. Here we extend the reach of exact diagonalization methods to handle model spaces with dimension exceeding $10^{10}$ on a single compute node. This allows us to perform no-core shell model (NCSM) calculations for 6Li in model spaces up to $N_\mathrm{max} = 22$ and to reveal the 4He+d halo structure of this nucleus. Still, the use of a finite harmonic-oscillator basis implies truncations in both infrared (IR) and ultraviolet (UV) length scales. These truncations impose finite-size corrections on observables computed in this basis. We perform IR extrapolations of energies and radii computed in the NCSM and with the coupled-cluster method at several fixed UV cutoffs. It is shown that this strategy enables information gain also from data that is not fully UV converged. IR extrapolations improve the accuracy of relevant bound-state observables for a range of UV cutoffs, thus making them profitable tools. We relate the momentum scale that governs the exponential IR convergence to the threshold energy for the first open decay channel. Using large-scale NCSM calculations we numerically verify this small-momentum scale of finite nuclei., Comment: Minor revisions.Accepted for publication in Physical Review C

We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($\chi$EFT) up to next-to-next-to-leading order (NNLO) and perform a rigorous uncertainty analysis of the results. We quantify the uncertainties due to (i) the computational method used to compute the $pp$ cross section in momentum space, (ii) the statistical uncertainties in the low-energy coupling constants of $\chi$EFT, (iii) the systematic uncertainty due to the $\chi$EFT cutoff, and (iv) systematic variations in the database used to calibrate the nucleon-nucleon interaction. We also examine the robustness of the polynomial extrapolation procedure, which is commonly used to extract the threshold $S$-factor and its energy-derivatives. By performing a statistical analysis of the polynomial fit of the energy-dependent $S$-factor at several different energy intervals, we eliminate a systematic uncertainty that can arise from the choice of the fit interval in our calculations. In addition, we explore the statistical correlations between the $S$-factor and few-nucleon observables such as the binding energies and point-proton radii of $^{2,3}$H and $^3$He as well as the $D$-state probability and quadrupole moment of $^2$H, and the $\beta$-decay of $^{3}$H. We find that, with the state-of-the-art optimization of the nuclear Hamiltonian, the statistical uncertainty in the threshold $S$-factor cannot be reduced beyond 0.7%., Comment: peer-reviewed version: some passages modified, typographic errors fixed

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.

Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to next-to-next-to-leading order. Our optimization protocol corresponds to a simultaneous fit to scattering and bound-state observables in the pion-nucleon, nucleon-nucleon, and few-nucleon sectors, thereby utilizing the full model capabilities of chiEFT. We study the effect on other observables by demonstrating error-propagation methods that can easily be adopted by future works. We employ mathematical optimization and implement automatic differentiation to attain efficient and machine-precise first- and second-order derivatives of the objective function with respect to the LECs. We use power-counting arguments to estimate the systematic uncertainty that is inherent to chiEFT and we construct chiral interactions at different orders with quantified uncertainties. Statistical error propagation is compared with Monte Carlo sampling showing that statistical errors are in general small compared to systematic ones. In conclusion, we find that a simultaneous fit to different sets of data is critical to (i) identify the optimal set of LECs, (ii) capture all relevant correlations, (iii) reduce the statistical uncertainty, and (iv) attain order-by-order convergence in chiEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector; in particlar when varying the cutoff in the chiral potentials. The methodology and results presented in this Paper open a new frontier for uncertainty quantification in ab initio nuclear theory., Comment: Accepted for publication in Physical Review X. This version includes Supplemental Material

We present an efficient implementation for the evaluation of Wigner 3j, 6j, and 9j symbols. These represent numerical transformation coefficients that are used in the quantum theory of angular momentum. They can be expressed as sums and square roots of ratios of integers. The integers can be very large due to factorials. We avoid numerical precision loss due to cancellation through the use of multi-word integer arithmetic for exact accumulation of all sums. A fixed relative accuracy is maintained as the limited number of floating-point operations in the final step only incur rounding errors in the least significant bits. Time spent to evaluate large multi-word integers is in turn reduced by using explicit prime factorisation of the ingoing factorials, thereby improving execution speed. Comparison with existing routines shows the efficiency of our approach and we therefore provide a computer code based on this work., Comment: 7 pages, 2 figures. Accepted for publication in SIAM Journal on Scientific Computing (SISC)

Few- and many-fermion systems on the verge of stability, and consisting of strongly interacting particles, appear in many areas of physics. The theoretical modeling of such systems is a very difficult problem. In this work we present a theoretical framework that is based on the rigged Hilbert space formulation. The few-body problem is solved by exact diagonalization using a basis in which bound, resonant, and non-resonant scattering states are included on an equal footing. Current experiments with ultracold atoms offer a fascinating opportunity to study universal properties of few-body systems with a high degree of control over parameters such as the external trap geometry, the number of particles, and even the interaction strength. In particular, particles can be allowed to tunnel out of the trap by applying a magnetic-field gradient that effectively lowers the potential barrier. The result is a tunable open quantum system that allows detailed studies of the tunneling mechanism. In this Contribution we introduce our method and present results for the decay rate of two distinguishable fermions in a one-dimensional trap as a function of the interaction strength. We also study the numerical convergence. Many of these results have been previously published (R. Lundmark, C. Forss\'en, and J. Rotureau, arXiv: 1412.7175). However, in this Contribution we present several technical and numerical details of our approach for the first time., Comment: 8 pages, 4 figures, contribution to the conference proceedings of the 7th International Workshop on the "Dynamics of Critically Stable Quantum Few-Body Systems"

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., Comment: 8 pages, 4 figures. Updated to address issues fixed during publication. We have prepared a python (v2.7) script that will compute the IR scales (both full scale and intrinsic scale). This script depends only on numpy and scipy and can be found at http://github.com/KyleWendt/ncsm_ho_ir_scale

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^\pi=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.

The creation of tunable open quantum systems is becoming feasible in current experiments with ultracold atoms in low-dimensional traps. In particular, the high degree of experimental control over these systems allows detailed studies of tunneling dynamics, e.g., as a function of the trapping geometry and the interparticle interaction strength. In order to address this exciting opportunity we present a theoretical framework for two-body tunneling based on the rigged Hilbert space formulation. In this approach, bound, resonant and scattering states are included on an equal footing, and we argue that the coupling of all these components is vital for a correct description of the relevant threshold phenomena. In particular, we study the tunneling mechanism for two-body systems in one-dimensional traps and different interaction regimes. We find a strong dominance of sequential tunneling of single particles for repulsive and weakly attractive systems, while there is a signature of correlated pair tunneling in the calculated many-particle flux for strongly attractive interparticle interaction., Comment: To be published in Phys. Rev. A (Rapid Communication)

Strongly interacting one-dimensional quantum systems often behave in a manner that is distinctly different from their higher-dimensional counterparts. When a particle attempts to move in a one-dimensional environment it will unavoidably have to interact and 'push' other particles in order to execute a pattern of motion, irrespective of whether the particles are fermions or bosons. A present frontier in both theory and experiment are mixed systems of different species and/or particles with multiple internal degrees of freedom. Here we consider trapped two-component bosons with short-range inter-species interactions much larger than their intra-species interactions and show that they have novel energetic and magnetic properties. In the strongly interacting regime, these systems have energies that are fractions of the basic harmonic oscillator trap quantum and have spatially separated ground states with manifestly ferromagnetic wave functions. Furthermore, we predict excited states that have perfect antiferromagnetic ordering. This holds for both balanced and imbalanced systems, and we show that it is a generic feature as one crosses from few- to many-body systems., Comment: 9 pages, 6 figures

We have quantified the statistical uncertainties of the low-energy coupling-constants (LECs) of an optimized nucleon-nucleon (NN) interaction from chiral effective field theory ($\chi$EFT) at next-to-next-to-leading order (NNLO). In addition, we have propagated the impact of the uncertainties of the LECs to two-nucleon scattering phase shifts, effective range parameters, and deuteron observables., Comment: Submitted to JPhysG Focus issue "Enhancing the interaction between nuclear experiment and theory through information and statistics" (ISNET) (updated author list)

The nature of strongly interacting Fermi gases and magnetism is one of the most important and studied topics in condensed-matter physics. Still, there are many open questions. A central issue is under what circumstances strong short-range repulsive interactions are enough to drive magnetic correlations. Recent progress in the field of cold atomic gases allows to address this question in very clean systems where both particle numbers, interactions and dimensionality can be tuned. Here we study fermionic few-body systems in a one dimensional harmonic trap using a new rapidly converging effective-interaction technique, plus a novel analytical approach. This allows us to calculate the properties of a single spin-down atom interacting with a number of spin-up particles, a case of much recent experimental interest. Our findings indicate that, in the strongly interacting limit, spin-up and spin-down particles want to separate in the trap, which we interpret as a microscopic precursor of one-dimensional ferromagnetism in imbalanced systems. Our predictions are directly addressable in current experiments on ultracold atomic few-body systems., Comment: 12 pages, 6 figures, published version including two appendices on our new numerical and analytical approach

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 \chi^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., Comment: 6 pages, 4 figures

A first-principles description of nuclear systems along the drip lines presents a substantial theoretical and computational challenge. In this paper, we discuss the nuclear theory roadmap, some of the key theoretical approaches, and present selected results with a focus on long isotopic chains. An important conclusion, which consistently emerges from these theoretical analyses, is that three-nucleon forces are crucial for both global nuclear properties and detailed nuclear structure, and that many-body correlations due to the coupling to the particle continuum are essential as one approaches particle drip lines. In the quest for a comprehensive nuclear theory, high performance computing plays a key role., Comment: Contribution to proceedings of Nobel Symposium 152: Physics with radioactive beams, June 2012, Gothenburg, Sweden

We report on recent progress of the implementation of the similarity renormalization group (SRG) for three-body interactions in a one-dimensional, bosonic model system using the plane wave basis. We discuss our implementation of the flow equations and show results that confirm that results in the three-body sector remain unchanged by the transformation of the Hamiltonian. We also show how the SRG transformation decouples low- from high-momentum nodes in the three-body sector and therefore simplifies the numerical calculation of observables., Comment: 7 pages, 5 figures, svjour style

Spectra of the 9Be(e,e') reaction have been measured at the S-DALINAC at an electron energy E_0 = 73 MeV and scattering angles of 93{\deg} and 141{\deg} with high energy resolution up to excitation energies E_x = 8 MeV. The astrophysically relevant resonance parameters of the first excited 1/2+ state of 9Be have been extracted in a one-level approximation of R-matrix theory resulting in a resonance energy E_R = 1.748(6) MeV and width Gamma_R = 274(8) keV in good agreement with the latest 9Be(gamma,n) experiment but with considerably improved uncertainties. However, the reduced B(E1) transition strength deduced from an extrapolation of the (e,e') data to the photon point is a factor of two smaller. Implications of the new results for a possible production of 12C in neutron-rich astrophysical scenarios are discussed., Comment: 8 pages, 7 figures, accepted for publication in Phys. Rev. C

The neutron capture cross section of 14C is of relevance for several nucleosynthesis scenarios such as inhomogeneous Big Bang models, neutron induced CNO cycles, and neutrino driven wind models for the r process. The 14C(n,g) reaction is also important for the validation of the Coulomb dissociation method, where the (n,g) cross section can be indirectly obtained via the time-reversed process. So far, the example of 14C is the only case with neutrons where both, direct measurement and indirect Coulomb dissociation, have been applied. Unfortunately, the interpretation is obscured by discrepancies between several experiments and theory. Therefore, we report on new direct measurements of the 14C(n,g) reaction with neutron energies ranging from 20 to 800 keV.

This contribution reviews a number of applications of the ab initio no-core shell model (NCSM) within nuclear physics and beyond. We will highlight a nuclear-structure study of the A = 12 isobar using a chiral NN + 3NF interaction. In the spirit of this workshop we will also mention the new development of the NCSM formalism to describe open channels and to approach the problem of nuclear reactions. Finally, we will illustrate the universality of the many-body problem by presenting the recent adaptation of the NCSM effective-interaction approach to study the many-boson problem in an external trapping potential with short-range interactions., Comment: 4 pages. Article based on the presentation by C. Forssen at the Fifth Workshop on Critical Stability, Erice, Sicily. Published in Few-Body Systems

Recently, charge radii and ground-state electromagnetic moments of Li and Be isotopes were measured precisely. We have performed large-scale ab initio no-core shell model calculations for these isotopes using high-precision nucleon-nucleon potentials. The isotopic trends of our computed charge radii and quadrupole and magnetic-dipole moments are in good agreement with experimental results with the exception of the 11Li charge radius. The magnetic moments are in particular well described, whereas the absolute magnitudes of the quadrupole moments are about 10% too small. The small magnitude of the 6Li quadrupole moment is reproduced, and with the CD-Bonn NN potential, also its correct sign., Comment: 6 pages, 3 figures, 3 tables (published manuscript)

We show that the convergence behavior of the many-body numerical diagonalization scheme for strongly interacting bosons in a trap can be significantly improved by the Lee-Suzuki method adapted from nuclear physics: One can construct an effective interaction that acts in a space much smaller than the original Hilbert space. In particular for short-ranged forces and strong correlations, the method offers a good estimate of the energy and the excitation spectrum, at a computational cost several orders of magnitude smaller than that required by the standard method., Comment: 5 pages, 4 figures

We demonstrate the existence of multiple converging sequences in the ab initio no-core shell model. By examining the underlying theory of effective operators, we expose the physical foundations for the alternative pathways to convergence. This leads us to propose a revised strategy for evaluating effective interactions for $A$-body calculations in restricted model spaces. We suggest that this strategy is particularly useful for applications to nuclear processes in which states of both parities are used simultaneously, such as for transition rates. We demonstrate the utility of our strategy with large-scale calculations in light nuclei.

Indirect methods play an important role in the determination of nuclear reaction cross sections that are hard to measure directly. In this paper we investigate the feasibility of using the so-called surrogate method to extract neutron-capture cross sections for low energy compound-nuclear reactions in spherical and near-spherical nuclei. We present the surrogate method and develop a statistical nuclear-reaction simulation to explore different approaches to utilize surrogate reaction data. We assess the success of each approach by comparing the extracted cross sections with a predetermined benchmark. In particular, we employ regional systematics of nuclear properties in the 34 <= Z <= 46 region to calculate (n,gamma) cross sections for a series of Zr isotopes, and to simulate a surrogate experiment and the extraction of the desired cross section. We identify one particular approach that may provide very useful estimates of the cross section, and we discuss some of the limitations of the method. General recommendations for future (surrogate) experiments are also given., Comment: 17 pages, 9 figures. Accepted for publication in Phys. Rev. C

We are presenting the first ab initio structure investigation of the loosely bound 11-Be nucleus, together with a study of the lighter isotope 9-Be. The nuclear structure of these isotopes is particularly interesting due to the appearance of a parity-inverted ground state in 11-Be. Our study is performed in the framework of the ab initio no-core shell model. Results obtained using four different, high-precision two-nucleon interactions, in model spaces up to 9\hbar\Omega, are shown. For both nuclei, and all potentials, we reach convergence in the level ordering of positive- and negative-parity spectra separately. Concerning their relative position, the positive-parity states are always too high in excitation energy, but a fast drop with respect to the negative-parity spectrum is observed when the model space is increased. This behavior is most dramatic for 11-Be. In the largest model space we were able to reach, the 1/2+ level has dropped down to become either the first or the second excited state, depending on which interaction we use. We also observe a contrasting behavior in the convergence patterns for different two-nucleon potentials, and argue that a three-nucleon interaction is needed to explain the parity inversion. Furthermore, large-basis calculations of 13-C and 11-B are performed. This allows us to study the systematics of the position of the first unnatural-parity state in the N=7 isotone and the A=11 isobar. The 11-B run in the 9\hbar\Omega model space involves a matrix with dimension exceeding 1.1 x 10^9, and is our largest calculation so far. We present results on binding energies, excitation spectra, level configurations, radii, electromagnetic observables, and 10-Be+n overlap functions., Comment: 17 pages, 12 figures To be published in Phys. Rev. C Resubmitted version. Minor changes

Electromagnetic processes in loosely bound nuclei are investigated using an analytical model. In particular, electromagnetic dissociation of $^8$B is studied and the results of our analytical model are compared to numerical calculations based on a three-body picture of the $^8$B bound state. The calculation of energy spectra is shown to be strongly model dependent. This is demonstrated by investigating the sensitivity to the rms intercluster distance, the few-body behavior, and the effects of final state interaction. In contrast, the fraction of the energy spectrum which can be attributed to E1 transitions is found to be almost model independent at small relative energies. This finding is of great importance for astrophysical applications as it provides us with a new tool to extract the E1 component from measured energy spectra. An additional, and independent, method is also proposed as it is demonstrated how two sets of experimental data, obtained with different beam energy and/or minimum impact parameter, can be used to extract the E1 component., Comment: Submitted to Phys. Rev. C. 10 pages, 7 figures

In this paper we develop an analytical model in order to study electromagnetic processes involving loosely bound neutron--rich and proton--rich nuclei. We construct a model wave function, to describe loosely bound few--body systems, having the correct behaviour both at large and small distances. The continuum states are approximated by regular Coulomb functions. As a test case we consider the two--body Coulomb dissociation of 8B and, the inverse, radiative capture reaction. The difference between using a pure two--body model and the results obtained when incorporating many--body effects, is investigated. We conclude that the interpretation of experimental data is highly model dependent and stress the importance of measuring few--body channels., Comment: Accepted for publication in Physics Letters B. Added a comparison with a potential model calculation in Fig. 2

An analytical model, recently developed to study the electromagnetic dissociation (EMD) energy spectra of two-neutron halo nuclei, is applied to 11-Li and 14-Be. We find that a reliable set of experimental data on EMD could help to resolve the problem concerning the structure of the 11-Li ground state. For 14-Be we find a mutual inconsistency between the existing experimental data on the binding energy, radius and EMD energy spectrum. We also conclude that the structure of 14-Be is essentially more complicated than the structure of other Borromean two--neutron halo nuclei., Comment: 17 pages, 2 figures See also previous paper in Nucl. Phys. A 697(2002)p.639 Revised version posted 020309 with 2 new references and a few revisions

An analytical model is developed to study the spectra of electromagnetic dissociation of two-neutron halo nuclei without precise knowledge about initial and final states. Phenomenological three-cluster bound state wave functions, reproducing the most relevant features of these nuclei, are used along with no interaction final states. The 6-He nucleus is considered as a test case, and a good agreement with experimental data concerning the shape of the spectrum and the magnitude of the strength function is found., Comment: 19 pages, 4 figures Accepted for publishing in Nuclear Physics A