Phenomenological and microscopic cluster models II. Phase transitions
Abstract
Based on the results of a previous paper (Paper I), by performing the geometrical mapping via coherent states, phase transitions are investigated and compared within two algebraic cluster models. The difference between the Semimicroscopic Algebraic Cluster Model (SACM) and the Phenomenological Algebraic Cluster Model (PACM) is that the former strictly observes the Pauli exclusion principle between the nucleons of the individual clusters, while the latter ignores it. From the technical point of view the SACM is more involved mathematically, while the formalism of the PACM is closer to that of other algebraic models with different physical content. First and secondorder phase transitions are identified in both models, while in the SACM a critical line also appears. Analytical results are complemented with numerical studies on cluster states of the Ne and Mg nuclei.
pacs:
21.60.n,21.60.Fw,21.60.GxI Introduction
In a former contribution pap1 (), called from now on Paper I, we investigated the geometric mapping of the Phenomenological Algebraic Cluster Model (PACM) and the Semimicroscopic Algebraic Cluster Model (SACM) sacm1 (); sacm2 (). The first does not observe the Pauli exclusion principle while the second one does. The PACM belongs to the same family as models like the vibron model vibron (). In both type of models we considered the same Hamiltonian, while the model space of each is quite different. In the SACM, observation of the Wildermuth condition wildermuth () means that the number of relative oscillation quanta is restricted from below while in the PACM the relative oscillation quanta start from zero.
The method of geometrical mapping is gleaned from geom () for the SACM, which reduces to the usual one roosmalen1 (); roosmalen2 (); kirson () when no Pauli exclusion principle is taken into account. We showed in pap1 () that the differences in the mapped geometrical potential are large. However, within the PACM one can reproduce the geometric potential of the SACM by including very complicated higherorder interactions.
We also included in pap1 () a discussion on how to define correctly the coherent state parameter, such that for a large total number of bosons (, which plays the role of a cutoff) the mapped potential will be independent of the cutoff .
In this paper we concentrate on the study of phase transitions and show that not only secondorder phase transitions may occur, but also firstorder phase transitions. Furthermore, in the SACM a critical line appears beyond which no phase transition occurs. Thus the structure of the phase diagram of the SACM will be much richer.
We will also apply the mapping to two kind of systems, one with two spherical clusters and the other one with a deformed and a spherical cluster. Numerical studies are added. We will show that the PACM leads to inconsistencies when applied to real nuclei, while the SACM performs well.
The paper is structured as follows: in section II a general discussion on phase transition, both in the PACM and SACM, is given, independent of a particular system. In Section III, results are illustrated with two particular cluster systems, the cluster states of Ne and Mg. In the former, both clusters are spherical, while in the latter, one of them is deformed, leading to a more complex physical situation. Finally, in Section IV conclusions are drawn and a discussion is presented on the differences between the PACM and the SACM and their importance in the study of nuclear clusters.
Ii Phase transitions in algebraic cluster models
In what follows we apply the formalism developed in Paper I pap1 () to discuss phase transitions in the SACM and PACM. Before that we present the general framework within which the discussion will be implemented.
ii.1 Definition of a phase transition
Phase transitions are investigated using the following steps and the recommendations of greinerthermo (); octavio (); octavio2 (). This method can be applied to any system and does not depend on the notion of symmetries. This presentation does not need the language of catastrophe theory gilmore (), which simplifies considerations.
a) In the first step the minima of the potential energy surface (PES) are determined in the space of the collective variables . In the present case there is only one relevant variable . This is due to the fact that the distance vector between the clusters can always be aligned along the intrinsic axis which connects the two clusters. The extrema are obtained from
(1) 
This determines the position of the extrema at
() and the values of the potential
there. The are the values of the variable
at the minimum, which is a function of the interaction
parameters , with as a shorthand notation for
the parameters , with being the number of
parameters.
b) Once the extrema are obtained, one determines at each
minimum the first and second derivatives of the potential with
respect to the parameters of the model, i.e.,
(2) 
When we discuss the PACM and SACM further below, we will give first a
general discussion on the properties of phase transitions, independent
of the values of a particular subset of interaction parameters. In
the subsequent concrete applications, however, we will fix most
interaction parameters and vary only and ,
which control the transition from one
effective symmetry huitz1 () to another.
This will give us
particular curves in the space of phase transitions.
c) The locations of phase transitions are determined by identifying the
borders in the parameter space where at least two minima are at
equal energy, i.e.,
(3) 
for some index values and of the minima. This results in a relation
(4) 
between the parameters, which allows one, in principle, to express one
parameter in terms of the others (see the particular examples
further below).
d) The phase transition is of order when, up to , the
derivatives of the potential with respect to the free parameters are
equal at the point of phase transition, while the derivative of
the potential, with respect to its parameters, is discontinuous at the
point of the phase transition, i.e.,
(5) 
Note again that this procedure for determining the order of a phase transition is quite general and does not depend on identifying the phase with a dynamical symmetry. Thus, this procedure can also be applied to systems that do not exhibit dynamical symmetries.
ii.2 A possible explanation for the differences in using coherent states and a numerical approach
In the numerical calculations further below, we will see that, though the use of coherent states shows a clear phase transition, this is not the case in the numerical study, because the total number of bosons is fixed. In what follows, we indicate why this is so.
The main reason for using coherent states is that the phase transition is given by the crossing of two potential minima. In the case of a firstorder phase transition, the two minima are well separated by a barrier and a simple crossing takes place. For the case of a secondorder phase transition, the spherical minimum coincides with the “deformed” one at at a certain point in the parameter space. Once one minimum is deeper than the other one, an immediate jump by hand to the global minimum is made, using coherent states. This appears to produce the well defined phase transition. However, in numerical calculations within a finite system, the states are a mixture of states in both minima. Around the point of phase transition, the eigenstates of the Hamiltonian are a mixture of states in both minima. Only far away the states are located clearly in one or the other minimum. This produces, for finite , a gradual change in the wave function and as a consequence a gradual change of the control parameters, like the energy or the expectation value of , which define the phase transition. For, the discontinuities in the derivatives of the potential will be clear cut.
As a consequence, the apparently clear phase transition for any number of is not a real one, because is finite. Nevertheless, the use of coherent states shows the points of phase transitions and of which order they will be for .
ii.3 Study of phase transitions in the SACM
Paper I pap1 () contains the general expression of the Hamiltonian (in subsection 2.2) and the geometrically mapped SACM potential (in subsection 4.1). Here we recall only the essential formula, Eq. (30) of Paper I, necessary for the present discussion:
(6)  
We also note that the case relevant to the to phase transition has to be discussed separately. (See Eqs. (35) to (38) in Paper I pap1 ().)
ii.3.1 Phase transition diagram for the SACM: General considerations
The complex structure of the geometrically mapped potential complicates an analytic treatment of the problem. It is therefore essential to formulate a set of criteria facilitating a straightforward way to determine the order of phase transitions.
Equation (6) demonstrates that the explicit dependence of the potential on the parameters , and is linear of the type
(7) 
with and being a shorthand notation for the parameter. The are given by ratios of the functions and are always greater than zero for . Only the function approaches zero for . (See Eq. (39) in Paper I pap1 ().)
According to the general discussion on phase transitions in Subsection
II.1, let us now turn to the potential minima located at
. We investigate their dependence on the parameters
, standing for , and . The structure of
(7) guarantees that is always an
extremum. Furthermore, not only its firstorder derivatives, but the
function value is also zero at . It is thus
sufficient to focus on the second, deformed minimum, with , for which the following consideration holds. There are
two possibilities:
a) , or
b) .
The consequences are seen by determining the first absolute derivative
of the potential with respect to the parameter . This first
derivative is given by
(8)  
because vanishes at the minimum. Taking into account (7), this is further expressed as
(9) 
For case a) this expression is different from zero (remember that ), while the derivative within the spherical minimum is equal to zero. In other words, there is a firstorder phase transition.
For case b), (9) is equal to zero, and therefore the phase transition must be of higher order. To determine which order, secondorder derivatives are also needed:
(10)  
Here we used the fact that in (9) there is no explicit dependence in left. Since in this case , (10) reduces to
(11) 
As both the functions and the partial derivatives are different from zero, so is (11) in general, i.e. the phase transition is of second order in case b).
In conclusion we have the simple identification of phase transitions:
a) When the deformed solution is different from zero,
then there is a firstorder phase transition.
b) When the deformed solution is equal zero,
then there is a secondorder phase transition.
ii.3.2 General discussion of phase transitions in the SACM
In Fig. 1 the surface where a phase transition takes place is plotted in the space of the independent parameters , and . The solid line marks a change from firstorder to secondorder transitions. For larger, positive the transition is of second order, while for smaller, negative it is of first order. For completeness we show in Fig. 2 the phase space diagram for the case when , as it is the case for the to transition.
A remarkable finding is that at approximately the surface of phase transition ceases to exist. Beyond that point, no phase transition can be observed, i.e., the straight dashed line in Fig. 1 at approximately represents a critical line. This is shown by a dashed line in Fig. 1. One then can trace a straight line from below the surface, around the critical line, ending up above the surface without passing through a phase transition, which is similar to the critical point in the twodimensional phase diagram of water. Fixing all interaction parameters except and (as we will do in the numerical applications), passing from one dynamical symmetry limit to another one will trace a line in this threedimensional space. Depending on the fixed parameters, this line will or will not cross the surface of phase transition.
Figure 3 displays corresponding to the deformed solution at the point of phase transition as a function in and . (Remember that the spherical solution corresponds to the always existing extremum at in the case it is a local minimum.) Each point of the surface represents also a given at which the phase transition occurs, i.e., is fixed by the requirement that there is a phase transition. In this figure the solid line also represents the change from one type of phase transition to the other one. For larger, positive holds, and according to the discussion of subsection II.3.1, it corresponds to a secondorder phase transition. For smaller, negative the holds at the point of phase transition, i.e., it corresponds to a firstorder phase transition.
The whole phase structure of the system is illustrated in Fig. 4, which corresponds to the fixed value . The solid line represents a cut through the phase transition surface. For the phase transition is of second order, while for is of first order. The Roman numbers indicate the following regions: i) Region I corresponds to the existence of two minima (one spherical and the other deformed), with the deformed one as the global minimum. ii) Region II also corresponds to two minima with the spherical as the global one. iii) In region III there is only one spherical minimum while iv) in region IV there is only one deformed minimum. The region denoted by a zero refers to potential with no minimum. The other dashed lines do not indicate phase transitions, rather they separate the areas where two minima exist and the ones where only one minimum exists. The line of phase transition ends at approximately , corresponding in this case to a point.
There is another point on the solid line, where the firstorder phase transition turns over into a secondorder phase transition. A distinction has to be made between this and the critical point of the phase diagram of water, where before the critical point the phase transition is of first order, while at the critical point it is of second order. In our case, however, the phase transition does not end at this point, rather it continues as secondorder one in the domain with . One can now draw a curve which encircles this line, going from a spherical phase to a deformed one. The value of at which this phenomenon happens, varies with and .
The horizontal dotted line represents a division (see discussion in II.1): above this line the potential approaches a negative infinite value for , i.e., the potential gets unstable and the cluster system dissolves, once getting to the “deformed” solution. Below that line the potential approaches plus infinity for , thus the deformed solution corresponds to a stable cluster system.
Fig. 5 gives a sample of energy functionals at a fixed ( and ) in different regions, the specific points being marked with lowercase letters in Fig. 4. The upper left panel shows functionals either side of the firstorder phase transitions, points a) and b). The upper right one shows graphs either side of the secondorder phase transitions, points c) and d). The lower left panel shows a functional either side of the dotted line in the Region III of greater , points e) and f), showing that while each have a barrier, one is bound as , and the other is unbound. These tails are dependent on , with the unbound approaching infinity and the bound negative infinity as . Thus, this region is unphysical. The Region III of lesser , shown in the second panel, has no barrier, and is bound as . The fourth panel shows a fully unbound potential with no minima from Region 0, point g).
For completeness, Fig. 6 traces the value of for at in Fig. 4, which thus represents a movement between several regions in that figure (and is therefore analogous to being a subset of Fig. 3). It shows that the critical point corresponds to asymptotically approaching at slightly greater than 22.
Note that the discussion of phase transitions is completely independent of the system considered, which is of great advantage. In the next section we will consider particular systems and study the properties of their phase transitions.
ii.4 The PACM case: Pauli principle not taken into account
Here we apply the results obtained in subsection 4.2 of Paper I pap1 (). We display only the expression of the normalized potential in Eq. (46) of Paper I:
(12) 
We also remind the reader that the to transition needs a separate discussion, because in that case the potential reduces to a quartic one. (See the discussion in subsection 4.2 of Paper I pap1 ().)
ii.4.1 Phase transitions
The extrema are found by setting the first derivative with respect to equal to zero. We obtain, in general, two solutions
(13) 
noting that in this case the potential value is always , and
(14) 
where the value of the potential at the extremum with the positive sign is given by
(15) 
The barrier maximum is obtained by choosing the negative sign in front of the square root expressions. The last equation is acceptable only when the solution in (14) is real and positive. The solution is real for . When and holds, no deformed solution exist. In fact, inspecting Eq. (12) for there is always a minimum at , while a negative always gives a positive contribution, . When and are negative, the potential is negative at small and is turned over by the positive contributions and . In this case, only a deformed minimum exists.
The structure explained above is summarized in Fig. 7. The horizontal line is the axis while the vertical one is the axis. On the left hand side () a spherical minimum exists for , while for the minimum is deformed. The phase transition takes place at the line . This will be discussed further below as case b). For and , a deformed minimum always exists. For and below the dashed curve a deformed minimum coexists with a spherical minimum until is reached. The solid curve is determined by requiring that the potential minima at and are degenerate, i.e.,
(16) 
which leads to the condition
(17) 
Between the upper dashed and the solid curves the spherical minimum is the global one, while below the solid curve the global minimum is the deformed one. Crossing the solid line a phase transition takes place, which will be called below case a). It is of no surprise that the phase diagram in Fig. 7 is similar to Fig. 2 of cejnar ().
Case a), :
The potential is exactly zero at , and thus
(18) 
for all . Contrary to this, the first derivative of the potential at , the deformed minimum, taken at where the height of the deformed minimum equals that of the spherical one, is given by
(19) 
which is clearly different from zero for . Thus the phase transition is of first order for .
We can repeat the calculation for the derivatives in , leading to the same conclusion. Also, any combination of curves given by the variable (with and , with as the angle between the tilted straight line and the axis) leads to the same conclusion.
Case b):
We will vary from positive to negative values for a fixed . The same steps as taken for case a) can be applied, setting . The first and second derivatives of the potential with respect to are given by
Evaluating these derivatives at and gives
(21) 
For this case we also have to determine the derivatives with respect to , given at the point of phase transition, i.e., setting afterward, and ,
(22)  
Changing corresponds only to moving on the line of phase transition, which does not make sense, because one has to cross the line of phase transition. However, a combination of changing and corresponds to an inclined straight line, crossing the line of phase transition. This means that one has to consider a general line given by the variable , with arbitrary and , with being an arbitrary rotation angle. Determining at the point of phase transition the first and secondorder derivatives of with respect to leads to
(23) 
Because the first derivative of the potential with respect to is zero, while the second is different from zero at the point of phase transition, we can claim also that the first derivative with respect to is zero, while the second is different from zero. We conclude that the phase transition is of second order.
Iii Numerical studies
This section deals with two widely known cluster systems: Ne as O+, where both clusters are spherical, and Mg as Ne + , where one of them is deformed.
Since our aim is studying transitions from one particular dynamical symmetry limit to another one, we first fix all the interaction parameters except and . The parameter was not fixed, but rather it was chosen according to the harmonic oscillator constant corresponding to the unified nucleus. The fixed interaction parameters are determined in such a way that in the dynamical symmetry limits the spectrum appears with the same scale as the physical measured one. The spectrum of the real nucleus would probably correspond to a single point in the parameter space. However, our aim is not reproducing the exact spectrum, rather to reach a conceptual understanding of phase transitions when going from one dynamical symmetry to another, and investigating the conjecture that a phase is defined by an effective symmetry. One alternative method would be to adjust several cluster systems and to try to find a series of systems which, for example, would cross the surface of phase transitions at one point. We do not follow this way, rather postponing it for later consideration.
iii.1 Two spherical clusters: O+ Ne
In this case the only degree of freedom is the radial motion, as the clusters do not have an internal structure apart from the fact that they are composed of fermions, which have to obey the Pauli exclusion principle. Therefore the cluster representation is ) = . As described in Paper I pap1 (), the limit does not exist as an independent limit in this case (the Hamiltonian is a reduced version of the the Hamiltonian), so the choice has to be made. The only transition to consider is thus between the and the limits.
Concerning the determination of the parameters, one has to take into account that some parameters appear in both dynamical symmetry limits, like preceding (see Eqs. (10) and (11) in Paper I pap1 ()). We first determine the parameters in the dynamical symmetry limit, which fixes , and then we determine the remaining parameter , which appears in the dynamical symmetry limit. The terms and do not contribute to the Hamiltonian in this case, so the corresponding parameters are kept zero. Note that in this case .
iii.1.1 The SACM
In the first step the parameters are adjusted within the limit, setting and and in the limit, setting and . The parameters are depicted in Table III.1.1. The SACM yields reasonable results, because the groundstate band belongs to , , where corresponds to the minimal number of relative oscillation quanta required by the Wildermuth condition. The first excited state corresponds to a 2 excitation and naturally lies at high energy as required by the experimental data. The spectra in the and limits are shown in the left and right extreme of right panel of Fig. 11, respectively. The spectrum of experimental Ne states, each corresponding to this clusterisation, is shown in the left panel. As already mentioned, the real nucleus will lie somewhere between and , though, the fit at is acceptable.



