Shape isomers in Pt, Hg and Pb isotopes with N≤126\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {N}}\le 126$$\end{document}

Deformation-energy surfaces of 54 even-even isotopes of Pt, Hg and Pb nuclei with neutron numbers up to 126 are investigated within a macroscopic-microscopic model based on the Lublin-Strasbourg-Drop macroscopic energy and shell plus pairing-energy corrections obtained from a Yukawa-folded mean-field potential at the desired deformation. A new, rapidly converging Fourier shape parametrization is used to describe nuclear shapes. The stability of shape isomeric states with respect to non-axial and higher-order deformations is investigated.


Introduction
Recent experimental investigations studying the structure of neutron deficient Lead and Mercury isotopes evidenced a shape coexistence phenomenon (see e.g. Refs. [1,2]) thus giving a strong impetus to theoreticians to explore this region of nuclei in more detail. Important advances in the study of neutron-deficient Z ≈ 82 nuclei have been realized using tagging techniques in the Accelerator Laboratory of the University of Jyväskylä, Finland [3], through Coulomb-excitation experiments undertaken at the REX-ISOLDE facility in CERN [4], or isomeric and β-decays of nuclei populated by relativistic energy-fragmentation experiments [5] performed at GSI within the RISING campaign. Several interesting theoretical studies on that subject, based either on a self-consistent approach (see e.g. Refs. [6][7][8][9]) or on macroscopic-microscopic models (see e.g. [10]), have been published in the last two decades.
The aim of the current paper is to present a systematic study of the potential energy surfaces of even-even Pt, Hg, and Pb isotopes with neutron number N ≤ 126 using a very efficient shape parametrization that we recently developed [11,12]. We have chosen to carry out our study using a e-mail: pomorski@kft.umcs.lublin.pl (corresponding author) the macroscopic-microscopic (mic-mac) approach with the Lublin-Strasbourg Drop (LSD) [13], able to reproduce accurately both nuclear masses and fission barriers, and a Yukawafolded single-particle potential [14,15] to account for shell and pairing-energy corrections. A similar kind of investigation of nuclear potential-energy surfaces for a wide range of nuclei has been carried out in Ref. [16] using the (ε 2 , γ, ε 4 ) deformation space. Our study is, however, different in several aspects: another expression is used for the macroscopic energy, the evaluation of pairing and shell-correction energy is different, and, most importantly, a better more general nuclear shape parametrization is used in our present approach. In addition, the range of nuclear quadrupole deformations in Ref. [16] is smaller than the one in our present study, with the consequence that the super-deformed minima found by us are not present in the energy landscapes presented in Ref. [16].
The present investigation is also similar to those described in our previous study on shape isomers in the same mass region [17], with, however, several differences: • the present calculations are carried out for a larger set of isotopes ( 170−204 Pt, 172−206 Hg, and 174−208 Pb), • the potential energy surfaces (PES) are studied here in two 4-dimensional deformation spaces: (q 2 , q 3 , q 4 , q 6 ) and (q 2 , q 3 , q 4 , η), where η describes a non-axial degree of freedom, and the q i parameters correspond roughly to multipole deformations, • a denser grid of deformation parameters than that of Ref. [17] is used, • an average pairing energy taking into account the effect of an approximate particle-number projection is included.

