Search

Giant Resonances are, with nuclear rotations, the most evident expression of collectivity in finite nuclei. These two categories of excitations, however, are traditionally described within different formal schemes, such that vibrational and rotational degrees of freedom are separately treated and coupling effects between those are often neglected. The present work puts forward an approach aiming at a consitent treatment of vibrations and rotations. Specifically, this paper is the last in a series of four dedicated to the investigation of the giant monopole resonance in doubly open-shell nuclei via the ab initio Projected Generator Coordinate Method (PGCM). The present focus is on the treatment and impact of angular momentum restoration within such calculations. The PGCM being based on the use of deformed mean-field states, the angular-momentum restoration is performed when solving the secular equation to extract vibrational excitations. In this context, it is shown that performing the angular momentum restoration only after solving the secular equation contaminates the monopole response with an unphysical coupling to the rotational motion, as was also shown recently for (quasi-particle) random phase approximation calculations based on a deformed reference state. Eventually, the present work based on the PGCM confirms that an a priori angular momentum restoration is necessary to handle consistently both collective motions at the same time. This further pleads in favor of implementing the full-fledged projected (quasi-particle) random phase approximation in the future., Comment: 9 pages, 4 figures

The use of quantality is discussed in the case of nuclei and other many-body systems such as atomic electrons. This dimensionless quantity is known to indicate when a many-body system behaves like a crystal or a quantum liquid. Its role is further analyzed by showing its relation to the scattering length. The emergence of a fundamental lengthscale, the limit radius, is also shown. It corresponds to the hard-core of the nucleon-nucleon interaction in the case of nucleons, and to a value close to the Bohr radius in the case of atomic electrons. The occurrence of a cluster phase in nuclei is analyzed using the quantality through its relation to the localization parameter, allowing for the identification of both the number of nucleons and the density as control parameters for the occurrence of this phase. The relation of the quantality to the magnitude of the interaction also exhibits a third dimensionless parameter, monitoring the magnitude of the spin-orbit effect in finite systems, through the realization of the pseudo-spin symmetry. The impact of quantality on the spin-orbit effect is compared in various many-body systems. The role of quantality in the relative effect of the binding energy and the shell one is also analyzed in nuclei. Finally, additional dimensionless quantities are proposed from the generalization of the quantality. Nuclei are found to be exceptional systems because all their dimensionless quantities are close to the order of unity, at variance with other many-body systems., Comment: 13 pages, 6 figures

The System for Measuring Overlap with Gas (SMOG2) at the LHCb detector enables the study of fixed-target ion-ion collisions at relativistic energies ($\sqrt{s_{\rm NN}}\sim100$ GeV in the centre-of-mass). With input from \textit{ab initio} calculations of the structure of $^{16}$O and $^{20}$Ne, we compute 3+1D hydrodynamic predictions for the anisotropic flow of Pb+Ne and Pb+O collisions, to be tested with upcoming LHCb data. This will allow the detailed study of quark-gluon plasma (QGP) formation as well as experimental tests of the predicted nuclear shapes. Elliptic flow ($v_2$) in Pb+Ne collisions is greatly enhanced compared to the Pb+O baseline due to the shape of $^{20}$Ne, which is deformed in a bowling-pin geometry. Owing to the large $^{208}$Pb radius, this effect is seen in a broad centrality range, a unique feature of this collision configuration. Larger elliptic flow further enhances the quadrangular flow ($v_4$) of Pb+Ne collisions via non-linear coupling, and impacts the sign of the kurtosis of the elliptic flow vector distribution ($c_2\{4\}$). Exploiting the shape of $^{20}$Ne proves thus an ideal method to investigate the formation of QGP in fixed-target experiments at LHCb, and demonstrates the power of SMOG2 as a tool to image nuclear ground states.

The paper is the third of a series dedicated to the ab initio description of monopole giant resonances in mid-mass closed- and open-shell nuclei via the so-called projected generator coordinate method. The present focus is on the computation of the moments $m_k$ of the monopole strength distribution, which are used to quantify its centroid energy and dispersion. First, the capacity to compute low-order moments via two different methods is developed and benchmarked for the $m_1$ moment. Second, the impact of the angular momentum projection on the centroid energy and dispersion of the monopole strength is analysed before comparing the results to those obtained from consistent quasi-particle random phase approximation calculations. Next, the so-called energy weighted sum rule (EWSR) is investigated. First, the appropriate ESWR in the center-of-mass frame is derived analytically. Second, the exhaustion of the intrinsic EWSR is tested in order to quantify the (unwanted) local-gauge symmetry breaking of the presently employed chiral effective field theory ($\chi$EFT) interactions. Finally, the infinite nuclear matter incompressibility associated with the employed $\chi$EFT interactions is extracted by extrapolating the finite-nucleus incompressibility computed from the monopole centroid energy., Comment: 19 pages, 8 figures

