Demonstration of the Universality of Molecular Structures in Prolate Deformed Nuclei

The relationship between the deformed harmonic oscillator and the formation of molecular cluster structures, whereby valence neutrons are exchanged between cluster cores, is examined. It is found that there is a strong connection between the properties of the valence orbitals associated with deformed structures in the deformed harmonic oscillator and the molecular orbitals created by linear combinations of single-centre orbitals around nuclear clusters. The conclusion is that in addition to the appearance of clustering in the deformed harmonic oscillator that \emph{every} prolate deformed cluster structure should have molecular orbitals built on that structure. This is demonstrated through a series of examples that range from $^{13}$C to $^{57}$Ni.


Introduction
The appearance of clusters in light nuclei is well documented [1,2]. The most prominent include the 8 Be ground-state and the 12 C excited Hoyle-state. Both have α-particle cluster structures and both play a key role in the synthesis of carbon in stars, through the triple-alpha process. There is good evidence for more complex cluster structures, where, for example, 20 Ne may be described in in terms of 16 O+α substructures and there is a long history of clustering in 24 Mg, where α+ 16 O+α clustering is one cluster configuration [1].
The formation of clusters presents an interesting possibility, which is well illustrated by the Beryllium isotopes. The 8 Be nucleus is unbound to decay to two α-particles in its ground-state, whereas 9 Be is stable. The additional neutron has been found to exist in a delocalised orbit, covalently exchanged between the two α-particle clusters. The existence of such structures was discussed by von Oertzen in a series of seminal papers [3,4]. These molecular orbitals in the Beryllium isotopes have both π and σ characteristics, where the notation follows from the atomic analogues. The experimental evidence for such structures comes from the associated rotational bands and decay properties of the states. Perhaps the best example is that of 10 Be, which provides evidence for a molecular structure of two α-particles and 2 σ covalent neutrons [5,6].
The extension of these ideas to three centres and carbon isotopes was made by Milin and von Oertzen [7], where the states in 13 C were analysed in terms of their potential molecular characteristics. Convincing evidence for the molecular structure of the carbon isotopes remains to be demonstrated, though there is significant contemporary interest in the molecular structure of 14 C [8]. Finally, the ideas of nuclear molecules were extended to asymmetric clusters by von Oertzen [9] and Kimura [10] through an analysis of the Neon isotopes.
These studies have provided a basis for the understanding of the link between the appearance of molecular orbitals and the underlying cluster structures. This was further cemented through the study of the relationship between the states of the deformed harmonic oscillator and molecular states created through the Hückle approach [11]. These molecular structures are found not only in such relatively schematic approaches, but also in more sophisticated methods that treat the nucleons individually, do not impose the clustering a priori and deal with realistic nucleon-nucleon interactions. There are many such approaches, see Ref. [2] for a review, but the Antisymmeterised Molecular Dynamics, AMD, method has been particularly successful in both reproducing the molecular structures and experimental data, providing a level of confidence that the molecular degrees of freedom are a good representation of the structure of a number of light nuclei.
The present work builds on Ref. [11], which examined the molecular states of α-like nuclei, and extends this understanding to the full spectrum of prolate deformed cluster states across the entire mass range of possible cluster structures. Figure 1. Energy levels of the deformed axially symmetric harmonic oscillator as a function of the quadrupole deformation (oblate and prolate, i.e. negative and positive values of 2 , respectively). Degeneracies (labelled) appear due to crossings of orbits at certain ratios of the length of the long axis (the symmetry axis) to the short perpendicular axis. The regions of high degeneracy define a shell closure for deformed shapes [12].

Theoretical Approach
The aim of the present paper is to explore the nature of molecular wave-functions across a wide range of nuclear systems. The framework for this is the deformed harmonic oscillator, DHO. The DHO solutions for molecular-like nuclei are compared with results of calculations where the molecular structures are created explicitly from the orbits of linear combinations of two-centre orbits created using the Hückle method.