Theoretical model and details of the calculations
In the macroscopic-microscopic method [18] the total energy of a nucleus at a given deformation can be calculated as the sum of the macroscopic (liquid-drop type) energy and quantum correction terms for protons and neutrons generated by shell and pairing effects: The LSD model [13] is used to evaluate the macroscopic part of the energy. Shell corrections are obtained, as usual, by subtracting the average energy E from the sum of the single-particle (s.p.) energies of occupied orbitals Here the s.p. energies e k are simply the eigenvalues of a mean-field Hamiltonian with a Yukawa-folded s.p. potential at the desired deformation. The average energy E is evaluated using the  with an 6th order correction polynomial. The pairing correction is determined as the difference between the BCS energy [23] and the singleparticle energy sum and the average pairing energy [22]: In the BCS approximation the ground-state energy of a system with an even number of particles and a monopole pairing force is given by where the sums run over the pairs of s.p. states contained in the pairing window defined below. The coefficients v k and is the pairing-energy correction due to the particle number projection done in the GCM+GOA approximation [24] where E k = (e k − λ) 2 + Δ 2 are the quasi-particle energies, and Δ and λ denote the pairing gap and the Fermi energy, respectively. The average projected pairing energy, for a pairing window of width 2Ω, which is symmetric in energy with respect to the Fermi energy, is equal to Hereg is the average single-particle level density andΔ the average paring gap corresponding to a pairing strength G The width of the pairing window for protons or neutrons is chosen so as to contain 2 √ 15N (N = N or Z) s.p. states around the Fermi surface. For such a window the pairing strength can be approximated by [25] where A is the mass number, and we have used the same value g p 0 = g n 0 = 0.28hω 0 for protons and neutrons with hω 0 = 41/A 1/3 MeV.
In the whole calculation the single-particle spectra are obtained by diagonalization of a Hamiltonian with a Yukawafolded mean-field potential [14,15] having the same parameters as those used in Ref. [26].
The macroscopic-microscopic energy landscape of each nucleus is evaluated as function of the different deformation degrees of freedom by using the Fourier parametrization of the nuclear shape that we recently developed [11]: where z 0 is the half-length of the shape and z sh locates the centre of mass of the nucleus at the origin of the coordinate system. It turns out that one can define what we will call "optimal coordinates", q n , through [12] in such a way that the liquid-drop energy as function of the elongation q 2 becomes minimal along a track that defines the liquid-drop path to fission. The a (0) 2n in Eq. (10) are the expansion coefficients of a spherical shape given by a (0) 2n = (−1) n−1 32 π 3 (2n−1) 3 . Non-axial shapes can easily be obtained assuming that, for a given value of the z-coordinate, the surface cross-section has the form of an ellipse with halfaxes a(z) and b(z) [12]: where the parameter η describes the non-axial deformation of the nuclear shapes. The volume conservation condition requires that ρ 2 s (z) = a(z)b(z). A similar kind of investigation has been undertaken in Ref. [16] for a much wider range of nuclei (7206 isotopes with masses between A = 31 and A = 290) but for a much more restricted deformation space (3 deformation parameters 2 , 4 and γ corresponding respectively to quadrupole, hexadecapole and axial asymmetry), whereas our deformation space contains the q n deformation parameters of Eq. 10) plus one axial asymmetry parameter η defined in Eq. (11). The study of Ref. [16] was also carried out using the macroscopicmicroscopic method, but with the finite-range liquid droplet model (FRLD), while in the present study we use the Lublin-Strasbourg Drop model (LSD) [13]. As it was shown in Ref. [13], the use of deformations of higher multipolarity, that are different in both models, induces a different stiffness of the PES, what has a non negligible effect on the energy landscape. At difference from Ref. [16] where a pairing treatment in the BCS approach was used together with the Lipkin-Nogami prescription to take into account the effect of the projection on the correct particle number, we use an approximate GCM+GOA particle-number projection in our approach [24,25], which yields results that are much closer to the exact ones, in particular in the weak-pairing limit.
The above description of non-axial shapes of deformed nuclei (q 2 , η) is more general than the commonly used (β, γ ) parametrization by Bohr [27,28]. For spheroidal shapes both descriptions are, however, equivalent. As one can see from Fig. 1, where the two parametrizations are compared, the periodicity of the nuclear shapes by an angle of 60 • is similar in the (q 2 , η) and in the (β, γ ) planes. One has, however, to keep in mind that this regularity is partially spoiled when higher multipolarity deformations q n (n > 2) come into play, which makes that our (η, q 2 , q 3 , q 4 , q 6 ) shape parametrization is substantially more general than the 3-dimensional (ε 2 , ε 4 , γ ) parametrization used in Ref. [16]. It is only for the very special case of spheroidal shapes that both parametrization coincide.
It is important, in this context, to stress that there are different points in the (β, γ ), as in the (q 2 , η) plane, that correspond, when higher q n , n > 2 degrees of freedom are neglected (and only under that assumption!), to precisely the same shape, with the only difference that the axes of the coordinate system have been interchanged (like y ↔ z).  When investigating potential energy landscapes, where the triaxial degree of freedom is taken into account, one therefore should be extremely careful, not to consider as two different configurations points in the (q 2 , η) deformation plane that are nothing but γ = 60 • rotation images of one another.

