Researcher profile

Vittorio Somà

Vittorio Somà contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2024arXiv

Deformed natural orbitals for ab initio calculations

The rapid development of ab initio nuclear structure methods towards doubly open-shell nuclei, heavy nuclei and greater accuracy occurs at the price of evermore increased computational costs, especially RAM and CPU time. While most of the numerical simulations are carried out by expanding relevant operators and wave functions on the spherical harmonic oscillator basis, alternative one-body bases offering advantages in terms of computational efficiency have recently been investigated. In particular, the so-called natural basis used in combination with symmetry-conserving methods applicable to doubly closed-shell nuclei has proven beneficial in this respect. The present work examines the performance of the natural basis in the context of symmetry-breaking many-body calculations enabling the description of superfluid and deformed open-shell nuclei at polynomial cost with system's size. First, it is demonstrated that the advantage observed for closed-shell nuclei carries over to open-shell ones. A detailed investigation of natural-orbital wave functions provides useful insight to support this finding and to explain the superiority of the natural basis over alternative ones. Second, it is shown that the use of natural orbitals combined with importance-truncation techniques leads to an even greater gain in terms of computational costs. The presents results pave the way for the systematic use of natural-orbital bases in future implementations of non-perturbative many-body methods.

preprint2024arXiv

Impact of correlations on nuclear binding energies

A strong effort will be dedicated in the coming years to extend the reach of ab initio nuclear-structure calculations to heavy doubly open-shell nuclei. In order to do so, the most efficient strategies to incorporate dominant many-body correlations at play in such nuclei must be identified. With this motivation in mind, the present work pedagogically analyses the inclusion of many-body correlations and their impact on binding energies of Calcium and Chromium isotopes. Employing an empirically-optimal Hamiltonian built from chiral effective field theory, binding energies along both isotopic chains are studied via a hierarchy of approximations based on polynomially-scaling expansion many-body methods. The corresponding results are compared to experimental data and to those obtained via valence-space in-medium similarity renormalization group calculations at the normal-ordered two-body level that act as a reference in the present study. The spherical mean-field approximation is shown to display specific shortcomings in Ca isotopes that can be understood analytically and that are efficiently corrected via the consistent addition of low-order dynamical correlations on top of it. While the same setting cannot appropriately reproduce binding energies in doubly open-shell Cr isotopes, allowing the unperturbed mean-field state to break rotational symmetry permits to efficiently capture the static correlations responsible for the phenomenological differences observed between the two isotopic chains. Eventually, the present work demonstrates in a pedagogical way that polynomially-scaling expansion methods based on unperturbed states that possibly break (and restore) symmetries constitute an optimal route to extend ab initio calculations to heavy closed- and open-shell nuclei.

preprint2022arXiv

Gorkov algebraic diagrammatic construction formalism at third order

Background. The Gorkov approach to self-consistent Green's function theory has been formulated in [V. Somà, T. Duguet, C. Barbieri, Phys. Rev. C 84, 064317 (2011)]. Over the past decade, it has become a method of reference for first-principle computations of semi-magic nuclear isotopes. The currently available implementation is limited to a second-order self-energy and neglects particle-number non-conserving terms arising from contracting three-particle forces with anomalous propagators. For nuclear physics applications, this is sufficient to address first-order energy differences, ground-state radii and moments on an accurate enough basis. However, addressing absolute binding energies, fine spectroscopic details of $N\pm1$ particle systems or delicate quantities such as second-order energy differences associated to pairing gaps, requires to go to higher truncation orders. Purpose. The formalism is extended to third order in the algebraic diagrammatic construction (ADC) expansion with two-body Hamiltonians. Methods. The expansion of Gorkov propagators in Feynman diagrams is combined with the algebraic diagrammatic construction up to the third order as an organization scheme to generate the Gorkov self-energy. Results. Algebraic expressions for the static and dynamic contributions to the self-energy, along with equations for the matrix elements of the Gorkov eigenvalue problem, are derived. It is first done for a general basis before specifying the set of equations to the case of spherical systems displaying rotational symmetry. Workable approximations to the full self-consistency problem are also elaborated on. The formalism at third order it thus complete for a general two-body Hamiltonian. Conclusion. Working equations for the full Gorkov-ADC(3) are now available for numerical implementation.

preprint2022arXiv

Multi-reference many-body perturbation theory for nuclei I -- Novel PGCM-PT formalism

Perturbative and non-perturbative expansion methods already constitute a tool of choice to perform ab initio calculations over a significant part of the nuclear chart. In this context, the categories of accessible nuclei directly reflect the class of unperturbed state employed in the formulation of the expansion. The present work generalizes to the nuclear many-body context the versatile method of Ref. \cite{burton20a} by formulating a perturbative expansion on top of a multi-reference unperturbed state mixing deformed non-orthogonal Bogoliubov vacua, i.e. a state obtained from the projected generator coordinate method (PGCM). Particular attention is paid to the part of the mixing taking care of the symmetry restoration, showing that it can be exactly contracted throughout the expansion, thus reducing significantly the dimensionality of the linear problem to be solved to extract perturbative corrections. While the novel expansion method, coined as PGCM-PT, reduces to the PGCM at lowest order, it reduces to single-reference perturbation theories in appropriate limits. Based on a PGCM unperturbed state capturing (strong) static correlations in a versatile and efficient fashion, PGCM-PT is indistinctly applicable to doubly closed-shell, singly open-shell and doubly open-shell nuclei. The remaining (weak) dynamical correlations are brought consistently through perturbative corrections. This symmetry-conserving multi-reference perturbation theory is state-specific and applies to both ground and excited PGCM unperturbed states, thus correcting each state belonging to the low-lying spectrum of the system under study.