Giant resonances (GRs) are a striking manifestation of collective motions in atomic nuclei. The present paper is the second in a series of four dedicated to the use of the projected generator coordinate method (PGCM) for the ab initio determination of the isoscalar giant monopole resonance (GMR) in closed- and open-shell mid-mass nuclei. While the first paper was dedicated to quantifying various uncertainty sources, the present paper focuses on the first applications to three doubly-open shell nuclei, namely $^{46}$Ti, $^{28}$Si and $^{24}$Mg. In particular, the goal is to investigate from an ab initio standpoint (i) the coupling of the GMR with the giant quadrupole resonance (GQR) in intrinsically-deformed nuclei, (ii) the possible impact of shape coexistence and shape mixing on the GMR, (iii) the GMR based on shape isomers and (iv) the impact of anharmonic effects on the monopole response. The latter is studied by comparing PGCM results to those obtained via the quasi-particle random phase approximation (QRPA), the traditional many-body approach to giant resonances, performed in a consistent setting. Eventually, PGCM results for sd-shell nuclei are in excellent agreement with experimental data, which is attributed to the capacity of the PGCM to capture the important fragmentation of the monopole response in light, intrinsically-deformed systems. Still, the comparison to data in $^{28}$Si and $^{24}$Mg illustrates the challenge (and the potential benefit) of extracting unambiguous experimental information., Comment: 19 pages, 27 figures

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of $^{16}$O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing $^{16}$O$^{16}$O data with collisions of $^{20}$Ne ions. We couple both NLEFT and PGCM ab initio descriptions of the structure of $^{20}$Ne and $^{16}$O to hydrodynamic simulations of $^{16}$O$^{16}$O and $^{20}$Ne$^{20}$Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of $^{20}$Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of $^{20}$Ne$^{20}$Ne collisions is enhanced by as much as 1.170(8)$_{\rm stat.}$(30)$_{\rm syst.}$ for NLEFT and 1.139(6)$_{\rm stat.}$(39)$_{\rm syst.}$ for PGCM relative to $^{16}$O$^{16}$O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system., Comment: 17 pages, 14 figures. The Trajectum code can be found at https://sites.google.com/view/govertnijs/trajectum and plotting routines can be found at http://wilkevanderschee.nl/trajectum

Giant resonances (GRs) are a striking manifestation of collective motions in mesoscopic systems such as atomic nuclei. Until recently, theoretical investigations have essentially relied on the (quasiparticle) random phase approximation ((Q)RPA), and extensions of it, based on phenomenological energy density functionals (EDFs). As part of a current effort to describe GRs within an ab initio theoretical scheme, the present work promotes the use of the projected generator coordinate method (PGCM). This method, which can handle anharmonic effects while satisfying symmetries of the nuclear Hamiltonian, displays a favorable (i.e. mean-field-like) scaling with system's size. Presently focusing on the isoscalar giant monopole resonance (GMR) of light- and medium-mass nuclei, PGCM's potential to deliver wide-range ab initio studies of GRs in closed- and open-shell nuclei encompassing pairing, deformation, and shape coexistence effects is demonstrated. The comparison with consistent QRPA calculations highlights PGCM's unique attributes and sheds light on the intricate interplay of nuclear collective excitations. The present paper is the first in a series of four and focuses on technical aspects and uncertainty quantification of ab initio PGCM calculations of GMR using the doubly open-shell $^{46}$Ti as an illustrative example. The second paper displays results for a set of nuclei of physical interest and proceeds to the comparison with consistent (deformed) ab initio QRPA calculations. While the third paper analyzes useful moments of the monopolar strength function and different ways to access them within PGCM calculations, the fourth paper focuses on the effect of the symmetry restoration on the monopole strength function., Comment: 23 pages, 26 figures

Alpha and cluster decays are analyzed for heavy nuclei located above $^{208}$Pb on the chart of nuclides: $^{216-220}$Rn and $^{220-224}$Ra, that are also candidates for observing the $2 \alpha$ decay mode. A microscopic theoretical approach based on relativistic Energy Density Functionals (EDF), is used to compute axially-symmetric deformation energy surfaces as functions of quadrupole, octupole and hexadecupole collective coordinates. Dynamical least-action paths for specific decay modes are calculated on the corresponding potential energy surfaces. The effective collective inertia is determined using the perturbative cranking approximation, and zero-point and rotational energy corrections are included in the model. The predicted half-lives for $\alpha$-decay are within one order of magnitude of the experimental values. In the case of single $\alpha$ emission, the nuclei considered in the present study exhibit least-action paths that differ significantly up to the scission point. The differences in alpha-decay lifetimes are not only driven by Q values, but also by variances of the least-action paths prior to scission. In contrast, the $2 \alpha$ decay mode presents very similar paths from equilibrium to scission, and the differences in lifetimes are mainly driven by the corresponding Q values. The predicted $^{14}$C cluster decay half-lives are within three orders of magnitudes of the empirical values, and point to a much more complex pattern compared to the alpha-decay mode., Comment: 9 pages, 13 figures