Results
The calculations of the potential-energy surfaces were performed for the isotopic chains of even-even Pt, Hg and Pb nuclei, with the same range of neutron numbers between N = 92 and N = 126. The macroscopic-microscopic energy of each isotope was evaluated in a four-dimensional (4D) deformation space spanned by the q 2 , q 3 , q 4 , and η coordinates [see Eqs. (10) and (11)]. In parallel, a calculation made in the 3D space of the q 2 , q 4 , and q 6 deformation parameters was carried out in order to test the effect of higher-order deformations on the PES in the considered nuclei.
Our calculations in the 4D deformation space spanned by the q 2 , q 3 , q 4 , and η coordinates, clearly show that the left-right asymmetry degree of freedom does not decrease the potential energy at any of the studied deformations, in any of the here considered nuclei, as one can see in the (q 2 , η), (q 2 , q 3 ) and (q 2 , q 4 ) PES presented in Fig. 2 for the 182 Hg nucleus, chosen here as a representative case for all the nuclei in the present study. The PES shown are generally defined as the difference of the total nuclear energy (1) of the deformed nucleus and the one of the corresponding spherical shape evaluated using the LSD macroscopic model. All energy points in the (q 2 , η) map (top of Fig. 2) Ref. [28]). In addition, red lines representing β = 0.1, 0.2, . . . values are drawn. The part of the (q 2 , η) map between the γ = 30 o and γ = 150 o lines is usually omitted when one projects a multidimensional macro-micro or self-consistent PES onto the (β, γ ) plane, loosing hereby, as we have explained above, some information on the impact of higher multipolarities. We have therefore decided to show in the following the full (q 2 , η) maps in order not to lose part of the results, which could be important when e.g. the deepest minimum would appear in this upper part of the map.
In the middle panel of Fig. 2 is shown the PES in the (q 2 , q 3 ) plane which has been obtained for the axial sym- Different line types indicate minimization of the energy with respect to different shape parameters: q 4 and q 6 set to zero (blue dotted line), minimization with respect to q 4 , but q 6 set to zero (red dashed line), minimization with respect to q 4 and q 6 (black solid line) metric case (η = 0) after minimization with respect to the q 4 degree of freedom. In this and in the rest of the energy landscapes shown in the following, solid lines correspond to layers separated by 1 MeV, while the distance between dashed and neighbouring solid lines is of 0.5 MeV. The labels on the layers denote energies in MeV.
The PES displayed in the bottom panel of Fig. 2 shows the energy landscape in the (q 2 , q 4 ) plane. The structure of the PES is seen to be rather complex and quite a few local minima can be observed in this cross-section of our 5D deformation space, two on the oblate (q 2 ≤ 0) and three on the prolate side. One has to bear in mind, however, that one is looking here only at a 2D cross-section of a 5D deformation space and that the stability of these local minima with respect to the other deformation degrees of freedom, especially the non-axial η, should be carefully studied in each case. Also the role of the higher multipole deformation q 6 needs to be analysed. We have therefore tried to investigate the impact of these higher-multipole degrees of freedom on the deformation energy of 182 Pt and 182 Hg. The results of this study is shown in Fig. 3, where the energy as function of the q 2 elongation parameter is shown in a 1D cross-section of the PES and where different lines indicate the importance of different higher-order deformations. These deformation energies which are typical for all the considered nuclei in the three here investigated isotopic chains, clearly show that the role of q 6 deformations is practically negligible, thus demonstrating that the "optimal coordinates" q n , defined in Eq. (10), have, indeed, been skilfully chosen.
An overview of the deformation energies as function of the elongation parameter q 2 is shown in Fig. 4 for the nuclei of the three Pt, Hg, and Pb isotopic chains with neutron numbers 92 ≤ N ≤ 126, where the potential energy of each isotope is minimized with respect to q 4 and q 6 . The left-right asymmetry degree of freedom q 3 is not included in this study, since for all the here considered nuclei it does not play any noticeable role as demonstrated above. Lines corresponding to different isotopes are drawn in different colours and the ground-state energy minimum of each isotope is shifted by 1 MeV with respect to the previous isotope in order to separate the curves. In each diagram the bottom and top lines are marked by the corresponding mass number A for better visibility. The experimental excitation energies of the superdeformed (SD) minima are also shown when available, as in the case of Hg and Pb nuclei.
In the upper parts of each diagram in Fig. 4 the corresponding values of the quadrupole deformation β 2 are also given. This deformation parameter is defined by the harmonic expansion of the radius of a deformed nucleus: where Y λμ (θ, ϕ) are the spherical harmonics. Such an expansion is frequently used in publications dealing with PES calculated in the mac-mic approach (see e.g. Ref. [29]) but turns out to be very slowly converging at large nuclear elongations (as shown in Refs. [11,12,30]). The experimental energies corresponding to the bottom of the SD bands, measured with respect the ground state, are marked by the red points in Fig. 4. These experimental data are taken from Ref. ([31], Fig. 46) and found in all cases to be in quite good agreement with our theoretical predictions. Quite pronounced oblate and prolate minima around q 2 ≈ 0.2 and q 2 ≈ −0.2 are found in all Pt isotopes, except that one finds, when looking just at such a kind of 1D deformation energy, that beyond 194 Pt (with neutron number larger than N ≈ 116) the prolate minimum at q 2 ≈ 0.2 gradually looses its importance at the profit of a super-deformed (SD) shape isomer located around q 2 ≈ 0.6. One has, however, to keep in mind that we have, for the moment, completely left out of our considerations the non-axial degree of freedom. As we will see in Fig. 7 in Sect. 3.1, that analysis will change, when we will take that additional degree of freedom tot ≈ 3.5 MeV is predicted by our calculation to occur in 178−184 Pt. A delicate structure appears on the prolate side in 172−178 Pt isotopes when two minima around q 2 = 0.25 and q 2 ≈ 0.35 begin to compete. At this point one cannot really judge about the stability of these local minima and an eventual shape coexistence by looking at such a 1D plot only. A more extended study which takes also the non-axial degree of freedom into account is necessary and will be carried out in Sect. 3.1.
Apart from the minima at deformations typical for the ground state, strongly oblate minima at q 2 ≈ −0.4 are found in 170−176 Pt isotopes, minima that are separated from the normal-deformed oblate one by a small barrier of only ≈ 0.5 MeV. In addition, a rather strongly deformed pronounced prolate minimum at q 2 ≈ 0.6 is observed both in the lightest here studied Pt isotope 170 Pt, as well as in Pt isotopes with A ≥ 186, with a barrier which separates them from the normal-deformed prolate minimum and thus guarantees their stability, a barrier that gradually grows with increasing neutron number to reach about 3 MeV in 204 Pt. In most cases these SD isomeric states are, however, rather high in energy above the ground state (E isomer tot − E g.s. tot ≈ 6 MeV), and will therefore be difficult to be observed experimentally. On the other hand, the moment of inertia of nuclei in such SD bands is several times larger as compared to a typical deformed ground-state band [32], which has a significant influence on the high-spin level structure since the SD states can appear at energies close to or even lower than the rotational states built on the ground-state.
The situation in the Hg isotopes (middle row of Fig. 4) is even more complex. Apart from the strongly deformed oblate minimum at q 2 ≈ −0.45 observed in the lighter isotopes 172−184 Hg, separated from the ground state minimum by a barrier of up to 1 MeV, a very elongated shape isomer is present at q 2 ≈ 0.6 in 172−176 Hg, disappears beyond 176 Hg, and reappears again in the 190−206 Hg isotopes. This strongly deformed shape isomer is protected from a decay towards to the normally deformed states by a fairly high barrier and is located between 4 MeV (in 190 Hg) and 16 MeV (in 206 Hg) above the ground state. In 176−186 Hg additional low-lying prolate or oblate minima appear, thus giving the chance to observe a shape-coexistence in these nuclei. One has of course to study the stability of these local minima with   [33].Thin dotted vertical line corresponds to the predicted equilibrium deformation of the SD isomer respect to non-axial deformations, an investigation that will be carried out in Sect. 3.2.
The ground state of all considered 174−208 Pb isotopes is found to be spherically symmetric, but in all of them deformed local minima can be observed. Similar as in the case of the Mercury isotopes, strongly deformed oblate minima with q 2 ≈ −0.45 appear in the isotopes with 174 ≤ A ≤ 186, located again rather high above the ground state, at about 5 MeV in 174 Pb and 2.5 MeV in 186 Pb. Very deformed prolate configurations with a main-axis ratio close to 2 (q 2 ≈ 0.6) are predicted by our calculations in the five lightest isotopes 174−182 Pb and in isotopes with N 104. It is interesting to notice that some additional local minima appear at smaller deformations q 2 ≈ −0.2 and 0.2 q 2 0.5 in 182−188 Pb isotopes, where even several oblate and prolate deformed minima are predicted. These results, already reported in our previous study [17], are in line with the experimental observations (see e.g. Ref. [1]). The stability of these minima with respect to the non-axial degree of freedom will be discussed in the subsequent Sect. 3.3.
In the three lightest Pb isotopes (A = 174-178) one observes a sudden wobbling of the energy at q 2 ≈ 0.5 and the appearance of a local SD minimum at q 2 ≈ 0.6. These are induced by the minimization process, when i.e. the energy minimum jumps with increasing elongation q 2 from one to another valley in the (q 4 , q 6 ) plane.
One could ask the question about the accuracy of our theoretical predictions. A comparison with available data, taken from Refs. [31,33,34], is presented in Figs. 5 and 6, where Fig. 5 shows the moment of inertia J x and the charge quadrupole moment Q 20 as function of the elongation param- Fig. 6 Predicted electric quadrupole moments of Pt, Hg, and Pb isotopes in the ground-state (black points) and in the super-deformed isomeric states (blue squares) are compared with the experimental data (stars) taken from Refs. [31,33,34] eter for 194 Pb, and Fig. 6 the charge quadrupole moment in the ground state and the SD shape isomeric state for the isotopic chains of Pt, Hg and Pb. The experimental value J exp x of the moment of inertia of the SD shape isomer in Fig. 5 is evaluated from the energies of two lowest experimentally observed members (4 + and 6 + ) of the SD rotational band [33]. One can see that the experimental values (J exp x and Q exp 20 ), given by the horizontal lines in Fig. 5 coincide quite well, at the predicted deformation of the SD minimum (vertical dotted line), with the theoretical values for both these quantities. Another measure of the accuracy of our predictions is presented in Fig. 6, where the theoretical absolute values of the charge quadrupole moment for the ground state (black points) and the SD minimum (blue squares) of Pt, Hg, and Pb are compared with the corresponding experimental data (red stars) taken from Refs. [31,33,34].
We would like to stress at this point that all results, obtained in models taking only axial symmetric shapes into account and giving an indication of a possible existence of shape isomers or shape coexistence, cannot seriously be considered as reliable without checking the stability of the identified minima with respect to other degrees of freedom, among which the non-axial deformation seems to be the most important one. That is precisely what will be investigated in the subsequent subsections.