Deformed Harmonic Oscillator
The DHO provides a good description of light nuclear systems, where the inclusion of the spin-orbit interaction and an l 2 dependence of the nuclear potential results in the Nilsson description of deformed nuclei [12]. The deformed harmonic oscillator has the well-known shell-like structure where at integer ratios of axial deformations the shell structure is enhanced, Fig. 1. This shell structure may be associated with stable deformed structures and indeed cluster partitions [13]. The total energy at each level for the DHO is given by: where ω x , ω y and ω z and n x , n y and n z are the oscillation frequencies and number of oscillation quanta in the x, y, and z directions respectively, and ω 0 is the average oscillation frequency. Furthermore, the oscillation frequencies must satisfy: Therefore, for a nucleus with, for example, a 2:1 deformation along the z-direction, the energy levels are given by: The orbitals in the DHO are labelled by [n x , n y , n z ], and the number of oscillator quanta in each Cartesian direction determines the orientation of that orbital. The association of cluster configurations with the shell structure [12,13] follows from the degeneracies illustrated in Fig. 1. Within the DHO the nucleus 8 Be, for example, is created through the population of the DHO orbitals, [n x , n y , n z ]=[0,0,0] 4 and [0,0,1] 4 . At a 2:1 deformation this is seen to be associated with the filling of the two orbitals with 2p+2n, i.e. 2α-particles, and hence 8 Be may be thought of as 2 α-clusters. The densities calculated from this DHO configuration also reveal the two-centre cluster structure [12]. For the next most complex case which is 20 Ne, this also corresponds to a 2:1 shell gap. Here the degeneracies are 2+(2+6) and may be linked to a 4 He+ 16 O cluster structure. The next shell gap at a deformation of 2:1 occurs for degeneracies 2+2+8+8 which would correspond to a 2:1 deformation of 32 S and a 16 O+ 16 O cluster structure. A similar situation occurs at a deformation of 3:1, for example the shell gap at cumulative degeneracies of 2+2+2+6+6+6 being linked to the 3:1 deformed state in 48 Cr and a 16 O+ 16 O+ 16 O cluster configuration.
To form molecular structures in these cluster-systems where the valence particle, typically a neutron, is exchanged between the cluster cores, one needs to examine the properties of the next available orbits beyond those associated with the α and 16 O cluster cores. These are then those which exist immediately above the shell gap. In the case of 8 Be+n then these are the three levels . In order to ascertain if these have a molecular character, their properties need to be compared with those created explicitly from the available two-centre orbitals. The present work has been framed in terms of neutron covalent orbitals as these are the most prevalent in terms of experimental counterparts. However, equally the ideas may be applied to valence protons or combinations of protons and neutrons.

Neutron Wave-functions in Two-Centred Systems
The neutron wave-functions may be created from linear combinations of the singlecentre obits of the valence neutron around the cluster core. In the case of the α-particle Hence, the nuclear molecular orbit (ψ) is then given by a sum of the neutron wavefunctions (φ) at each cluster centre, with the appropriate normalisation constants, c i , (determined by the number of cluster centres in the molecule): To find the relative contribution which each single-centre wave-function makes to the overall molecular one, the Hückel method was followed. The appropriate coefficients are found by solving the secular determinant equation which arises from the Schrödinger equation for the nuclear molecule (Section 2.3).
When making linear combinations of the orbitals in the multiple centre model, the quantum numbers defining the orbital obey certain rules which allowed the combination to be associated with one of the single-centre orbitals [14,11]: • The n z quantum numbers of the single centre orbits should be summed over the number of centres in the molecule, C, i.e. n z = C i n zi , to form the lowest energy molecular configuration.
• The n x and n y quantum numbers, perpendicular to the direction of deformation, are conserved; the value for n x and n y remains unchanged in the composite system for symmetric clusters (n x1 = n x2 = n x and n y1 = n y2 = n y ).
• When forming the orbitals for the excited molecular states, the quantum number along the direction of deformation can be found using: n z = C i n zi + C − 1 where n z is the quantum number for the orbital in the molecule, n zi are the quantum numbers at each centre in the multiple centre molecule.