Background: The study of additional terms in functionals is relevant to better describe nuclear structure phenomenology. Among these terms, the tensor one is known to impact nuclear structure properties, especially in neutron-rich nuclei. However, its effect has not been studied on the whole nuclear chart yet. Purpose: The impact of terms corresponding to the tensor at the Hartree level, is studied for infinite nuclear matter as well as deformed nuclei, by developing new density-dependent functionals including these terms. In particular, we study in details the improvement such a term can bring to the description of specific nuclear observables. Methods: The framework of covariant energy density functional is used at the Hartree-Bogoliubov level. The free parameters of covariant functionals are optimized by combining Markov-Chain-Monte-Carlo and simplex algorithms. Results: An improvement of the RMS binding energies, spin-orbit splittings and gaps is obtained over the nuclear chart, including axially deformed ones, when including tensors terms. Small modifications of the potential energy surface and densities are also found. In infinite matter, the Dirac mass is shifted to a larger value, in better agreement with experiments. Conclusions: Taking into account additional terms corresponding to the tensor terms in the vector-isoscalar channel at the Hartree level, improves the description of nuclear properties, both in nuclei and in nuclear matter.

The present paper is the second of a series of publications that aim at investigating relevant directions to turn the nuclear energy density functional (EDF) method as an effective field theory (EFT). The EDF approach has known numerous successes in nuclear theory over the past decades and is currently the only microscopic technique that can be applied to all atomic nuclei. However, the phenomenological character of the EDF method also comes with important limitations, such as the lack of an explicit connection with quantum chromodynamics (QCD). As was argued in the first paper of this series, reformulating the EDF framework as an EFT would enable us to overcome these limitations. In particular, path-integral (PI) techniques are suited to achieve such a purpose as they allow to design numerous non-perturbative approximations and can take Lagrangians possibly derived from EFTs of QCD as inputs. In our previous paper, we have illustrated such technical features for diagrammatic PI techniques in a study of the (0+0)-D $O(N)$-symmetric $\varphi^{4}$-theory. In the present work, we consider another class of PI techniques, i.e. functional renormalization group (FRG) approaches, that we apply to the same toy model. Despite our explicit interest for the nuclear many-body problem, the presented study is also directed towards FRG practitioners from various fields: technical details are provided for FRG techniques based on 1-particle-irreducible (1PI), 2-particle-irreducible (2PI) and 2-particle-point-irreducible (2PPI) effective actions (EAs), coined respectively as 1PI-, 2PI- and 2PPI-FRGs, and the treatment of the $O(N)$ symmetry is also addressed thoroughly. Connections between these various FRG methods are identified as well., Comment: This article is part of a PhD project. The corresponding manuscript can be found at: arXiv:2210.16676

The nucleonic localization function has been used for a decade to study the formation of alpha-particles in nuclei, by providing a measure of having nucleons of a given spin in a single place. However, differences in interpretation remain, compared to the nucleonic density of the nucleus. In order to better understand the respective role of the nucleonic localization function and the densities in the alpha-particle formation in cluster states or in alpha-decay mechanism, both an analytic approximation and microscopic calculations, using energy density functionals, are undertaken. The nucleonic localization function is shown to measure the anti-centrifugal effect, and is not sensitive to the level of compactness of the alpha-particle itself. It probes the purity of the spatial overlap of four nucleons in the four possible (spin, isospin) states. The density provides, in addition, information on the compactness of an alpha-particle cluster. The respective roles of the nucleonic localization function and the density are also analyzed in the case of alpha-particle emission. More generally, criteria to assess the prediction of alpha-cluster in nuclear states are provided., Comment: 9 pages, 6 figures

Recently, ab initio techniques have been successfully connected to the traditional valence-space shell model. In doing so, they can either explicitly provide ab initio shell-model effective Hamiltonians or constrain the construction of empirical ones. In the present work, the possibility to follow a similar path for the nuclear energy density functional (EDF) method is analyzed. For this connection to be actualized, two theoretical techniques are instrumental: the recently proposed ab initio PGCM-PT many-body formalism and the MR-IMSRG pre-processing of the nuclear Hamiltonian. Based on both formal arguments and numerical results, possible new lines of research are briefly discussed, namely to compute ab initio EDF effective Hamiltonians at low computational cost, to constrain empirical ones or to produce them directly via an effective field theory that remains to be invented., Comment: 20 pages, 7 figures