Pt nuclei
The potential energy surfaces of 170−204 Pt isotopes on the (q 2 , η) plane are displayed in Fig. 7. The energy in each of the deformation points on these maps is minimized with respect to q 4 deformations only, since we have seen, as explained above, that q 3 and q 6 degrees of freedom are practically ineffective in these nuclei at the investigated deformations. A very elongated local minimum appears here at q 2 ≈ 0.6 in the two lightest 170−172 Pt isotopes, when a non-axial deformation is taken into account. It turns out that this local minimum is deeper here than in Fig. 4 where that degree of freedom is absent. Both these isotopes also appear to be unstable in their ground state (|q 2 | ≈ 0.2) with respect to the non-axial degree of freedom. They can therefore be considered as good examples of the so-called γ -instability. Please recall in this respect the discussion related to the two different descriptions, (q 2 , η) and (β, γ ), of the non-axial degree of freedom at the end of Sect. 2 and visualized in Fig. 1.
When looking at the deformation energies of the 174−180 Pt nuclei in Fig. 4, where axial symmetry had been imposed, one has the impression that there is a coexistence of a prolate minimum at q 2 ≈ 0.28 with an oblate minimum at q 2 ≈ −0.20 separated by a barrier of about 1.5-2.5 MeV. When, however, taking the non-axial degree of freedom into account, as this has been done in the calculations represented in Fig. 7, the picture changes. When looking e.g. at the case of 178 Pt, one again finds a coexistence of two minima, one prolate at (q 2 = 0.27, η = 0), the other triaxial at (q 2 = −0.15, η = 0.09), which are separated by a much lower barrier now, that now the non-axial degree of freedom is taken into account Remembering however the discussion at the end of Section 2, one concludes that these two local minima just represent the same prolate shape, since the triaxial minimum with negative q 2 value is nothing but the symmetry image of the axially symmetric prolate ground-state minimum, as can easily be verified by taking into account the 60 • γ symmetry represented in Fig. 1.
The same situation is observed in all of the 174−180 Pt isotopes. Let us recall at this point that the above sextant symmetry, which is exact in the case of spheroid shapes, is broken when higher multipolarity deformations are taken into account. In Ref. [16] the special choice of the axial and non-axial hexadecapole deformations, originally proposed by Rohoziński and Sobiczewski [35], was used to keep this 60 • regularity. The shape parametrization used in the present paper does not posses exactly this symmetry. One could in principle introduce similar constraints as in Ref. [35] on the q i deformation parameters to restore such a sextant symmetry, but such constraints would cause a loss of generality of our deformation space which we want to avoid. Oblate-prolate shape coexistence is calculated in the 182−186 Pt nuclei, but the barrier between the two minima is small. Also in 188 Pt two minima are observed, one oblate at (q 2 ≈ −0.2, η = 0) and another one triaxial at (q 2 ≈ 0.2, η≈ 0.02). In the heavier Pt isotopes with 190 ≤ A ≤ 200, oblate shapes are preferred in the ground state, whereas, approaching the N = 126 shell closure, the 202−204 Pt nuclei turn out to be spherical. The very deformed shape isomers around q 2 ≈ 0.6 observed for the 188−204 Pt isotopes in Fig. 4 persist when non-axial deformations are taken into account, with a slight reduction (of up to 0.5 MeV) of the barrier height between these isomers and the ground state for the four heaviest isotopes.