In order to see if a phase transition appears, and of which order it is, we added a curve to the phase space in Fig. 8, depicting the transition from to . The figure shows only the relevant part of the phase space, i.e. the one where the curve crosses the surface associated to the secondorder phase transition. This occurs at approximately the parameter value. The actual values of , and associated with this can be seen in Fig. 9. In summary, the situation in this example corresponds to a phase transition of second order.
Figure 10 displays how the expectation value of changes as the function of . The lighter (yellow online) and darker (red online) curves depict the result of the geometrical mapping and the numerical diagonalization, respectively. From up to the point of phase transition, the effective limit is realized and the expectation value is equal to the minimal number of bosons, i.e., . Below the darker (red) curve begins to rise slightly, indicating that the structure of the system is changing.
iii.1.2 The PACM
This type of model was considered in zhang (), where a pure schematic investigation on possible phase transitions was presented. The model space was restricted to low excitations. No parameter fit was applied to a physical system. We show here that the relation to a physical system is of utmost importance and can discriminate between physical and nonphysical models.
As a first step, we tried to adjust the parameters of the model in the
limit ( and ) and in the limit ( and
). However, we already encountered severe problems in the
limit: The model space for even angular momentum starts with
, thus the lowest states are comprised by , , which contains one state, and by
, , which contains a and a
state. The problems encountered are:
i) Since the factor in front of has to be positive and
supposing that the ground state belongs to and the next
excited positive parity state belongs to , the first
excited state will always be lower than the first excited
state. This contradicts the experimental spectrum with
MeV and MeV.
ii) Since the state belongs to a 2 excitation,
with MeV, the quadrupolequadrupole interaction
has to be unnaturally strong in order to shift the energy to
1.634 MeV. This, in turn, will move very high excitations to
low energy, even below the supposed ground state with .
iii) Due to the completely different structure of the states
within the ground state “band”, one can not talk about a rotation
band.
Restricting to the limit, the eigenvalues of the Hamiltonian are given by
(24) 
Here we already see that for and a fixed , higher spin states are higher in energy. In order to adjust the and the states to the experimental energies, the fitting routine assigns to both the ground state and the state a different . Using , as a result we obtain (the total number of bosons was set to be 20) for the and the state, while the state belongs to . When we change the total number of bosons, we get similar results. The mere fact that we have to involve states with indicates that no convergence is achieved, considering that represents a cutoff value.
Similar results are also obtained when the limit is considered. In the limit, the energy is given by
(25) 
where refers now to the quantum number with = , , …, 0 or 1. The lowest state is normally taken as , which contains . Thus, choosing, the groundstate band () is composed of the angular momentum states = 0, 1, …, 20 and the first excited band with even spin ( = 18) is given by the states = 0, 1, …, 18. Thus, the first excited state can be set at higher energies than the first excited state, adjusting the parameter . Everything seems to be in order, except for the problems which the following discussion demonstrates.
The difference with respect to the limit is that a state is a mixture of many basis states in . We adjusted the spectrum of Ne to the limit and confirmed that for the expectation value of the operator is given by frank ()
(26) 
The problem here is that when the cutoff is increased, the structure of the states change: increases, implying that no convergence has been reached. This also implies high shell excitations, if we assume that the two clusters are moving in a shell model mean field, which has been proven in many microscopic calculations TE1 (). In the PACM, however, the mean field has no specific meaning, i.e., the parameter can be very small. As shown above, the spectrum can be easily adjusted within the limit. The fact that it adjusts the spectrum suggests that there must be some truth in it. Nevertheless, the basic degrees of freedom (clusters plus relative motion) cannot be interpreted as real clusters or relative motion, as is done in microscopic cluster studies, but must be in a complicated relation with them. Just what relations these are remains a big problem.
This result demonstrates that the model shows inconsistencies and the reasons are exposed in the arguments i) and ii) above. The interpretational problems of the PACM have already been indicated in huitz1 ().
A possible solution to this problem is to redefine the pairing operator as , i.e., the introduction of a new parameter which can be set proportional to . In this way, the dependence can be eliminated and physical results can be expected. This procedure was adopted in bijker (), where the vibron model was extended to three clusters describing the C nucleus as an oblate symmetric top. Though, in bijker () the limit in a algebraic model, which describes three clusters, is considered, the same ideas can be translated here.
iii.2 One spherical and one deformed cluster: Ne+ Mg
This is the first example where the cluster part has a structure owing to the deformed Ne.
iii.2.1 The SACM
In order to analyze the transitions between the , and limits, first the parameters were fixed in the three limits. The results are displayed in Table III.2.1.