The energy density functional (EDF) method is currently the only microscopic theoretical approach able to tackle the entire nuclear chart. Nevertheless, it suffers from limitations resulting from its empirical character and deteriorating its reliability. This paper is part of a larger program that aims at formulating the EDF approach as an effective field theory (EFT) in order to overcome these limitations. A relevant framework to achieve this is the path-integral (PI) formulation of quantum field theory (QFT). The latter indeed provides a wide variety of treatments of the many-body problem well suited to deal with non-perturbative interactions and to exploit a Lagrangian resulting from an EFT as a starting point. While developing the formalism in a general setting, we present a comparative study of such techniques applied to a toy model, i.e. the (0+0)-D $O(N)$-symmetric $\varphi^{4}$-theory. More specifically, our focus will be on the following diagrammatic techniques: loop expansion (LE), optimized perturbation theory (OPT) and self-consistent perturbation theory (SCPT). With these methods, we notably address the spontaneous breakdown of the $O(N)$ symmetry with care especially since spontaneous symmetry breakings (SSBs) play a paramount role in current implementations of the EDF approach., Comment: This article is part of a PhD project. The corresponding manuscript can be found at: arXiv:2210.16676

We investigate the localization and clustering features in $^{20}$Ne ($N=Z$) and neutron-rich $^{32}$Ne nuclei at zero and finite temperatures. The finite temperature Hartree-Bogoliubov theory is used with the relativistic density-dependent meson-nucleon coupling functional DD-ME2. It is shown that clustering features gradually weaken with increasing temperature and disappear when the shape phase transition occurs. Considering thermal fluctuations in the density profiles, the clustering features vanish at lower temperatures, compared to the case without thermal fluctuations. The effect of the pairing correlations on the nucleon localization and the formation of cluster structures are also studied at finite temperatures. Due to the inclusion of pairing in the calculations, cluster structures are preserved until the critical temperatures for the shape phase transition are reached. Above the critical temperature of the shape phase transition, the clustering features suddenly disappear, which differs from the results without pairing., Comment: 11 pages, 10 figures

Ab initio approaches to the nuclear many-body problem have seen their reach considerably extended over the past decade. However, collective excitations have been scarcely addressed so far due to the prohibitive cost of solving the corresponding equations of motion. Here, a numerically efficient method to compute electromagnetic response functions at zero- and finite-temperature in superfluid and deformed nuclei from an ab initio standpoint is presented and applied to $^{16}$O, $^{28}$Si, $^{46}$Ti and $^{56}$Fe. This work opens the path to systematic ab initio calculations of nuclear responses to electroweak probes across a significant portion of the nuclear chart., Comment: 7 pages, 6 figures

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.

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.

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., Comment: The present paper is the first in a series of three and discusses the PGCM-PT formalism in detail. The second paper displays numerical zeroth-order results, i.e. the outcome of PGCM calculations. Second-order, i.e. PGCM-PT(2), calculations performed in both closed- and open-shell nuclei are the object of the third paper

Background : Several recent experiments report significant low-energy isoscalar monopole strength, below the giant resonance, in various nuclei. In light $\alpha$-conjugate nuclei, these low-energy resonances were recently interpreted as cluster vibration modes. However, the nature of these excitations in neutron-rich nuclei remain elusive. Purpose : The present work provides a systematic analysis of the low-energy monopole strength in isotopic chains, from Neon to Germanium, in order to monitor and understand its nature and conditions of emergence. Methods : We perform covariant quasiparticle random phase approximation (QRPA) calculations, formulated within the finite amplitude method (FAM), on top of constrained relativistic Hartree-Bogoliubov (RHB) reference states. Results : Neutron excess leads to the appearance of low-energy excitations according to a systematic pattern reflecting the single-particle features of the underlying RHB reference state. With the onset of deformation, these low-energy resonances get split and give rise to more complex patterns, with possible mixing with the giant resonance. At lower energy, cluster-like excitations found in $N=Z$ systems survive in neutron-rich nuclei, with valence neutrons arranging in molecular-like orbitals. Finally, at very low energy, pair excitations are also found in superfluid nuclei, but remain negligible in most of the cases. Conclusions : The low-energy part of the monopole strength exhibits various modes, from cluster vibrations ($\sim$ 5-10 MeV) to components of the giant resonance downshifted by the onset of deformation, including soft modes ($\sim$ 10-15 MeV) as well as pair excitation ($<$ 5 MeV), with possible mixing, depending on neutron-excess, deformation, and pairing energy., Comment: 15 pages, 17 figures