Hg nuclei
Just in the same way as done in Fig. 7 for the Pt isotopes, we show in Fig. 8 the deformation energies of the Mercury isotopes 172−206 Hg on the (q 2 , η) plane. One observes that the shape isomeric states corresponding to both the very oblate deformations (q 2 ≈−0.45) in 172−182 Hg, and the very prolate deformations (q 2 ≈ 0.6) in 172−176 Hg and in 190−206 Hg all survive when the non-axial degree of freedom is taken into account with the exception of the prolate isomer in 174 Hg and 176 Hg which becomes slightly triaxial with η ≈ 0.03. A significant reduction of the height of the barrier between the shape-isomeric states with q 2 ≈ 0.6 and the ground state minimum can, however, be observed in the 200−206 Hg isotopes due to the non-axial degree of freedom. The ground state of 172 Hg and of 200−206 Hg are spherically symmetric, while in the 174−176 Hg isotopes it is the prolate minimum that corresponds to the ground state. The additional triaxial minima observed at (q 2 ≈ −0.05, η ≈ 0.03) correspond again to the 60 • mirror image of that ground-state deformation, as explained at the end of Section 2. For the 178−188 Hg isotopes shape-coexistence is calculated at normal deformations, while for isotopes between 190 Hg and 198 Hg the ground-state minimum corresponds to an oblate shape.