Figures 12 and 13 display the expectation value of and the lowest states in energy, respectively, for all three transitions. The darker (red online) curve in Fig. 12 shows the results of the numerical diagonalization, while the lighter (yellow online) curve shows the one of the geometrical mapping. The results are qualitatively similar to those encountered in the former example. The expectation value of starts in the limit at 8 and increases towards the limit. The transition is smooth in the numerical calculation, while in the geometrical mapping the transition is pronounced well. The transition is indicated by a sudden change in the slope at above in the central panel of Fig. 12. The energy spectrum in Figure 13 does not show a particular structure at points of phase transition. Therein, the leftmost panel is again the experimental spectrum used, and the right three are theoretical results as and are adjusted. In the to transition no phase change appears, because the global minimum of the potential is always at =0. One can observe a distinct behavior as a function of and between the states with positive and negative parity, marked as solid and dashed curves, respectively. The latter are more sensitive to the change in and .
iii.2.2 The PACM
Here we find similar inconsistencies with respect to the model space as in the case of Ne in III.1.2. The cluster irrep of Ne is (8,0), while the relative oscillation irreps are , with = 0,1,2,… Restricting to positiveparity states only, the lowest energy model space, in the limit, consists of at 0, and (10,0) and (6,2) for the 2 excitation. The groundstate band is a (8,0) irrep and the lowest band is the (10,0) irrep at 2. Even if we change to an arbitrary small parameter , in order to bring the 2 states down in energy, the internal structure of the ground state band is not what we expected, namely (8,4). Also there is no B transition between the ground state band and the band, because they belong to different irreps. (See the discussion in Subsection 2.1 in Paper I pap1 ().)
For the limit one obtains a satisfactory fit, but again with the fact that the expectation value of the number operator depends on the cutoff .
For these reasons, we do not present figures of the spectra and expectation values.
Iv Conclusions
Phase transitions were investigated in two algebraic cluster models, one of which observed the Pauli exclusion principle between the nucleons of the individual clusters (SACM), while the other (PACM) did not. This analysis was based on the results of a previous work pap1 (), in which the geometric mapping of the two models had been performed using the coherent state formalism, leading to appropriate potential energy surfaces. The dynamical symmetries of the SACM and PACM had also been identified in pap1 (), and in the present analysis special attention was paid to transitions between phases associated with the , and dynamical symmetries. The potential energy surfaces depended on the parameters appearing in the Hamiltonian shared by the SACM and PACM, including also the and variables controlling the transitions between the three dynamical symmetries. The phase space was reparameterized in terms of three parameters , and .
In the case of the PACM, the potential energy surface was a sextic oscillator in the intercluster distance variable, while for the SACM the potential shape was more complex due to the restrictions enforced by the Pauli principle. The potential energy surface typically contained up to two minima, one spherical and one deformed. The analysis identified both first and secondorder phase transitions for the PACM and the SACM, while in the latter case a critical line was also found.
The results were illustrated with numerical studies on the O+ and Ne+ systems, which correspond to two spherical clusters and to one spherical and one deformed cluster, respectively. The limit was found to be the most appropriate one in reproducing the data of the cluster systems. Clear phase transitions were identified in the parameter controlling the transition between the and limits. It was found that the PACM led to energy spectra that are rather different from the observed physical ones.
Acknowledgments
We gratefully acknowledge financial help from DGAPA, from the National Research Council of Mexico (CONACyT), OTKA (grant No. K72357), and from the MTACONACyT joint project. POH acknowledges very useful discussions with Octavio Castaños (ICNUNAM), related to the differences in phase transitions in finite systems to the use of coherent states. The authors are also thankful to József Cseh for illuminating discussions on the subject.
References
 (1) H. YépezMartínez, P. R. Fraser, P. O. Hess, and G. Lévai, Paper I.
 (2) J. Cseh, Phys. Lett. B 281 (1992), 173.
 (3) J. Cseh and G. Lévai, Ann. Phys. (N.Y.) 230 (1994), 165.
 (4) F. Iachello, Phys. Rev. C 23 (1981), 2778.
 (5) K. Wildermuth and Y. C. Tang, A Unified Theory of the Nucleus (Friedr. Vieweg & Sohn Verlagsgesselschaft mbH, Braunschweig, 1977).
 (6) P. O. Hess, G. Lévai and J. Cseh, Phys. Rev. C54 (1996), 2345.
 (7) O. S. Roosmalen, Algebraic Description of Nuclear and Molecular RotationVibration Spectra, (Ph.D. Thesis, Groningen, Netherland, 1982).
 (8) O. S. Roosmalen and A. E. L. Dieperink, Ann. Phys. (N.Y.) 139 (1982), 198.
 (9) A. Leviatan and M. W. Kirson, Ann. Phys. (N.Y.) 188 (1988), 142.
 (10) W. Greiner, L. Neise and H. Stöcker, Thermodynamcis and Statistical Mechanics, (Springer, Heidelberg, 1995).
 (11) E. LópezMoreno and O. Castaños, Phys. Rev. C 54 (1996), 2374.
 (12) O. Castaños, R. LópezPeña, J. G. Hirsch and E. LópezMoreno, Phys. Rev. B 74 (2006), 104118.
 (13) R. Gilmore, Catastrophe Theory for Scientists and Engineers (Wiley, New York, 1981).
 (14) H. YépezMartínez, J. Cseh and P. O. Hess, Phys. Rev. C 74 (2006), 024319.
 (15) L. Parra Rodríguez, Master Thesis, UNAM, 2011.
 (16) P. R. Fraser, H. YépezMartínez, P. O. Hess and L. ParraRodríguez, J. Phys. Conf. Proc. (2011), accepted.
 (17) Cejnar P and Iachello F, 2007 J. Phys. A 40 581
 (18) Y. Zhang, Z.F. Hou, H. Chen, H. Wei and Y.X. Liu, Phys. Rev. C 78 (2008), 024314.
 (19) A. Frank, P. Van Isacker, Symmetry Methods in Molecules and Nuclei, (SyG editores, México D.F., 2005).
 (20) G. Rosensteel and J. P. Draayer, Nucl. Phys. A 436 (1985), 445.
 (21) R. Bijker and F. Iachello, Ann. Phys. 298 (2002), 334.