The computational cost of ab initio nuclear structure calculations is rendered particularly acute by the presence of (at least) three-nucleon interactions. This feature becomes especially critical now that many-body methods aim at extending their reach beyond mid-mass nuclei. Consequently, state-of-the-art ab initio calculations are typically performed while approximating three-nucleon interactions in terms of effective, i.e. system-dependent, zero-, one- and two-nucleon operators. While straightforward in doubly closed-shell nuclei, existing approximation methods based on normal-ordering techniques involve either two- and three-body density matrices or a symmetry-breaking one-body density matrix in open-shell systems. In order to avoid such complications, a simple, flexible, universal and accurate approximation technique involving the convolution of the initial operator with a sole symmetry-invariant one-body matrix is presently formulated and tested numerically. Employed with a low-resolution Hamiltonian, the novel approximation method is shown to induce errors below $2-3\%$ across a large range of nuclei, observables and many-body methods., Comment: 31 pages, 14 figures

Significant transition strength in light $\alpha$-conjugate nuclei at low energy, typically below 10 MeV, has been observed in many experiments. In this work the isoscalar low-energy response of N=Z nuclei is explored using the Finite Amplitude Method (FAM) based on the microscopic framework of nuclear energy density functionals. Depending on the multipolarity of the excitation and the equilibrium deformation of a particular isotope, the low-energy strength functions display prominent peaks that can be attributed to vibration of cluster structures: $\alpha$+$^{12}$C+$\alpha$ and $\alpha$+$^{16}$O in $^{20}$Ne, $^{12}$C+$^{12}$C in $^{24}$Mg, 4$\alpha$+$^{12}$C in $^{28}$Si, etc. Such cluster excitations are favored in light nuclei with large deformation.

The reduction of the neutron spin-orbit splitting $2p_{3/2} - 2p_{1/2}$ between the $^{41}$Ca and $^{35}$Si isotones is a unique feature throughout the chart of nuclides, as the spin-orbit splitting usually increases with $A$. Moreover, its way of decrease, gradual between $^{41}$Ca and $^{35}$Si or abrupt between $^{37}$S and $^{35}$Si, as well as its origin, caused by the weak binding energy of the $p$ states or by the sudden central proton density depletion in $^{35}$Si, are subject of debate. The results reported here using the self-consistent Covariant Energy Density Functional calculations with the DD-ME2 parametrization rather point to an abrupt, local decrease in $^{35}$Si, and to the large dominance of the central density depletion effect. It is concluded that weak binding, central density depletion as well as correlations must be taken into account to fully evaluate the amplitude and causes of this spin-orbit reduction., Comment: 5 pages, 4 figures, 1 table, to be submitted to Phys.Lett.B

A microscopic calculation of half-lives for the recently observed $^{108}$Xe $\to$ $^{104}$Te $\to$ $^{100}$Sn $\alpha$-decay chain is performed using a self-consistent framework based on energy density functionals. The relativistic density functional DD-PC1 and a separable pairing interaction of finite range are used to compute axially-symmetric deformation energy surfaces of $^{104}$Te and $^{108}$Xe as functions of quadrupole, octupole and hexadecupole collective coordinates. Dynamic least-action paths are determined that trace the $\alpha$-particle emission from the equilibrium deformation to the point of scission. The calculated half-lives: 197 ns for $^{104}$Te and 50 $\mu$s for $^{108}$Xe, are compared to recent experimental values of the half-lives of superallowed $\alpha$-decay of $^{104}$Te: $< 18$ ns, and $^{108}$Xe: 58$^{+106}_{-23}$ $\mu$s., Comment: 4 pages, 4 figures

Relativistic energy density functional approaches are known to well describe nuclear states which involve alpha clusters. Here, alpha emitting nuclei are analysed through the behavior of the spatial localisation of nucleonic states, calculated with an axially deformed RHB approach over the nuclear chart. The systematic occurrence of more localised valence states, having a n = 1 radial quantum number, allows to pinpoint nuclei in agreement with experimentally known alpha-emitters. The cases of 212Po and 104Te are investigated, showing the concomitant contributions of the pseudospin symmetry and the presence of n = 1 states, on the alpha preformation probability. The impact of the localisation of valence states, on alpha preformation probability, is then analysed. It allows to study shell effects on this probability, over isotopic and isotonic chains. Finally, a phenomenological law is also provided, relating this probability to the radial quantum number of the valence states., Comment: 8 pages, 5 figures

