BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

C. Forssén

C. Forssén appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thnucl-excond-mat.quant-gasastro-ph.SRphysics.comp-phphysics.data-anquant-phcond-mat.str-elhep-phmath-phmath.MPphysics.atom-ph

Catalog footprint

What is connected

17works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fast & rigorous predictions for $A=6$ nuclei with Bayesian posterior sampling

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 ($χ$EFT). We utilize eigenvector continuation and Bayesian methods to quantify uncertainties stemming from the many-body method, the $χ$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 $χ$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_{β^-}({}^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Uncertainty analysis and order-by-order optimization of chiral nuclear interactions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Uncertainty quantification for proton-proton fusion in chiral effective field theory

We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($χ$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 $χ$EFT, (iii) the systematic uncertainty due to the $χ$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 $β$-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%.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast and accurate evaluation of Wigner 3j, 6j, and 9j symbols using prime factorisation and multi-word integer arithmetic

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strongly interacting few-fermion systems in a trap

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én, and J. Rotureau, arXiv: 1412.7175). However, in this Contribution we present several technical and numerical details of our approach for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tunneling Theory for Tunable Open Quantum Systems of Ultracold Atoms in One-Dimensional Traps

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fermionization of two-component few-fermion systems in a one-dimensional harmonic trap

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum magnetism in strongly interacting one-dimensional spinor Bose systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical uncertainties of a chiral interaction at next-to-next-to leading order

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 ($χ$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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Similarity Renormalization Group for Three-Body Interactions in One Dimension

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Resonance parameters of the first 1/2+ state in 9Be and astrophysical implications

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° and 141° 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.

0 citations0 reviewsNo code linkedOpen access