BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
L. M. Robledo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The existing formalism used to compute the operator overlaps necessary to carry out generator coordinate method calculations using a set of Hartree- Fock- Bogoliubov wave functions, is generalized to the case where each of the HFB states are expanded in different arbitrary bases spanning different sub-space of the Hilbert space.
We apply a formalism recently developed to carry out Generator Coordinate Method calculations using a set of Hartree- Fock- Bogoliubov wave functions, where each of the members of the set can be expanded in an arbitrary basis. In this paper it is assumed that the HFB wave functions are expanded in Harmonic Oscillator (HO) bases with different oscillator lengths. General expressions to compute the required matrix elements of arbitrary operators are given. The application of the present formalism to the case of fission is illustrated with an example.
The prediction of single particle level crossing phenomenon between $2p_{3/2}$ and $1f_{5/2}$ orbitals in $Ni$- and $Cu$-isotopic chains by the finite range simple effective interaction without requiring the tensor part is discussed. In this case the experimentally observed crossing could be studied as a function of nuclear matter incompressibility, $K(ρ_0)$. The estimated crossing for the neutron number $N$=46 could be reproduced by the equation of state corresponding to $K(ρ_0)$=240 MeV. However, the observed proton gaps between the $1h_{11/2}$ and $1g_{7/2}$ shells in $Sn$ and $Sb$ isotopic chain, and the neutron gaps between the $1i_{13/2}$ and $1h_{9/2}$ shells in $N$=82 isotones, as well as the shell closure properties at $N$=28 require explicit consideration of a tensor part as the central contribution is not enough to initiate the required level splittings.
We compute the $β$-decay $ft$-values within the frameworks of the energy density functional (EDF) and the interacting boson model (IBM). Based on the constrained mean-field calculation with the Gogny-D1M EDF, the IBM Hamiltonian for an even-even nucleus and essential ingredients of the interacting boson-fermion-fermion model (IBFFM) for describing the neighboring odd-odd nucleus are determined in a microscopic way. Only the boson-fermion and residual neutron-proton interaction strengths are determined empirically. The Gamow-Teller (GT) and Fermi (F) transition rates needed to compute the $β$-decay $ft$-values are obtained without any additional parameter or quenching of the $g_A$ factor. The observed $\ft$ values for the $β^+$ decays of the even-even Ba into odd-odd Cs nuclei, and of the odd-odd Cs to the even-even Xe nuclei, with mass $A\approx 130$ are reasonably well described. The predicted GT and F transition rates represent a sensitive test of the quality of the IBM and IBFFM wave functions.
The low-energy excitations and $β$ decays of odd-A nuclei are studied within the interacting boson-fermion model (IBFM), based on the Gogny-D1M nuclear energy density functional (EDF). The constrained Hartree-Fock-Bogoliubov (HFB) approximation is employed to compute potential energy surfaces in terms of triaxial quadrupole degrees of freedom for even-even Xe and Ba nuclei in the mass $A\approx 130$ region. The mean field approximation also provides spherical single-particle energies and occupation probabilities for the neighboring odd-A nuclei. Those quantities represent a microscopic input for spectroscopic calculations in odd-A Xe and Ba, Cs and La isotopes. The Gamow-Teller (GT) and Fermi (F) transition matrix elements, needed to compute $β$-decay $\log{ft}$ values are obtained without any phenomenological fitting. It is shown that both the low-lying states and $β$ decays of the studied odd-A systems are described reasonably well within the employed theoretical framework.
The validity of the rotational formula used to compute E1 and E3 transition strengths in even-even nuclei is analyzed within the Generator Coordinate Method framework based on mean field wave functions. It turns out that those nuclei with spherical or near spherical shapes the E1 and E3 strengths computed with this formula are strongly underestimated and a sound evaluation of them requires angular momentum projected wave functions. Results for several isotopic chains with proton number equal or near magic numbers are analyzed and compared with experimental data. The use of angular momentum projected wave functions greatly improves the agreement with the scarce experimental data.
The constrained Hartree-Fock-Bogoliubov approximation, based on the recent parametrization D1M$^{*}$ of the Gogny energy density functional, is used to describe fission in 435 superheavy nuclei. The Gogny-D1M$^{*}$ parametrization is benchmarked against available experimental data on inner and second barrier heights, excitation energies of the fission isomers and half-lives in a selected set of Pu, Cm, Cf, Fm, No, Rf, Sg, Hs and Fl nuclei. Results are also compared with those obtained with the Gogny-D1M energy density functional. A detailed study of the minimal energy fission paths is carried out for isotopic chains with atomic numbers 100 $\le$ Z $\le$ 126 including very neutron-rich sectors up to around 4 MeV from the two-neutron driplines. Single-particle energies, ground state deformations, pairing correlations, two-nucleon separation energies and barrier heights are also discussed. In addition to fission paths, the constrained Hartree-Fock-Bogoliubov framework provides collective masses and zero-point quantum rotational and vibrational energies. Those quantities are building blocks within the Wentzel-Kramer-Brillouin formalism employed to evaluate the systematic of the spontaneous fission half-lives t$_\mathrm{SF}$. The competition between spontaneous fission and $α$-decay is studied, through the computation of the $α$-decay half-lives t$_\mathrmα$ using a parametrization of the Viola-Seaborg formula. From the comparison with the available experimental data and the results obtained with other theoretical approaches, it is concluded that D1M$^{*}$ represents a reasonable starting point to describe fission in heavy and superheavy nuclei.
The quadrupole-octupole coupling and the related spectroscopic properties have been studied for the even-even light actinides $^{218-238}$Ra and $^{220-240}$Th. The Hartree-Fock-Bogoliubov approximation, based on the Gogny-D1M energy density functional, has been employed as a microscopic input, i.e., to obtain (axially symmetric) mean-field potential energy surfaces as functions of the quadrupole and octupole deformation parameters. The mean-field potential energy surfaces have been mapped onto the corresponding bosonic potential energy surfaces using the expectation value of the $sdf$ Interacting Boson Model (IBM) Hamiltonian in the boson condensate state. The strength parameters of the $sdf$-IBM Hamiltonian have been determined via this mapping procedure. The diagonalization of the mapped IBM Hamiltonian provides energies for positive- and negative-parity states as well as wave functions which are employed to obtain transitional strengths. The results of the calculations compare well with available data from Coulomb excitation experiments and point towards a pronounced octupole collectivity around $^{224}$Ra and $^{226}$Th.
The spectroscopic properties of the odd-odd isotopes $^{124-132}$Cs have been studied within the interacting boson-fermion-fermion model based on the Gogny-D1M energy density functional framework. Major ingredients to build the interacting boson-fermion-fermion Hamiltonian, such as the ($β,γ$)-deformation energy surfaces for the even-even core nuclei $^{124-132}$Xe as well as single-particle energies and occupation probabilities of the odd nucleons, have been computed microscopically with the constrained Hartree-Fock-Bogoliubov method. A few coupling constants of the boson-fermion and residual neutron-proton interactions are fitted to reproduce with a reasonable accuracy the experimental excitation energy of the low-lying levels of the odd-mass and odd-odd nuclei. The method is applied to describe the low-energy low-spin spectra of the odd-odd Cs nuclei and the band structures of higher-spin higher-energy states, mainly based on the $(νh_{11/2})^{-1}\otimesπh_{11/2}$ configuration. Many of those odd-odd Cs nuclei have been identified as candidates for exhibiting chiral doublet bands.
Constrained mean-field calculations, based on the Gogny-D1M energy density functional, have been carried out to describe fission in the isotopes $^{250-260}$No. The even-even isotopes have been considered within the standard Hartree-Fock-Bogoliobov (HFB) framework while for the odd-mass ones the Equal Filling Approximation (HFB-EFA) has been employed. Ground state quantum numbers and deformations, pairing energies, one-neutron separation energies, inner and outer barrier heights as well as fission isomer excitation energies are given. Fission paths, collective masses and zero-point quantum vibrational and rotational corrections are used to compute the systematic of the spontaneous fission half-lives t$_\mathrm{SF}$ both for even-even and odd-mass nuclei. Though there exists a strong variance of the predicted fission rates with respect to the details involved in their computation, it is shown that both the specialization energy and the pairing quenching effects, taken into account within the self-consistent HFB-EFA blocking procedure, lead to larger t$_\mathrm{SF}$ values in odd-mass nuclei as compared with their even-even neighbors. Alpha decay lifetimes have also been computed using a parametrization of the Viola-Seaborg formula. The high quality of the Gogny-D1M functional regarding nuclear masses leads to a very good reproduction of $Q_α$ values and consequently of lifetimes.
The structure of even-even neutron-rich Ru, Mo, Zr and Sr nuclei in the $A\approx 100$ mass region is studied within the interacting boson model (IBM) with microscopic input from the self-consistent mean-field approximation based on the Gogny-D1M energy density functional. The deformation energy surface in the quadrupole deformation space $(β,γ)$, computed within the constrained Hartree-Fock-Bogoliubov framework, is mapped onto the expectation value of the appropriately chosen IBM Hamiltonian with configuration mixing in the boson condensate state. The mapped IBM Hamiltonian is used to study the spectroscopic properties of $^{98-114}$Ru, $^{96-112}$Mo, $^{94-110}$Zr and $^{92-108}$Sr. Several cases of $γ$-soft behavior are predicted in Ru and Mo nuclei while a pronounced coexistence between strongly-prolate and weakly-oblate deformed shapes is found for Zr and Sr nuclei. The method describes well the evolution of experimental yrast and non-yrast states as well as selected $B$(E2) transition probabilities.
The available data for E2 transition strengths in the region between neutron-deficient Hf and Pt isotopes are far from complete. More and precise data are needed to enhance the picture of structure evolution in this region and to test state-of-the-art nuclear models. In a simple model, the maximum collectivity is expected at the middle of the major shell. However, for actual nuclei, this picture may no longer be the case, and one should use a more realistic nuclear-structure model. We address this point by studying the spectroscopy of Hf. We remeasure the 2^+_1 half-lives of 172,174,176Hf, for which there is some disagreement in the literature. The main goal is to measure, for the first time, the half-lives of higher-lying states of the rotational band. The new results are compared to a theoretical calculation for absolute transition strengths. The half-lives were measured using γ-γand conversion-electron-γdelayed coincidences with the fast timing method. For the determination of half-lives in the picosecond region, the generalized centroid difference method was applied. For the theoretical calculation of the spectroscopic properties, the interacting boson model is employed, whose Hamiltonian is determined based on microscopic energy-density functional calculations. The measured 2^+_1 half-lives disagree with results from earlier γ-γfast timing measurements, but are in agreement with data from Coulomb excitation experiments and other methods. Half-lives of the 4^+_1 and 6^+_1 states were measured, as well as a lower limit for the 8^+_1 states. We show the importance of the mass-dependence of effective boson charge in the description of E2 transition rates in chains of nuclei. It encourages further studies of the microscopic origin of this mass dependence. New data on transition rates in nuclei from neighboring isotopic chains could support these studies.
This is a short review of methods and results of calculations of fission barriers and fission half-lives of even-even superheavy nuclei. An approvable agreement of the following approaches is shown and discussed: The macroscopic-microscopic approach based on the stratagem of the shell correction to the liquid drop model and a vantage point of microscopic energy density functionals of Skyrme and Gogny type selfconsistently calculated within Hartree-Fock-Bogoliubov method. Mass parameters are calculated in the Hartree-Fock-Bogoliubov cranking approximation. A short part of the paper is devoted to the nuclear fission dynamics. We also discuss the predictive power of Skyrme functionals applied to key properties of the fission path of $^{266}$Hs. It applies the standard techniques of error estimates in the framework of a $χ^2$ analysis.
Mean field calculations, based on the D1S, D1N and D1M parametrizations of the Gogny energy density functional, have been carried out to obtain the potential energy surfaces relevant to fission in several Ra isotopes with the neutron number 144 $\le$ N $\le$ 176. Inner and outer barrier heights as well as first and second isomer excitation energies are given. The existence of a well developed third minimum along the fission paths of Ra nuclei, is analyzed in terms of the energetics of the "fragments" defining such elongated configuration. The masses and charges of the fission fragments are studied as functions of the neutron number in the parent Ra isotope. The comparison between fission and $α$-decay half-lives, reveals that the former becomes faster for increasing neutron numbers. Though there exists a strong variance of the results with respect to the parameters used in the computation of the spontaneous fission rate, a change in tendency is observed at N=164 with a steady increase that makes heavier neutron-rich Ra isotopes stable against fission, diminishing the importance of fission recycling in the r-process.
We calculate properties of the ground and excited states of nuclei in the nobelium region for proton and neutron numbers of 92 <= Z <= 104 and 144 <= N <= 156, respectively. We use three different energy-density-functional (EDF) approaches, based on covariant, Skyrme, and Gogny functionals, each within two different parameter sets. A comparative analysis of the results obtained for odd-even mass staggerings, quasiparticle spectra, and moments of inertia allows us to identify single-particle and shell effects that are characteristic to these different models and to illustrate possible systematic uncertainties related to using the EDF modelling
Collective quadrupole and octupole states are described in a series of Sm and Gd isotopes within the framework of the interacting boson model (IBM), whose Hamiltonian parameters are deduced from mean field calculations with the Gogny energy density functional. The link between both frameworks is the ($β_2β_3$) potential energy surface computed within the Hartree-Fock-Bogoliubov framework in the case of the Gogny force. The diagonalization of the IBM Hamiltonian provides excitation energies and transition strengths of an assorted set of states including both positive and negative parity states. The resultant spectroscopic properties are compared with the available experimental data and also with the results of the configuration mixing calculations with the Gogny force within the generator coordinate method (GCM). The structure of excited $0^{+}$ states and its connection with double octupole phonons is also addressed. The model is shown to describe the empirical trend of the low-energy quadrupole and octupole collective structure fairly well, and turns out to be consistent with GCM results obtained with the Gogny force.
Given the promising features of the recently proposed Barcelona-Catania-Paris (BCP) functional \cite{Baldo.08}, it is the purpose of this paper to still improve on it. It is, for instance, shown that the number of open parameters can be reduced from 4-5 to 2-3, i.e. by practically a factor of two. One parameter is tightly fixed by a fine-tuning of the bulk, a second by the surface energy. The third is the strength of the spin-orbit potential on which the final result does not depend within the scatter of the values used in Skyrme and Gogny like functionals. An energy rms value of 1.58 MeV is obtained from a fit of these three parameters to the 579 measured masses reported in the Audi and Waspra 2003 compilation. This rms value compares favorably with the one obtained using other successful mean field theories. Charge radii are also well reproduced when compared with experiment. The energies of some excited states, mostly the isoscalar giant monopole resonances, are studied within this model as well.
Levels densities of independent-particle Hamiltonians can be calculated easily by using the real-time representation of the evolution operator together with the fast Fourier transform. We describe the method and implement it with a set of Python programs. Examples are provided for the total and partial levels densities of a heavy deformed nucleus (Dy-164). The partial level densities that may be calculated are the projected ones on neutron number, proton number, azimuthal angular momentum, and parity.
The most recent parametrization D1M of the Gogny energy density functional is used to describe fission in the isotopes $^{232-280}$ Pu. We resort to the methodology introduced in our previous studies [Phys. Rev. C \textbf{88}, 054325 (2013) and Phys. Rev. C \textbf {89}, 054310 (2014)] to compute the fission paths, collective masses and zero point quantum corrections within the Hartree-Fock-Bogoliubov framework. The systematics of the spontaneous fission half-lives t$_{SF}$, masses and charges of the fragments in Plutonium isotopes is analyzed and compared with available experimental data. We also pay attention to isomeric states, the deformation properties of the fragments as well as to the competition between the spontaneous fission and $α$-decay modes. The impact of pairing correlations on the predicted t$_{SF}$ values is demonstrated with the help of calculations for $^{232-280}$Pu in which the pairing strengths of the Gogny-D1M energy density functional are modified by 5 $\%$ and 10 $\%$, respectively. We further validate the use of the D1M parametrization through the discussion of the half-lives in $^{242-262}$Fm. Our calculations corroborate that, though the uncertainties in the absolute values of physical observables are large, the Gogny-D1M Hartree-Fock-Bogoliubov framework still reproduces the trends with mass and/or neutron numbers and therefore represents a reasonable starting point to describe fission in heavy nuclear systems from a microscopic point of view.
The most recent parametrizations D1S, D1N and D1M of the Gogny energy density functional are used to describe fission in the isotopes $^{232-280}$ U. Fission paths, collective masses and zero point quantum corrections, obtained within the constrained Hartree-Fock-Bogoliubov approximation, are used to compute the systematics of the spontaneous fission half-lives $t_\mathrm{SF}$, the masses and charges of the fission fragments as well as their intrinsic shapes. The Gogny-D1M parametrization has been benchmarked against available experimental data on inner and second barrier heights, excitation energies of the fission isomers and half-lives in a selected set of Pu, Cm, Cf, Fm, No, Rf, Sg, Hs and Fl nuclei. It is concluded that D1M represents a reasonable starting point to describe fission in heavy and superheavy nuclei. Special attention is also paid to understand the uncertainties in the predicted $t_\mathrm{SF}$ values arising from the different building blocks entering the standard semi-classical Wentzel-Kramers-Brillouin formula. Although the uncertainties are large, the trend with mass or neutron numbers are well reproduced and therefore the theory still has predictive power. In this respect, it is also shown that modifications of a few per cent in the pairing strength can have a significant impact on the collective masses leading to uncertainties in the $t_\mathrm{SF}$ values of several orders of magnitude.
The quadrupole collectivity in Nd, Sm, Gd, Dy, Er, Yb, Hf and W nuclei with neutron numbers 122 $\le$ N $\le$ 156 is studied, both at the mean field level and beyond, using the Gogny energy density functional. Besides the robustness of the N=126 neutron shell closure, it is shown that the onset of static deformations in those isotopic chains with increasing neutron number leads to an enhanced stability and further extends the corresponding two-neutron driplines far beyond what could be expected from spherical calculations. Independence of the mean field predictions with respect to the particular version of the Gogny energy density functional employed is demonstrated by comparing results based on the D1S and D1M parameter sets. Correlations beyond mean field are taken into account in the framework of the angular momentum projected generator coordinate method calculation. It is shown that N=126 remains a robust neutron magic number when dynamical effects are included. The analysis of the collective wave functions, average deformations and excitation energies indicate that, with increasing neutron number, the zero-point quantum corrections lead to dominant prolate configurations in the 0$_{1}^{+}$, 0$_{2}^{+}$, 2$_{1}^{+}$ and 2$_{2}^{+}$ states of the studied nuclei. Moreover, those dynamical deformation effects provide an enhanced stability that further supports the mean field predictions, corroborating a shift of the r-process path to higher neutron numbers. Beyond mean field calculations provide a smaller shell gap at N=126 than the mean field one in good agreement with previous theoretical studies. However, the shell gap still remains strong enough in the two-neutron driplines.
The evolution of the total energy surface and the nuclear shape in the isotopic chain $^{172-194}$Pt are studied in the framework of the interacting boson model, including configuration mixing. The results are compared with a self-consistent Hartree-Fock-Bogoliubov calculation using the Gogny-D1S interaction and a good agreement between both approaches shows up. The evolution of the deformation parameters points towards the presence of two different coexisting configurations in the region 176 $\leq$ A $\leq$ 186.
The interacting boson model with configuration mixing, with parameters derived from the self-consistent mean-field calculation employing the microscopic Gogny energy density functional, is applied to the systematic analysis of the low-lying structure in Hg isotopes. Excitation energies, electromagnetic transition rates, deformation properties, and ground-state properties of the $^{172-204}$Hg nuclei are obtained by mapping the microscopic deformation energy surface onto the equivalent IBM Hamiltonian in the boson condensate. These results point to the overall systematic trend of the transition from the near spherical vibrational state in lower-mass Hg nuclei close to $^{172}$Hg, onset of intruder prolate configuration as well as the manifest prolate-oblate shape coexistence around the mid-shell nucleus $^{184}$Hg, weakly oblate deformed structure beyond $^{190}$Hg up to the spherical vibrational structure toward the near semi-magic nucleus $^{204}$Hg, as observed experimentally. The quality of the present method in the description of the complex shape dynamics in Hg isotopes is examined.
The Gallagher-Moszkowski rule in the spectroscopy of odd-odd nuclei imposes a new spin constraint on the energy functionals for self-consistent mean field theory. The commonly used parameterization of the effective three-body interaction in the Gogny and Skyrme families of energy functionals is ill-suited to satisfy the spin constraint. In particular, the Gogny parameterization of the three-body interaction has the opposite spin dependence to that required by the observed spectra. The two-body part has a correct sign, but in combination the rule is violated as often as not. We conclude that a new functional form is needed for the effective three-body interaction that can take into better account the different spin-isospin channels of the interaction.
Spectroscopic observables such as electromagnetic transitions strengths can be related to the properties of the intrinsic mean-field wave function when the latter are strongly deformed, but the standard rotational formulas break down when the deformation decreases. Nevertheless there is a well-defined, non-zero, spherical limit that can be evaluated in terms of overlaps of mean-field intrinsic deformed wave functions. We examine the transition between the spherical limit and strongly deformed one for a range of nuclei comparing the two limiting formulas with exact projection results. We find a simple criterion for the validity of the rotational formula depending on $<Δ\vec{J}^2>$, the mean square fluctuation in the angular momentum of the intrinsic state. We also propose an interpolation formula which describes the transition strengths over the entire range of deformations, reducing to the two simple expressions in the appropriate limits.
The interplay between the collective dynamics of the quadrupole and octupole deformation degree of freedom is discussed in a series of Sm and Gd isotopes both at the mean field level and beyond, including parity symmetry restoration and configuration mixing. Physical properties like negative parity excitation energies, E1 and E3 transition probabilities are discussed and compared to experimental data. Other relevant intrinsic quantities like dipole moments, ground state quadrupole moments or correlation energies associated to symmetry restoration and configuration mixing are discussed. For the considered isotopes, the quadrupole-octupole coupling is found to be weak and most of the properties of negative parity states can be described in terms of the octupole degree of freedom alone.
We discuss the octupole deformation properties of many even-even actinide isotopes. The analysis is carried out in the mean field framework (Hartree-Fock-Bogoliubov approximation) by using the axially symmetric octupole moment as a constraint. A one-dimensional octupole collective Hamiltonian is used to obtain properties like excitation energies or $E1$ and $E3$ transition probabilities of the negative parity band-heads associated to the lowest lying $1^{-}$ and $3^{-}$ states. The evolution of these values with neutron number is discussed and a comparison with available experimental data is made. In order to minimize the uncertainties associated to the energy density functional used, the calculations have been carried out for an assorted set ranging from the BCP1 functional to the finite range Gogny interaction with the D1S, D1N and D1M parametrization.
As part of a program to study odd-A nuclei in the Hartree-Fock-Bogoliubov (HFB) theory, we have developed a new calculational tool to find the HFB minima of odd-A nuclei based on the gradient method and using interactions of Gogny's form. The HFB minimization includes both time-even and time-odd fields in the energy functional, avoiding the commonly used "filling approximation". Here we apply the method to calculate neutron pairing gaps in some representative isotope chains of spherical and deformed nuclei, namely the Z=8,50 and 82 spherical chains and the Z=62 and 92 deformed chains. We find that the gradient method is quite robust, permitting us to carry out systematic surveys involving many nuclei. We find that the time-odd field does not have large effect on the pairing gaps calculated with the Gogny D1S interaction. Typically, adding the T-odd field as a perturbation increases the pairing gap by ~100 keV, but the re-minimization brings the gap back down. This outcome is very similar to results reported for the Skyrme family of nuclear energy density functionals. Comparing the calculated gaps with the experimental ones, we find that the theoretical errors have both signs implying that the D1S interaction has a reasonable overall strength. However, we find some systematic deficiencies comparing spherical and deformed chains and comparing the lighter chains with the heavier ones. The gaps for heavy spherical nuclei are too high, while those for deformed nuclei tend to be too low. The calculated gaps of spherical nuclei show hardly any A-dependence, contrary to the data. Inclusion of the T-odd component of the interaction does not change these qualitative findings.
We investigate the emergence and evolution of shape coexistence in the neutron-deficient Lead isotopes within the interacting boson model (IBM) plus configuration mixing with microscopic input based on the Gogny energy density functional (EDF). The microscopic potential energy surface obtained from the constrained self-consistent Hartree-Fock-Bogoliubov method employing the Gogny-D1M EDF is mapped onto the coherent-state expectation value of the configuration-mixing IBM Hamiltonian. In this way, the parameters of the IBM Hamiltonian are fixed for each of the three relevant configurations (spherical, prolate and oblate) associated to the mean field minima. Subsequent diagonalization of the Hamiltonian provides the excitation energy of the low-lying states and transition strengths among them. The model predictions for the $0^{+}$ level energies and evolving shape coexistence in the considered Lead chain are consistent both with experiment and with the indications of the Gogny-EDF energy surfaces.
A computer code is presented for solving the equations of Hartree-Fock-Bogoliubov (HFB) theory by the gradient method, motivated by the need for efficient and robust codes to calculate the configurations required by extensions of HFB such as the generator coordinate method. The code is organized with a separation between the parts that are specific to the details of the Hamiltonian and the parts that are generic to the gradient method. This permits total flexibility in choosing the symmetries to be imposed on the HFB solutions. The code solves for both even and odd particle number ground states, the choice determined by the input data stream. Application is made to the nuclei in the $sd$-shell using the USDB shell-model Hamiltonian.
An interacting boson model Hamiltonian determined from Hartree-Fock-Bogoliubov calculations with the new microscopic Gogny energy density functional D1M, is applied to the spectroscopic analysis of neutron-rich Yb, Hf, W, Os and Pt isotopes with mass $A\sim 180-200$. Excitation energies and transition rates for the relevant low-lying quadrupole collective states are calculated by this method. Transitions from prolate to oblate ground-state shapes are analyzed as a function of neutron number $N$ in a given isotopic chain by calculating excitation energies, $B$(E2) ratios, and correlation energies in the ground state. It is shown that such transitions tend to occur more rapidly for the isotopes with lower proton number $Z$, when departing from the proton shell closure Z=82. The triaxial degrees of freedom turn out to play an important role in describing the considered mass region. Predicted low-lying spectra for the neutron-rich exotic Hf and Yb isotopes are presented. The approximations used in the model and the possibilities to refine its predictive power are addressed.
We present a new exactly solvable Hamiltonian with a separable pairing interaction and non-degenerate single-particle energies. It is derived from the hyperbolic family of Richardson-Gaudin models and possesses two free parameters, one related to an interaction cutoff and the other to the pairing strength. These two parameters can be adjusted to give an excellent reproduction of Gogny self-consistent mean-field calculations in the canonical basis.
We present a computational methodology for a theory of the lowest octupole excitations applicable to all even-even nuclei beyond the lightest. The theory is the well-known generator-coordinate extension (GCM) of the Hartree-Fock-Bogoliubov self-consistent mean field theory (HFB). We use the discrete-basis Hill-Wheel method (HW) to compute the wave functions with an interaction from the Gogny family of Hamiltonians. Comparing to the compiled experimental data on octupole excitations, we find that the performance of the theory depends on the deformation characteristics of the nucleus. For nondeformed nuclei, the theory reproduces the energies to about 20 % apart from an overall scale factor of about 1.6. The performance is somewhat poorer for (quadrupole) deformed nuclei, and for both together the dispersion of the scaled energies about the experimental values is about 25 %. This compares favorably with the performance of similar theories of the quadrupole excitations. Nuclei having static octupole deformations in HFB form a special category. These nuclei have the smallest measured octupole excitation energies as well as the smallest predicted energies. However, in these cases the energies are seriously underpredicted by the theory. We find that a simple two-configuration approximation, the Minimization After Projection method, (MAP) is almost as accurate as the full HW treatment, provided that the octupole-deformed nuclei are omitted from the comparison. This article is accompanied by a tabulation of the predicted octupole excitations for 818 nuclei extending from dripline to dripline, computed with several variants of the Gogny interaction.
Cluster radioactivity is the emission of a fragment heavier than $α$ particle and lighter than mass 50. The range of clusters observed in experiments goes from $^{14}$C to $^{32}$Si while the heavy mass residue is always a nucleus in the neighborhood of the doubly-magic $^{208}$Pb nucleus. Cluster radioactivity is described in this paper as a very asymmetric nuclear fission. A new fission valley leading to a decay with large fragment mass asymmetry matching the cluster radioactivity products is found. The mass octupole moment is found to be more convenient than the standard quadrupole moment as the parameter driving the system to fission. The mean-field HFB theory with the phenomenological Gogny interaction has been used to compute the cluster emission properties of a wide range of even-even actinide nuclei from $^{222}$Ra to $^{242}$Cm, where emission of the clusters has been experimentally observed. Computed half-lives for cluster emission are compared with experimental results. The noticeable agreement obtained between the predicted properties of cluster emission (namely, clusters masses and emission half-lives) and the measured data confirms the validity of the proposed methodology in the analysis of the phenomenon of cluster radioactivity. A continuous fission path through the scission point has been described using the neck parameter constraint.
The isotopic evolution of the ground-state nuclear shapes and the systematics of one-quasiproton configurations are studied in neutron-rich odd-A Yttrium and Niobium isotopes. We use a selfconsistent Hartree-Fock-Bogoliubov formalism based on the Gogny energy density functional with two parametrizations, D1S and D1M. The equal filling approximation is used to describe odd-A nuclei preserving both axial and time reversal symmetries. Shape-transition signatures are identified in the N=60 isotopes in both charge radii and spin-parities of the ground states. These signatures are a common characteristic for nuclei in the whole mass region. The nuclear deformation and shape coexistence inherent to this mass region are shown to play a relevant role in the understanding of the spectroscopic features of the ground and low-lying one-quasiproton states. Finally, a global picture of the neutron-rich A=100 mass region from Krypton up to Molybdenum isotopes is illustrated with the systematics of the nuclear charge radii isotopic shifts.
Structural evolution in neutron-rich Os and W isotopes is investigated in terms of the Interacting Boson Model (IBM) Hamiltonian determined by (constrained) Hartree-Fock-Bogoliubov (HFB) calculations with the Gogny-D1S Energy Density Functional (EDF). The interaction strengths of the IBM Hamiltonian are produced by mapping the potential energy surface (PES) of the Gogny-EDF with quadrupole degrees of freedom onto the corresponding PES of the IBM system. We examine the prolate-to-oblate shape/phase transition which is predicted to take place in this region as a function of neutron number $N$ within the considered Os and W isotopic chains. The onset of this transition is found to be more rapid compared to the neighboring Pt isotopes. The calculations also allow prediction of spectroscopic variables (excited state energies and reduced transition probabilities) which are presented for the neutron-rich $^{192,194,196}$W nuclei, for which there is only very limited experimental data available to date.
We present a pfaffian formula for projection and symmetry restoration for wave functions of the general Bogoliubov form, including quasiparticle excited states and linear combinations of them. This solves a long-standing problem in calculating states of good symmetry, arising from the sign ambiguity of the commonly used determinant formula. A simple example is given of projecting good particle number and angular momentum from a Bogoliubov wave function in the Fock space of a single j-shell.
Several technical aspects concerning the evaluation of the overlap between two mean field wave functions of the Hartree Fock Bogoliubov type, are discussed. The limit when several orbitals become fully occupied is derived as well as the formula to reduce the dimensionality of the problem when exactly empty orbitals are present. The formalism is also extended to deal with the case where the bases of each of the wave functions are different. Several practical results concerning the evaluation of pfaffians as well as the canonical decomposition of norm overlaps are also discussed in the appendices.
The evolution of the ground-state nuclear shapes in neutron-rich Sr, Zr, and Mo isotopes, including both even-even and odd-A nuclei, is studied within a self-consistent mean-field approximation based on the D1S Gogny interaction. Neutron separation energies and charge radii are calculated and compared with available data. A correlation between a shape transition and a discontinuity in those observables is found microscopically. While in Sr and Zr isotopes the steep behavior observed in the isotopic dependence of the charge radii is a consequence of a sharp prolate-oblate transition, the smooth behavior found in Mo isotopes has its origin in an emergent region of triaxiality.
The evolution of the nuclear shapes along the triaxial landscape is studied in the Pt isotopic chain using the selfconsistent Hartree-Fock-Bogoliubov approximation based on the Gogny interaction. In addition to the parametrization D1S, the new incarnations D1N and D1M of this force are also included in our analysis to asses to which extent the predictions are independent of details of the effective interaction. The considered range of neutron numbers 88<N<26 includes prolate, triaxial, oblate and spherical ground state shapes and serves for a detailed comparison of the predictions obtained with the new sets D1N and D1M against the ones provided by the standard parametrization Gogny-D1S in a region of the nuclear landscape for which experimental and theoretical fingerprints of shape transitions have been found. Structural evolution along the Pt chain is discussed in terms of the deformation dependence of single particle energies.
Evaluation of pfaffians arises in a number of physics applications, and for some of them a direct method is preferable to using the determinantal formula. We discuss two methods for the numerical evaluation of pfaffians. The first is tridiagonalization based on Householder transformations. The main advantage of this method is its numerical stability that makes unnecessary the implementation of a pivoting strategy. The second method considered is based on Aitken's block diagonalization formula. It yields to a kind of LU (similar to Cholesky's factorization) decomposition (under congruence) of arbitrary skew-symmetric matrices that is well suited both for the numeric and symbolic evaluations of the pfaffian. Fortran subroutines (FORTRAN 77 and 90) implementing both methods are given. We also provide simple implementations in Python and Mathematica for purpose of testing, or for exploratory studies of methods that make use of pfaffians.
We discuss the octupole deformation properties of the recently proposed Barcelona-Catania-Paris (BCP) energy density functionals for two sets of isotopes, those of radium and barium, where it is believed that octupole deformation plays a role in the description of the ground state. The analysis is carried out in the mean field framework (Hartree- Fock- Bogoliubov approximation) by using the axially symmetric octupole moment as a constraint. The main ingredients entering the octupole collective Hamiltonian are evaluated and the lowest lying octupole eigenstates are obtained. In this way we restore, in an approximate way, the parity symmetry spontaneously broken by the mean field and also incorporate octupole fluctuations around the ground state solution. For each isotope the energy of the lowest lying $1^{-}$state and the $B(E1)$ and $B(E3)$ transition probabilities have been computed and compared to both the experimental data and the results obtained in the same framework with the Gogny D1S interaction, which are used here as a well established benchmark. Finally, the octupolarity of the configurations involved in the way down to fission of $^{240}$Pu, which is strongly connected to the asymmetric fragment mass distribution, is studied. We confirm with this thorough study the suitability of the BCP functionals to describe octupole related phenomena.
Two-body matrix elements of arbitrary local interactions are written as the sum of separable terms in a way that is well suited for the exchange and pairing channels present in mean-field calculations. The expansion relies on the transformation to center of mass and relative coordinate (in the spirit of Talmi's method) and therefore it is only useful (finite number of expansion terms) for harmonic oscillator single particle states. The converge of the expansion with the number of terms retained is studied for a Gaussian two body interaction. The limit of a contact (delta) force is also considered. Ways to handle the general case are also discussed.
The isotopic evolution of the ground-state nuclear shapes and the systematics of one-quasiproton configurations are studied in odd-A Rubidium isotopes. We use a selfconsistent Hartree-Fock-Bogoliubov formalism based on the Gogny energy density functional with two parametrizations, D1S and D1M, and implemented with the equal filling approximation. We find clear signatures of a sharp shape transition at N=60 in both charge radii and spin-parity of the ground states, which are robust, consistent to each other, and in agreement with experiment. We point out that the combined analysis of these two observables could be used to predict unambiguously new regions where shape transitions might develop.
Spectroscopic calculations are carried out, for the description of the shape/phase transition in Pt nuclei in terms of the Interacting Boson Model (IBM) Hamiltonian derived from (constrained) Hartree-Fock-Bogoliubov (HFB) calculations with the finite range and density dependent Gogny-D1S Energy Density Functional. Assuming that the many-nucleon driven dynamics of nuclear surface deformation can be simulated by effective bosonic degrees of freedom, the Gogny-D1S potential energy surface (PES) with quadrupole degrees of freedom is mapped onto the corresponding PES of the IBM. Using this mapping procedure, the parameters of the IBM Hamiltonian, relevant to the low-lying quadrupole collective states, are derived as functions of the number of valence nucleons. Merits of both Gogny-HFB and IBM approaches are utilized so that the spectra and the wave functions in the laboratory system are calculated precisely. The experimental low-lying spectra of both ground-state and side-band levels are well reproduced. From the systematics of the calculated spectra and the reduced E2 transition probabilities $B$(E2), the prolate-to-oblate shape/phase transition is shown to take place quite smoothly as a function of neutron number $N$ in the considered Pt isotopic chain, for which the $γ$-softness plays an essential role. All these spectroscopic observables behave consistently with the relevant PESs and the derived parameters of the IBM Hamiltonian as functions of $N$. Spectroscopic predictions are also made for those nuclei which do not have enough experimental E2 data.
The systematics of one-quasiparticle configurations in neutron-rich Sr, Zr, and Mo odd isotopes is studied within the Hartree-Fock-Bogoliubov plus Equal Filling Approximation method preserving both axial and time reversal symmetries. Calculations based on the Gogny energy density functional with both the standard D1S parametrization and the new D1M incarnation of this functional are included in our analysis. The nuclear deformation and shape coexistence inherent to this mass region are shown to play a relevant role in the understanding of the spectroscopic features of the ground and low-lying one-quasineutron states.
The behavior of several nuclear properties with temperature is analyzed within the framework of the Finite Temperature Hartree-Fock-Bogoliubov (FTHFB) theory with the Gogny force and large configuration spaces. Thermal shape fluctuations in the quadrupole degree of freedom, around the mean field solution, are taken into account with the Landau prescription. As representative examples the nuclei $^{164}$Er, $^{152}$Dy and $^{192}$Hg are studied. Numerical results for the superfluid to normal and deformed to spherical shape transitions are presented. We found a substantial effect of the fluctuations on the average value of several observables. In particular, we get a decrease in the critical temperature ($T_c$) for the shape transition as compared with the plain FTHFB prediction as well as a washing out of the shape transition signatures. The new values of $T_c$ are closer to the ones found in Strutinsky calculations and with the Pairing Plus Quadrupole model Hamiltonian.
The particle number projection method is formulated for density dependent forces and in particular for the finite range Gogny force. Detailed formula for the projected energy and its gradient are provided. The problems arising from the neglection of any exchange term, which may lead to divergences, are throughly discussed and the possible inaccuracies estimated. Numericala results for the projection after variation method are shown for the nucleus 164Er and for the projection before variation approach for the nuclei 48-50Cr. We also confirm the Coulomb antipairing effect found in mean field theories.
We derive the equations for approximate particle number projection based on mean field wave functions with finite range density dependent forces. As an application ground bands of even-A superdeformed nuclei in the A=150 and A=190 regions are calculated with the Gogny force. We discuss nuclear properties such as quadrupole moments, moments of inertia and quasiparticle spectra, among others, as a function of the angular momentum. We obtain a good overall description.