When the density of a nuclear system is decreased, homogeneous states undergo the so-called Mott transition towards clusterised states, e.g. alpha clustering, both in nuclei and in nuclear matter. Here we investigate such a quantum phase transition (QPT) by using microscopic energy density functional (EDF) calculations both with the relativistic and the Gogny approaches on the diluted $^{16}$O nucleus. The evolution of the corresponding single-particle spectrum under dilution is studied, and a Mott-like transition is predicted at about 1/3 of the saturation density. Complementary approaches are used in order to understand this QPT. A study of spatial localisation properties as a function of the density allows to derive a value of the Mott density in agreement with the one obtained by fully microscopic calculations in $^{16}$O and in nuclear matter. Moreover a study of the spontaneous symmetry breaking of the rotational group in $^{16}$O, down to the discrete tetrahedral one, provides further insight on the features displayed by the single-particle spectrum obtained within the EDF approach.The content of the tetrahedrally deformed A-nucleon product state in terms of spherical particle-hole configurations is investigated. Finally a study of quartet condensation and the corresponding macroscopic QPT is undertaken in infinite matter., Comment: 13 pages, 13 figures

We demonstrate that a committee of deep neural networks is capable of predicting the ground-state and excited energies of more than 1800 atomic nuclei with an accuracy akin to the one achieved by state-of-the-art nuclear energy density functionals (EDFs) and a major speed-up. An active learning strategy is proposed to train this algorithm with a minimal set of 210 nuclei. This approach enables future fast studies of the influence of EDFs parametrizations on structure properties over the whole nuclear chart and suggests that for the first time an artificial intelligence successfully encoded the laws of nuclear deformation.

Spectroscopic properties of low-lying states and cluster structures in $^{12}$C are analyzed in a "beyond mean-field framework" based on global energy density functionals (EDFs). To build symmetry-conserving collective states, axially-symmetric and reflection-asymmetric solutions of the relativistic Hartree-Bogoliubov equations are first projected onto good values of angular momentum, particle number, and parity. Configuration mixing is implemented using the generator coordinate method formalism. It is shown that such a global microscopic approach, based on a relativistic EDF, can account for the main spectroscopic features of $^{12}$C, including the ground-state and linear-chain bands as well as, to a certain approximation, the excitation energy of the Hoyle state. The calculated form factors reproduce reasonably well the available experimental values, and display an accuracy comparable to that of dedicated microscopic cluster models., Comment: 10 pages, 6 figures

We describe the first version (v1.0.0) of the code ADG that automatically (1) generates all valid Bogoliubov many-body perturbation theory (BMBPT) diagrams and (2) evaluates their algebraic expression to be implemented for numerical applications. This is achieved at any perturbative order $p$ for a Hamiltonian containing both two-body (four-legs) and three-body (six-legs) interactions (vertices). The automated generation of BMBPT diagrams of order $p$ relies on elements of graph theory, i.e., it is achieved by producing all oriented adjacency matrices of size $(p+1) \times (p+1)$ satisfying topological Feynman's rules. The automated evaluation of BMBPT diagrams of order $p$ relies both on the application of algebraic Feynman's rules and on the identification of a powerful diagrammatic rule providing the result of the remaining $p$-tuple time integral. The diagrammatic rule in question constitutes a novel finding allowing for the straight summation of large classes of time-ordered diagrams at play in the time-independent formulation of BMBPT. Correspondingly, the traditional resolvent rule employed to compute time-ordered diagrams happens to be a particular case of the general rule presently identified. The code ADG is written in Python2.7 and uses the graph manipulation package NetworkX. The code is also able to generate and evaluate Hartree-Fock-MBPT (HF-MBPT) diagrams and is made flexible enough to be expanded throughout the years to tackle the diagrammatics at play in various many-body formalisms that already exist or are yet to be formulated., Comment: 32 pages, 22 figures, 6 tables

The spatial dispersion of the single-nucleon wave functions is analyzed using the self-consistent mean-field framework based on nuclear energy density functionals, and with the harmonic oscillator approximation for the nuclear potential. It is shown that the dispersion depends on the radial quantum number n, but displays only a very weak dependence on the orbital angular momentum. An analytic expression is derived for the localization parameter that explicitly takes into account the radial quantum number of occupied single-nucleon states. The conditions for single-nucleon localization and formation of cluster structures are fulfilled in relatively light nuclei with $A \leq 30$ and $n=1$ states occupied. Heavier nuclei exhibit the quantum liquid phase of nucleonic matter because occupied levels that originate from $n > 1$ spherical states are largely delocalized. Nevertheless, individual $\alpha$-like clusters can be formed from valence nucleons filling single-particle levels originating from $n=1$ spherical mean-field states., Comment: 5 pages, 4 figures

