-
Spectroscopy of deeply bound orbitals in neutron-rich Ca isotopes
Authors:
P. J. Li,
J. Lee,
P. Doornenbal,
S. Chen,
S. Wang,
A. Obertelli,
Y. Chazono,
J. D. Holt,
B. S. Hu,
K. Ogata,
Y. Utsuno,
K. Yoshida,
N. L. Achouri,
H. Baba,
F. Browne,
D. Calvet,
F. Château,
N. Chiga,
A. Corsi,
M. L. Cortés,
A. Delbart,
J-M. Gheller,
A. Giganon,
A. Gillibert,
C. Hilaire
, et al. (63 additional authors not shown)
Abstract:
The calcium isotopes are an ideal system to investigate the evolution of shell structure and magic numbers. Although the properties of surface nucleons in calcium have been well studied, probing the structure of deeply bound nucleons remains a challenge. Here, we report on the first measurement of unbound states in $^{53}$Ca and $^{55}$Ca, populated from \ts{54,56}Ca($p,pn$) reactions at a beam en…
▽ More
The calcium isotopes are an ideal system to investigate the evolution of shell structure and magic numbers. Although the properties of surface nucleons in calcium have been well studied, probing the structure of deeply bound nucleons remains a challenge. Here, we report on the first measurement of unbound states in $^{53}$Ca and $^{55}$Ca, populated from \ts{54,56}Ca($p,pn$) reactions at a beam energy of around 216 MeV/nucleon at the RIKEN Radioactive Isotopes Beam Factory. The resonance properties, partial cross sections, and momentum distributions of these unbound states were analyzed. Orbital angular momentum $l$ assignments were extracted from momentum distributions based on calculations using the distorted wave impulse approximation (DWIA) reaction model. The resonances at excitation energies of 5516(41)\,keV in $^{53}$Ca and 6000(250)\,keV in $^{55}$Ca indicate a significant $l$\, =\,3 component, providing the first experimental evidence for the $ν0f_{7/2}$ single-particle strength of unbound hole states in the neutron-rich Ca isotopes. The observed excitation energies and cross-sections point towards extremely localized and well separated strength distributions, with some fragmentation for the $ν0f_{7/2}$ orbital in $^{55}$Ca. These results are in good agreement with predictions from shell-model calculations using the effective GXPF1Bs interaction and \textit{ab initio} calculations and diverge markedly from the experimental distributions in the nickel isotones at $Z=28$.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Ab initio computations of strongly deformed nuclei around $^{80}$Zr
Authors:
B. S. Hu,
Z. H. Sun,
G. Hagen,
T. Papenbrock
Abstract:
Nuclei around $N\approx Z\approx 40$ are strongly deformed and exhibit coexistence of shapes. These phenomena have challenged nuclear models. Here we perform ab initio coupled-cluster computations of low-lying collective states and electromagnetic quadrupole transitions of the even-even nuclei $^{72}$Kr, $^{76,78}$Sr, $^{78,80}$Zr and $^{84}$Mo starting from chiral nucleon-nucleon and three-nucleo…
▽ More
Nuclei around $N\approx Z\approx 40$ are strongly deformed and exhibit coexistence of shapes. These phenomena have challenged nuclear models. Here we perform ab initio coupled-cluster computations of low-lying collective states and electromagnetic quadrupole transitions of the even-even nuclei $^{72}$Kr, $^{76,78}$Sr, $^{78,80}$Zr and $^{84}$Mo starting from chiral nucleon-nucleon and three-nucleon forces. Our calculations reproduce the coexistence of oblate and prolate shapes in these nuclei, yield rotational bands and strong electromagnetic transitions, but are not accurate for some observables and nuclei. These results highlight the advances and challenges of ab initio computations of heavy deformed nuclei.
△ Less
Submitted 8 May, 2024;
originally announced May 2024.
-
Surface Normal Estimation with Transformers
Authors:
Barry Shichen Hu,
Siyun Liang,
Johannes Paetzold,
Huy H. Nguyen,
Isao Echizen,
Jiapeng Tang
Abstract:
We propose the use of a Transformer to accurately predict normals from point clouds with noise and density variations. Previous learning-based methods utilize PointNet variants to explicitly extract multi-scale features at different input scales, then focus on a surface fitting method by which local point cloud neighborhoods are fitted to a geometric surface approximated by either a polynomial fun…
▽ More
We propose the use of a Transformer to accurately predict normals from point clouds with noise and density variations. Previous learning-based methods utilize PointNet variants to explicitly extract multi-scale features at different input scales, then focus on a surface fitting method by which local point cloud neighborhoods are fitted to a geometric surface approximated by either a polynomial function or a multi-layer perceptron (MLP). However, fitting surfaces to fixed-order polynomial functions can suffer from overfitting or underfitting, and learning MLP-represented hyper-surfaces requires pre-generated per-point weights. To avoid these limitations, we first unify the design choices in previous works and then propose a simplified Transformer-based model to extract richer and more robust geometric features for the surface normal estimation task. Through extensive experiments, we demonstrate that our Transformer-based method achieves state-of-the-art performance on both the synthetic shape dataset PCPNet, and the real-world indoor scene dataset SceneNN, exhibiting more noise-resilient behavior and significantly faster inference. Most importantly, we demonstrate that the sophisticated hand-designed modules in existing works are not necessary to excel at the task of surface normal estimation.
△ Less
Submitted 11 January, 2024;
originally announced January 2024.
-
Magnetic dipole transition in $^{48}$Ca
Authors:
B. Acharya,
B. S. Hu,
S. Bacca,
G. Hagen,
P. Navrátil,
T. Papenbrock
Abstract:
The magnetic dipole transition strength $B(M1)$ of $^{48}$Ca is dominated by a single resonant state at an excitation energy of 10.23 MeV. Experiments disagree about $B(M1)$ and this impacts our understanding of spin flips in nuclei. We performed ab initio computations based on chiral effective field theory and found that $B(M1:0^+\rightarrow1^+)$ lies in the range from $7.0$ to $10.2~μ_N^2$. This…
▽ More
The magnetic dipole transition strength $B(M1)$ of $^{48}$Ca is dominated by a single resonant state at an excitation energy of 10.23 MeV. Experiments disagree about $B(M1)$ and this impacts our understanding of spin flips in nuclei. We performed ab initio computations based on chiral effective field theory and found that $B(M1:0^+\rightarrow1^+)$ lies in the range from $7.0$ to $10.2~μ_N^2$. This is consistent with a $(γ,n)$ experiment but larger than results from $(e,e^\prime)$ and $(p,p')$ scattering. Two-body currents yield no quenching of the $B(M1)$ strength and continuum effects reduce it by about 10%. For a validation of our approach, we computed magnetic moments in $^{47,49}$Ca and performed benchmark calculations in light nuclei.
△ Less
Submitted 7 June, 2024; v1 submitted 19 November, 2023;
originally announced November 2023.
-
Isomeric excitation energy for $^{99}$In$^{m}$ from mass spectrometry reveals constant trend next to doubly magic $^{100}$Sn
Authors:
L. Nies,
D. Atanasov,
M. Athanasakis-Kaklamanakis,
M. Au,
K. Blaum,
J. Dobaczewski,
B. S. Hu,
J. D. Holt,
J. Karthein,
I. Kulikov,
Yu. A. Litvinov,
D. Lunney,
V. Manea,
T. Miyagi,
M. Mougeot,
L. Schweikhard,
A. Schwenk,
K. Sieja,
F. Wienholtz
Abstract:
The excitation energy of the 1/2$^-$ isomer in $^{99}$In at ${N=50}$ is measured to be 671(37) keV and the mass uncertainty of the 9/2$^+$ ground state is significantly reduced using the ISOLTRAP mass spectrometer at ISOLDE/CERN. The measurements exploit a major improvement in the resolution of the multi-reflection time-of-flight mass spectrometer. The results reveal an intriguing constancy of the…
▽ More
The excitation energy of the 1/2$^-$ isomer in $^{99}$In at ${N=50}$ is measured to be 671(37) keV and the mass uncertainty of the 9/2$^+$ ground state is significantly reduced using the ISOLTRAP mass spectrometer at ISOLDE/CERN. The measurements exploit a major improvement in the resolution of the multi-reflection time-of-flight mass spectrometer. The results reveal an intriguing constancy of the $1/2^-$ isomer excitation energies in neutron-deficient indium that persists down to the $N = 50$ shell closure, even when all neutrons are removed from the valence shell. This trend is used to test large-scale shell model, \textit{ab initio}, and density functional theory calculations. The models have difficulties describing both the isomer excitation energies and ground-state electromagnetic moments along the indium chain.
△ Less
Submitted 3 June, 2023;
originally announced June 2023.
-
Ab initio descriptions of $A=16$ mirror nuclei with resonance and continuum coupling
Authors:
S. Zhang,
F. R. Xu,
J. G. Li,
B. S. Hu,
Z. H. Cheng,
N. Michel,
Y. Z. Ma,
Q. Yuan,
Y. H. Zhang
Abstract:
We have used an {\it ab initio} Gamow shell model to study the isospin symmetry breaking in the $A=16$ mirror nuclei of $^{16}$F, $^{16}$N, $^{16}$Ne and $^{16}$C. Starting from a chiral interaction with two-nucleon force (2NF) at N$^3$LO and three-nucleon force (3NF) at N$^2$LO, a complex-momentum ${\it psd}$-shell Hamiltonian was constructed by employing the many-body perturbation theory in the…
▽ More
We have used an {\it ab initio} Gamow shell model to study the isospin symmetry breaking in the $A=16$ mirror nuclei of $^{16}$F, $^{16}$N, $^{16}$Ne and $^{16}$C. Starting from a chiral interaction with two-nucleon force (2NF) at N$^3$LO and three-nucleon force (3NF) at N$^2$LO, a complex-momentum ${\it psd}$-shell Hamiltonian was constructed by employing the many-body perturbation theory in the Gamow Hartree-Fock basis which includes bound, resonant and continuum states self-consistently. Such an elaborated {\it ab initio} Gamow shell model with both continuum coupling and 3NF included can properly treat the many-body correlations of weakly bound and unbound nuclei. The mirror partners of $^{16}$F and $^{16}$N exhibit different level orders in their excitation spectra, which can be well explained by the inclusion of 3NF in the calculation. The isospin asymmetry between the mirror partners $^{16}$Ne and $^{16}$C was studied in detail by insight into their configuration structures. The interplay between 3NF and the continuum coupling is discussed in the weakly bound and unbound nuclear states.
△ Less
Submitted 27 December, 2023; v1 submitted 5 January, 2023;
originally announced January 2023.
-
Proton distribution radii of $^{16-24}$O: signatures of new shell closures and neutron skin
Authors:
S. Kaur,
R. Kanungo,
W. Horiuchi,
G. Hagen,
J. D. Holt,
B. S. Hu,
T. Miyagi,
T. Suzuki,
F. Ameil,
J. Atkinson,
Y. Ayyad,
S. Bagchi,
D. Cortina-Gil,
I. Dillmann,
A. Estradé,
A. Evdokimov,
F. Farinon,
H. Geissel,
G. Guastalla,
R. Janik,
R. Knöbel,
J. Kurcewicz,
Yu. A. Litvinov,
M. Marta,
M. Mostazo
, et al. (16 additional authors not shown)
Abstract:
The root mean square radii of the proton density distribution in $^{16-24}$O derived from measurements of charge changing cross sections with a carbon target at $\sim$900$A$ MeV together with the matter radii portray thick neutron skin for $^{22 - 24}$O despite $^{22,24}$O being doubly magic. Imprints of the shell closures at $N$ = 14 and 16 are reflected in local minima of their proton radii that…
▽ More
The root mean square radii of the proton density distribution in $^{16-24}$O derived from measurements of charge changing cross sections with a carbon target at $\sim$900$A$ MeV together with the matter radii portray thick neutron skin for $^{22 - 24}$O despite $^{22,24}$O being doubly magic. Imprints of the shell closures at $N$ = 14 and 16 are reflected in local minima of their proton radii that provide evidence for the tensor interaction causing them. The radii agree with it ab initio calculations employing the chiral NNLO$_{\mathrm{sat}}$ interaction, though skin thickness predictions are challenged. Shell model predictions agree well with the data.
△ Less
Submitted 1 September, 2022;
originally announced September 2022.
-
Deformed in-medium similarity renormalization group
Authors:
Q. Yuan,
S. Q. Fan,
B. S. Hu,
J. G. Li,
S. Zhang,
S. M. Wang,
Z. H. Sun,
Y. Z. Ma,
F. R. Xu
Abstract:
We have developed an {\it ab initio} deformed in-medium similarity renormalization group (IMSRG) for open-shell nuclei. This is a single-reference IMSRG in deformed Hartree-Fock (HF) basis. Deformed wave functions are more efficient in describing deformed nuclei. The broken spherical symmetry needs to be restored by angular momentum projection, which is computational expensive. The angular momentu…
▽ More
We have developed an {\it ab initio} deformed in-medium similarity renormalization group (IMSRG) for open-shell nuclei. This is a single-reference IMSRG in deformed Hartree-Fock (HF) basis. Deformed wave functions are more efficient in describing deformed nuclei. The broken spherical symmetry needs to be restored by angular momentum projection, which is computational expensive. The angular momentum mainly capture the static correlations and can be estimated by the projection of the HF state. In this work, we do deformed IMSRG calculation and add the correlation energy from projected HF as a leading order approximation. As the test ground, we have calculated the deformed $^{8,10}\rm Be$ isotopes with the optimized chiral interaction NNLO$_{\rm opt}$. The results are benchmarked with the no-core shell model and valence space IMSRG calculations. Then we systematically investigated the ground-state energies and charge radii of even-even isotopes from light beryllium to medium-mass magnesium. The calculated energies are extrapolated to infinite basis space by an exponential form, and compared with the extrapolated valence-space IMSRG results and experimental data available.
△ Less
Submitted 14 April, 2022;
originally announced April 2022.
-
Electromagnetic moments of scandium isotopes and $N=28$ isotones in the distinctive $0f_{7/2}$ orbit
Authors:
S. W. Bai,
Á. Koszorús,
B. S. Hu,
X. F. Yang,
J. Billowes,
C. L. Binnersley,
M. L. Bissell,
K. Blaum,
P. Campbell,
B. Cheal,
T. E. Cocolios,
R. P. de Groote,
C. S. Devlin,
K. T. Flanagan,
R. F. Garcia Ruiz,
H. Heylen,
J. D. Holt,
A. Kanellakopoulos,
J. Krämer,
V. Lagaki,
B. Maaß,
S. Malbrunot-Ettenauer,
T. Miyagi,
R. Neugart,
G. Neyens
, et al. (9 additional authors not shown)
Abstract:
The electric quadrupole moment of $^{49}$Sc was measured by collinear laser spectroscopy at CERN-ISOLDE to be $Q_{\rm s}=-0.159(8)$ $e$b, and a nearly tenfold improvement in precision was reached for the electromagnetic moments of $^{47,49}$Sc. The single-particle behavior and nucleon-nucleon correlations are investigated with the electromagnetic moments of $Z=21$ isotopes and $N=28$ isotones as v…
▽ More
The electric quadrupole moment of $^{49}$Sc was measured by collinear laser spectroscopy at CERN-ISOLDE to be $Q_{\rm s}=-0.159(8)$ $e$b, and a nearly tenfold improvement in precision was reached for the electromagnetic moments of $^{47,49}$Sc. The single-particle behavior and nucleon-nucleon correlations are investigated with the electromagnetic moments of $Z=21$ isotopes and $N=28$ isotones as valence neutrons and protons fill the distinctive $0f_{7/2}$ orbit, respectively, located between magic numbers, 20 and 28. The experimental data are interpreted with shell-model calculations using an effective interaction, and ab-initio valence-space in-medium similarity renormalization group calculations based on chiral interactions. These results highlight the sensitivity of nuclear electromagnetic moments to different types of nucleon-nucleon correlations, and establish an important benchmark for further developments of theoretical calculations.
△ Less
Submitted 22 March, 2022; v1 submitted 17 March, 2022;
originally announced March 2022.
-
The roles of three-nucleon force and continuum coupling in mirror symmetry breaking of oxygen mass region
Authors:
S. Zhang,
Y. Z. Ma,
J. G. Li,
B. S. Hu,
Q. Yuan,
Z. H. Cheng,
F. R. Xu
Abstract:
With both three-nucleon force and continuum coupling included, we have developed a self-consistent {\it ab initio} Gamow shell model within the Gamow Hartree-Fock (GHF) basis obtained by the realistic interaction itself. With the chiral two-nucleon N$^3$LO and three-nucleon N$^2$LO interactions, the Gamow shell model has been applied to the mirror systems of $Z=8$ neutron-rich isotopes and $N=8$ p…
▽ More
With both three-nucleon force and continuum coupling included, we have developed a self-consistent {\it ab initio} Gamow shell model within the Gamow Hartree-Fock (GHF) basis obtained by the realistic interaction itself. With the chiral two-nucleon N$^3$LO and three-nucleon N$^2$LO interactions, the Gamow shell model has been applied to the mirror systems of $Z=8$ neutron-rich isotopes and $N=8$ proton-rich isotones, giving good agreements with data in binding energies, dripline positions and excitation spectra. The GHF calculated that the $0d_{3/2}$, $1s_{1/2}$ and $1p_{3/2}$ orbitals are resonances. The resonance states and their interplay with nonresonant continua play a crucial role in the descriptions of nuclei around driplines. Excitation spectra and Thomas-Ehrman shifts observed can be better described when both three-nucleon force and continuum coupling are considered in calculations. The three-nucleon force and continuum coupling produce a combined effect on the Thomas-Ehrman shift, e.g., for the ${1/2}^+$ resonance level of $^{19}$Na. The calculations help the understandings of related nuclear astrophysical processes.
△ Less
Submitted 15 February, 2022; v1 submitted 6 December, 2021;
originally announced December 2021.
-
Tensor force role in $β$ decays analyzed within the Gogny-interaction shell model
Authors:
B. Dai,
B. S. Hu,
Y. Z. Ma,
J. G. Li,
S. M. Wang,
C. W. Johnson,
F. R. Xu
Abstract:
Background: The half-life of the famous $^{14}$C $β$ decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. Purpose: We revisit and analyze the role of the tensor force for the $β$ decay of $^{14}$C as well as of neighboring isotopes. Methods: We add a…
▽ More
Background: The half-life of the famous $^{14}$C $β$ decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. Purpose: We revisit and analyze the role of the tensor force for the $β$ decay of $^{14}$C as well as of neighboring isotopes. Methods: We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a $p$-$sd$ model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum $L$ in order to analyze the effects of the tensor force and cross-shell mixing.
Results: The inclusion of the tensor force significantly improves the shell-model calculations of the $β$-decay properties of carbon isotopes. In particular, the anomalously slow $β$ decay of $^{14}$C can be explained by the isospin $T=0$ part of the tensor force, which changes the components of $^{14}$N with the orbital angular momentum $L=0,1$, and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby $β$ decays are improved.
Conclusions: Decomposition of wave function into $L$ components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of $β$-decay matrix elements. Cross-shell mixing also has a visible impact on the $β$-decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.
△ Less
Submitted 5 September, 2021;
originally announced September 2021.
-
Ab initio structure factors for spin-dependent dark matter direct detection
Authors:
B. S. Hu,
J. Padua-Argüelles,
S. Leutheusser,
T. Miyagi,
S. R. Stroberg,
J. D. Holt
Abstract:
We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: $^{19}$F, $^{23}$Na, $^{27}$Al, $^{29}$Si, $^{73}$Ge, $^{127}$I, and $^{129,131}$Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon…
▽ More
We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: $^{19}$F, $^{23}$Na, $^{27}$Al, $^{29}$Si, $^{73}$Ge, $^{127}$I, and $^{129,131}$Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon currents at the one-body level, including effects from axial-vector two-body currents. We then apply the in-medium similarity renormalization group to construct effective valence-space Hamiltonians and consistently transformed operators of nuclear responses. Combining the recent advances of natural orbitals with three-nucleon forces expressed in large spaces, we obtain basis-space converged structure factors even in heavy nuclei. Generally results are consistent with previous calculations, but large uncertainties in $^{127}$I highlight the need for further study.
△ Less
Submitted 22 December, 2021; v1 submitted 1 September, 2021;
originally announced September 2021.
-
Neutron-rich calcium isotopes within realistic Gamow shell model calculations with continuum coupling
Authors:
J. G. Li,
B. S. Hu,
Q. Wu,
Y. Gao,
S. J. Dai,
F. R. Xu
Abstract:
Based on the realistic nuclear force of the high-precision CD-Bonn potential, we have performed comprehensive calculations for neutron-rich calcium isotopes using the Gamow shell model (GSM) which includes resonance and continuum. The realistic GSM calculations produce well binding energies, one- and two-neutron separation energies, predicting that $^{57}$Ca is the heaviest bound odd isotope and…
▽ More
Based on the realistic nuclear force of the high-precision CD-Bonn potential, we have performed comprehensive calculations for neutron-rich calcium isotopes using the Gamow shell model (GSM) which includes resonance and continuum. The realistic GSM calculations produce well binding energies, one- and two-neutron separation energies, predicting that $^{57}$Ca is the heaviest bound odd isotope and $^{70}$Ca is the dripline nucleus. Resonant states are predicted, which provides useful information for future experiments on particle emissions in neutron-rich calcium isotopes. Shell evolutions in the calcium chain around neutron numbers \textit{N} = 32, 34 and 40 are understood by calculating effective single-particle energies, the excitation energies of the first $2^+$ states and two-neutron separation energies. The calculations support shell closures at $^{52}$Ca (\textit{N} = 32) and $^{54}$Ca (\textit{N} = 34) but show a weakening of shell closure at $^{60}$Ca (\textit{N} = 40). The possible shell closure at $^{70}$Ca (\textit{N} = 50) is predicted.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
Continuum and three-nucleon force in Borromean system: The 17Ne case
Authors:
Y. Z. Ma,
F. R. Xu,
N. Michel,
S. Zhang,
J. G. Li,
B. S. Hu,
L. Coraggio,
N. Itaco,
A. Gargano
Abstract:
Starting from chiral two-nucleon (2NF) and chiral three-nucleon (3NF) potentials, we present a detailed study of 17Ne, a Borromean system, with the Gamow shell model which can capture continuum effects. More precisely, we take advantage of the normal-ordering approach to include the 3NF and the Berggren representation to treat bound, resonant and continuum states on equal footing in a complex-mome…
▽ More
Starting from chiral two-nucleon (2NF) and chiral three-nucleon (3NF) potentials, we present a detailed study of 17Ne, a Borromean system, with the Gamow shell model which can capture continuum effects. More precisely, we take advantage of the normal-ordering approach to include the 3NF and the Berggren representation to treat bound, resonant and continuum states on equal footing in a complex-momentum plane. In our framework, 3NF is essential to reproduce the Borromean structure of 17Ne, while the continuum is more crucial for the halo property of the nucleus. The two-proton halo structure is demonstrated by calculating the valence proton density and correlation density. The astrophysically interesting $3/2^-$ excited state has its energy above the threshold of the proton emission, and therefore the two-proton decay should be expected from the state.
△ Less
Submitted 4 August, 2020; v1 submitted 4 August, 2020;
originally announced August 2020.
-
Chiral three-nucleon force and continuum for dripline nuclei and beyond
Authors:
Y. Z. Ma,
F. R. Xu,
L. Coraggio,
B. S. Hu,
J. G. Li,
T. Fukui,
L. De Angelis,
N. Itaco,
A. Gargano
Abstract:
Three-nucleon force and continuum play important roles in reproducing the properties of atomic nuclei around driplines. Therefore it is valuable to build up a theoretical framework where both effects can be taken into account to solve the nuclear Schrödinger equation. To this end, in this letter, we have expressed the chiral three-nucleon force within the continuum Berggren representation, so that…
▽ More
Three-nucleon force and continuum play important roles in reproducing the properties of atomic nuclei around driplines. Therefore it is valuable to build up a theoretical framework where both effects can be taken into account to solve the nuclear Schrödinger equation. To this end, in this letter, we have expressed the chiral three-nucleon force within the continuum Berggren representation, so that bound, resonant and continuum states can be treated on an equal footing in the complex-momentum space. To reduce the model dimension and computational cost, the three-nucleon force is truncated at the normal-ordered two-body level and limited in the $sd$-shell model space, with the residual three-body term being neglected. We choose neutron-rich oxygen isotopes as the test ground because they have been well studied experimentally, with the neutron dripline determined. The calculations have been carried out within the Gamow shell model. The quality of our results in reproducing the properties of oxygen isotopes around the neutron dripline shows the relevance of the interplay between three-nucleon force and the coupling to continuum states. We also analyze the role played by the chiral three-nucleon force, by dissecting the contributions of the $2π$ exchange, $1π$ exchange and contact terms.
△ Less
Submitted 2 August, 2020;
originally announced August 2020.
-
Nuclear multipole responses from chiral effective field theory interaction
Authors:
B. S. Hu,
Q. Wu,
Q. Yuan,
Y. Z. Ma,
X. Q. Yan,
F. R. Xu
Abstract:
We probe nuclear multipole resonances in the framework of the random-phase approximation by using the interaction obtained from the chiral effective field theory. The three-nucleon force is included in a form of the in-medium two-nucleon interaction which was derived from the chiral three-nucleon force. The isoscalar monopole, isoscalar dipole, isovector dipole and isoscalar quadrupole resonances…
▽ More
We probe nuclear multipole resonances in the framework of the random-phase approximation by using the interaction obtained from the chiral effective field theory. The three-nucleon force is included in a form of the in-medium two-nucleon interaction which was derived from the chiral three-nucleon force. The isoscalar monopole, isoscalar dipole, isovector dipole and isoscalar quadrupole resonances of the closed-shell $^{56,68,78}$Ni have been investigated. The calculations reasonably reproduce the experimental multipole resonances of $^{56,68}$Ni, and well describe the pygmy dipole resonance and dipole polarizability measured in $^{68}$Ni. The multipole resonances of $^{78}$Ni, including pygmy dipole resonance and dipole polarizability, are predicted. The detailed effects of the tensor force and three-body force are analyzed by dissecting the chiral interaction. We find that in general the tensor force effect on electric giant resonances is not as significant as the effect from the three-body force, although the tensor force provides more than half of the binding energy. The effect from three-body force is strong in light nuclei. Particularly, three-body force is crucial for the formation of the pygmy resonance in calculations.
△ Less
Submitted 23 February, 2020;
originally announced February 2020.
-
An {\it ab-initio} Gamow shell model approach with a core
Authors:
B. S. Hu,
Q. Wu,
J. G. Li,
Y. Z. Ma,
Z. H. Sun,
N. Michel,
F. R. Xu
Abstract:
Gamow shell model (GSM) is usually performed within the Woods-Saxon (WS) basis in which the WS parameters need to be determined by fitting experimental single-particle energies including their resonance widths. In the multi-shell case, such a fit is difficult due to the lack of experimental data of cross-shell single-particle energies and widths. In this paper, we develop an {\it ab-initio} GSM by…
▽ More
Gamow shell model (GSM) is usually performed within the Woods-Saxon (WS) basis in which the WS parameters need to be determined by fitting experimental single-particle energies including their resonance widths. In the multi-shell case, such a fit is difficult due to the lack of experimental data of cross-shell single-particle energies and widths. In this paper, we develop an {\it ab-initio} GSM by introducing the Gamow Hartree-Fock (GHF) basis that is obtained using the same interaction as the one used in the construction of the shell-model Hamiltonian. GSM makes use of the complex-momentum Berggren representation, then including resonance and continuum components. Hence, GSM gives a good description of weakly bound and unbound nuclei. Starting from chiral effective field theory and employing many-body perturbation theory (MBPT) (called nondegenerate $\hat Q$-box folded-diagram renormalization) in the GHF basis, a multi-shell Hamiltonian ({\it sd-pf} shells in this work) can be constructed. The single-particle energies and their resonance widths can also been obtained using MBPT. We investigated $^{23-28}$O and $^{23-31}$F isotopes, for which multi-shell calculations are necessary. Calculations show that continuum effects and the inclusion of the {\it pf} shell are important elements to understand the structure of nuclei close to and beyond driplines.
△ Less
Submitted 8 January, 2020;
originally announced January 2020.
-
Ab-initio no-core Gamow shell model calculations of multi-neutron systems
Authors:
Jian Guo Li,
Nicolas Michel,
Bai Shan Hu,
Wei Zuo,
Fu Rong Xu
Abstract:
The existence of multi-neutron systems has always been a debatable question. Indeed, both inter-nucleon correlations and a large continuum coupling occur in these states. We then employ the ab-initio no-core Gamow shell model to calculate the resonant energies and widths of the trineutron and tetraneutron systems with realistic interactions. Our results indicate that trineutron and tetraneutron ar…
▽ More
The existence of multi-neutron systems has always been a debatable question. Indeed, both inter-nucleon correlations and a large continuum coupling occur in these states. We then employ the ab-initio no-core Gamow shell model to calculate the resonant energies and widths of the trineutron and tetraneutron systems with realistic interactions. Our results indicate that trineutron and tetraneutron are both unbound and bear broad widths. The calculated energy and width of tetraneutron are also comparable with recent experimental data. Moreover, our calculations suggest that the energy of trineutron is lower than that of tetraneutron, while its resonance width is also narrower. This strongly suggests that trineutron is more likely to be experimentally observed than tetraneutron. We thus suggest experimentalists to search for trineutron at low energy.
△ Less
Submitted 18 November, 2019; v1 submitted 15 November, 2019;
originally announced November 2019.
-
High-Resolution, Respiratory-Resolved Coronary MRA Using a Phyllotaxis-Reordered Variable-Density 3D Cones Trajectory
Authors:
Srivathsan P. Koundinyan,
Corey A. Baron,
Mario O. Malave,
Frank Ong,
Nii Okai Addy,
Joseph Y. Cheng,
Phillip C. Yang,
Bob S. Hu,
Dwight G. Nishimura
Abstract:
Purpose: To develop a respiratory-resolved motion-compensation method for free-breathing, high-resolution coronary magnetic resonance angiography using a 3D cones trajectory.
Methods: To achieve respiratory-resolved 0.98 mm resolution images in a clinically relevant scan time, we undersample the imaging data with a variable-density 3D cones trajectory. For retrospective motion compensation, tran…
▽ More
Purpose: To develop a respiratory-resolved motion-compensation method for free-breathing, high-resolution coronary magnetic resonance angiography using a 3D cones trajectory.
Methods: To achieve respiratory-resolved 0.98 mm resolution images in a clinically relevant scan time, we undersample the imaging data with a variable-density 3D cones trajectory. For retrospective motion compensation, translational estimates from 3D image-based navigators (3D iNAVs) are used to bin the imaging data into four phases from end-expiration to end-inspiration. To ensure pseudo-random undersampling within each respiratory phase, we devise a phyllotaxis readout ordering scheme mindful of eddy current artifacts in steady state free precession imaging. Following binning, residual 3D translational motion within each phase is computed using the 3D iNAVs and corrected for in the imaging data. The noise-like aliasing characteristic of the combined phyllotaxis and cones sampling pattern is leveraged in a compressed sensing reconstruction with spatial and temporal regularization to reduce aliasing in each of the respiratory phases.
Results: In a volunteer and 5 patients, respiratory motion compensation using the proposed method yields improved image quality compared to non-respiratory-resolved approaches with no motion correction and with 3D translational correction. Qualitative assessment by two cardiologists indicates the superior sharpness of coronary segments reconstructed with the proposed method (P < 0.01).
Conclusion: The proposed method better mitigates motion artifacts in free-breathing, high-resolution coronary angiography exams compared to translational correction.
△ Less
Submitted 27 October, 2019;
originally announced October 2019.
-
Unraveling the Effect of Spatial Resolution and Scan Acceleration on 3D Image-Based Navigators for Respiratory Motion Tracking in Coronary MR Angiography
Authors:
Srivathsan P. Koundinyan,
Joseph Y. Cheng,
Mario O. Malave,
Phillip C. Yang,
Bob S. Hu,
Dwight G. Nishimura,
Corey A. Baron
Abstract:
Purpose: To study the accuracy of motion information extracted from beat-to-beat 3D image-based navigators (3D iNAVs) collected using a variable-density cones trajectory with different combinations of spatial resolutions and scan acceleration factors.
Methods: Fully sampled, breath-held 4.4 mm 3D iNAV datasets for six respiratory phases are acquired in a volunteer. Ground truth translational and…
▽ More
Purpose: To study the accuracy of motion information extracted from beat-to-beat 3D image-based navigators (3D iNAVs) collected using a variable-density cones trajectory with different combinations of spatial resolutions and scan acceleration factors.
Methods: Fully sampled, breath-held 4.4 mm 3D iNAV datasets for six respiratory phases are acquired in a volunteer. Ground truth translational and nonrigid motion information is derived from these datasets. Subsequently, the motion estimates from synthesized undersampled 3D iNAVs with isotropic spatial resolutions of 4.4 mm (acceleration factor = 10.9), 5.4 mm (acceleration factor = 7.2), 6.4 mm (acceleration factor = 4.2), and 7.8 mm (acceleration factor = 2.9) are assessed against the ground truth information. The undersampled 3D iNAV configuration with the highest accuracy motion estimates in simulation is then compared with the originally proposed 4.4 mm undersampled 3D iNAV in six volunteer studies.
Results: The simulations indicate that for navigators beyond certain scan acceleration factors, the accuracy of motion estimates is compromised due to errors from residual aliasing and blurring/smoothening effects following compressed sensing reconstruction. The 6.4 mm 3D iNAV achieves an acceptable spatial resolution with a small acceleration factor, resulting in the highest accuracy motion information among all assessed undersampled 3D iNAVs. Reader scores for six volunteer studies demonstrate superior coronary vessel sharpness when applying an autofocusing nonrigid correction technique using the 6.4 mm 3D iNAVs in place of 4.4 mm 3D iNAVs.
Conclusion: Undersampled 6.4 mm 3D iNAVs enable motion tracking with improved accuracy relative to previously proposed undersampled 4.4 mm 3D iNAVs.
△ Less
Submitted 27 October, 2019;
originally announced October 2019.
-
Ab initio Gamow in-medium similarity renormalization group with resonance and continuum
Authors:
B. S. Hu,
Q. Wu,
Z. H. Sun,
F. R. Xu
Abstract:
We have developed a novel ab initio Gamow in-medium similarity renormalization group (Gamow IMSRG) in the complex-energy Berggren framework. The advanced Gamow IMSRG is capable of describing the resonance and nonresonant continuum properties of weakly bound and unbound nuclear many-body systems. As test grounds, carbon and oxygen isotopes have been calculated with chiral two- and three-nucleon for…
▽ More
We have developed a novel ab initio Gamow in-medium similarity renormalization group (Gamow IMSRG) in the complex-energy Berggren framework. The advanced Gamow IMSRG is capable of describing the resonance and nonresonant continuum properties of weakly bound and unbound nuclear many-body systems. As test grounds, carbon and oxygen isotopes have been calculated with chiral two- and three-nucleon forces from the effective field theory. Resonant states observed in the neutron-dripline 24O are well reproduced. The halo structure of the known heaviest Borromean nucleus 22C is clearly seen by calculating the density distribution in which the continuum s channel plays a crucial role. Furthermore, we predict low-lying resonant excited states in 22C. The Gamow IMSRG provides tractable ab initio calculations of weakly bound and unbound open quantum systems.
△ Less
Submitted 23 June, 2019;
originally announced June 2019.
-
$p$-shell hypernuclear structures using the Gogny-interaction shell model
Authors:
X. Y. Chen,
Z. Zhou,
W. G. Jiang,
B. S. Hu,
F. R. Xu
Abstract:
We have systematically investigated the excitation spectra of $p$-shell hypernuclei within the shell model based on the nucleon-nucleon and hyperon-nucleon interactions. For the effective nucleon-nucleon interaction, we adopt the Gogny force instead of the widely-used empirical $p$-shell Cohen-Kurath interaction, while the hyperon-nucleon interaction takes the $ΛN$ interaction including the $ΛN$-…
▽ More
We have systematically investigated the excitation spectra of $p$-shell hypernuclei within the shell model based on the nucleon-nucleon and hyperon-nucleon interactions. For the effective nucleon-nucleon interaction, we adopt the Gogny force instead of the widely-used empirical $p$-shell Cohen-Kurath interaction, while the hyperon-nucleon interaction takes the $ΛN$ interaction including the $ΛN$-$ΣN$ coupling effect. We find that the shell model with the Gogny force can give reasonable descriptions of both spectra and binding energies of the $p$-shell nuclei. With this confidence, combined with the $ΛN$ interaction, we have performed shell-model calculations for the $p$-shell hypernuclei. We compare our results with $γ$-ray data as well as various theoretical calculations, and explain recent experimental hypernuclear excitation spectra observed at JLab.
△ Less
Submitted 20 October, 2018;
originally announced October 2018.
-
Ab initio Rayleigh-Schrödinger perturbation calculation including three-body force
Authors:
B. S. Hu,
T. Li,
F. R. Xu
Abstract:
We first derive the Rayleigh-Schrödinger many-body perturbation theory up to third order (RSPT3) for Hamiltonians with three-body interaction. The structure of closed-shell nuclei in a wide mass range from 4He to 48Ca has been investigated by the RSPT3 with explicit NN+3N Hamiltonian. The RSPT3 calculations are performed within Hartree-Fock bases. The perturbative contribution of antisymmetrized G…
▽ More
We first derive the Rayleigh-Schrödinger many-body perturbation theory up to third order (RSPT3) for Hamiltonians with three-body interaction. The structure of closed-shell nuclei in a wide mass range from 4He to 48Ca has been investigated by the RSPT3 with explicit NN+3N Hamiltonian. The RSPT3 calculations are performed within Hartree-Fock bases. The perturbative contribution of antisymmetrized Goldstone diagrams (diagrammatic expansion for RSPT) with normal-ordered interaction has been analyzed. We demonstrate that the normal-ordered two-body level (NO2B) approximation which neglects the residual three-body term can catch the main effect of full three-body force in RSPT calculation. We also present rigorous benchmarks for RSPT3 with non-perturbative coupled cluster and in-medium similarity renormalization group using the same chiral NNLOsat NO2B interaction. The three methods are in good agreement with each other for binding energies and radii. However, the RSPT3 provides an alternative that is computationally inexpensive but comparable in accuracy to state-of-the-art non-perturbative many-body methods.
△ Less
Submitted 20 October, 2018;
originally announced October 2018.
-
Resonances of unbound quantum many-body systems of nuclei
Authors:
B. S. Hu,
Q. Wu,
Z. H. Sun,
F. R. Xu
Abstract:
Resonance is a general phenomenon which can happen in classic or quantum systems. An unbound many-body quantum system can undergo a self-resonant process. It has long been a challenge how to describe unbound many-body quantum systems in resonances. In this paper, we develop a novel first-principles method that is capable of describing resonant quantum systems. We exploit, for the first time, the a…
▽ More
Resonance is a general phenomenon which can happen in classic or quantum systems. An unbound many-body quantum system can undergo a self-resonant process. It has long been a challenge how to describe unbound many-body quantum systems in resonances. In this paper, we develop a novel first-principles method that is capable of describing resonant quantum systems. We exploit, for the first time, the advanced in-medium similarity renormalization group (IMSRG) in the complex-energy Gamow-Berggren representation with resonance and non-resonant continuum. The ab initio Gamow IMSRG has broad applications, such as, to the electromagnetic-interaction systems of atoms, molecules or quantum dots, and strong-interaction atomic nuclei. In the present work, we apply the method to loosely bound or unbound nuclear systems. Carbon and oxygen isotopes have been investigated with an optimized chiral effective field theory potential. The resonant states observed in neutron-rich 22O and 24O are well reproduced. The loose halo structure of the Borromean nucleus 22C is clearly seen by the density calculation, in which the continuum s-waves play a crucial role. Further, we predict low-lying resonant excited states in 22C. This method provides rigorous and tractable theoretical calculations for weakly-bound or unbound open quantum systems.
△ Less
Submitted 22 September, 2018;
originally announced September 2018.
-
{\it Ab initio} perturbation calculations of realistic effective interactions in the Hartree--Fock basis
Authors:
Q. Wu,
F. R. Xu,
B. S. Hu,
J. G. Li
Abstract:
We perform two types of {\it ab initio} perturbation calculations of effective interactions in the Hartree--Fock (HF) basis instead of the harmonic-oscillator basis: one is called the Brillouin--Wigner (BW) perturbation and another is called the Rayleigh--Schrödinger (RS) perturbation. It is shown that the HF basis provides good convergences. We also benchmark the perturbation calculations with th…
▽ More
We perform two types of {\it ab initio} perturbation calculations of effective interactions in the Hartree--Fock (HF) basis instead of the harmonic-oscillator basis: one is called the Brillouin--Wigner (BW) perturbation and another is called the Rayleigh--Schrödinger (RS) perturbation. It is shown that the HF basis provides good convergences. We also benchmark the perturbation calculations with the in-medium similarity renormalization group (IM-SRG) which is a nonperturbative method. In the HF basis some-type perturbation diagrams can be cancelled out, while the cancellation does not happen in the harmonic-oscillator basis. We have investigated the {\it sd} shell using the chiral N$^3$LO potential softened by $V_{\text{low-}k}$. With the low-momentum N$^3$LO potential, we first perform the spherical HF calculations for the $^{16}$O core of the {\it sd} shell, and the realistic effective two-body interactions for the HF {\it sd}-shell space can be derived by the perturbation calculations. The calculations give simultaneously effective single-particle energies and excitation spectra of two valence particle systems (i.e., $^{18}$O, $^{18}$F, and $^{18}$Ne in the {\it sd} shell). Convergences have been analyzed order by order. The perturbation calculations are in fairly good agreement with nonperturbative method IM-SRG. We find that the HF RS perturbation gives even better results compared with the BW perturbation. The HF realistic effective interactions derived by the perturbations can be used for further shell-model calculations.
△ Less
Submitted 24 June, 2018;
originally announced June 2018.
-
Chiral NNLO$_{\text{sat}}$ descriptions of nuclear multipole resonances within the random phase approximation
Authors:
Q. Wu,
B. S. Hu,
F. R. Xu,
Y. Z. Ma,
S. J. Dai,
Z. H. Sun,
G. R. Jansen
Abstract:
We study nuclear multipole resonances in the framework of the random phase approximation using the chiral potential NNLO$_{\text{sat}}$. This potential includes two- and three-body terms that has been simultaneously optimized to low-energy nucleon-nucleon scattering data and selected nuclear structure data. Our main foci have been the isoscalar monopole, isovector dipole, and isoscalar quadrupole…
▽ More
We study nuclear multipole resonances in the framework of the random phase approximation using the chiral potential NNLO$_{\text{sat}}$. This potential includes two- and three-body terms that has been simultaneously optimized to low-energy nucleon-nucleon scattering data and selected nuclear structure data. Our main foci have been the isoscalar monopole, isovector dipole, and isoscalar quadrupole resonances of the closed-shell nuclei, $^4$He, $^{16,22,24}$O, and $^{40,48}$Ca. These resonance modes have been widely observed in experiment. In addition, we use a renormalized chiral potential $V_{\text{low-}k}$, based on the N$^3$LO two-body potential by Entem and Machleidt [Phys. Rev. \textbf{C68}, 041001 (2011)]. This introduces a dependency on the cutoff parameter used in the normalization procedure as reported in previous works by other groups. While NNLO$_{\text{sat}}$ can reasonably reproduce observed multipole resonances, it is not possible to find a single cutoff parameter for the $V_{\text{low-}k}$ potential that simultaneously describe the different types of resonance modes. The sensitivity to the cutoff parameter can be explained by missing induced three-body forces in the calculations. Our results for neutron-rich $^{22,24}$O show a mixing nature of isoscalar and isovector resonances in the dipole channel at low energies. We predict that $^{22}$O and $^{24}$O have low-energy isoscalar quadrupole resonances at energies lower than 5 MeV.
△ Less
Submitted 13 January, 2018;
originally announced January 2018.
-
Brueckner-Hartree-Fock and its renormalized calculations for finite nuclei
Authors:
B. S. Hu,
F. R. Xu,
Y. Z. Ma,
Q. Wu,
Z. H. Sun
Abstract:
We have performed self-consistent Brueckner-Hartree-Fock (BHF) and its renormalized theory to the structure calculations of finite nuclei. The $G$-matrix is calculated within the BHF basis, and the exact Pauli exclusion operator is determined by the BHF spectrum. Self-consistent occupation probabilities are included in the renormalized Brueckner-Hartree-Fock (RBHF). Various systematics and converg…
▽ More
We have performed self-consistent Brueckner-Hartree-Fock (BHF) and its renormalized theory to the structure calculations of finite nuclei. The $G$-matrix is calculated within the BHF basis, and the exact Pauli exclusion operator is determined by the BHF spectrum. Self-consistent occupation probabilities are included in the renormalized Brueckner-Hartree-Fock (RBHF). Various systematics and convergences are studies. Good results are obtained for the ground-state energy and radius. RBHF can give a more reasonable single-particle spectrum and radius. We present a first benchmark calculation with other {\it ab initio} methods using the same effective Hamiltonian. We find that the BHF and RBHF results are in good agreement with other $\it{ab}$ $\it{initio}$ methods.
△ Less
Submitted 9 September, 2016;
originally announced September 2016.