Pb nuclei
The 174−208 Pb isotopes are analysed in the same way in Fig. 9 through PES landscapes in the (q 2 , η) plane. The situation turns out, however, to be slightly different here, since the ground state is found to be spherically symmetric in all the considered isotopes and no competition between prolate and oblate ground-state deformation is present. The issue of a possible shape-coexistence can, in fact, not be limited to prolate versus oblate shapes, but is much more general, as can be seen through the (q 2 , η) maps at small deformations in the 178−186 Pb isotopes, where the landscapes in η directions are very soft, thus allowing for local minima, particularly visible in the maps of 182−184 Pb. Similarly as in the case of the Pt and Hg isotopes, a well deformed prolate shape isomer at q 2 ≈ 0.6 is predicted in practically all Pb isotopes except 182−188 Pb. This shape isomer is particularly visible in  Fig. 9 The same as in Fig. 7 but for the Pb isotopes the heavier isotopes, as already noticed on the lower part of Fig. 4, even though we now realize that the barrier between the ground state and that prolate shape isomeric state is lowered by as much as 1 MeV by the inclusion of deformations that break the axial symmetry, in particular in the heavier isotopes 204−208 Pb. Rather shallow oblate minima appear at an elongation of q 2 ≈ −0.45 in the 176−186 Pb isotopes, with a small barrier of about 0.5 MeV that separates them from the ground-state minimum. Since these minima turn out to be axially symmetric they already appear clearly in Fig. 4. An additional axially symmetric prolate minimum is predicted at q 2 ≈ 0.35 by our calculations in 182−188 Pb. It is located at about 1 MeV above the ground state well, separated from it by a small barrier of about 0.5 MeV. Similarly, small oblate minima are visible in Fig. 4 for elongations in the range −0.3 q 2 −0.1 in 184−194 Pb, as was found previously in some mainly self-consistent approaches (see e.g. Refs. [6][7][8][9]). These turn out, however, to disappear when the non-axial degree of freedom is taken into account, as this is clearly seen in Fig. 9.
Let us refer the interested reader to a recent publication [36] of our group, where shape isomers in nuclei of the three here studied isotopic chains have been extensively analysed and listed in a table in the appendix together with their energies with respect the ground-state, axial Q 20 and non-axial Q 22 quadrupole moments, and moments of inertia.