We analyze the localization properties of two-body correlations induced by pairing in the framework of relativistic mean field (RMF) models. The spatial properties of two-body correlations are studied for the pairing tensor in coordinate space and for the Cooper pair wave function. The calculations are performed both with Relativistic-Hatree-Bogoliubov (RHB) and RMF+Projected-BCS (PBCS) models and taking as examples the nuclei $^{66}$Ni, $^{124}$Sn and $^{200}$Pb. It is shown that the coherence length have the same pattern as in previous non-relativistic HFB calculations, i.e., it is maximum in the interior of the nucleus and drops to a minimum in the surface region. In the framework of RMF+PBCS we have also analysed, for the particular case of $^{120}$Sn, the dependence of the coherence length on the intensity of the pairing force. This analysis indicates that pairing is reducing the coherence length by about 25-30 $\%$ compared to the RMF limit., Comment: 8 pages, 7 figures

The lowest positive- and negative-parity bands of $^{20}$Ne and neutron-rich even-even Ne isotopes are investigated using a theoretical framework based on energy density functionals. Starting from a self-consistent relativistic Hartree-Bogoliubov calculation of axially-symmetric and reflection-asymmetric deformation energy surfaces, the collective symmetry-conserving states are built using projection techniques and the generator coordinate method. Overall a good agreement with the experimental excitation energies and transition rates is obtained. In particular, the model provides an accurate description of the excitation spectra and transition probabilities in $^{20}$Ne. The contribution of cluster configurations to the low-energy states is discussed, as well as the transitional character of the ground state. The analysis is extended to $^{22}$Ne and the shape-coexisting isotope $^{24}$Ne, and to the drip-line nuclei $^{32}$Ne and $^{34}$Ne. The role of valence neutrons in the formation of molecular-type bonds between clusters is discussed., Comment: 16 pages, 15 figures, 2 tables, Accepted for Publication in Physical Review C

Background: Ab initio many-body methods whose numerical cost scales polynomially with the number of particles have been developed over the past fifteen years to tackle closed-shell mid-mass nuclei. Open-shell nuclei have been further addressed by implementing variants based on the concept of spontaneous symmetry breaking (and restoration). Purpose: In order to access the spectroscopy of open-shell nuclei more systematically while controlling the numerical cost, we design a novel many-body method that combines the merit of breaking and restoring symmetries with those brought about by low-rank individual excitations. Methods: The recently proposed truncated configuration-interaction method based on optimized symmetry-broken and -restored states is extended to the z-signature symmetry associated with a discrete subgroup of SU(2). The highly-truncated N-body Hilbert subspace within which the Hamiltonian is diagonalized is spanned by a z-signature broken and restored Slater determinant vacuum and associated low-rank excitations. Results: The proposed method provides an excellent reproduction of the ground-state energy and of low-lying excitation energies of various z-signatures and total angular momenta. In doing so, the successive benefits of (i) breaking the symmetry, (ii) restoring the symmetry, (iii) including low-rank particle-hole excitations and (iv) optimizing the amount by which the underlying vacuum breaks the symmetry are illustrated. Conclusions: The numerical cost of the newly designed variational method is polynomial with respect to the system size. The present study confirms the results obtained previously for the attractive pairing Hamiltonian in connection with the breaking and restoration of U(1) global gauge symmetry. These two studies constitute a strong motivation to apply this method to realistic nuclear Hamiltonians., Comment: 14 pages, 10 figures, full abstract available in text

Many properties of the atomic nucleus, such as vibrations, rotations and incompressibility, can be interpreted as due to a two component quantum liquid of protons and neutrons. Electron scattering measurements on stable nuclei demonstrate that their central densities are saturated, as for liquid drops. In exotic nuclei near the limits of mass and charge, with large imbalances in their proton and neutron numbers, the possibility of a depleted central density, or a 'bubble' structure, has been discussed in a recurrent manner since the 1970s. Here we report first experimental evidence that points to a depletion of the central density of protons in the short-lived nucleus 34Si. The proton-to-neutron density asymmetry in 34Si offers the possibility to place constraints on the density and isospin dependence of the spin--orbit force-on which nuclear models have disagreed for decades-and on its stabilizing effect towards limits of nuclear existence., Comment: 6 pages, 3 figures

Dimensionless ratios characterizing many-body systems are a powerful tool to reveal the main universal quantities involved. The recently-introduced localisation parameter allow to study the occurrence of crystal, clusterisation, and quantum liquid states in many-body systems such as nuclei. Its concomitant use with other dimensionless quantities such as the quantality, allows to pinpoint fundamental lengthscales and dimensionless ratios at work in nuclei. Within the present approach, the impact of delocalisation on nuclear saturation is discussed. The transition from a homogeneous to a clusterised state in low density systems is studied. The spin-orbit effect in nuclei and fermionic systems is also reviewed in the light of this analysis. A generalised localisation parameter is derived, showing that the delocalisation properties are not always driven by the quantality. The concepts of quantum fluidity and mobility are finally introduced., Comment: 16 pages, 8 figures