The Hückel Method
The wave-functions describing all quantum particles can be found by solving the Schrödinger equation, for the molecular system whereĤ is the Hamiltonian of the nucleus, E is the energy eigenvalue associated with the particular molecular state, and ψ is the overall wave-function of the nuclear molecule as given above. To demonstrate the method the case of 33 S will be considered. Here there are two 16 O clusters and a single valence neutron. The wave-function describing the valence neutron 33 S is then: Substituting the molecular wave-function of 33 S into its eigenvalue equation yields: Where the equation has been simplified using H ij = ψ * iĤ ψ j dτ and S ij = ψ * i ψ j dτ . The actual wave-function of the nuclear molecule, given by the values of c 1 and c 2 will yield the lowest energy state (rather than an excited state); therefore the coefficients can be found by minimising the energy eigenvalue with respect to c 1 and c 2 . The following secular equations are obtained by finding the partial derivatives and setting them equal to zero: These can be equivalently written in matrix form: The energy eigenvalue, E, which satisfies this equation is found by solving the secular determinant, because for a non-trivial solution the matrix must be singular: The secular determinant equation may be simplified by making assumptions about the interactions within the nuclear molecule. H 11 and H 22 represent the self-interaction energies of the free 17 O nuclei at each centre (i.e. the n+ 16 O binding energy), which are identical due to the symmetry of the molecule, this value is set equal to α. The interaction energy between the two clusters is given by H 12 = H 21 = β if the 16 O clusters are neighbouring (as is the case for 32 S) or set to zero otherwise. Additionally, due to the neutron wave-functions φ 1 and φ 2 being located at different cluster centres, their separation enables the simplifying assumption of small overlap to be made, so S 12 = S 21 0. Finally, if the neutron wave-functions at each cluster centre have been normalised then S 11 = S 22 = 1. Using all of these assumptions gives the secular determinant equation: The solution to this equation gives E = α ± β, and the solutions for which are c 1 = c 2 and c 1 = − c 2 . As the molecular wave-function must be normalised, for 33 S the solutions are: The same method for calculating the molecular wave-function for nuclei comprising more complex arrangements of clusters, i.e. for 3 or more centre systems. The wave-function for each neutron orbital can be constructed in the same way (by solving the secular determinant equation), and linear combinations were made by using the appropriate coefficients. The coefficients for collinear molecules (i.e. when the clusters form a 1d linear structure) of up to six clusters are given in Table 1 [11]. For the two-centre system molecular wave-function is given by: Here, the HO scale factor mω/h is set to unity and n denotes the number of oscillator quanta, r 1 and r 2 give the locations of the two clusters with respect to the origin and H n is the appropriate Hermite polynomial. For this investigation, the wave-functions were projected onto the y − z plane and integrated up over the x dimension, and the z-direction was chosen to be the deformation axis, or the axis of separation of the clusters. Table 1. The coefficients of the single particle wave-functions in the molecular orbits of prolate nuclei. The magnitude of the coefficients determine the relative contribution of the wave-function to the particular molecular orbit, taken from [11]. The increasing orbit number corresponds to higher number of nodes in the composite wave-function and hence higher energy levels.

Results
Figures 2 to 5 show the densities of the cluster cores built from both the DHO (singlecentred nuclear core) and also by explicitly creating a superposition of the densities from the cluster components in a multi-centre approach. These cores correspond to the underlying structures around which the molecular neutrons would orbit. What is clear is that the symmetries that are found in the case of the degeneracies of the DHO (Fig. 1) are broadly reproduced in terms of the densities, revealed in the multi-centre calculations. The multi-centre calculations were performed for separations of the cluster components which are approximately equivalent to the peaks in the density distributions indicated in the DHO densities. The oscillator scale parameter was adjusted from unity for the individual clusters according to the relationship mω/h = 0.9A 1/3 + 0.7 fm 2 . It can be observed that where the densities of the clusters in the multicentre calculations overlap the densities increase and the details of the multicentre structure are diminished. This can be seen, for example, in the case of the 4α structure of 16 O, in Fig. 2, where the 4α clusters are most evident in the low density region. However, in the deformed harmonic oscillator the cluster structure is often more evident, demonstrating the importance of the Pauli Exclusion Principle, PEP, in enhancing the cluster structure. Here the  overlapping densities from the separated cluster cores does not obey the PEP. This is explored further in the latter part of this paper. Fig. 2 shows the two sets of calculations for the three isotopes 12 C, 16 O and 20 Ne, which demonstrate the expected cluster structure associated with the 3:1, 4:1 and 2:1 shell closures, respectively. In the case of 20 Ne the multi-centred calculations are for 16 O+α which is reflection antisymmetric. Similarly Fig. 3 shows the α-16 O-α, 16 Fig. 5, demonstrating that the DHO produces the cluster symmetries. There is no significant reduction of the cluster structure in the deformed harmonic oscillator even in the heaviest systems. As noted, it is these clustered structures which then create the cores for the valence neutrons, covalently exchanged between the clusters.