Summary and conclusions
In our study of isotopic chains of neutron deficient Pt, Hg and Pb nuclei with neutron numbers in the range 92 ≤ N ≤ 126, a certain number of general remarks can be made.
First of all we have found that the left-right asymmetry degree of freedom does not play any important role in any of the here considered nuclei, and this not only around the ground state deformation, but throughout the here considered range of deformations. However, this conclusion does not necessarily hold when going to very large elongations (beyond q 2 ≈ 1) as encountered in the fission process.
Another important conclusion is that it appears absolutely essential to take into account deformations that break the axial symmetry. It indeed appears in all three of the here investigated isotopic chains that, at least for small neutron numbers 92 ≤ N ≤ 106, i.e. away from the N = 126 shell closure, the energy landscape in the (q 2 , η) plane is rather flat around the ground-state deformation, which always stays rather close to the spherical configuration, while a very pronounced spherical ground-state minimum gradually develops when approaching the magic number N = 126. For the more neutron deficient nuclei that present a very rich (q 2 , η) landscape in the vicinity of the spherical configuration, small differences can decide whether the ground-state deformation will turn out to be prolate or oblate, as this is demonstrated when going from 186 Pt, which is prolate deformed, to the neighbouring 188 Hg nucleus, by just adding a pair of protons, which then turns out to have an oblate ground-state deformation. A similar, even more pronounced effect was already observed in in the Polonium isotopes [17], on the other side of the Z = 82 shell closure, when going from 182 Po with a prolate ground-state at q 2 ≈ 0.4 with an ≈ 1.5 MeV higher oblate isomer at q 2 ≈−0. 25, to 192 Po where the oblate minimum becomes the ground-state, before both these merge at 198 Po with a ground-state which is spherically symmetric.
As far as the shape coexistence between prolate and oblate minima is concerned, several nuclei in all three of the here studied isomeric chains, clearly emerge: it is often away from the magic numbers, i.e. in regions where the shell corrections are not that dominant, that such a phenomenon comes into play. In the Platinum nuclei, we have found that the region 182−188 Pt exhibit minima at q 2 ≈ ±0.2 separated by a small barrier, while in the Mercury isotopic chain it is in the 178−188 Hg isotopes the shape coexistence is present at small deformation. But the question of a possible shapecoexistence can, in fact, not be limited to prolate versus oblate shapes, being more general, as this has been illustrated through the (q 2 , η) maps in the Lead isotopes, in particular in 178−186 Pb, where the landscapes in η directions are very soft, thus allowing for local minima, particularly visible in the landscapes of 182−184 Pb.
Future investigations of some additional isotopic chains in the Z ≈ 82 mass region, but also in other quite different mass regions are anticipated, in order to probe, from a comparison with the experimental data, as that has been done here, the predictive power of our theoretical approach.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: No further data are associated with this paper.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.