We introduce an improved model that links the frequency shift of the $^{133}\text{Cs}$ hyperfine Zeeman transitions $\vert F = 3, m_F> \longleftrightarrow \vert F = 4, m_F >$ to the Lorentz-violating Standard-Model Extension (SME) coefficients of the proton and neutron. The new model uses Lorentz transformations developed to second order in boost and additionally takes the nuclear structure into account, beyond the simple Schmidt model used previously in SME analyses, thereby providing access to both proton and neutron SME coefficients including the isotropic coefficient $\tilde{c}_{TT}$. Using this new model in a second analysis of the data delivered by the FO2 dual Cs/Rb fountain at Paris Observatory and previously analysed in arXiv:hep-ph/0601024v1, we improve by up to 12 orders of magnitude the present maximum sensitivities (see arXiv:0801.0287v9) on the $\tilde{c}_{Q}$, $\tilde{c}_{TJ}$ and $\tilde{c}_{TT}$ coefficients for the neutron and on the $\tilde{c}_{TT}$ coefficient for the proton, reaching respectively $10^{-20}$, $10^{-17}$, $10^{-13}$ and $10^{-15}$ GeV.

Background: Ab initio many-body methods have been developed over the past ten years to address mid-mass nuclei... As progress in the design of inter-nucleon interactions is made, further efforts must be made to tailor many-body methods. Methods: We formulate a truncated configuration interaction method that consists of diagonalizing the Hamiltonian in a highly truncated subspace of the total N-body Hilbert space. The reduced Hilbert space is generated via the particle-number projected BCS state along with projected seniority-zero two and four quasi-particle excitations. Furthermore, the extent by which the underlying BCS state breaks U(1) symmetry is optimized in presence of the projected two and four quasi-particle excitations... The quality of the newly designed method is tested against exact solutions of the so-called attractive pairing Hamiltonian problem. Results: By construction, the method reproduce exact results for N=2 and N=4. For N=(8,16,20) the error on the ground-state correlation energy is less than (0.006, 0.1, 0.15) % across the entire range of inter-nucleon coupling defining the pairing Hamiltonian and driving the normal-to-superfluid quantum phase transition. The presently proposed method offers the advantage to automatically access the low-lying spectroscopy, which it does with high accuracy. Conclusions: The numerical cost of the newly designed variational method is polynomial (N$^6$) in system size. It achieves an unprecedented accuracy on the ground-state correlation energy, effective pairing gap and one-body entropy as well as on the excitation energy of low-lying states of the attractive pairing Hamiltonian. This constitutes a strong enough motivation to envision its application to realistic nuclear Hamiltonians in view of providing a complementary, accurate and versatile ab initio description of mid-mass open-shell nuclei in the future., Comment: 15 pages, 11 figures, full abstract available in text

Relativistic self-consistent mean-field (SCMF) models naturally account for the coupling of the nucleon spin to its orbital motion, whereas non-relativistic SCMF methods necessitate a phenomenological ansatz for the effective spin-orbit potential. Recent experimental studies aim to explore the isospin properties of the effective spin-orbit interaction in nuclei. SCMF models are very useful in the interpretation of the corresponding data, however standard relativistic mean-field and non-relativistic Hartree-Fock models use effective spin-orbit potentials with different isovector properties, mainly because exchange contributions are not treated explicitly in the former. The impact of exchange terms on the effective spin-orbit potential in relativistic mean-field models is analysed, and it is shown that it leads to an isovector structure similar to the one used in standard non-relativistic Hartree-Fock. Data on the isospin dependence of spin-orbit splittings in spherical nuclei could be used to constrain the isovector-scalar channel of relativistic mean-field models. The reproduction of the empirical kink in the isotope shifts of even Pb nuclei by relativistic effective interactions points to the occurrence of pseudospin symmetry in the single-neutron spectra in these nuclei.

Spin-orbit coupling characterizes quantum systems such as atoms, nuclei, hypernuclei, quarkonia, etc., and is essential for understanding their spectroscopic properties. Depending on the system, the effect of spin-orbit coupling on shell structure is large in nuclei, small in quarkonia, perturbative in atoms. In the standard non-relativistic reduction of the single-particle Dirac equation, we derive a universal rule for the relative magnitude of the spin-orbit effect that applies to very different quantum systems, regardless of whether the spin-orbit coupling originates from the strong or electromagnetic interaction. It is shown that in nuclei the near equality of the mass of the nucleon and the difference between the large repulsive and attractive potentials explains the fact that spin-orbit splittings are comparable to the energy spacing between major shells. For a specific ratio between the particle mass and the effective potential whose gradient determines the spin-orbit force, we predict the occurrence of giant spin-orbit energy splittings that dominate the single-particle excitation spectrum.