Symmetric Two Centre Molecules
The 9 Be nucleus has been shown previously to contain two alpha clusters plus a covalent neutron [11]. The molecular characteristics are also found in the AMD calculations of the Beryllium isotopes, described at length in Ref. [2]. The next two-centre, symmetric, nuclear molecule is 33 S. 33 S can be considered as a 2:1 deformed nucleus with two 16  Given that the neutron could exist at either cluster centre with equal probability, the coefficients associated with the linear combination of the orbitals are 1/ √ 2 (Eqn. 15). Using the 'rules' described earlier, the resulting two-centre molecular wave-functions can Here the n z quantum number in the DHO is the sum of the two-centre n z quantum numbers and the n x and n y values remain unchanged. All of these linear combinations correspond to the addition (+) of the two-centre orbitals. The negative (−) possibility produces the DHO n z quantum number 2n z +1, i.e. an additional node in the wavefunction and hence higher energy orbitals, which can be matched to DHO counterparts in the next oscillator shell at this deformation; but not shown here. Fig. 6 shows the corresponding wave-functions. In each case there is a correspondence between the patterns of the wave-function in terms of shape and nodes, though of cause there is not a precise agreement. For example, in the case of the [0,0,4] DHO wave-function, there is good agreement with the molecular wave-function calculated from the linear combination of the two axially displaced [0,0,2] orbitals. Here the final n z quantum number, 4, is the sum of the two separated n z quantum numbers, 2+2. This agreement across all molecular wavefunctions demonstrates the molecular character of the valence orbitals available for the neutron in 32 S above the shell gap associated with the 16 O+ 16 O cluster structure.
The next symmetric two-centred molecule to be examined is 81 Zr, which can be thought of as two 40 Ca clusters with a single delocalised neutron. 40 Ca is the cluster formed in the spherical HO model when orbitals are filled to the next shell closure, after that of the 16  As the delocalised neutron in 81 Zr exists at a much higher energy than 33 S, the degeneracy is much larger. This means that there are a greater number of molecular orbitals which must be constructed using linear combinations, but there is still only one possible combination of wave-functions from each cluster centre which give the correct orbital shape, only a subset of these orbitals are shown in Fig. 7, but all combinations have been checked. Once again, it can be seen there is a good match between the properties of the 2:1 deformed orbitals and the explicit two-centre molecular counterparts. These two examples demonstrate that at a deformation of 2:1 the deformed harmonic oscillator orbitals above the deformed shell closure all have a molecular structure and thus the deformed harmonic oscillator explicitly includes the 2-centre molecular behaviour.

Three, and More, Centre Nuclear Molecules
In this section we will use selected examples to show that the behaviour found at 2:1 is reproduced at a deformation of 3:1 and that the inclusion of orbitals with a molecular character is more widely embedded within the deformed harmonic oscillator. The most simple nuclear molecule with three centres is the collinear arrangement of 12 C, comprising three alpha clusters. The nucleus can also be described with a 3:1 deformation in the DHO model, where the oscillation frequency along the z-direction is a third of that along the x or y directions. Following the approach in the previous section, the orbitals filled, by the core nucleons, when the nucleus is considered to exist at a single centre are When the 13 C nucleus is considered as three separate alpha particles, bound by a delocalised neutron, each alpha particle's nucleons occupy orbital [0,0,0]. So, the neutron could be found in a p-state above any of the three centres, and correspond to the degenerate [1,0,0], [0,1,0] and [0,0,1] spherical HO orbits. The coefficients found using the Hückel method, when linear combinations of these neutron wave-functions were made, the overall molecular wave-functions possess the form: The quantum numbers parallel to the deformation axis, n z , are summed from each centre to give the overall number; the quantum numbers perpendicular to the deformation direction, n x and n y , were conserved. The results are shown in Fig. 8.
The lowest energy configuration for the [0,0,3] orbital comes from adding the three oscillation quanta situated at each cluster centre in the z-direction. Here the phases of the three wave-functions are chosen to match in the overlap region maintaining the number of total nodes in the composite wave-function as 3×n z . There is another linear combination in which the central wave-function gives no contribution, Eqn. 19, which produces a wave-function which resembles that of the DHO [0,0,3] orbital. Here an additional node is created in the composite wave-function by the opposite phase of the two wave-functions, generating an additional node and an additional HO oscillator quantum; 2n z + 1. For the [0,1,0] orbital the combination of three oscillator quanta in the y-direction yields one overall oscillator quanta perpendicular to the deformation axis. This corresponds to the 3 single-centre wave-functions being added such that the phases combine constructively with no additional nodes in the z-direction.
The next three-centre (or 3:1 deformed) nucleus with 3 16   The structures at a deformation of 4:1 have also been explored in exactly the same way as at 2:1 and 3:1 and the same conclusions emerge; namely there is a 1 to 1 match between the nature of the wave-functions associated with the next available energy levels beyond closing the cluster cores and those predicted in the multi-centred approach.
In other words, for all structures with identical cluster cores the DHO contains the molecular symmetries for the valence neutrons.