preprint2022arXiv

Multi-reference many-body perturbation theory for nuclei II -- Ab initio study of neon isotopes via PGCM and IM-NCSM calculations

The neon isotopic chain displays a rich phenomenology, ranging from clustering in the ground-state of the self-conjugate doubly open-shell stable $^{20}$Ne isotope to the physics of the island of inversion around the neutron-rich $^{30}$Ne isotope. This second (i.e. Paper II) of the present series proposes an extensive ab initio study of neon isotopes based on two complementary many-body methods, i.e. the quasi-exact in-medium no-core shell model (IM-NCSM) and the projected generator coordinate method (PGCM) that is ideally suited to capturing strong static correlations associated with shape deformation and fluctuations. Calculations employ a state-of-the-art generation of chiral effective field theory Hamiltonians and evaluate the associated systematic uncertainties. In spite of missing so-called dynamical correlations, which can be added via the multi-reference perturbation theory proposed in the first paper (i.e. Paper I) of the present series, the PGCM is shown to be a suitable method to tackle the low-lying spectroscopy of complex nuclei. Still, describing the physics of the island of inversion constitutes a challenge that seems to require the inclusion of dynamical correlations. This is addressed in the third paper (i.e. Paper III) of the present series.

preprint2022arXiv

Multi-reference many-body perturbation theory for nuclei III -- Ab initio calculations at second order in PGCM-PT

In spite of missing dynamical correlations, the projected generator coordinate method (PGCM) was recently shown to be a suitable method to tackle the low-lying spectroscopy of complex nuclei. Still, describing absolute binding energies and reaching high accuracy eventually requires the inclusion of dynamical correlations on top of the PGCM. In this context, the present work discusses the first realistic results of a novel multi-reference perturbation theory (PGCM-PT) that can do so within a symmetry-conserving scheme for both ground and low-lying excited states. First, proof-of-principle calculations in a small ($e_{\mathrm{max}}=4$) model space demonstrate that exact binding energies of closed- (\nucl{O}{16}) and open-shell (\nucl{O}{18}, \nucl{Ne}{20}) nuclei are reproduced within $0.5-1.5\%$ at second order, i.e. through PGCM-PT(2). Moreover, profiting from the pre-processing of the Hamiltonian via multi-reference in-medium similarity renormalization group transformations, PGCM-PT(2) can reach converged values within smaller model spaces than with an unevolved Hamiltonian. Doing so, dynamical correlations captured by PGCM-PT(2) are shown to bring essential corrections to low-lying excitation energies that become too dilated at leading order, i.e., at the strict PGCM level. The present work is laying the foundations for a better understanding of the optimal way to grasp static and dynamical correlations in a consistent fashion, with the aim of accurately describing ground and excited states of complex nuclei via ab initio many-body methods.

preprint2021arXiv

Evidence of the triaxial structure of $\boldsymbol{^{129}}$Xe at the Large Hadron Collider

The interpretation of the emergent collective behaviour of atomic nuclei in terms of deformed intrinsic shapes [1] is at the heart of our understanding of the rich phenomenology of their structure, ranging from nuclear energy to astrophysical applications across a vast spectrum of energy scales. A new window onto the deformation of nuclei has been recently opened with the realization that nuclear collision experiments performed at high-energy colliders, such as the CERN Large Hadron Collider (LHC), enable experimenters to identify the relative orientation of the colliding ions in a way that magnifies the manifestations of their intrinsic deformation [2]. Here we apply this technique to LHC data on collisions of $^{129}$Xe nuclei [3-5] to exhibit the first evidence of non-axiality in the ground state of ions collided at high energy. We predict that the low-energy structure of $^{129}$Xe is triaxial (a spheroid with three unequal axes), and show that such deformation can be determined from high-energy data. This result demonstrates the unique capabilities of precision collider machines such as the LHC as new means to perform imaging of the collective structure of atomic nuclei.

preprint2020arXiv

Examining the $N$ = 28 shell closure through high-precision mass measurements of $^{46-48}$Ar

The strength of the $N$ = 28 magic number in neutron-rich argon isotopes is examined through high-precision mass measurements of $^{46-48}$Ar, performed with the ISOLTRAP mass spectrometer at ISOLDE/CERN. The new mass values are up to 90 times more precise than previous measurements. While they suggest the persistence of the $N$ = 28 shell closure for argon, we show that this conclusion has to be nuanced in light of the wealth of spectroscopic data and theoretical investigations performed with the \emph{SDPF-U} phenomenological shell model interaction. Our results are also compared with \emph{ab initio} calculations using the Valence Space In-Medium Similarity Renormalization Group and the Self-Consistent Green's Function approaches. Both calculations provide a very good account of mass systematics at and around $Z$ = 18 and, generally, a consistent description of the physics in this region. This combined analysis indicates that $^{46}$Ar is the transition between the closed-shell $^{48}$Ca and collective $^{44}$S.

preprint2019arXiv

Improved many-body expansions from eigenvector continuation

Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei however, effects such as strong short-range repulsion between nucleons can spoil the convergence of the expansion and make the reliability of perturbation theory unclear. Mathematicians have devised an extensive machinery to overcome the problem of divergent expansions by making use of so-called resummation methods. In large-scale many-body applications such schemes are often of limited use since no a priori analytical knowledge of the expansion is available. We present here eigenvector continuation as an alternative resummation tool that is both efficient and reliable because it is based on robust and simple mathematical principles.