Asymmetric Nuclear Molecules
Up to this stage, the nuclear molecules considered all have clusters of the same type, e.g. 32 S= 16 O+ 16 O and 80 Zr= 40 Ca+ 40 Ca clusters. Within the DHO energy level diagram (Fig. 1), this means that the nuclear core will fill orbitals at each centre with the same degeneracies, e.g. 2, 6, 12, ... , and same energies; the valence neutron must then exist in a shell closure which has a new value of degeneracy. However, nuclear molecules which have a different types of cluster may appear in many systems, e.g. 20 Ne≡ 16 O+α. As pointed out in [1], it is also possible to construct molecular orbitals from the valence neutron orbital around the two cores, e.g. from 5 He and 17 O. In this instance, the simplest nuclear molecule of this kind is 21 Ne.
For the nucleus 21 Ne at a deformation of 2:1, the orbitals that are available to the valence neutron are [0,0,3], [0,1,1] and [1,0,1]. In the multi-centre picture, there are different possible neutron orbitals which could be combined. At one centre the nuclear cluster is an alpha particle, and at the other centre the nuclear cluster is a 16  The linear combinations made to produce the single-centre wave-function follow the same rules in terms of quantum numbers as in the case of identical cluster cores. The quantum number parallel to the direction of deformation (z-axis) is found by summing the n z values of the orbitals at each centre, and the quantum numbers perpendicular to the deformation are conserved/unchanged. However, to create wave-functions that have counterparts in the DHO, the transverse quantum numbers, in the x and y directions must match for the multi-centred wave-functions. So for example, as illustrated in Fig. 10 the linear combination of the [0,1,0] neutron orbit around the α-cluster and [0,1,1] neutron orbit about the 16 O-core results in the two centred orbit which has the representation of [0,1,1] within the DHO. If the α-particle lies to the right-hand-side of the 16 O then the [0,1,0] wave-function is superimposed on the right-hand component of the clover-like wave-function of the [0,1,1] orbit (see Fig. 10). The 16 O+α system is mass-asymmetric and correspondingly the resulting molecular wave-functions do not possess good parity. Parity can be restored by creating linear combinations of the molecular wave-functions formed from having the α-particle on either the left or right side of the 16 O core. The resulting wave-functions then more closely resemble the symmetric wave-functions of the deformed harmonic oscillator.
Another asymmetric two-centre nuclear molecule that was investigated was 57 Ni, which could be thought of as 40 Ca and 16 O clusters, with a delocalised neutron, at a deformation of 2:1 in the deformed harmonic oscillator. Again, it is possible to construct DHO orbitals from combinations of the single-centre orbitals around each cluster. For    Table 2.
It is also possible to build nuclear molecules from unequal clusters but into symmetric systems. For example, 25 Mg can be constructed from the clusters α+ 16 O+α with a valence neutron. Such systems were also explored and found that again the wave-functions formed from the linear combination of those at the three centres created the available valence orbital for the neutron in the single centre. Indeed, the ordering of the clusters did not need to be symmetric to achieve this result. Examples explored include 73 Kr (two 16 O clusters and a 40 Ca cluster) and 97 Cd (two 40 Ca clusters and a 16 O cluster). The orbitals for the 73 Kr case are mapped in Table 3. The result is very similar to that in the case of 57 Ni.

Discussion and Conclusion
This extensive study of the mapping of molecular orbitals for symmetric and asymmetric cluster systems has found that in every case the molecular orbitals have a one-to-one correspondence with the valence orbitals for neutrons at the deformation given by the underlying cluster structure. In other words, molecular states are not just found for a few select systems but are a ubiquitous feature of nuclear systems where there is an underlying cluster symmetry.
The connection between the deformed harmonic oscillator and the appearance of molecular structures can be traced to the symmetries of the deformed harmonic oscillator [15] and the ideas of Harvey [16]. The spherical harmonic oscillator possesses a SU(3) symmetry and it has been demonstrated [15] that at a deformation of 2:1 and 3:1, etc..., that the deformed harmonic oscillator can be described by 2 or 3 coupled SU(3) groups and hence at a deformation of N:1 the appropriate symmetry is given by N⊗SU (3). In other words, the cluster structures are explicitly encoded in the symmetries. This is seen in the degeneracy patterns in Fig. 1 and the densities calculated in Figs. 2 to 5 and Ref. [12].
Harvey [16] noted the rules for combining two oscillator potentials which have been widely used to understand the nature and decay of cluster structures [14]. These rules relate to the combination of wave-functions from the two centres to the single centre and recognised it is possible to either preserve the total number of nodes in the wavefunction or create an additional node at the point where the two wave-functions connect. This means for two wave-functions with harmonic oscillator quantum numbers n z then two new wave-functions are produced with quantum numbers 2n z and 2n z + 1. These rules are depicted in the Harvey Diagram in Fig 12. For three centres, merging three wave-functions with the same n z values, produces three possible values of 3n z , 3n z + 1 and 3n z + 2. This can be further extended to N-centres.
For symmetric clusters, the mapping of the clusters observed in the deformed harmonic oscillator to their multi-centre counterparts fully accounts for all of the levels in the separated components and the levels up to the associated shell closure in the

Zero Separation
Infinite Separation Figure 12. The Harvey-diagram which demonstrates the connection between the orbitals of the single centre and two-centre system. Infinite separation is when the two clusters are far apart and zero separation is when they are fused.
deformed harmonic oscillator. The valence orbits above the deformed shell closure map to those outside the closed spherical core. The merging of the two centres into a single-centre oscillator is a process which respects the Pauli Exclusion Principle, PEP, inasmuch that one of the two merging levels with identical quantum number has an additional node generated. This is also the case for the merging of multi-centres. This effect results in the molecular orbitals not being the low-lying orbitals, which tend to be occupied by the core nucleons, and hence are at higher single-particle energy. The fact that the molecular, and indeed cluster, structures are consistent with the PEP is confirmed through the AMD calculations which explicitly respect the PEP.
For two identical clusters, e.g. two 16 O nuclei, the set of degenerate levels outside the core (e.g. 6 levels in 16 O) map two sets of levels at a deformation of 2:1 split in energy by 3/5hω (see Eqn. 3). These set of levels are identical to each other aside from the addition of one unit of n z (2n z and 2n z +1). These are the two degeneracies of 12 observed at a deformation of 2:1 in Fig. 1, and are the original two-centre levels which have undergone a 2n z and 2n z +1 transformation. Similarly, the set of levels associated with the deformed harmonic cluster structure, up to and including the shell closure, at the deformation of 2:1 are the levels of the separated 16 O clusters that have also undergone the 2n z and 2n z +1 transformation. So for the original 16  . These are again identical aside from a difference of n z =1. Thus, the 2n z and 2n z +1 transformation applies equally for the core and valence nucleons. This process of merging from multicentre to single centre preserves the number of nodes in the final merged wave-function plus creates a new set of wave-functions with an additional node along the separation axis. These principles apply to more than two centres, e.g. 3 and 4-centre cluster structures, but the transformation will be 3n z , 3n z +1 and 3n z +2.
For asymmetric cores the principle is the same. The merging of the two different valence orbitals around the two different cores preserves the number of total nodes in the final wave-functions and creates a second wave-function with one additional node. These resulting orbitals are the next available orbitals in the composite system, which then have a molecular structure. Linear combinations of orbitals at the two centres should have the same number of transverse nodes to create the molecular wave-functions.
The conclusion from this work is, beyond the observation that all the shell structures in the prolate deformed harmonic oscillator can be associated with cluster structures, that the next available orbitals beyond the shell-gap have molecular structure which is associated with the linear combination of the valence nucleon orbitals around the separated cluster cores. In this way molecular structures should be a wide spread phenomenon and that the molecular behaviour in nuclei is a ubiquitous phenomenon.
The spin-orbit interaction is known to play an important role in the breakdown of clustering, and though it is included in the effective interaction of the AMD approach [2] where the molecular structures of light nuclei are reproduced, there is no guarantee that the molecular structure persists to much heavier systems described in the present work. The confirmation of the robustness of the states predicted here needs to be tested with more microscopic approaches such as those described in Ref. [2]. Moreover, there are experimental challenges for the observation of such states, given that they exist at high excitation energy and will be embedded in a sea of states of other character. This will fragment the molecular states and experimental methods such as that described in Ref. [17] may be required.
This work is presently being extended to oblate deformations.