1 Introduction
Nucleons are easy to assemble and disassemble because of the saturation property of binding energy and density, which is one of the important properties that characterize the nuclear many-body system. Nuclear clustering [
1-
3] is the physics of dynamical assembling and disassembling of nucleons. Formation of clusters is a fundamental aspect of nuclear many-body dynamics together with the formation of mean field.
In the case of liquid-gas transition of infinite nuclear matter, the nuclear matter in the transition region is the mixture of liquid nuclear matter and nucleon-gas. But in the case of finite nuclei, the nucleus in the transition region of liquid-gas transition is composed of not simple mixture of liquid and nucleon-gas but mixture of various clusters, namely hot gas of clusters. Figure 1 shows the constant-pressure caloric curve of the finite nucleus
36Ar (
) obtained by AMD (antisymmetrized molecular dynamics) in Ref. [
4]. We can see that the figure shows clearly the existence of nuclear liquid-gas phase transition in finite nucleus. We see the existence of the region of negative heat capacity which is a characteristic feature of finite systems and which corresponds to the plateau region of the constant-pressure caloric curve in the case of infinite nuclear matter. Figure 2 shows fragment mass distribution for the ensembles along the
MeV/fm
3 line. When the energy is low (
MeV), the distribution shows
U-shape with two peaks so that the typical configuration at this energy is a large nucleus coexisting with a few gaseous nucleons. When the energy is increased (
MeV), the peak at the large fragment becomes smaller and the distribution changes into shoulder-like and power-law-like shapes. Thus complex configurations with many intermediate and light mass fragments are typical at these energies, and the proportion of light fragments increases as the energy increases (from
MeV).
Recently, study of (cold) cluster gas has started in the field of nuclear structure. Since the breakup thresholds into clusters are much lower than the threshold of nucleon gas, cluster-gas-like structure can be studied spectroscopically. Figure 3 shows excitation energies of
α-cluster gas-like state and nucleon gas-like state in the case of
12C. It is now widely accepted that the second
state which is known as the Hoyle state is the
-cluster gas-like state having the
-condensate-like structure [
5,
6]. The excitation energy
MeV of this
state is slightly above the 3
α and
8Be
thresholds.
For the study of the Hoyle state, a new type of cluster-model wave function was introduced [
8] which was named THSR (Tohsaki–Horiuchi–Schuck–Röpke) wave function. The THSR wave function in the case of
system has the form
Here
Xk and
XCM stand for the center-of-mass coordinates of the
cluster and the total system, respectively. The intrinsic wave function of the
kth cluster is represented by the
. For the Hoyle state the size parameter
B is much larger than the single-nucleon H.O. (harmonic oscillator) size parameter
b. It was discovered that the
cluster model wave functions of the Hoyle state obtained by solving
RGM (resonating group method) and GCM (generator coordinate method) equations in 1970’s are almost 100% equivalent to the single THSR wave functions of the Hoyle state [
9]:
In spite of the fact that the THSR wave function was devised for describing gas-like cluster states, it was found that the THSR wave function can describe very well also non-gaslike ordinary cluster states with spatial localization of clusters [
10]. It was reported that the traditional RGM/GCM wave functions of the inversion-doublet band states of
20Ne are almost 100% equivalent to the
16O+
α single THSR wave functions [
10];
This discovery clarified that the spatial localization of clusters in two-cluster systems is due to the kinematic effect of inter-cluster Pauli principle [
11], which means that the dynamical motion of clusters is of nonlocalized nature. This result teaches us that the cluster dynamics is of nolocalized nature not only in the gas-like cluster states but also in non-gaslike cluster states.
The nonlocalized cluster dynamics can be well formulated by a newly introduced model which is named container model of cluster dynamics [
11]. The container model gives us both gas-like state and non-gaslike states as the solutions of the Hill–Wheeler equation with respect to the size parameter of THSR wave function
which is just the size parameter of the container,
Considering the fact that, in the limit of the smallest value of the size parameter B, the THSR wave function is just equivalent to the shell-model wave function, the Hill–Wheeler equation with respect to the size parameter of the container describes the evolution of the cluster structure from the ground state with shell-model like nature to the ordinary cluster states and further to the gas-like states of clusters.
In order to describe variety of cluster states which show up in the process of evolution of cluster states, the THSR wave function should be extended from the simple form with one kind of size parameter to the form of two or more kinds of size parameters. For the sake of the development of the study of cluster dynamics, it is necessary to devise better extension of the THSR wave function.
The content of this review article is as follows. In the next section (Section 2), we discuss gas-like states of clusters. Since the Hoyle state is now widely accepted as a gas-like state having the
-condensate-like structure, we first review the research of the Hoyle state. Early in 1970’s, the Hoyle state was recognized as being a gas-like state with very large radius composed of
clusters interacting weakly via
relative
S waves [
12-
14]. Then, more than 20 years later, the Hoyle state was reinterpreted as having Bose-condensate-like structure of
clusters [
8]. The THSR wave function was devised for studying
α condensate-like structure and it was found the THSR wave function is almost 100% equivalent to the old
cluster-model wave functions [
9]. It is discussed that the recent quantum Monte Carlo calculation gives the wave function of the Hoyle state [
15] whose density distribution is quite similar to that of the
THSR wave function namely that of the old
cluster-model wave function [
7]. We also review recent studies of the breathing-like excited state [
16-
18] of the Hoyle state which is assigned to the
state of
12C. Next we review the research of the
condensate-like state in
16O which is assigned to the observed 6th
state at 15.1 MeV. We explain how this assignment of the
condensate-like state is supported by theoretical calculations and their comparison with experiments. Theoretical calculations are those by
OCM (orthogonality condition model) [
19] and by the GCM with
THSR wave functions [
20] with respect to the size parameter in Eq. (4).
The THSR wave function is a very novel cluster wave function and it has some different features compared with the traditional cluster model. In Section 3, the characters of the THSR wave function will be discussed. Firstly, the THSR wave function has the shell-model limit when the size parameter becomes very small while it also has the free-cluster limit when the size parameter becomes very large. Secondly, it was found that the prolate and oblate THSR wave functions are almost equivalent after angular-momentum projection. This is because the rotation-average of a prolate THSR wave function is almost equivalent to an oblate THSR wave function. Thirdly, although the original THSR wave function has the positive-parity character, by introducing another shift parameter, the THSR wave function can also handle the negative-parity cluster states. At last, we also introduce the single-particle property of the THSR wave function of 13C.
In Section 4, we discuss that the THSR wave function can describe very well also non-gaslike ordinary cluster states with spatial localization of clusters. It is explained that the spatial localization of clusters in two-cluster systems is due to the kinematic effect of inter-cluster Pauli principle. Namely the inter-cluster Pauli principle gives rise to the Pauli-forbidden states of the inter-cluster relative motion which are against to the close approach of two clusters. With this understanding of the spatial localization of clusters, it is explained that the dynamical motion of clusters is of nonlocalized nature just like in the case of gas-like states. We review that this understanding was obtained from the study of the inversion doublet bands in
20Ne [
21] by the use of the THSR wave function for
16O+
α clustering [
10,
22]. We explain that the inversion doublet bands in
20Ne have been regarded as a convincing evidence of the spatial localization of clusters since the intrinsic structure with
16O+
α clustering has parity-violating deformation. The equivalence of the RGM/GCM wave function with the THSR wave function shown in Eq. (3) cast a strong doubt on the use of the energy curve by the Brink wave function for determining the spatial localization of clusters. The
16O+
α Brink wave function has the inter-cluster distance parameter
D and the optimum value of |
D| is obtained from the minimum energy point of the energy curve. It is explained that in Ref. [
10] it was shown that this way of determining the inter-cluster distance is misleading and cannot be adopted. Since we prove that the THSR wave function is equivalent to the shell-model wave function in the limit of the smallest value of the size parameter, we see that the THSR wave function can describe, in addition to the gas-like cluster states, the ordinary non-gaslike cluster states and also shell-model-like states.
In Section 5, we explain the container model of cluster dynamics. We discuss the solutions of the Hill–Wheeler equation with respect to the single size parameter B in the case of and systems. The solutions of the Hill–Wheeler equation especially in the case of the system show us that, for describing variety of observed cluster states, the THSR wave function should be extended from the simple form with one kind of size parameter to the forms of two or more kinds of size parameters. We review some recent results obtained by using extended THSR wave functions having more than one size parameters. We see that in the case of 12C the extended THSR calculations predict rich spectra of states above the Hoyle state which are in good accordance with experiments when observed data are available. We also report some recent results of the extended THSR calculations in the case of 16O. The extended THSR calculations are necessary for the study of neutron-rich nuclei because the container for only clusters needs to be different from the containers accommodating valence neutrons. We review the studies of neutron-rich Be isotopes. We also introduce the new development of analytical treatment of the THSR wave function.
In Section 6, we discuss the study of evolution of cluster structure by the container model in general. We discuss the duality property of the shell-model ground state which is well described by the THSR wave function as is seen in the fact that the THSR wave function in the limit of the smallest value of the size parameters is equivalent to the shell-model wave function. The important ingredients of the container model is firstly the nonlocalized motion of clusters which is described by the Hill–Wheeler equation of motion with THSR wave functions and secondly the inter-cluster Pauli principle coming from the Fermi statistics of nucleons. We explain that the actions of the inter-cluster Pauli principle are well studied by calculating the Pauli-allowed states of cluster systems, not only two-cluster systems but also many-cluster systems. Based on these studies we discuss the evolution of cluster structure from the ground state with shell-model structure to excited states composed of two clusters and further to excited states composed of many clusters including the gas-like states of clusters. Finally, we give a summary in Section 7.
In this review paper, we are mainly focusing on the concept of nonlocalized clustering originating from the THSR wave function and its applications in nuclear cluster physics. It is worth mentioning that the clustering plays a very important role in understanding nuclear structure and there are many related studies in this field recent years. For example, some
ab initio methods [
23,
24] and also mean-filed approaches like density functional theory [
25-
27] have been applied to study formation of clustering in nuclei. The monopole [
28] and dipole transitions [
29,
30] are proposed as two effective probes for identifying clustering structure from experiments. More information about recent progress of nuclear cluster physics can be referred to Ref. [
3] and the included references.
2 Hoyle state and gas-like structure of clusters
In this section, we discuss gas-like-structure states with Bose–Einstein condensate character of
α clusters. According to the Ikeda threshold rule [
31], gas-like
cluster states in self-conjugate
nuclei are expected to appear around the
decaying threshold energy. These energies are in general low enough for the spectroscopic study to be possible, whereas when we consider nucleon gas, its excitation energy is expected to be much higher. The typical example is the Hoyle state in
, whose excitation energy is
MeV, located at
MeV above the
threshold energy. The structure study of the Hoyle state is done by introducing the so-called THSR wave function, and its gas-like and Bose–Einstein condensate characters are investigated.
2.1 condensate-like structure of the Hoyle state
2.1.1 Reduced α-decay width and S-wave dominantstructure of the Hoyle state
The Hoyle state, which is predicted by Fred Hoyle, is known to play a crucial role in the helium burning in red giant stars [
32,
33]. This state is very closely located at the
breakup threshold and is well known as one of the mysterious
states in light nuclei. The understanding of its structure has been actually one of the most difficult and challenging problems of nuclear physics. Its small excitation energy of
MeV has been regarded to be difficult to be explained by the shell model. The no-core shell model, which is the most advanced modern shell-model approach, has so far not succeeded to reproduce it [
23,
34]. About 50 years ago, Morinaga proposed a famous idea that this state was to be assigned as a linear-chain state of
clusters [
35,
36]. However, a simple theoretical assumption of linear-chain structure turned out to give the reduced
α-decay width of at most one-third of the Wigner limit value, while from observation it is much larger and exceeds the Wigner limit value [
37]. This is because the linear-chain state necessarily includes high partial waves between
and
α clusters, whereas only the
S wave contributes to the
α decay, due to the small
Q value for the Hoyle state.
The observed large
α-decay reduced width of the Hoyle state in the
channel was successfully reproduced by semi-microscopic three-body calculation [
12], called the Orthogonality Condition Model (OCM) [
38-
40]. This
OCM calculation gave the
S-wave dominant structure between
and
α. Since
is composed of loosely coupled two
α clusters in relative
S wave, the Hoyle state was concluded to have a loosely coupled
structure in relative
S waves with large spatial extension, namely a gas-like structure of the
clusters, which is quite different from the
linear-chain structure. This understanding of the structure of the Hoyle state was confirmed, a few years later, by full microscopic
calculations, i.e., the
RGM calculation [
14] and the
Brink-GCM calculations [
13], where the Brink basis function is used [
41]. The
Brink-GCM calculation with resonance boundary condition was also performed by Brussels group [
42].
Both methods, the
RGM and
Brink-GCM, are to fully solve the relative motions of the
α clusters, and are proved to be mutually equivalent [
43,
44]. For the former method, the
RGM wave function is obtained by solving the following equation of motion,
where and is the internal wave function of the α cluster with the four nucleons inside it being constrained to harmonic oscillator configuration, and is the relative wave function between the clusters to be obtained.
These full microscopic calculations nicely reproduced not only the excitation energy of the Hoyle state but also other experimental properties including inelastic electron form factor and and transition properties. In Table 1 we show that the experimental data available in are nicely reproduced by the RGM and GCM calculations. The S-wave dominant nature of the Hoyle state mentioned above was again clarified from the microscopic point of view. Figure 4 shows the reduced width amplitude (RWA) of in calculated with the GCM, which is defined below:
For the Hoyle state ( state), the S-wave component, which is denoted as , is shown to be dominant and to have the largest amplitude in outer region at around 5 fm, while the inner oscillation is suppressed, which typically characterizes a cluster structure.
2.1.2 Alpha-condensate-like character and the THSR wave function of Hoyle state
The THSR wave function is originally introduced to represent the Bose–Einstein condensate states in self-conjugate
nuclei in 2001 [
8], and later it was extended so as to include deformation [
45], in which the following operator form is assumed:
where is a creation operator of one nucleon with spin σ, isospin τ at a spatial point r, and is a harmonic oscillator (H.O.) wave function centered at R with the form of . The coordinate representation of this THSR state is expressed as
where
is the Brink wave function in the
cluster system [
41]. The
Brink wave function can be written in a form of Slater determinant and another equivalent form as shown below:
where Xi and are the center-of-mass (c.o.m.) coordinates and the internal wave function of the i-th α cluster, respectively, and is described by the H.O. shell-model wave function as follows:
and is the antisymmetrizer exchanging the nucleons belonging to different α clusters. The integration over the position parameters R1,…,Rn in Eq. (9) can simply be performed by substituting Eq. (11) into Eq. (9), which leads to
with
and , where and for , and . and are the total spurious c.o.m. wave function to be eliminated and the c.o.m. free form of the THSR wave function, respectively. The above coordinate representation of the THSR wave function shows that all the constituent α clusters take an identical motion, occupying the same deformed orbit, with the size parameter Bk , under the condition that the effect of the antisymmetrizer is negligible. This condition is satisfied when the size parameter Bk is much larger than the size parameter of the α cluster, b, or equivalently . In this situation, the α clusters can move rather freely inside a mean-field-like potential that is produced self-consistently from a nucleon-nucleon interaction. The parameters Bk or then parameterize the width of the mean-field-like potential.
On the contrary, if the width parameters Bk are small enough, the α clusters are resolved by the stronger effect of the antisymmetrizer , so that there only remains the α-cluster-like correlation. Furthermore in the limit of , or , the THSR wave function Eq. (13) coincides with the shell-model Slater determinant. This can be easily proved by considering the well known fact that the Brink wave function coincides with the shell-model Slater determinant in the limit of . By using an alternative expression of the delta function, , this limiting operation in Eq. (9) gives the Brink wave function in the limit of , namely the shell-model Slater determinant.
We should note that this wave function has a parity symmetric form, as is clearly understood from this coordinate representation in Eq. (13). We later extend this THSR wave function to parity-violating forms, to discuss inversion doublet bands of and Be isotopes in the following sections. While this form of Eq. (13) also allows us to grasp easily its physical picture such as mentioned above, it is very difficult to handle it as it is in practical calculations, since Eq. (13) is written in a general form of -body microscopic wave function so that the nucleons should be explicitly antisymmetrized by the operator . Instead for practical calculations it is much better to start from the form of the Brink wave function, Eq. (9), since its Slater determinant form, Eq. (10), makes the calculations of matrix elements be much more tractable.
Starting from the form of Eq. (9), we can use the following transformation, to factorize the spurious c.o.m. motion, like Eq. (13),
where is the set of Jacobi coordinates transformed from the position-parameter set , with the relation, for and , and is the internal part of the Brink wave function defined by
While in Eq. (15) the first half which is in the first line is nothing but the spurious total c.o.m. wave function, in Eq. (13), the latter half in the second line is the c.o.m. free THSR wave function defined in Eq. (14), i.e.,
This is written in the form of superposing c.o.m. free Brink wave functions, which allows us to easily perform any calculations. However, one might be afraid that the c.o.m. free Brink wave function
in the above equation deviates from the pure Slater determinant,
, and hence the above non-Slater-determinant form is difficult to be handled in practical calculations. However, once we know an analytical formula for the matrix element of an operator
, with a pure Slater determinant, like
, the matrix element to be desired,
, is also easily derived, since the transformation matrix
is also analytically given [
7].
Since now the parameter in the THSR wave function can be taken in an anisotropic way, the THSR wave function in Eq. (17) is not rotationally symmetric. The recovery of the rotational symmetry can be done by employing the following well known angular-momentum projection method:
with , Ω, , , and , are Euler angle, a rotational operator, Wigner D-function, and a normalization constant, respectively. Then in order to calculate a matrix element of any operator for the THSR wave function , i.e., , one should first calculate the matrix element for the intrinsic Brink wave function , i.e., . This is basically represented as a form of Gaussians. Thus integrating over the set of variables is shown to result in multiple Gaussian integral, and therefore essentially it can be performed in an analytical way, while the integration over the Euler angle Ω can be numerically performed. This aspect will be discussed in more detail in Section 3.
All calculations in this section are performed within axially symmetric deformation, where the z-axis is taken as the symmetry axis: . The above relation of Eq. (18) is simplified and can be redefined as
with a redefined normalization constant .
The energy eigenstates can be represented by superposing the above wave functions, Eq. (19), as follows:
The coefficients can then be determined by solving the following Hill–Wheeler equation,
where the microscopic Hamiltonian
is composed of the kinetic energy with the c.o.m. kinetic energy subtracted,
T-
TG, nucleon-nucleon effective interaction
, and Coulomb interaction
. The important criterion to choose an effective nucleon-nucleon interaction is to reasonably reproduce the
α+
α scattering phase shift [
46].
2.1.3 Equivalence of the 3α THSR wave function to 3α RGM/GCM wave functions
As we discussed in Section 2.1.1, the fully microscopic cluster model calculations of
RGM and
GCM were performed for
around 40 years ago [
2,
13,
14,
42]. The wave function for the Hoyle state reproduces the experimental data available, such as the binding energy,
α-decay width, monopole transition strength to the ground state, and inelastic electron scattering [
2]. We showed in Section 2.1.1 that the Brink-GCM wave function for the Hoyle state dominantly has
component, together with a large r.m.s. radius (3.4 fm), compared with the one of the ground state (2.4 fm). Thus the Hoyle state has been considered to have a well-developed
-cluster structure, which are loosely coupled in relative
S-waves, rather than the
-linear-chain structure. In the
RGM, or equivalently in the
Brink-GCM, relative motions between the
α clusters are solved without any model assumption, as shown in Eq. (5) in Section 2.1.1, in which the
RGM wave function can be expressed in the following form:
where the Jacobi coordinates for . On the other hand, as we can easily see from the r-space representation of the THSR wave function in Eq. (13) for , the relative wave function between the clusters in the THSR wave function is written like
Obviously the special form is assumed for the relative wave function between the
clusters in the THSR wave function, compared with the RGM formalism. What we try to discuss here is then how well the THSR wave function describes the
RGM/GCM wave function that is the full solution of the microscopic
-cluster problem. And as we emphasized repeatedly in the preceding sections, we show that the THSR wave function can perfectly describe the RGM/GCM wave function for the Hoyle state [
7,
9,
47,
48].
The energy spectrum and the corresponding wave functions for
in the THSR ansatz explained in the previous section are obtained by solving the Hill–Wheeler equation Eq. (21). Table 2 shows comparison of the calculated energies, r.m.s. radii, and monopole matrix elements with those of the
RGM/GCM calculation and the experimental data. In the GCM [
13] and RGM [
14] different effective nucleon-nucleon forces are adopted. In the former Volkov No. 1 with Majorana parameter
is adopted and in the latter a slightly modified version of Volkov No. 2 with
[
2] is adopted. Hence we adopt the same effective forces as adopted in the GCM and RGM to give precise comparison of our results with their results. We denote the former as Force 1 and the latter Force 2 in Table 2. One can see that the THSR wave function succeeds in reproducing the properties of the ground state and the Hoyle state. The large value of the r.m.s. radius for the Hoyle state indicates that the state has a volume 3 to 4 times larger than that of the ground state of
. It should be noted that for the Hoyle state the energies calculated with the THSR ansatz are lower than those calculated with the RGM/GCM. In particular, the THSR ansatz gives considerably lower energy, by
MeV, than the GCM calculation. This is surprising since, according to the framework of the formalism, no model is assumed with respect to the relative motions between the
clusters for the RGM/GCM, and therefore, due to variational principle, binding energies obtained by RGM/GCM must be lower than those of the THSR ansatz. Also the calculated r.m.s. radius for the Hoyle state is larger than those calculated with the RGM/GCM. In particular in comparison with that of the GCM, THSR ansatz gives much larger r.m.s. radius, by 1 fm. This is because the THSR wave function can much better describe a long tail nature which the Hoyle state located slightly above the
threshold is expected to have, than the GCM wave function, in which the finite number of spatially localized
-cluster basis functions, (Brink wave functions), are superposed in the practical calculation.
For the ground state the THSR wave function gives slightly lower binding energy than the RGM wave function obtained, though in comparison with the GCM wave function, it gives slightly higher energy. This suggests that even for very compact cluster state like the ground state the THSR ansatz works very well. We will discuss later an extended version of the THSR wave function, in which a
-like correlation is taken into account following a “container picture” discussed in Section 5, which is shown to give identical value to the one of the GCM [
50].
The inelastic form factor of
from the ground state to the Hoyle state is another important experimental information available. Rather recently those at forward angle are also measured in Darmstadt [
51,
52]. Figure 5 shows the comparison between the data and the result calculated by the THSR wave function. We can see that the curve obtained by the THSR wave function is almost exactly on the experimental data, except for higher momentum region,
, in which the size of nucleon is probed and hence ordinary microscopic nuclear models cannot be applied.
The THSR wave function obtained by solving the Hill–Wheeler equation Eq. (21) is expressed as a sum, like in Eq. (20), of the original form of the THSR wave function with a single parameter in Eq. (19). The concept of the α condensation originates from the THSR wave function with a single parameter value, which we hereafter call a single THSR wave function. It is important to know how much this α-condensation picture kept by the single THSR wave function survives in the superposed wave function in Eq. (20). We then study how the superposed THSR wave function, which is shown to be equivalent to or better than the RGM/GCM wave function, can be approximated by the single THSR wave function.
We first calculate the following energy expectation value by the single THSR wave function in the two-parameter space , with the Force 2,
In Fig. 6(a), we show the energy contour map in the two-parameter space and for the above quantity. A minimum appears at a spherical position, . This corresponds to the ground state which has a compact density structure. The squared overlap between this single wave function and the ground state wave function obtained by solving the Hill–Wheeler equation, is calculated to be . This means that the single THSR wave function with a rather small value of the parameter very well approximates the ground state wave function , though the minimum energy MeV is about MeV higher than the energy MeV obtained by solving the Hill–Wheeler equation.
In this energy surface, however, there only appears one energy minimum, which corresponds to the ground state, and there is no other minimum corresponding to the Hoyle state. This is because the orthogonality to the ground state plays an important role in the formation of the Hoyle state. We then calculate the energy surface given in a configuration space orthogonal to the ground state, which can be defined below:
, is the projection operator onto the functional space which is orthogonal to the ground state. In Fig. 6(b), we show the contour map of this energy surface in the two parameter space . We see an energy minimum at in the oblate region of the map and a second energy minimum at in the prolate region. The minimum energy value is MeV. This value is almost the same as the energy of MeV given by solving the Hill–Wheeler equation. The minimum energy MeV is close to the second minimum energy of MeV, and there is a valley with an almost flat bottom connecting these two minima. This means that the energy gain due to the deformation is small. The wave functions at the two minima are almost identical to each other, and indeed, the squared overlap between the wave functions giving the two minima is calculated to be more than . The fact that the minimum energy is almost the same as the eigenenergy given by solving the Hill–Wheeler equation means that the single THSR wave function at the minimum energy position is also the same as the superposed THSR wave function given by solving the Hill–Wheeler equation.
The single THSR wave function corresponding to this minimum energy, i.e., corresponding to the Hoyle state, can be defined, in the orthogonal configuration space, as
We then show in Fig. 6(c) the contour map of the squared overlap surface defined below:
where is the Hoyle state wave function obtained by solving the Hill–Wheeler equation. Two maxima appear corresponding to the two minima in the previous figure. The largest value is at . This means that the Hoyle state wave function obtained by solving the Hill–Wheeler equation is almost identical to the single THSR wave function . Since we keep in mind that the wave function is practically better than the solutions of the RGM/GCM in Table 2 since the former gives deeper binding energy than those of the RGM/GCM, we can say that the Hoyle state is the best described by the single THSR wave function in the form of Eq. (26). This result gives us a crucial evidence that the Hoyle state is the condensate state such as described by the single THSR wave function in Eq. (26).
2.1.4 Quantum Monte Carlo calculation of the Hoyle state
One of the most remarkable recent progresses in nuclear theory is the development of
approach in light nuclei. The so-called lattice Quantum Monte Carlo (QMC) approach [
54,
55], which is based on chiral perturbation of an effective field theory, reasonably reproduces the low-lying spectrum of
including the Hoyle state [
56,
57], where the authors find an obtuse triangular configuration of the
α clusters for the Hoyle state. However this interpretation is quite a lot contradictory to the finding of many other theoretical investigations discussed up to now, that the Hoyle state has a gas-like configuration of loosely coupled
clusters in relative
S-wave.
On the other hand, the no-core symplectic model (NCSpM), which is based on a special symmetry to use only a fraction of the model space extended beyond current no-core shell-model limits, also has a great progress in the description of the cluster states including the Hoyle state [
58-
60]. By the help of its symmetry guided framework, which allows for selecting physically relevant configurations, the authors in Ref. [
59] adjust the energy positions of the Hoyle state and the
state in
, and discuss the physical observables of the states. The calculated r.m.s. radius of the Hoyle state is
, which is sizably larger than that of the ground state
. The monopole transition strength is also calculated, which is
, compared with the experimental value,
(see below and Table 2).
The pioneering work on the
ab-
initio theory has obviously been done by the Argonne group. In 2000, Wiringa
et al. published the outstanding result that the
α+
α cluster structure in
is verified [
61]. They take a variational approach, consisting of two steps, Variational Monte Carlo (VMC) and Green’s Function Monte Carlo (GFMC) with a constrained path algorithm, in which the Argonne
two-body force (AV18) and Urbana IX (UIX) three-body force, and eventually a new class of three-body potential, Illinois model (IL7), in
p-shell nuclei, are adopted as realistic nuclear force. Their calculations are now proved to work very well for low-lying states up to
, in spite of the fact that the computational cost increases exponentially with the number of particles [
15].
In particular, their result about the Hoyle state is quite amazing. We show in Fig. 7 the inelastic electron scattering form factor, where most of the experimental points are reproduced very accurately [
15] as much as the THSR wave function (see Fig. 5). In the inset the calculated
transition form factor is compared to the experimental data, which leads to
at
, given in [
51,
52]. We can see that the GFMC result gives excellent agreement with the data. We should keep in mind that the THSR wave function gives
and
in the case of adopting Volkov No. 1 and No. 2 forces, respectively (see Table 2). The calculated excitation energy of the Hoyle state is still rather high,
versus
MeV. In Fig. 8, we compare the density distributions of the Hoyle state and the ground state (weighted with
r2) obtained with the THSR wave function and in Ref. [
15]. For the ground state both curves are in quite good agreement with each other. For the Hoyle state both curves characteristically give a plateau region between 1 and 4 fm, and both are in good agreement, though in the GFMC calculation it is more pronounced around 1-2 fm than in the THSR ansatz. Beyond 4 fm, the density in the GFMC calculation falls off more rapidly, which is the reason why the r.m.s. radius of the Hoyle state is calculated to be smaller in the GFMC calculation than in the THSR ansatz. The former calculation gives around
fm and the latter gives 3.8 fm and 4.4 fm when Volkov No. 1 and No. 2 forces are adopted, respectively (see Table 2).
2.1.5 Breathing-like excitation of the Hoyle state
There has been known a broad
state above the Hoyle state, at
MeV with a width of
MeV. This state has triggered much interest of nuclear physicists, since there has been long controversy between experimentalists and theorists about the spin assignment of this state. Theorists predicted the existence of the second
state around 10 MeV almost 40 years ago by using the microscopic cluster models [
13,
14,
42]. Amazingly the
state has been observed recently, as a different state from the broad
MeV state. The pioneering work on the observation was originally done by Itoh
et al. [
63], and later it was confirmed in several experiments [
64-
68]. Furthermore, Itoh
et al. pointed out that the broad
peak is decomposed into two peaks, giving the
and
states at
MeV and
MeV above the
threshold, with the widths of
MeV and
MeV, respectively [
69]. Theoretically it is very difficult to obtain both
states several MeV above the Hoyle state, due to their broad widths. However, with proper treatments of resonances, some theoretical works succeeded in reproducing both states [
16-
18,
70,
71]. For example, the OCM plus complex scaling method (CSM) [
72-
76] properly gives two energy poles corresponding to the two
states [
70,
71]. Though the interpretation of these states is not yet converged, the result of this study suggests that the
state has a higher nodal structure between
in relative
S-wave [
70]. Also the AMD and FMD calculations suggest that the
state has a bent-armed shape of the
clusters, resembling the
linear-chain configuration, though in these calculations the
state is missing [
51,
77].
The THSR ansatz could also reproduce the
and
states by using a so-called radius constraint method [
78] that improves treatment of excited states in resonance region [
16-
18]. The translationally invariant THSR wave function discussed in this section is what is extended to contain the
correlation [
50], as follows:
where and with and , and is the angular-momentum projection operator onto . The degree of freedom of the cluster motions characterized by is divided into and , corresponding to () and α, respectively. Thus the case comes back to the original THSR wave function, and this natural extension allows for including the correlation and its asymptotic configuration. More detail of this extended THSR wave function is discussed in Section 5.
In Fig. 9, the result of solving the Hill–Wheeler equation of the type of Eq. (21), in which the four kinds of parameters,
and
are taken as the generator coordinates, denoted as GCM, is compared with the corresponding experimental data. The so-called radius constraint method [
78] mentioned above is also applied to obtain the
and
states. We can see that the
and
states are reasonably reproduced. Amazing feature of the
state is the large monopole transition to the Hoyle state, whose strength is around
. The r.m.s. radius of this state is calculated to be much larger than the Hoyle state. These results suggest that the
state is excited from the Hoyle state by monopole-like dilatation.
We should also emphasize that for the formation of the state, the orthogonality to the lower-lying states, i.e., the most compact ground state and much more dilute Hoyle state, is crucial. This can be demonstrated as follows: First we introduce the projection operator onto the orthogonal space to the ground state and the Hoyle state, like . Second we construct the THSR wave function in the orthogonal space, like . We then calculate the energy of this THSR state, which is defined as
In Fig. 10, the contour map of the energy surface in two-parameter space
and
is shown. In construction of the projection operator
, optimal single configurations in the THSR states are taken as the
and
states (see Ref. [
18] for more detail). We can see that two local minima appear, both of which, however, do not indicate two different modes of excitation, but are, in principle, identical state. In fact both minimum states have very large squared overlap of 84%. Furthermore the calculated squared overlap of these local minimum states with the wave function of the
state obtained by solving the Hill–Wheeler equation, which is shown in the spectra of Fig. 9, is more than 90%. This means that the
state appears as the minimum state, obviously indicating that the orthogonalization to the lower states is essential, since otherwise it never appears.
In order to investigate in more detail how the state is excited from the Hoyle state with the strong monopole strength, the following quantities are calculated:
where and , with and normalization factors, and and . and are the Jacobi coordinates. The upper equation is nothing but the RWA of channel, where Bx, By, and Bz in is best parameterized to describe nucleus. The lower equation is defined as representing the correlation of the clusters associated with the coordinate .
In Fig. 11(a), the RWAs in the channel for the – states are shown. The RWAs for the ground state and the state have two and four nodes, respectively. However, for the Hoyle state, the nodal behavior almost disappears and only a remnant of three nodes can be seen as an oscillatory behavior. Since the outermost nodal position corresponds to a radius of repulsive core between the core and the α cluster, due to the effect of the Pauli principle, the disappearance of the nodes for the Hoyle state indicates a dissolution of the core, and hence formation of the loosely coupled cluster state.
On the contrary, the
state recovers the distinct nodal behavior and, with one additional node, forms a higher nodal
structure. Obviously the relative motion between the
and
α is excited from the Hoyle state by the monopole transition. This is consistent with the argument in Ref. [
70], as mentioned above. On the other hand, in Fig. 11(b), the correlation function of the
part associated with
is shown, where the relative motion between the
and
α is chosen to best simulate the
state. In this calculation, this part is taken as
with
and
, since the
state obtained by solving the Hill–Wheeler equation has the largest squared overlap (around
%) with the single THSR wave function for this
value. While for the
state the two nodes in inner region remain, which come from the repulsive core between the
clusters, an oscillatory behavior is also clearly seen in outer region more than 4 fm. This indicates that the
relative motion associated with
is also excited by the monopole transition from the Hoyle state, although further analyses are necessary to clarify this point of view. This is, however, natural considering the fact that the monopole transition operator is written as
[
79]. The strong monopole transition strength from the Hoyle state can be from the two relative motions associated with the coordinates
and
, giving the picture that the
state is featured as monopole dilatation from the Hoyle state with a vibrational mode.
2.2 condensate-like state in
2.2.1 4α OCM calculation
We discuss in this section about the Hoyle-analog state in
calculated by the
OCM, in which the
problem is fully solved without any model assumption with respect to the
α+
α relative motions within a semi-microscopic treatment [
19]. This approach allows us to obtain a bosonic wave function of the
α clusters in an easier way than the fully microscopic approach like the RGM [
80]. The RGM equation of motion and the RGM wave function in the
-cluster system are written in the following way, respectively,
and
where , , and are the Jacobi coordinates between the α clusters and H is the microscopic Hamiltonian composed of 16 nucleons. Eq. (32) is rewritten in the following form,
with
where , , and denote the sets , , and , , respectively, which are all the Jacobi coordinates between the clusters. In the above equation, the dynamical coordinates of 16 nucleons except for the total c.o.m. coordinate, i.e., internal coordinates of the clusters and Jacobi coordinates between the clusters, are integrated out. Therefore, due to the delta functions in the bra and ket, the Jacobi coordinates between the α clusters, and in the bra and ket, respectively, only remain in and . This Eq. (35) can further be transformed into the following form:
where can be defined as a bosonic wave function of the clusters with only the Jacobi coordinates left, as follows:
In the OCM, the following approximation with respect to the treatment of the antisymmetrization in these kernels is taken,
where Λ is a projection operator onto Pauli allowed space [
40,
44] and
is the effective local Hamiltonian of the OCM written below:
where
Ti,
,
,
and
stand for the operators of kinetic energy for the
ith
α cluster, two-body, Coulomb, three-body and four-body forces between the
α clusters, respectively. The c.o.m. kinetic energy
is subtracted from the Hamiltonian.
is the Pauli exclusion operator, by which the Pauli forbidden states between the
clusters in
,
, and
states are eliminated, so that the ground state with the shell-model-like configuration can be described correctly. The two-body
α+
α force is constructed by folding procedure from effective nucleon-nucleon interaction, for which we adopt Modified–Hasegawa–Nagata [
81,
82] (MHN) force, which is known to reproduce the
α+
α scattering phase shift well. The three-body and four-body forces are phenomenologically introduced to reproduce the ground states of
and
at the same time. Readers can refer to Ref. [
19] for more details about the choice of the force parameters and refer to Refs. [
38,
40,
44,
80] for details about the OCM formalism.
In Fig. 12, we present the comparison of the experimental
spectrum with the calculated ones by using the
OCM [
19], together with the results calculated by using the THSR ansatz [
20], which will be discussed in the next section. For clarification we denote the six
states calculated by the
OCM and the four
states calculated by the
THSR ansatz as
and the
, respectively. We can see that the
OCM calculation gives a satisfactory one-to-one correspondence with the experimental spectrum. The low-lying states (
–
) of this spectrum are actually in good agreement with earlier
OCM calculations [
85,
86] (see also Ref. [
28]). We should emphasize that the orthogonality to the low-lying
states and the ground state is indispensable for the formation of the
-cluster gas-like states. This respect is particularly discussed in Section 5.1.
In Table 3 the comparison of the energy
E, r.m.s. radius
, monopole matrix element to the ground state
, and
α-decay width Γ between the experimental data and the results of the
OCM and THSR ansatz is shown. The
α-decay width is calculated based on the
R-matrix theory [
28,
87]. The
state has the largest r.m.s. radii among the calculated states and this large value
fm indicates that this state is the
gas-like state and considered to be the Hoyle-analog state. We can see that the observed widths of the experimental
,
, and
states are nicely reproduced by the
OCM calculations. In particular, the observed width of the
state at
MeV is very well reproduced by
state, and hence the Hoyle-analog state corresponds to the
state at
MeV.
We then investigate the Hoyle-analog state by calculating the RWAs for the channels, which are defined by
with and being the OCM states of and , respectively. We show in Fig. 13 the RWAs of the state. We can see that the Hoyle state component is the most largely included, which indicates that the state is the Hoyle-analog state because the condensate state must have a large overlap with the structure.
The α-condensate character of the is also clearly seen in the momentum distribution of the c.o.m. motion of α-clusters, which is defined by
where
is the one-body density matrix constructed by the bosonic OCM wave functions. How to calculate the matrix element and how to choose the internal coordinates as integral variables are in detail discussed in Refs. [
88,
89]. We show in Fig. 14(a) the momentum distributions of the
states. Only the one of the
state has a markedly sharp peak near zero momentum compared with the other
states. In Fig. 14(b), we show for comparison the momentum distribution of the Hoyle state by
OCM (dotted curve), which is also sharply peaked near zero momentum.
In addition to the
α condensate state, we also discussed the other
states. The structures of the
and
states are well established as having the
+alpha
and
+alpha
cluster structures, respectively. These structures of the
and the
states are confirmed in the present
OCM calculation. The present calculations also show that the
and
states mainly have
structure with higher nodal behavior and
structure, respectively [
28].
2.2.2 4α THSR calculation
The Hoyle-analog state in
is also investigated with the
THSR wave function. While in the
OCM discussed in the previous subsection the
problem is fully solved without any model assumption with respect to the
α+
α relative motions within the semi-microscopic treatment, in the
THSR the THSR-type constraint is imposed on the
α+
α relative motions. In this calculation, we only consider the spherical case, i.e.,
, within a fully microscopic treatment, and as a microscopic Hamiltonian, we adopt the nucleon-nucleon effective interaction named F1 force [
90], which includes a three-body term representing a density dependent force.
In Fig. 12, we presented the spectrum calculated with the THSR ansatz, denoted as the
–
[
20]. We can see that in the THSR ansatz the two excited states are missing. This can be understood by the fact that in the OCM the relative motions of
α particles are solved in a sufficiently large model space but the THSR ansatz does not possess nonzero spin waves with respect to the
α+
α relative motions. In particular, due to the OCM calculation, the
and
states are shown to have a large fraction of
S-factor from the
and
channels, respectively (see Ref. [
28]). These states can thus hardly be represented by the present spherical-type THSR wave function. In order to properly describe the
cluster states, it is necessary to further extend the THSR ansatz in Eq. (14) so as to separate the degree of freedom between the
core and the remaining
α cluster, by assigning different parameter sets of
to the first
’s and the fourth
α. This corresponds to the extended THSR wave function, which is in detail discussed in Section 5.
Thus in the THSR ansatz we consider the and states to be completely missing and the and states to be averaged representations of the and states, respectively. We can therefore say that the state corresponds to the state with the condensate character. In fact, in Table 3 we can see that the state has the largest r.m.s. radius of the four states, indicating that it has a dilute gas-like structure. The monopole matrix element and decay width well correspond to those calculated by the OCM.
We also calculate the RWAs for the channels in the THSR ansatz, which are defined by
with
and
being the wave functions of the
and
states, respectively, calculated with the (spherical) THSR ansatz in Eqs. (20) and (21). We show in Figs. 15(a)–(d) the amplitudes for the
, respectively, in the
and
channels. The ones of the
well correspond to the amplitudes for the
shown in Fig. 13. This again indicates that the
state is the counterpart of the
state and hence the Hoyle-analog state. We can also see that the
and
states have large amplitudes from the
channel and their peak positions are in outer region, beyond 5 fm, which are also consistent with the results of the
OCM calculation (see Ref. [
28]).
In the same way as the
OCM, we calculate the momentum distribution of the
α clusters in the THSR ansatz. The
THSR wave function is, however, the fermionic wave function of the 16 nucleons. We therefore have to extract a bosonic degree of freedom from
in Eq. (14), like in Eq. (37) for the RGM wave function. Since in this extraction harmonic oscillator expansion is necessary, which is, however, very difficult for dilute-density states in especially heavier system, the following approximation is taken in Ref. [
20],
with
where
can be determined by solving the Hill–Wheeler equation Eq. (21) for the spherical case. This approximation was first examined in Ref. [
91]. Thus we can calculate the momentum distributions of the c.o.m. motion of
α-clusters for the THSR ansatz, which is a microscopic treatment. In Fig. 16 those for the
states with
are shown. We can see that the
state has a sharp peak near zero momentum and its shape is in fairly good agreement with the one of the
state in Fig. 14, where both curves have a peak height around 2.0 fm and width around
.
3 Characters of the THSR wave function
3.1 Shell-model limit of the THSR wave function
The original THSR wave function [
8] is designed for the study of the
α condensation state. With in-depth study for this wave function, it is found that the THSR wave function can describe not only
gas-like cluster states but also the normal and compact cluster states [
11,
22]. One important reason why the THSR wave function can also describes the compact cluster states in nuclei is the THSR wave function has a shell-model limit when its size parameter
(or
). This is quite similar with the well-known shell-model limit of the Brink wave function when the inter-cluster distance approaches 0 as discussed in Section 2.
To illustrate the shell-model limit of the THSR wave function, let us take the THSR wave function as an example,
where
ξ1=
X2-
X1 and
ξ2=
X3-(
X1+
X2)/2.
and
c0 are normalization constants.
is the radial harmonic oscillator function with size parameter
with and
N stands for the number of harmonic oscillator quanta,
. The obtained Eq. (47) shows that the
THSR wave function at the limit of
becomes a shell-model wave function. This is why the THSR wave function not only describes well the gas-like Hoyle state but also gives a good description of the shell-model-like ground state of
[
9,
11].
Furthermore, if the composed clusters of a nuclear system are closed-shell nuclei like
α and
, the THSR wave function for this cluster system also has the simple shell-model limit. In the
+
α system, after angular-momentum projection with even
L, the THSR wave function can be written,
. The shell-model limit of this THSR wave function can be expressed as follows [
11]:
where
and
cL are normalization constants. The derived Eq. (50) shows that at the limit of
, the THSR wave function becomes the most important
shell-model wave function with spatial symmetry [
4] and SU(3) symmetry
[
92]. This point is very important for the description of shell-model-like states in
with the THSR wave function.
What is interesting is the proposed THSR wave function for the α condensation or gas-like states has a shell-model limit when the size parameter tends to zero and this character is more than expected. It should be noted that, on the other limit, that is , the correlations of clusters disappear and the THSR wave function becomes the simple product of the cluster wave functions.
3.2 Equivalence of projected prolate and oblate THSR wave functions
The THSR is a novel microscopic cluster model and it has some special characteristics, in particular, the equivalence of prolate and oblate THSR wave functions after angular-momentum projections. Figure 17 shows the squared overlaps between a prolate
THSR wave function with
THSR wave functions with various deformations in
[
11]. The deformed THSR wave function
of
has the form
[
(
α)
φ(
)], where
is
and
. We see in this figure that the prolate THSR wave function with
=
= 0.9 fm,
= 2.5 fm is almost 100% equivalent to oblate THSR wave functions with
fm,
fm after angular-momentum projections onto the
state. The equivalence of prolate and oblate THSR wave functions after angular-momentum projection is also true for other the spin-parity states of
.
In spite of the equivalence of prolate and oblate THSR wave functions after angular-momentum projection, we can say that the states expressed by THSR wave functions have all prolate deformation as the actual deformation. This conclusion is obtained from the fact that the expectation values of the quadrupole moments of all the states expressed by THSR wave functions have negative signs. From the well-known formula (intrinsic), we know that when the expectation value of the quadrupole moment by the wave function with good spin J is negative, the quadrupole moment of the intrinsic state Q (intrinsic) is positive which means that the deformation of the state is prolate. The THSR wave function after angular-momentum projection has the form and we can prove that this type of wave function gives us the following formula for ,
This formula shows that has negative value and explains why the calculated values of by THSR wave functions have all negative signs. Of course the negative sign of by THSR wave functions is in accordance with the prolate distribution of nucleon density shown in Fig. 26.
The reason why projected prolate and oblate THSR wave functions are almost equivalent is because the rotation-average of a prolate THSR wave function is almost equivalent to an oblate THSR wave function. The rotation-average wave function generated from a prolate THSR wave function is defined as follows:
If a prolate THSR wave function is rotated around an axis (x-axis) perpendicular to the symmetry axis (z-axis) of the prolate deformation and we can construct a wave function by taking an average of this rotation, the density distribution of the rotation-average wave function will be oblate (see Fig. 18). In the case of the state, when we construct the rotation-average wave function from the prolate THSR wave function with which gives the minimum energy for , it is almost 100% equivalent to the oblate THSR wave function with ,
When we construct the good-spin wave function from , it is clear that the resulting wave function is the same as the good-spin wave function constructed from . Therefore in the above example, the oblate THSR wave function is almost 100% equivalent to the prolate THSR wave function after angular momentum projection to the state.
3.3 Negative-parity THSR wave function
As we know, the original THSR has a positive-parity nature. In this subsection, we will show that the THSR wave function with a simple extension can be applied to study the negative-parity cluster states without any difficulty.
The first example for studying the negative-parity states in the framework of the THSR is the
case. We proposed a new type of microscopic cluster wave function [
10], which we call hybrid-Brink–THSR wave function,
Here
and
Xi and
are the center-of-mass coordinate and the internal wave function of the cluster
Ci, respectively. Different clusters
Ci can have different mass numbers
Ai and variational parameters
. The oscillator parameter of the cluster
Ci is called
b, which is taken as the same size value for different clusters.
is the corresponding general Brink model wave function [
41].
In Eq. (56), another generator coordinate Si is introduced to the original THSR wave function. It can be seen from Eq. (57) that this hybrid wave function combines the important characters of the Brink model and the THSR wave function in a very simple way. When Si=0, Eq. (56) corresponds to the THSR wave function and or Bi becomes the size parameter. When , i.e., , this equation is nothing more than the Brink wave function Eq. (58) and Si is the position parameter of the cluster Ci. Thus, due to the violation of the parity symmetry, we can easily projected out the negative-parity components from this hybrid wave function.
As we know, the THSR wave function is a nonlocalized clustering wave function for the cluster structure rather than the localized clustering represented by the Brink model [
48]. Since these two different kinds of pictures for clustering are both included in the hybrid-Brink–THSR wave function as two limits, this hybrid wave function provides a very nice way for verifying which picture is more adequate for understanding the relative motions of the cluster structures in nuclei. The details will be discussed in Section 4.
To deal with the negative-parity states of
9Be, another parity-violating THSR wave function was developed [
93], which is much more suitable for the calculations of Monte Carlo techniques. Details will be discussed in Section 5.5.1. In the constructed THSR wave function of
9Be, the spatial wave function of the extra neutron Φ
n(
rn) can be written as [
93]
When we change the integration variables to , in the integral representation of by Eq. (59), we obtain . It is because the azimuthal angle of is and . When , Eq. (59) is a standard THSR wave function and has a positive parity as already known. So the total parity of 9Be is now determined by m,
From Eq. (59) we can also see that the z component of orbital angular momentum is a good quantum number for the extra neutron. Because of the rotational symmetry of the two-α-cluster subsystem about z-axis, we have for α clusters. Thus the z-component of the orbital angular momentum lz of total system is
3.4 Single-particle property of the THSR wave function of 13C
The THSR wave function for the
13C system employed in Ref. [
94] is described as
with
where is the spherical THSR wave function for with in Eq. (14) and , with and being the coordinates of the valence neutron and the c.o.m. coordinate of the system, respectively. The coefficient of the superposition in Eq. (62), , and the corresponding energy eigenvalue can be obtained as a function of the -n relative distance parameter R by solving the following Hill–Wheeler-type equation,
The adiabatic energy curves with for the lowest and the second lowest eigenstates are shown in Fig. 19. It is noted that the cases of and correspond to the ground state and the Hoyle state for the core system, respectively, and hence demonstrate, respectively, that the valence neutron couples with the ground state and the Hoyle state in or . When the distance parameter R is large enough for the neutron to be far away from the core, we can see that the neutron with is energetically favored for both cases. On the contrary, as the S-wave neutron approaches the core, it starts to feel the repulsive force from the core, due to the Pauli blocking effect. This effect is considered to be specially important when the core is the ground state. In fact, in Fig. 19, the eigenenergy becomes higher than inside the surface region of the ground state of , i.e., for R less than around 3 fm. At the center of the core, the neutron sits in the orbit rather than in the orbit.
However, a more interesting result is for the case, where the core is in the Hoyle state. When the S-wave neutron comes closer to the core, it is also seen to feel repulsive core and its energy goes higher than that of the P-wave neutron, , at a larger value of R, according to the larger size for the Hoyle state. Again at the center of the core, the orbit for the valence neutron becomes higher than the orbit, though the energy spacing between them is much smaller than in the case. This result means that even for the Hoyle state, which is a dilute object with about of the normal density, the effect of the antisymmetrization is not still non-negligible, though the effect is much more smoothed than when the core is in the ground state.
4 Localized vs. nonlocalized clustering
4.1 Inversion-doublet bands in 20Ne and their description by the THSR wave function
4.1.1 Energy curve by the hybrid-Brink–THSR wave function
As one of the most typical cluster nuclei, the
α+
cluster system of
has been extensively studied by many nuclear models [
95-
99]. The observed ground-state band and
=
band in
were regarded as being an inversion doublet band arising from the hetero-polar di-nucleus configuration [
21] of
α+
cluster structure. See the schematic diagram in Fig. 20 This kind of inversion doublet band in
can be considered as the direct manifestation of cluster structures in
. Furthermore, based on the traditional understanding of the cluster structures, in asymmetry two-cluster systems, the inversion doublet band can also be regarded as a clear indication of the existence of the localized cluster structure together with the observation of large cluster decay widths. Actually, in heavier nuclei, many core+
α cluster structures [
100-
102] were identified based on the observed inversion doublet bands. In this situation, it seems that we have to regard states of the inversion doublet band of
20Ne as having a
α+
16O localized clustering. However, the study from the THSR wave function updates this kind of traditional understanding.
To study the inversion doublet bands in
, a hybrid-Brink–THSR wave function wave function was proposed as follows [
10]:
where
,
,
.
and
represent the center-of-mass coordinates of the
α cluster and the
cluster, respectively. All calculations are performed with restriction to axially symmetric deformation, that is,
and
. The spin and parity eigenfunctions can be obtained by the angular-momentum projection technique. As the nuclear interaction, we adopt the effective nucleon-nucleon force, Volkov No. 1 [
49] with the Majorana parameter
, and the oscillator parameter
fm.
Figure 21 shows the contour map of the energy surface of the intrinsic wave function of in the two-parameter space, Sz and . We can see the minimum energy MeV appears at and fm. The value means, after variational calculations, the obtained optimum intrinsic hybrid-Brink–THSR wave function becomes the pure THSR wave function in describing the ground state of . This result indicates that the intrinsic THSR wave function based on the concept of nonlocalized clustering is more suitable for describing the ground state of than the Brink wave function.
After angular-momentum projection on the intrinsic hybrid-Brink–THSR wave function with the spherical case , the optimal projected wave function can be obtained by variational calculations. Figures 22 and 23 show the contour maps for the , positive-parity states and , negative-parity states of in the two-parameter space, Sz and , respectively. It is surprising to find that the obtained minimum energies with respect to the projected states all appear at . For instance, the obtained minimum energy MeV for state appears at and fm. For the state, the minimum energy MeV appears at and fm in the contour map. For the other states of inversion doublet bands of , they have the same character after variational calculations. The inter-cluster distance parameter means that this hybrid-Brink–THSR wave function tends to become a pure THSR wave function in describing the cluster states of the inversion doublet bands in . One should note that although Sz does not give any contribution to the energy gain, it still plays an important role in providing negative-parity states. This is because the parameter Sz is the only component which breaks the parity symmetry, as is clearly seen in the form of the intrinsic wave function in Eq. (65).
The above variational calculations with the hybrid-Brink–THSR wave function lead the parameter Sz to zero for the inversion doublet bands. This means that, in spite of the fact that a pure Brink wave function gives a distinct energy minimum point with non-zero Sz, the localized clustering picture cannot be supported. We can realize this situation by looking at Fig. 24, which shows the energy curves of the lower excited states of the inversion doublet bands with different widths of the Gaussian relative wave functions in the hybrid model. If is fixed at 0, the hybrid-Brink–THSR wave function becomes the Brink wave function. In this case, Sz is the inter-cluster distance parameter and it is usually regarded as a dynamics parameter for describing the cluster system. For instance, the minimum energy of the ground state of appears at fm. For the state, the optimum position appears at fm. The non-zero values of Sz seem to indicate that the α+ structure of favours localized clustering. This is just the traditional concept of localized clustering. Now, we believe that this argument is misleading. The non-zero minimum point Sz simply occurs since the width of the Gaussian wave function of the relative motion in the Brink model is fixed to a narrow wave packet, characterized by the parameter b. If we take non-zero values for , namely, = fm, 1.5 fm, 2.4 fm, and 1.9 fm for , , , and states, respectively, according to their minimum positions in the contour maps, then we find that the minima appear at in Fig. 24. This indicates that the separation distance parameter Sz does not play any physical role in describing the α+ cluster structure, even for the negative-parity states. Instead of that, the new parameterization by , which characterizes nonlocalized clustering, is more appropriate for describing the cluster structure in .
4.1.2 Equivalence between the 16O+α THSR wave function and 16O+α RGM/GCM wave functions
Now that the parameter Sz tends to be zero in the obtained hybrid wave functions for the inversion doublet bands in , we can make further variational calculations in the two-parameter space, and using the projected hybrid-Brink–THSR wave function with the parameter (In practical calculations Sz is fixed at a very small value close to zero for negative-parity states). The obtained minimum energies and the corresponding values of and using the THSR-type wave functions are listed in Table 4.
On the other hand, the exact solution of the α+ cluster system can be obtained by superposing the Brink wave functions, that is the Brink-GCM wave function,
Here, can be obtained directly from the projected Brink wave function with . By solving the Hill–Wheeler equation Eq. (66), we can obtain the following Brink-GCM wave function,
Thus, we can compare the single THSR-type wave function with the exact Brink-GCM wave function by calculating their squared overlap . In Table 4, we can find that the obtained single THSR-type wave functions have 99.29%, 98.79%, 97.75%, 99.98%, and 99.87% squared overlaps for =, , , , and states of , respectively, with the corresponding Brink-GCM solutions. These very high squared overlaps mean that the single THSR-type wave functions are almost 100% equivalent to the corresponding RGM/GCM wave functions, thus, these obtained single angular-momentum projected THSR-type wave functions can accurately describe the states of the inversion doublet bands in . Moreover, the concept of nonlocalized clustering indicated by the THSR-type wave function obtained from the hybrid-Brink–THSR wave function is essential to correctly understand the cluster structure in .
Since we have shown that the single THSR wave functions are almost equivalent to the GCM wave function, it is natural that their relative wave functions, or any physical quantities calculated by the wave functions should be also equivalent. In Ref. [
103], the author analyzed the
α-cluster wave functions in cluster states of
and
by comparing the exact relative wave function obtained by GCM wave function with various types of trial functions. The adopted trial functions are the Brink wave function (BB), the spherical Gaussian THSR wave function (sTHSR), the deformed Gaussian THSR wave function (dTHSR), and a little more mathematical Yukawa tail (YT) function. By comparing with these trial functions and the GCM wave function, it was confirmed again that the THSR wave function can give a very good description of the positive-parity states and negative-parity states in
. To be more clear for the behaviors of the localized wave function, nonlocalized wave function, and the exact GCM wave function, we can focus on their relative wave functions, which are essential parts for the description of the motion of clusters. Figure 25 shows the relative wave functions of the optimized trial wave functions and the corresponding GCM wave function, it can be seen clearly whether it is the compact shell-model-like ground state of
or the developed negative-parity state, the THSR wave function characterized the nonlocalized clustering is much better than the Brink wave function. Furthermore, in the tail parts of the relative wave functions, the single deformed THSR wave functions can almost give exact description while this is only can be done by superposed many localized Brink wave functions in GCM calculations. The evidence presented thus far supports the idea that the nonlocalized clustering should be one essential feature of cluster motion.
4.2 Localized clustering from the inter-cluster Pauli principle
Now we explain how the idea of the parity-violating deformation of localized 16O+α clustering for the inversion-doublet bands of 20Ne can be justified within the THSR wave function which assumes nonlocalized clusters. The parity-violating deformation is a property of the intrinsic state which is the instantaneous (or adiabatic) quantum state of the rotation of the nucleus. Since the instantaneous configuration of two clusters is of prolate shape, the prolate THSR wave function is the intrinsic state of the system and the oblate THSR wave function is not the intrinsic state but rather a mathematical object which expresses the rotation-average of the intrinsic state as discussed in Section 3.2. The spherical THSR wave function expresses the time average of the fully three-dimensional rotational motion, namely the angular-momentum projected state of the intrinsic state (the prolate THSR wave function). We, however, also need to notice the fact that two clusters can not come close to each other because, as just mentioned, of the inter-cluster Pauli repulsion, which implies that two clusters in the intrinsic state (the prolate THSR wave function) are effectively localized in space. Thus, the prolate THSR wave function has the parity-violating deformation of localized 16O+α clustering. We can say that dynamics prefers nonlocalized clustering but kinematics makes the system look like localized clustering. Figure 26 shows density distribution of the 16O+α hybrid-Brink–THSR intrinsic wave function and this kind of localized clustering appears. Of course, this localization is most pronounced in the necessarily strongly prolate two cluster systems. In systems with low density α clusters in number more than two have more space to move independently and are, therefore, less localized in spherical containers.
The effective localization of clusters in the prolate THSR wave function of the two-cluster system is clearly seen in the density distribution of the prolate THSR wave function. In Fig. 27, it is shown the density distribution of a 2α prolate THSR wave function with . Since the THSR wave function before the antisymmetrization operation is obviously composed of nonlocalized clusters, it is evident that the clear spatial localization of clusters shown in this figure is attributed to the inter-cluster Pauli principle.
Recently it has been reported [
104] that the density distribution of a 3
α THSR wave function with strong prolate deformation with
shows clear spatial localization of three
α clusters aligned linearly, which is displayed in Fig. 29 in Section 5.2. It is to be noted that because of the almost zero values of
, three
α clusters are not allowed to expand into the
x and
y directions, which means that three
α clusters are only allowed to make one-dimensional motion along the
z direction. Therefore the inter-cluster Pauli principle acts only along the
z direction, which is the reason of the spatial localization of the three
α clusters. In Ref. [
104] it is reported that the
α linear-chain Brink-GCM wave function is almost 100% equivalent to a single 3
α THSR wave function with strong prolate deformation which is just the 3
α THSR wave function in Fig. 29 with
.
All these results provide important insights into the concept of nonlocalized clustering. In two-cluster systems, cluster states generally have effective localization of clusters because of the inter-cluster Pauli repulsion. However, in three or more cluster systems, the spatial arrangement of clusters are not necessarily geometrical, namely clusters can be nonlocalized, although the inter-cluster separations are non-zero simultaneously because of the inter-cluster Pauli repulsion. However, as is discussed in Ref. [
104], if a cluster state is forced to have strongly-prolate deformation, the state can have effective localization of clusters like in the case of 3
α linear-chain structure. When the inter-cluster separations are large, the spatial arrangement of clusters can be non-rigid and gas-like. We will discuss this case in the following section.
5 Container model of cluster dynamics
In the previous sections, we showed that a single configuration of the THSR wave function does not only describe very accurately gas-like cluster states like the Hoyle state in but also non-gaslike localized cluster structure states like the inversion doublet band states in . This means that the THSR ansatz, which was originally introduced to represent the α condensates with a dilute density, can also be applied to more general cluster states. As was discussed in Section 4, the important ingredient of the THSR wave function is to be characterized by a nuclear size parameter, . This is completely different from conventional cluster model descriptions like the Brink wave function, including the AMD and FMD models, in which the relative distance between clusters, or between nucleons, are the key parameters. Thus the THSR model does not provide geometrical configurations of clusters but is closer to a cluster-mean-field approach, where constituent clusters make independent motions in the cluster-mean-field self-consistent potential, under the constraint of the antisymmetrization of nucleons. Furthermore in this model, a shape and size of the nucleus (characterized by ) can be varied flexibly, according to dynamical excitation of the nucleus. This model therefore gives the picture that the individual clusters move rather freely in a shape- and size-changeable container. In this section, we investigate the container dynamics of clusters with respect to its size (Section 5.1) and shape (Section 5.2), rich spectra of 3α (Section 5.3) and (Section 5.4), and discuss more extension and application of the THSR ansatz following this picture (Section 5.5).
5.1 Nuclear size as the generator coordinate for the evolution of cluster structure and extension of the THSR wave function
In Section 2.2, we argued that the , , and states in the THSR ansatz correspond to the , , and states in the OCM calculation, which have the structures of and α in an S-wave, of and α in a higher nodal S-wave, and of the condensate, respectively. The THSR wave functions of these states, with of Eq. (20), are obtained by solving the Hill–Wheeler equation Eq. (21), with a constraint of , which are represented as follows: . In this subsection we clarify that the cluster excitation given by solving the Hill–Wheeler equation are associated with an evolution of the size of container, which is characterized by .
We first construct orthonormal wave functions from the THSR wave functions as follows:
where are normalization constants and are projection operators defined by
Here the case gives . The wave function depends on the parameter and is orthogonal to the wave functions of the states satisfying . We then calculate the following squared overlap amplitudes between the wave functions and ,
In Fig. 28 we show of Eq. (70) as a function of . For any wave functions , take values close to 100% when is optimally chosen. This result means that all the states are almost represented by single THSR wave functions parametrized by single values in the spaces orthogonal to the lower-lying states. We discussed this point of view in Section 2.1.3 (see Fig. 6), to show that the Hoyle state can be represented by a single THSR configuration with an optimal parameter value, and hence can be regarded as the condensate state. In particular, in this case for the state the maximum of amounts to 96% at fm, for which this large value of again supports the picture that this state is considered to be the condensate state, in which the clusters are trapped into a very large (spherical) container characterized by the nuclear size parameter .
However, this figure shows that not only the dilute gas-like state but also the lower states with , in which cluster dynamics prevails, are expressed as single THSR wave functions in the proper orthogonal spaces. This corresponds to the case of discussed in the previous section, Section 4, in which the inversion doublet band states are all expressed with almost 100% precision by single THSR configurations. We further show in Table 5 the largest squared overlaps and the corresponding values for the four states. We can see that the optimal values increase from the to the states. The r.m.s. radii are also proportional to the parameter values. Therefore we can say that the clustering excitations, which are usually described by solving the Hill–Wheeler equation, are very accurately characterized by the variation of the value, i.e., dilatation of the container, from the most compact state to the most dilute gas-like state with the condensate structure, through the intermediate density and more dilute states with the cluster structures. The evolution of the cluster excitations can be described by an evolution of the container, from a smaller value of to a larger one. This also supports the correctness of the container picture.
We should also note that the projection operator , which removes the lower-lying-state components with more compact structure, plays a role as a repulsive force to prevent the α clusters from being resolved or from coming closer to form the more compact structure, due to its orthogonality condition. As a result, higher-lying states are constructed one by one. This is essentially the same as the situation of , where the antisymmetrization operator removes the Pauli forbidden states, to construct a structural repulsive core between the two α clusters.
The more sophisticated treatment of configuration in the container picture is to use an extended version of the THSR wave function, which we will discuss in Section 5.4.
5.2 One-dimensional gas of α clusters – α-linear-chain structure
In this subsection, we consider the so-called
α-linear-chain state, which has been studied for many years [
35,
105] and it is one of the most typical representatives of the localized cluster-structure states. A special interest is whether the
α-linear-chain states, which are ordinary considered as the
α clusters geometrically arranged in a line, can be reinterpreted by the container picture provided by the THSR wave function, as was clarified in the case of the
in Section 4 and even of lower-lying states of
in Section 5.1. The ideal situation of the
linear-chain structure states are considered to be obtained by the Brink-GCM wave function, which was discussed repeatedly in the precedent sections but is in this case considered in only one dimension. If that is
z-axis, the (ideal) linear-chain states (LCS) can be obtained below:
where
is the Brink wave function in Eq. (16) with
, and the coefficients
are determined by solving the Hill–Wheeler equation like Eq. (21) taking the variables,
, as the generator coordinates instead of
parameter in this equation. In Ref. [
104], we calculated the
and
LCS in this way and efficiently obtained convergent results by superposing the basis functions in Eq. (71) in a stochastic way, whose numbers amount to 100 and 300 for the
and
cases, respectively. Then we calculated the squared overlap of these LCS wave functions and the single configuration of the THSR wave function, i.e., the quantity defined as:
. The detailed parameter set chosen in this calculation can be found in Ref. [
104]. Most amazingly the squared overlap values are found to be very large, more than 98%, and 93%, for all
,
, and
states in
and
nuclei, respectively. The optimal
parameter values in this situation all indicate extremely-prolate deformation, namely,
and
.
In Table 6, we show the maximum squared overlap values and the corresponding values (denoted as and LCS). These large overlap values are to be compared with the maximum square overlap values with and single Brink wave functions, instead of the single THSR wave functions, which are only 78% and 48%, respectively. Therefore for the LCS, the localized α-cluster picture aligned in a line, which is provided by the Brink wave function, is less appropriate while the one-dimensional container picture, newly provided by the THSR ansatz, which gives much larger values close to almost 100%, is much more appropriate to understand the LCS. More precisely, the LCS can be represented by a container with extremely prolate shape, in which it seems the α clusters again move rather freely without mutual dynamical correlation.
Nevertheless, just like the case of inversion doublet, the nucleon-density distributions of the intrinsic THSR wave functions of the LCS show the localized α clusters, as shown in Fig. 29 for the and cases. This is because of inter-α Pauli repulsion as discussed in the previous section, Section 4. Therefore, this again indicates that the localization of nuclear clustering comes from the kinematic effect, namely the effect of the antisymmetrization, since if we believe the present results, there is no way to make α clusters localized except for the operator , considering the form of the single THSR wave function of Eq. (18), which essentially represents a free motion of the α clusters in a container with a certain shape and size characterized by parameter.
We also summarize in Table 6 the squared overlap values between the single THSR configurations and the corresponding Brink-GCM wave functions in
[
45,
48],
[
9,
48,
50],
[[
10,
11,
22], and
[
106], some of which are already discussed up to now. All these values are very large and close to 100%. These results suggest that the single THSR wave functions can describe gas-like and non-gaslike cluster states in a unified way, allowing us to consider that cluster states in general retain a container structure, where the container corresponds to the cluster-mean-field potential characterized by the size parameter
in the THSR wave function. The constituent clusters are confined into the container in a nonlocalized way and occupy an identical orbit, under the effect of the antisymmetrization. This is very different from the conventional cluster model wave function, like the Brink model wave function [
41], in which relative motions of clusters are characterized by inter-cluster separation distance parameters, in a localized way.
5.3 Rich spectra of 3α states above the Hoyle state in 12C
5.3.1 Positive-parity states in 12C
In Section 2.1.5, the nature of the state in was discussed, which was identified as the monopole excitation from the Hoyle state. In this subsection, we argue that there exist besides it a number of other α cluster states above the Hoyle state, which we can qualify as the excited states of the Hoyle state. It is then indispensable to generalize the THSR ansatz, particularly to take into account -type cluster structure. On the other hand, in Section 4, we showed that the single THSR wave functions excellently describe the cluster states in , and they are much better than the single Brink wave functions, whose variational parameter is the inter-cluster separation distance. This result suggests that using the parameter of the THSR ansatz is much better than the Brink-type parameter choice to describe the structure in , or even structure in .
We thus introduced in Section 2.1.5 the extended version of the THSR wave function Eq. (28), in which the two different width parameters and are associated to the two Jacobi coordinates and , corresponding to α+α and -α motions, respectively. We show it again below,
This generalization is a natural extension of the original THSR wave function, since taking
results in it. This also follows the container concept that is discussed in the previous sections, since the case of
corresponds to the
cluster structure, in which the
core and the remaining
α cluster are confined in a small-size container, while in the core
the
clusters are confined in a large-size container. The shapes and sizes of the two containers are then characterized by the corresponding two parameters,
and
. We should note that the fact that the THSR wave function includes the
asymptotic feature for
allows for the sophisticated approaches to resonances like the CSM [
72-
76] and the analytic continuation of the coupling constant (ACCC) method [
107].
For the excited states above the
threshold, it is well known that the application of the bound state approximation gives accidental mixing between spurious continuum states and resonances. By using the fact that the root mean square (r.m.s.) radii of spurious continuum states are calculated to be extremely large within the bound state approximation, we developed a new method to remove the spurious continuum components [
78]. First we diagonalize the operator of mean square radius as follows:
We then remove out of the present model space the eigenstates belonging to unphysically large eigenvalues. By taking the following bases,
with γ satisfying , we diagonalize Hamiltonian as follows:
Since the present extended THSR wave function can include
asymptotic form by taking the large values of the two width parameters
and
, the
continuum components, as well as the
continuum components, can be successfully removed by imposing the cutoff for the mean square radius
. Although we could not obtain the excited states except for the
and
states by using the original THSR wave function [
47], we can now obtain the other observed
,
, and
states by using the present extended THSR wave function with a treatment of resonances. This method is also used in Section 2.1.5 to obtain the
state and is there referred to as the radius constraint method.
With this type of generalized THSR wave function, one can get a much richer spectrum of
. In solving the Hill–Wheeler equation with the radius constraint method Eqs. (73)–(75), axial symmetry is assumed and the four
β parameters are taken as generator coordinates. In Fig. 30(a), the calculated energy spectrum is shown. One can see that besides the ground state band, there are many excited states obtained above the Hoyle state. All these states turn out to have large rms radii (3.7–4.7 fm), and therefore can be considered as excitation states of the Hoyle state. The Hoyle state can, thus be considered as the “ground state” of a new class of excited states in
. In particular, the nature of the series of states (
) and the
and
states have recently been much discussed from the experimental side. As was mentioned in Section 2.1.5, the
state that theoretically has been predicted at a few MeV above the Hoyle state already in the early works of
Brink-GCM [
13,
42] and
RGM [
14] was recently confirmed by several experiments [
64-
68]. A strong candidate for a member of the Hoyle family of states with
was also reported by Freer
et al. [
108]. Itoh
et al. recently pointed out that the broad
resonance at
MeV should be decomposed into two states:
and
[
69]. This finding is consistent with theoretical predictions where the
state is considered as a breathing excitation of the Hoyle state [
16-
18,
70] and the
state as the bent arm or linear chain configuration [
17,
51,
77].
In Fig. 30(a), the
transition strengths between
J and
states and monopole transitions between
states are also shown with corresponding arrows. We can note the very strong
transitions inside the Hoyle band,
and
. The transition between the
and
states is also very large,
. In Fig. 30(b), the calculated energy levels are plotted as a function of
, together with the experimental data. There have been attempts to interpret this as a rotational band of a spinning triangle as this was successfully done for the ground state band [
109,
110]. However, the situation may not be as straightforward as it seems. This is because the two transitions
and
are of similar magnitude, and hence no clear band head can be identified. It was also pointed out in Refs. [
51,
77] that the states
form a rotational band not with the
but with the
state. The line which connects the two other hypothetical members of the rotational band, in Fig. 30(b), has a slope, which points to somewhere in between of the
and
states. The similar effect is also argued in the study of the
condensate and
rotational band in
[
111,
112]. Furthermore, we show in Fig. 31 the
S2 factors of the
components, which are defined below:
where
is the RWAs of the
states in the
channel. Except for the
state, all the states have dominant components from the channels
, which is consistent with
and
α rotation. However, In the
and
states, the components from the channels [0, 2]
2 and [0,
4]
4, respectively, are dominant. However, in the
and
state, the mixtures of the other components are also found to be large. These suggest that to conclude from there this gives rise to a simple rotational band, is premature [
16].
Finally we briefly discuss the structure of the
state in
. In the AMD and FMD calculations [
51,
77], the dominant intrinsic configuration of the
state is a bent-armed structure of the
clusters, resembling the linear-chain structure. However, any firm conclusions cannot still be drawn, since in the AMD, FMD, or even the
Brink-GCM calculation, the four
states are not consistently reproduced. As we mentioned above, in these calculations, the
state is missing and only the
state was obtained. Indeed, it is reported that the orthogonality to the lower states plays an essential role in the emergence of the linear-chain structure in the
state [
113]. On the other hand, in Fig. 31 the state is shown to have the largest component from [2, 2]
0 channel, i.e., the
α cluster orbiting in a
D-wave in the
core. This is consistent with the former calculations in Refs. [
16,
70]. We then examine whether the interpretation of the linear-chain-like structure is feasible for the
state.
In Fig. 32, we show the contour map of the squared overlap between the
state and the single THSR configuration in two-parameter space,
. The contour lines giving more than 0.4 are denoted by dot in red, in a step of 0.01. This contour map has a characteristic feature, where the strongly prolate deformation is only allowed to have a non-negligible squared overlap amplitude. Except for this prolate deformed region, the squared overlap is less than 0.1. The largest value is 0.47, which is not so much large but clearly indicates the
linear-chain structure. The parameter values giving the maximal squared overlap is calculated to be
, which is close to
that was obtained in a rather ideal one-dimension situation in Ref. [
104]. In Ref. [
104], it is discussed that largely prolate-deformed THSR wave function shows one-dimensional
α condensate of
clusters, which is fairly different from the ordinary picture of the linear-chain state with rigid-body
-cluster configuration arranged in a line in a spatially localized way. Thus we can say that the present
state has the one-dimensional
α condensate structure by around 50%, where the
clusters are loosely trapped into a prolate-deformed potential like a one-dimensional gas.
5.3.2 Container picture for the 3- and 4- states in 12C
As we know, the negative-parity states in
also have many well-developed cluster states [
13,
14,
114]. Quite recently, the
state was reported to have the enhanced dipole transition strengths from the ground state and it probably has the dilute cluster structure [
115]. On the other hand, the negative-parity states of
are closely related with the intrinsic structure of 3
α clusters. Bijker
et al. [
109,
116] proposed a description of cluster states in nuclei in terms of representations of unitary algebras. The
was assumed to have
3
α cluster structure and many cluster states
were predicted. Particularly, the
state at 22.4(2) MeV in
was observed from experiments [
110]. Within our microscopic container model, the nonlocalized motion is the essential feature of 3
α and it seems rigid geometrical arrangements of
were not favored, especially the Hoyle state and other developed positive-parity cluster states. The positive-parity states of
have been well described within the THSR wave function as discussed in Section 5.3.1. As a first step to study the negative-parity states in
, we focus on the 3
- and 4
- states in
in this subsection.
To deal with the negative-parity states, we can introduce the shift parameter to violate the parity of the THSR wave function, which has been discussed for dealing with the negative-parity states of two-cluster structure in in Section 4. We start from the general two-β THSR wave function,
where , . , . and . The ,, and are the center-of-mass coordinates of the three clusters in . In the above equation, is the 3α Brink wave function with the corresponding generator coordinates and . The represents the specified distance parameter of two clusters and these two clusters and the third α cluster. The and are introduced shift parameters.
At present, in order to simplify calculations and facilitate the analysis, we adopt the one-deformed THSR wave function by letting
and
in the above equations. And the introduced shift parameter
and
can arrange arbitrary triangle shapes on
xy-plane. Also, we have known that these shift parameters play an important role for separating the 2
α clusters and a third
α cluster in the wave function, which can be seen in Fig. 33. Thus, by the angular momentum and parity projection technique, we can obtain the negative-parity components from the extended THSR wave function. As for the Hamiltonian, the Volkov No. 2 [
49] (modified version) with Majorana parameter
M=0.59 and harmonic-oscillator size parameter
b=1.35 fm, which was used by Kamimura et al. for 3
α RGM calculation [
14].
The first 3- and 4- states in 12C have been known to have the developed cluster structure for a long time. In the one-β container picture, we performed the variational calculations after the angular-momentum projections for the 3- and 4- states in two-parameter βx=βy and βz space. The shift parameter is adopted as S1=(1/2, 0, 0), .
Figure 34 shows the variational result for the 3
- state. It can be seen that [
117] there are two local minimum points appeared in the contour plot. There is a long and narrow valley connecting the two points. The deeper local minimum point appear at
βx=
βy=1.5 fm and
βz=3 fm. The obtained energy for the 3
- state is about -80.9 MeV. Figure 35 shows the contour plot for the 4
- state. We found two local minimum points and the deeper one appears at
βx=
βy=1.0 fm and
βz=3.4 fm. The obtained energy for the 4
- state is -76.8 MeV. It should be noted that the THSR wave functions with respect to the two local minima for 3
- and 4
- states are also quite similar, whose squared overlap values are about 98%. This means the equivalence of the oblate and prolate THSR wave functions after angular-momentum projections, which is similar to the case of
20Ne in Section 4.1.2.
Following the study of 20Ne system in Section 4, we further compare the results obtained by the THSR wave functions with those from the Brink cluster model shown in Fig. 36. Firstly, we assume the equilateral triangle shape with the length of the side is S for 3α clusters in the Brink model. At the same time, S is also the shift distance on xy plane in our extended THSR wave function. It can be seen that, in the Brink cluster mode, there are two distinct pockets around 3 fm for the 3- and 4- states. It seems that this is a support for the rigid geometrical structure from the inter-cluster distance parameter S. However, if we consider the width variables of the relative wave function, that is the spirit of the container picture, the inter-cluster distance parameter is becoming a shift parameter instead of the dynamic variable. It can be seen that the minimum energy points appeared around S≈0. Furthermore, the obtained energies in the container picture for 3- and 4- states are much more deeper than that obtained from the Brink model. And actually these energies are very closed to the GCM solutions. Moreover, both energy curves are very flat especially in the range S<2 fm. The obtained energies are almost degenerate without depending the inter-cluster distance parameter S.
Figure 34 and Fig. 35 show that the 3- and 4- states can be considered to have quite a similar intrinsic oblate shape, i.e., βx=βy=2.0 fm and βz=0.5 fm. By fixing the size parameter in the intrinsic wave function, various triangular shapes can be constructed taking the shift parameters S1 and S2 on xy plane. It should be noted that this oblate wave function (βz→0) actually can be treated as one two-dimensional container. As we have discussed the effective localized clustering in 20Ne in Section 4 and one-dimensional linear chain localized clustering in Section 5.2, this kind of localized clustering mainly originates from the Pauli principle. How about the three clusters released from one-dimensional space in the container picture? Figure 37 shows the densities of various intrinsic wave functions with different values of shift parameters. It can be seen that in Fig. 37(a) the length of the shift triangle is too small to violate the positive-parity character of the THSR wave function. To break the symmetry, we take different lengths and orientations of the shift parameters as shown in Figs. 37(b)–(d). It can be found that, in spite of various triangle geometrical shapes, their corresponding single projected wave functions are quite similar and they are more than 92% equivalent to the corresponding 3- and 4- GCM wave functions. It seems that the α clusters can almost move freely without much influence from Pauli principle in our two-dimensional container and they show the real nonlocalized clustering feature.
5.4 Variety of cluster states in 16O
In this subsection, we discuss various cluster states in
16O, such as
α+
12C and 4
α gas-like cluster states, which are obtained via the extended THSR ansatz. The special interest is in how those various cluster states are formed with the increase of excitation energy. As we will discuss in detail in Section 6, the
α+
12C cluster structure in
16O is formed by the activation of cluster degree of freedom in the ground state having a dual property [
28,
79], i.e., by the excitation of relative motion between the
α and
12C clusters. The gas-like 4
α cluster state is then produced as a result of further excitation of the
12C core, to the 3
α cluster state, i.e., to the Hoyle state. The path of this cluster evolution is shown in the Ikeda diagram, together with many other paths in many other nuclei (see Section 6). We show in Fig. 38 a schematic picture for the path of cluster evolution along the excitation energy in
16O, which can be regarded as the path of a size and shape evolution of a container.
In order to describe the above mentioned cluster evolution in 16O, the following extended THSR wave function with a double-container (eTHSR) is employed. This is an extension of Eq. (28), which is adopted for the 3α system in the previous subsection, to the 4α system, and is expressed in the following form:
where ξi is again the Jacobi coordinates between the α particles, and , for i=1, 2, 3. While the parameter b characterizes the size of the constituent α particle, the parameters β1 and β2 characterize the size and shape of a container, in which the α clusters are confined. We also assume the axial symmetry , so as to deal with the four parameters, , in the practical calculations.
A schematic picture representing the eTHSR wave function in the present 16O system is shown in Fig. 39, in which the 3α clusters and another α cluster are confined in different containers characterized by the parameters β1 and β2, respectively. This is contrasted with the Brink wave function, in which the cluster configurations are described by their relative distance parameters.
In Fig. 40, the calculated energy spectrum for
states is shown. The corresponding experimental data and results by the 4
α OCM calculation [
19] discussed in Section 2.2.1 are also shown. The solution of Hill–Wheeler equation with the
r2 constraint method is shown. The
state is actually the seventh 0
+ state obtained by solving the Hill–Wheeler equation, i.e., the fifth and sixth eigenstates are kicked out from the present consideration, since they have larger r.m.s. radii and are regarded as spurious continuum states accidentally mixed with the physical states.
As we discussed in Section 2.2.1, the 4α OCM calculation gives six 0+ states. The state has the 4α condensate character and the - states all have α+12C cluster structures, i.e., , , , and cluster structures, respectively. The difference between the and states are that in the latter the α and 12C relative motion is further excited and has a higher nodal S-wave, to have a larger r.m.s. radius than the former.
Since in the present eTHSR wave function of Eq. (79) the
α clusters occupy positive parity orbits, such a state as having the
cluster structure, like the
state in the OCM calculation, is missing. The treatment of negative parity orbit in the THSR ansatz has been discussed in the subsequent Section 3.3. However, for the other states, a one-to-one correspondence to the experimental data as well as to the 4
α OCM calculation is consistently obtained. It should be noted that in the OCM calculations [
19,
86], the binding energies of the ground states of
16O and the
12C are phenomenologically fitted to the corresponding experimental values. In the present calculation, however, there is no adjustable parameter in the microscopic Hamiltonian.
In Table 8, r.m.s. radii and monopole matrix elements with the ground state are shown. The experimental data available are reasonably reproduced. We can also see that from the to the states, i.e., as the states are excited, the r.m.s. radius becomes larger and the monopole matrix element becomes smaller. This indicates that the higher the excitation energy is, the more evolved the clustering is. The evolution of the clustering can be described by solving the Hill–Wheeler equation concerning the model parameters β1 and β2.
This respect is made much clearer by calculating the following squared overlap:
with the λth eigenfunction and normalized single eTHSR wave function in a space orthogonal to the lower eigenstates, i.e.,
where with , and , and are the corresponding normalization constants.
This quantity indicates how these five states are expressed by single configurations of the eTHSR wave functions, and therefore, gives direct information of whether the container structure is realized or not in these states, and if so, what kind of containers represent the states. In Fig. 41, the contour maps of squared overlap of the states – with single configurations in the and parameter space in Eq. (80) are shown. Here the parameter, i.e., is fixed at the position denoted by in these figures, so that the maximum value of the squared overlap in the four parameter space appears at the position denoted by ×.
As we have discussed in Section 2.1, the single 3α THSR wave function can very precisely describe the ground state and excited states of 12C. It is here useful to show in Table 9 the maximum values of the squared overlap of the , and states in 12C with the single 3α THSR configuration.
In Fig. 41, the squared overlaps of the five eigenstates with the single eTHSR configuration defined in Eq. (80) are shown. In the ground state, shown in Fig. 41(a), 3
α clusters are put into an oblate-deformed and very compact container with
, while the remaining
α cluster is put into a prolate-deformed and very compact container with
. This means that the first 3
α clusters move in a
xy-plane and the last
α cluster moves in
z-direction. This supports the idea that the ground state has a tetrahedral shape of the 4
α clusters proposed by several authors [
118,
119]. Our calculation indicates that this configuration is contained in the
state by 98%.
In the state, shown in Fig. 41(b), the 3α clusters are in a spherical container with . The fourth α cluster is put into a larger size container with spherical shape, i.e., . In particular, the parameter set is almost the same as that for 12C in Table 9, i.e., . This means that the first 3α clusters are confined in a compact container to form the ground state of 12C, since the state can be very precisely described by the single configuration with these parameter values. The fourth α cluster moves in a larger spherical container, because of , which gives the largest squared overlap 94%. This is the new interpretation of the α+12C cluster structure, whose traditional understanding is that the α cluster orbits in an S-wave around the state.
The state, shown in Fig. 41(c), is similar to the state but both containers are not spherical but deformed. The β1 parameter takes almost the same value as that of the isolated state, as shown in Table 9, which means that the first 3α clusters form the state, since the state is described by the single parameter value of . The configuration of the remaining α cluster , giving the largest value 76%, means that the α cluster moves in a deformed and larger container. This is present understanding of the state, which is conventionally considered to have the structure.
In the state, shown in Fig. 41(d), one can see that the 3α clusters are put in slightly larger container than that for the state, which is slightly deformed in a oblate shape. The fourth α cluster, however, moves in a much larger and almost spherical container, like a satellite. This configuration expresses the state dominantly by 84%. This means that the second container characterized by is further evolved from that in the state. We can say that this state corresponds to the state in the former 4α OCM calculation, which predicts the higher nodal structure for the state.
The state, shown in Fig. 41(e), is the most striking. All the α clusters occupy an identical orbit, with . This is qualified to call the α condensation. This configuration is contained in this state by 78%, which is still very large. Furthermore, this container is very close to the one of the Hoyle state, with in Table 9. This means that the state is regarded as the Hoyle analog state, in which the fourth α cluster is also put into the container occupied with the 3α clusters in the Hoyle state. The large size of this container indicates that the 4α clusters are loosely coupled with each other and configured like a gas. Note that the 4α condensate state is also predicted by the 4α OCM calculation slightly above the 4α threshold, as the state.
These results tell us that the evolution of cluster structures is described by the container evolution with respect to its size and shape. The reason why the container evolution arises is the orthogonality to the lower states, which is explicitly taken into account in the definition of the single configuration in Eq. (80). The orthogonality condition prevents a higher state configuration from overlapping with the lower-states more compact configurations. It thus plays a role as a repulsive core and is considered to give the container evolution.
5.5 Application to other nuclear systems
5.5.1 Neutron-rich isotopes of beryllium
Recent years, the study for the cluster structures in neutron-rich nuclei is becoming a very interesting and important subject [
120-
122]. Compared with the typical cluster structures in
nuclei, neutron-rich nuclei could have some novel kinds of cluster structures due to the dynamics of the excess neutrons. This kind of clustering, which is sometimes enhanced by valence neutrons, can have some molecular cluster structures with some kind of valence neutron orbit and this never appears in stable nuclei. For example,
14C is reported to have a 3
α equilateral-triangular shape [
63] and linear-chain structure [
123] with attractive interaction between the excess neutron and
α particles.
As typical examples, the Be isotopes have been studied intensively with cluster models, molecular-orbit models, and AMD [
124] these years. In AMD calculations, the existence of clusters are not a priori assumption like the traditional cluster models. However, AMD calculations show the appearance of 2
α cluster structures surrounded by valence neutrons with some kind of orbit in Be isotopes, which clearly supports the existence of clustering structures in neutron-rich nuclei. Figure 42 shows the density distributions of the intrinsic states of the band heads of
10Be and
12Be obtained by AMD. It can be seen that various cluster structures appear in different states in the Be isotopes. For example, in
10Be, there is a larger clustering deformation for the
state than the ground state. From AMD calculations, we know that in the ground state of
10Be two valence neutrons occupy
π orbits while those of the
state occupy
σ orbits.
Now, the THSR wave function has been extended to describe the general cluster states from the α condensation states or gas-like states to the compact shell-model-like states. The concept of nonlocalized clustering has been proved to be very successful in the nα nuclei. Considering the different features of the clustering in stable nuclei and neutron-rich nuclei, a natural question is, is it possible to describe the molecular cluster structures of neutron-rich nuclei using the THSR wave function?
Recently, the Be isotopes have been investigated using the THSR wave functions [
93,
125]. The constructed THSR wave function for the
9Be can be written as [
93]
where and are creation operators of α-particle and neutron, respectively. represents the coordinate of the extra neutron and its corresponding generate coordinate is . is the creation operator of the extra neutron with spin σ. Here, spin is assumed to be up. The is the azimuthal angle of , which is introduced here for obtaining the negative-parity states of 9Be as discussed in Section 3.3.
As for the Hamiltonian, the effective nucleon-nucleon Volkov No. 2 force [
49] is employed and G3RS [
126,
127] is used as the spin–orbit interaction. In practical calculations, the Monte Carlo method is used for the multi-integration calculations from the Gaussian integration and angular-momentum projection.
In Table 10, it is shown the calculated results of ground-state rotational band by the single THSR wave functions and GCM Brink wave functions. By superposing different configurations for the 2α+n cluster structures, the GCM Brink wave function usually can provide us with a very exact solution for this kind of cluster system. Now, from Table 1, it can be seen that the single THSR wave functions give a very similar energies compared with the results from the GCM Brink wave functions. Especially, there are very high-percentage squared overlap between the single THSR wave functions and GCM Brink wave functions, e.g., as for the ground state 9Be, their squared overlap is as high as 0.96. That means the constructed THSR wave functions in Eq. (82) give a very good description of the ground state, 5/2-, and 7/2- states and the correlations between α clusters and the extra nucleon are well grasped in the constructed THSR wave function.
Figure 43 shows the density distribution of the extra neutron. It can be seen that the extra neutron in 9Be can spread more than 6 fm along the z direction and it has the overlap with the distribution of two α clusters. This is just the well-known π orbit, which has also been obtained by other cluster models. Different with the traditional molecular orbit models, in the THSR wave function, the extra nucleon is assumed to make a nonlocalized motion inside the nucleus rather than occupy some kind of molecular orbit. Finally, the π-orbit emerges from the antisymmetrization which cancels out nonphysical distributions in a very natural way. The reproduction of nuclear molecular orbit structure provides another support for the extension of the nonlocalized clustering concept to 9Be.
Next, we show some calculated results about the
10Be using the THSR wave function [
125]. One constructed THSR wave function for the
10Be can be written as
Here,
is the creation operator of the
α cluster.
and
are the creation operators of the two valence neutrons with spin up and down, respectively. See details in Ref. [
125].
As we know,
state in
10Be was reported to be a developed molecular cluster state from AMD calculations [
124]. The calculated energy spectrum of the
rotational band with the THSR wave function are shown in Fig. 44. It can be seen that the obtained results agree well with those of other theoretical models [
128,
129]. It should be noted that again the only one single THSR wave function can describe well the
state of
10Be while they need a large number of basis for other cluster models. This shows that the constructed THSR wave function is also very suitable for the description of the molecular cluster structures in the neutron-rich nuclei.
By using the single and also simple THSR wave functions, some low-lying cluster states of the 9Be and 10Be are described very well. This means the THSR wave function can also be extended to the non- nuclei and the concept of nonlocalized clustering is also suitable for the neutron-rich nuclei. In the future, to give a unified description of the various molecular cluster structures in neutron-rich nuclei, how to construct some kind of more general microscopic cluster wave function according to the THSR spirit is still a challenging task.
5.5.2 Λ-hypernuclei
Light Λ hypernuclei have been studied by using the microscopic and semi-microscopic cluster models for a long time [
130]. Since a Λ particle is not affected by antisymmetrization on nucleons, it plays a role as an impurity when it is injected into a core nucleus. Many studies in the past made an emphasis on structural change of the core nucleus with the impurity, while the pursuit of knowledge of Λ-nucleon interaction through the structure study of Λ hypernuclei has also been one of the most important issues. A well-known feature in Λ hypernuclei is that the additional Λ particle contracts the nucleon spatial distribution. In the following subsections, the THSR ansatz is shown to be easily applicable to Λ hypernuclei and to better demonstrate the so-called shrinkage effect as contraction of a container of the core.
· cluster structure in . We propose here a new-type microscopic cluster model wave function, which is applicable to hypernuclei and we hereafter refer to as the Hyper–THSR (H-THSR) wave function [
106,
131]. This is based on the above mentioned deformed
THSR wave function.
The H-THSR wave function describing the hypernucleus with good angular momentum is then introduced as follows:
where is the angular-momentum projection operator and the Λ particle simply couples to the core nucleus in an S wave. Its radial part is expanded in terms of Gaussian basis functions, with respect to the width parameter κ, where with a position vector of the Λ particle. In the practical calculations, we assume an axially-symmetric deformation for , .
The coefficients of the expansion are then determined by solving the following Hill–Wheeler-type equation of motion,
This H-THSR wave function is characterized only by the parameters and , which correspond to a spatial extension of the whole nucleus. It is well known that the Λ particle invokes spatial core shrinkage in many hypernuclei. It should be mentioned that such a core shrinkage effect is expected to be taken into account very naturally by this parameterization of , in general hypernuclei, since it specifies the dilatation of the whole nucleus. In this section, only the S-wave component of the Λ particle is considered, for simplicity. Thus the monopole-like shrinkage is expected to be described nicely by this wave function. The extension to inclusion of the other angular-momentum channels is discussed in the next section for the study of .
In order to compare the single H-THSR wave function with the Brink-GCM wave function, we calculate the following squared overlap:
The above wave function is the Brink-GCM wave function projected onto the model space with the angular-momentum channel , in which L, λ are the angular momentum between the clusters and the one between the Λ particle and , respectively.
In Fig. 45, we show the contour maps of this quantity in the two-parameter space and , calculated for the (a), (b), and (c) states. For these states, we can see two maxima in prolate and oblate regions (denoted as × and +), as in the case of the energy surfaces. The values of the maxima are extremely large and close to unity. The values for the states are 99.5% and 99.4%, respectively. For the state, the maximum value is slightly down to 97.7%. This reduction may originate from a mixture of spurious scattering-state components. These practically % squared overlap values of course mean that, at least on the subspace with and , the Brink-GCM wave function obtained by solving the Hill–Wheeler equation is equivalent to the single configuration of the H-THSR wave function. We should note that the Brink-GCM wave function in this subspace is shown to be almost the same as the one in the full angular-momentum-channel space. This result also means that the size parameter , which specifies the monopole-like dilatation of the whole nucleus, quite well takes into account the effect of the spatial core shrinkage by the additional Λ particle. We can thus conclude that in the clusters are trapped into a container, which is specified by the optimal value of the size parameter , under the influence of the antisymmetrizer acting on the nucleons, such as realized in the form of the wave function Eq. (85).
We here discuss the intrinsic structure of , together with that of . Using the intrinsic wave function defined as the THSR wave function before the angular-momentum projection in Eq. (85), we can calculate the following intrinsic density of nucleons,
This density is normalized as usual to the total number of nucleons, .
In Fig. 46, we show the intrinsic density profiles of defined by Eq. (88) and of . The single optimal values giving the minimum energies are adopted, i.e., for and for . While both ones clearly show the 2α-cluster structure with the prolate-deformed shape, has a gas-like tail of the 2α clusters and not. In the Λ particle gives rise to strong shrinkage and the gas-like tail in disappears. The r.m.s. radius of the core is accordingly changed from fm for to fm for . Nevertheless the 2α-cluster structure definitely remains in this very compact object, which is produced by the competition between the quite strong effect of inter-α Pauli repulsion originating from the antisymmetrizer and the fairly strong attractive effect among the 2α and Λ particles, and therefore the 2α clusters in this intrinsic state are effectively localized in space. This means that even for the states which are described by nonlocalized-type wave function with container structure, localized nature of clustering can appear in density distribution due to the Pauli principle.
We next simulate the shrinkage effect by varying interaction artificially, with overall factor δ multiplied with . The values to give the energy minima after projection onto space, which are listed in Table 11, are adopted. As the interaction is strengthened, from to 0.8, the α+α distance is shortened along z-direction. The corresponding values, Λ binding energies, the maximum squared overlaps in Eq. (87) with the use of the same artificial potential, rms radii of the core, are also shown in Table 11. and 1 correspond to the cases of and , respectively. Starting from the case of with very large value of and small value of , i.e., fm and fm, only the value drastically gets smaller as the increase of interaction, and eventually for , becomes much smaller while is almost unchanged, i.e., fm and fm.
We should note that at each step of the change of
interaction, from dilute 2
α cluster structure to compact localized 2
α cluster structure, all are described precisely by single H-THSR wave functions, which have optimal sizes of deformed container, i.e., those with large-size container to small-size container. We should emphasize that at every step of
δ, the squared overlaps with the full Brink-GCM solutions are more than 99.3%. This may imply that the THSR-type container picture is essentially important in understanding every type of cluster structures, from the compact to dilute ones. We should also emphasize that for
a single optimal Brink wave function can only approximate the Brink-GCM wave function by 0.722 [
48], while for
the former can approximate the latter by 0.940 [
106].
· Multi-cluster dynamics in . We discuss
in this subsection by using the H-THSR ansatz.
is theoretically investigated and the second
state is found as a bound state, which is considered to correspond to the
state [
133,
134], although experimentally the second
state in
has not yet been observed [
135]. The H-THSR wave function introduced in this study is slightly different from the one used in the previous subsection for the study of
, as follows:
with
where the is the wave function of the Λ particle with its width parameter, , which can be out of the antisymmetrization, with mN and being nucleon and Λ-particle masses, respectively, and . In this wave function, the Λ particle is also confined in a deformed container, so that the present H-THSR wave function is more extended than the previous one, which only contains the S-wave component.
Also in this calculation, the axially symmetric deformation is taken into account, i.e., and . After performing the angular-momentum projection for the H-THSR wave function in Eq. (89), the energy eigenstates can be obtained by solving the Hill–Wheeler equation by taking the six parameters as generator coordinates.
In Fig. 47, the calculated
spectra of
12C and
and the observed
spectrum of
12C are shown. We discussed the 0
+ states of
12C in detail in Sections 2.1.5 and 5.3 (see Figs. 9 and 30). Let us keep in mind that the
state is concluded to be the excited state from the Hoyle state by monopole vibration and the
state dominantly has a 3
α linear-chain configuration. While the
state is considered to be a family of the Hoyle state, the
state has a quite different structure from the
and
states. In
, first we should mention that we simply denote the spectrum of
as
states, neglecting the intrinsic spin 1/2 of the Λ particle, for simplicity. From this figure, one can see that
gives a much larger number of decay channels than
12C, i.e.,
,
,
,
,
,
, and
. Up to the
threshold, we have the four
states, corresponding to the four
states in
12C. We can see that the
in
is obtained as a bound state which is located below the lowest threshold,
, by 0.8 MeV. The
and
states are connected to the
and
states of
12C, respectively, as you can also see from the figure. This assignment is because the
and
states in
largely include
and
components, respectively [
131].
In fact, the state in is found to keep essentially the same features as the state in 12C, which are shown in Section 5.3. For example, we obtain very similar result as in Fig. 32, also for the state in , where the squared overlap with the single THSR configuration is calculated. The maximum squared overlap value is calculated to be 47%, which is almost the same as state in 12C, but the parameter values to give the maximum is shifted from in the state of 12C, to in the state of . The linear-chain configuration in then shrinks along z-direction, compared with the case of 12C, due to the additional Λ particle.
In Table 12, we show Λ binding energy
, the root mean square (r.m.s.) radius of the core
, the r.m.s. distance between the core and Λ particle
, and the r.m.s. radius of
,
, together with the r.m.s. radii of isolated
calculated in Ref. [
16]. Here the Λ binding energy is defined as the binding energy measured from the corresponding
threshold. From the calculated r.m.s. radii for these states, we can expect that the
state is qualified to be a gas-like cluster state, since only this state has a larger r.m.s. radius
fm than the one of the Hoyle state
fm. From the RWA analysis, we find that the
state has the largest component of
, in which
has a dilute
structure [
131]. Therefore this state is a loosely coupled
state, i.e., a Hoyle-analog state. In this table, we can also see that the
state suffers from much less shrinkage effect from the additional Λ particle than the other
and
states, except for the ground state that has saturation density and gives almost no shrinkage effect. The
value for this state is also smaller than the one of the other states. This is because for the
state the Λ particle couples with only one
α cluster, to form
. The little overlap with nucleons reduces the shrinkage as well as the gain of the Λ binding energy. In fact, this state appears around the
threshold energy, for the threshold rule. The slightly higher energy position,
MeV above the
threshold, allows for the diluteness of this state.
On the other hand, the dominant component included in the state is configuration. The Λ particle then overlaps with the clusters and therefore gives much stronger shrinkage effect and larger Λ binding energy. The binding energy of the is so large that it is energetically favored that the injected Λ particle shrinks the core and forms cluster configuration. In fact, this state is located at around the threshold as well as around the threshold (see Fig. 47), which follows the so-called threshold rule.
We further show in Fig. 48 the RWAs in the
channel of the
(solid curve in red) and the
(solid curve in blue) states in
. In comparison, the corresponding RWAs of the Hoyle state (dotted curve in red) and the
state (dotted curve in blue) in
for the
channel are also shown [see also Fig. 11(a)]. As we discussed in Section 2.1.5, in the Hoyle state (dotted curve in red) the nodal behavior appears as oscillation because of dissolution of
core, while in the
state in
(dotted curve in blue) there exist four nodes. On the other hand, in
, the most prominent feature is the disappearance of nodes for the
state, though the innermost node only remains. This means that the
core is dissolved in the
state, completely unlike the case of the
state in
. We should also be aware that for the
state a long tail is still very much developed, which is similar to the Hoyle state. This is also completely different from the
state, in which the whole amplitude is pushed inside with no long tail any more. All these results give a strong support that the
state has a gas-like structure of
clusters, as a Hoyle analog state. It is important to mention that this Hoyle analog state uniquely appears in
and we cannot expect in
gas, since
is not bound [
136], whereas
is bound by 3.1 MeV.
5.6 More generalization and manipulation of the THSR wave function
5.6.1 Two-center THSR wave function
Recently, to investigate the correlation between the valence proton and 2
α clusters in
9B, a two-center THSR wave function was proposed [
137]. Based on the designed THSR wave function for the description of the
9Be(2
α+
n) cluster system [
93], one similar THSR wave function for
9B(2
α+
p) can be written as
where is the size parameter for the nonlocalized motion of 2α clusters and is the size parameter for the valence proton around the 2α clusters. The is the phase factor for the treatment of the negative-parity states. is the Brink wave function, in which the (), and are corresponding generator coordinates for the α clusters and the valence proton, respectively.
In the
system, many configurations should be considered for some developed cluster states due to the nonlocalized motion of the valence proton. To describe these cluster states and also discuss the correlations of valence proton and
α clusters in a more accuracy way, the following THSR wave functions are constructed [
137],
It should be noted that, compared with the THSR wave function in Eq. (91), the correlation of the valence proton and α clusters are considered in the by slightly changing the generator coordinate of the valence proton in the Brink wave function in the construction.
Finally, the THSR wave function of 9B can be constructed as a superposition of both cluster-correlated configuration and large spreading configuration
Here and can include cluster-correlated configurations of 9B describing the correlated motion of valence proton around each α-cluster. The wave function is corresponding to the large spreading configuration of 9B in which the valence proton orbits around the 8Be core. The c and d are coefficient parameters of the two-center THSR wave function.
Figure 49 shows the obtained energy spectrum of the
9B using the two-center THSR wave function and the comparisons with those obtained by the THSR wave function
, GCM Brink wave function, and the experimental data. It can be seen that the significant improvements for the binding energy of
9B using the two-center THSR wave function which are in much better agreements with the experimental values. Furthermore, by analyzing the weight (c,d) in the two-center THSR wave function, it can be seen that the cluster-correlated and the large spreading configurations both have significant contributions to the dynamics of the ground state of
9B. Details discussions can be found in Ref. [
137].
5.6.2 Cluster systems with different size variables
In nuclear cluster physics, the widths of Gaussian packets describing clusters are very important dynamics variables especially for the asymmetrical cluster systems of nuclei, which have been discussed in many works for many years [
138-
141]. However, the introduced width variables will lead to the serious center-of-mass problem, i.e., the intrinsic wave function and the center-of-mass wave function cannot be separated exactly, which is a long-standing problem in nuclear cluster physics.
In RGM [
141], the trial cluster wave function is constructed in term of the translational invariant coordinates and the width variables can be naturally included. Nevertheless, the practical calculation is becoming very tedious. A more practical method for the description of cluster structure is GCM combined with in the traditional Brink wave function [
41] and other microscopic cluster wave functions [
142], in which the cluster wave function can be expressed as the superposition of Slater determinants. However, the center-of-mass problem still remains. Quite recently, we found that, in the THSR framework, cluster systems with different size variables can be described in a simple way, in which an exact solution to the center-of-mass problem was given. The key to separate exactly the center-of-mass part and the translation-invariant intrinsic wave function part is the following constructed integral operator [
143],
Here, the is a shift operator. With the auxiliary generator coordinate , this operator can perform the integral transformation for exactly separating the center-of-mass part from a many-body wave function. This handled wave function can be antisymmetrizd or not, nuclear cluster or nucleon wave function.
By means of the integral operator, a new THSR-type wave function can be naturally obtained, in which the center-of-mass part can be exactly separately. Take a two-cluster system (C=C1+C2) as an example. The THSR-type wave function can be constructed as follows:
Here the clusters C1 and C2 have the atomic mass number A1 and A2, respectively. and are the corresponding center-of-mass coordinate, respectively. and . and are the width parameters or the size parameters of cluster C1 and C2. is the general Brink wave function with different widths of clusters. It can be seen that compared with the previous THSR wave function in Section 4 for , the main change of Eq. (95) is just one more integral for the generator coordinate R. The computation is very straightforward and simple, in particularly, compared with the other RGM-based approaches.
Following the two-cluster wave function in Eq. (95), the wave function can be written directly,
Here
and
.
and
are the center-of-mass coordinate of
α cluster and
cluster, respectively.
and
are the width parameters or the size parameters of
α cluster and
cluster, respectively. The
is a generalized Brink wave function [
41,
144] for
system with different width variables. As we know, in the traditional Brink wave function, the widths of different clusters usually are fixed at the same values, i.e.,
in Eq. (97), to avoid the center-of-mass problem.
By making variational calculations for widths , , and variables, the energy curve of and the variation of width variables can be obtained. Figure 50 shows the energy curve of along with the obtained optimum widths and using the intrinsic wave function in Eq. (98). Clearly, it can be seen that with the increase of the size parameter indicating the cluster system expands, the gap between and is becoming larger and larger. Another saying, imagine a situation where the stable double-closed-shell nuclei α and are approaching each other to form a stable nuclear state (the ground state of ), in the process Fig. 50 tells us the is slightly shrinking while the α nucleus is becoming more loosely bound. It has been known that the ground state of has a very compact cluster state. In this case, the obtained minimum energy for the ground state is -159.74 MeV. The corresponding wave function is characterized by a rather small value of the size parameter fm and slightly different values of widths fm and fm. As a comparison, by using the traditional angular-momentum projected Brink wave function in a variational calculation with the common width variable and inter-cluster distance parameter R, the minimum energy fm, fm) MeV is obtained. This is more than 1 MeV higher than the optimum energy obtained by our new intrinsic wave function.
In the new THSR-type wave function, not only the size parameter for the relative motion of clusters but also the width of the clusters are also included. Actually, without considering the width variables, many physical quantities like radius, differential cross section, phase shift actually cannot be well reproduced and it has been discussed for a long time in the RGM approach [
1,
145,
146]. Therefore, the gas-like cluster structure, halo nuclei, the distortion effects of clusters in light nuclei are very promising to be investigated in a more realistic way in this new framework.
5.6.3 Monte Carlo integration in the THSR wave function
As we discuss in Section 2, the original THSR wave function has a very simple mathematics form and it is also related with the Slater determinant-based Brink wave function. In the practical calculations, the THSR matrix element of Hamiltonian can be directly expressed as
To compute the above THSR matrix element, the first step is to obtain the analytical kernels of the Brink wave function, which has been studied for many years and there are many well-developed techniques and methods for this kind of analytical calculations [
44]. The second step is to perform the integrals on the obtained analytical expressions. Because the weight function has a simple Gaussian form, the integrals for the generator coordinates can also be done the analytical way, which does not lose any accuracy. For example, the calculations for the 2
α, 3
α, 4
α, and
(
α+
) systems discussed in previous sections are all performed in the exact analytical way.
However, in some heavier multi-cluster systems, e.g., 6α cluster system, this kind of analytical calculations for the THSR kernels are becoming very difficult. Firstly, even the calculations for analytical expressions for the Brink kernel is not so easy due to its huge size. Moreover, analytically, the computations for the high-dimensional integral is very expensive. For example, based on the Eq. (99), to study the 6α cluster system, we need to perform Gaussian integrals for the generator coordinates to the extremely large analytical expressions from the Brink kernels, which will lead to time-consuming computations. To overcome the high-dimensional integral problem, a natural method is to use the Monte Carlo technique for the integrals. As we know, Monte Carlo integration does not depend on the number of integral dimensions, which can be considered as a very economic advantage of Monte Carlo integration against the analytical integrals. Besides, once we use the Monte Carlo technique for the integral, the Brink kernel part is not necessary to be obtained analytically. Thus, the complexity of the calculations can be reduced.
To study the possible
α condensation states around the heavy core nucleus, Itagaki and his collaborators introduce the Monte Carlo techniques [
147,
148] for the treatment of the THSR wave function. The original THSR wave function can be reformulated as
Here is the wave function for the i α cluster centered at , which has a Gaussian form Exp with the width ν. To simplify the calculation of this THSR wave function, the integral over the generator coordinate parameters in the original THSR wave function can be replaced by the summation of many Slater determinants as follows:
The generator coordinate parameter can be randomly produced, however, the random numbers are generated by the weight function with the same Gaussian shape:
With the increasing ensemble number, the distribution of will be close to the Gaussian with the width parameter . Therefore, this kind of approximation method can be considered that the integration of the THSR wave function is performed by using a Monte Carlo technique. Clearly, when the number of the superposed Slater determinants is becoming a very large number, the obtained wave function will agree with the exact THSR wave function. Most importantly, this numerical treatment for the integral is much easier than the analytical integral, especially the high-dimensional integral problem.
To check the reliability of the Monte Carlo technique for the integral of the THSR wave function, Itagaki et al. studied from the well-studied
8Be and
[
147]. It was found that as for the
8Be(
α+
α) cluster system, except for the line with small
β (=2 fm) value, all other lines have a good converge character to the 2
α threshold energy even with a small number of superposed wave functions. By using the Monte Carlo, the Hoyle state was also studied and the results are consistent with the original THSR calculations. This means, this kind of Monte Carlo technique is reliable for the calculations of the THSR wave function.
By using the Monte Carlo technique, the THSR wave function can be extended to study the heavier cluster system. Itagaki et al., studied the α-condensed state with a heavy core. One simple case is the α-condensed state around the core in . Figure 51 shows the energy convergence of measured from the +2α threshold, where two α-clusters occupy the α-condensed state around the core. The solid and dotted lines (, 3 fm) converge to the energies below the threshold, however, the dashed- and dash-dotted-lines (, 5 fm) converge almost around the threshold energy. Here, the converged energies of the these lines show that they are rather insensitive to the β values. Thus, the α-condensed state is considered to be realized in around the +α +α threshold energy.
Recently, the Monte Carlo technique is also applied to study the non-n
α nuclei [
93]. As for the non-n
α nuclei, the constructed THSR wave function has a different weight function compared with the original n
α THSR wave function, which is not suitable for the analytical integral. However, using the Monte Carlo technique, the weight function is not necessary to be limited the Gaussian form and the calculation is becoming much easier. The Monte Carlo technique provides us with a flexible way to extend the THSR wave function to heavier multi-cluster systems.
5.6.4 New development of analytical treatment of the THSR wave function
We are extending the THSR framework for dealing with various nuclear systems, which is one of the key ansatz for the nuclear cluster aggregate. Therefor, let us consider the practical analytical treatments of the THSR ansatz beyond the original container model which has been restricted to one container by Gaussian confinement including deformation. Here we present three kinds of possible extension as follows: Firstly, the introduction of a distant parameter between two aggregates of nuclear cluster in the integrand of the wave function. This makes it possible to estimate the parity dependence including the odd state in the container. The kind of wave function for
α+
has been discussed in Section 4. Secondly, the introduction of the same distant parameter as the previous one in the integral operator of the THSR ansatz. This can give two containers with different centers each other. It is very interesting to clarify the interaction between two cluster-aggregates through two containers, which has been applied for the calculations for the
n–3
α system in Section 3.4. Thirdly, the partial condensation of the nuclear cluster aggregate. For instance, let us see 3
α clusters around
16O core. The integral operator acts only on outer 3
α clusters. Itagaki et al. applied this ansatz to
24Mg [
147] which is represented by 2
α–
16O core. They carried out multiple integral by the Monte Carlo method. Here we devote our efforts to analytical derivation of exchange kernel coming from the Pauli principle, for 2
α–
α system as an example for three kinds of extension. We will concentrate on the norm kernel for three kinds of extensions for simplicity.
· Original THSR ansatz for 3α system. At first, we should give a brief explanation of treatment of the original 3α THSR ansatz in order to compare possible extensions of the THSR ansatz. The THSR ansatz consists of two parts; that is, the integrand which is fully antisymmetrized wave function on the parameter coordinate space, and the integral operator which acts on each α cluster on the space. The integrand is described as follows:
where suffix p means the spin and the isospin part. and represents the coordinate of the α-cluster in the parameter space R. The spatial part of a single nucleon wave function is written by
where . The variable stands for the real coordinate for nucleon with spin-isospin state. We call the integrand the Brink-type wave function.
The integral operator, which is a key stone in the THSR ansatz, and which acts on the parameter coordinates of clusters, is given by
where the integral is carried out over the infinite parameter space. The THSR wave function is represented by
The size parameter of the container is denoted as β.
Let us obtain the norm kernel by the analytical expression:
The overlap of the wave function (105) is given after carrying the integral with respect to real variable :
where the overlap between single nucleon wave functions is written by
where the suffix p is for all the spin-isospin states, where . Therefor, all the terms, which come from the permutation, in the integrand of Eq. (109) are written by the Gaussian function. For instance, the first term in the integrand corresponding to has the exponential form:
where matrices are given by
where only the matrix is varied depending on the exchange property, but the matrices are unchanged during all the permutations.
The multiple integral with respect to the parameters, is easily obtained because the integrand is written by the summation of the Gaussian form. The result of the first term including the coefficients is as follows:
where
We concentrate our thought on a computing procedure of the analytical derivation of the determinant; that is, the most tedious and primitive approach is the best way for programming code. We rewrite the determinant part in Eq. (114):
At first, we have to divide all the columns into two parts; that is, for instance, the first column is written by
As consequences, the determinant can be estimated by the summation of 26 separated determinants which are represented by the polynomial with respect to x and . Almost all of the separated determinants are exactly 0 except for the power of due to the physical background. Thus, we can imagine that the derivation of the determinant of each term in Eq. (114) is easy task, and all the summation of the terms is obtained with analytical form of the norm kernel. Of course, we have to elegantly arrange them because there are a lot of the exactly same terms. The analytical formula of the norm kernel is written by the following equation:
where
The weight Wk, which comes from the same exchange character of
and the coefficients
are listed in Table 13. All the 9 terms with respect to the polynomial have a common factor
because of the inclusion of c.o.m parameter. Unfortunately, we cannot find the common factor automatically by the programming code above mentioned. We have already presented the more advanced method in Ref. [
150], in which we pointed out that the permutation property of nucleons among the 3
α clusters is naturally found out. Note that the important parameters
relating to the
α condensation strength are defined in the range of
. In the case of
, the system goes to shell-model state, on the other hand, for
the system shows a free
α gas aggregate.
· The integrand with distance parameter. Let us assume –α system. The distance parameter d between and α clusters is introduced; namely, In Eq. (110) the matrix elements of overlapping are explicitly written by
Concerning Eq. (113) related to , we write the power part in exponential as
where matrices are given by
and the other scalars are
The D matrices and the scalars are independent of the permutation property. The matrices , and are common to those in Eq. (113).
The multiple integral with respect to the parameters, where the parameter d and d' remain being not carried out the integral is easily obtained because the integrand is written by the summation of the Gaussian form. The result of the first term is as follows:
where
The analytic form of the numerator should be considered as well as that of the denominator, where the same procedure is applied to how to analytically derive the numerator. Using the representation of double tilde followed by Eq. (114), we rewrite the matrix element as
where
The explicit form of the numerator in Eq. (123) is
The same procedure as that of the original THSR derivation is carried out. As a consequence, the analytical form is written by
where
The coefficients are listed in Table 14, where is satisfied for each type of k in the original THSR ansatz. The property of the function with respect to the polynomial of is as follows:
As for this primitive and tedious method, it is difficult to carry out the unique factorization, but easy to do an exact angular-momentum projection because of analytical expression of .
· The integral operator with the distant parameter. Now we consider the interaction between two aggregates of α clusters. Here, we take interacting system of 2α as one aggregate and α as another, which is the simplest aggregate. It is reasonable to introduce two kinds of condensation strengths. Eq. (107) has more complex forms as
In the same way as the previous discussions, the Eq. (123) is
where matrices are
The scalars are as follows:
First of all, we should scrutinize two-strength THSR ansatz as the case of . Namely, for the first term of the deteminant is given by
where
The function of should be also exchanged by as
After arranging all the terms, the number of the independent terms is 25, which is exactly the same as that of the previous section. This is conserved in the case of . However, the function of is written by the combination of 4 kinds of x parameter. As for the case of , the function for each term is analytically written by the same form as that of Eq. (122). The numerator for the term of is given by dividing the final column as
Thus, we can easily write the numerator such as
where we use the relation of
for the final row. Thus, we can give the analytical expression of the case of two-center THSR ansatz after obtaining the independent terms. The property is exactly same as that of the introduction of distant parameter in the single THSR ansatz. This is because the special role is given for two α-cluster aggregates. Nevertheless, the Pauli principle through all the nucleons is completely carried out.
· Partial integral for the THSR ansatz. According to the previous discussions, the partial integral can be easily done in the THSR framework. Fixing to the origin , we carry out the integral with respect to . For the first term of , the determinant is written by
The resultant norm kernel is as follows:
where
The coefficients in this ansatz are listed in Table 15. The number of the independent terms is 20, which is larger than that of the original THSR ansatz because of the break of the symmetry. However, comparing the case of the introduction of the distant parameter d, the number of the independent terms is degenerated.
In summary, on the basis of the complete consideration of the Pauli principle, the THSR ansatz has been successfully applied to various cluster systems. The positions of clusters are not fixed, then they can freely move in the container. The introduction of the dimensionless parameters in the THSR ansatz plays an decisive role in clarifying the property of the cluster aggregates.
Here, we mainly present two kinds of extensions including distant parameter d. The first one is related to the inner part in the single container. Its introduction enables us to find out the even-odd difference of two-cluster system in the container. The other is introduced to consider the interaction between two containers. Analytical derivation of the exchange kernel needs to arrange a tremendous number of the permutation coming prom the Pauli principle. In addition, the multi-integral concerning the number of clusters looks difficult task. In order to avoid such a cumbersomeness, numerical integral has been proposed by some authors. We cannot deny such efforts because, for instance, even 10 α system has terms. Nevertheless, the analytical derivation of the exchange kernel makes an important base to scrutinize the physical property of nuclear cluster system.
Naturally, as the next step, we should present a new method to treat resonant states or scattering states within the framework of the THSR ansatz. One of the key notes is to express explicitly the real variables, with nucleon number n. We are now going to prepare a report on this subject.
6 Evolution of cluster structure
6.1 Duality property of mean-field states and formation of cluster states
6.1.1 Duality property of shell-model states with good SU(3) symmetry
As we discussed in Section 3.1, the wave functions of the ground states of the self-conjugate
nuclei (“
” nuclei) have duality property. The duality property means that the wave function has both the character of mean-field-type structure and that of cluster structure. The ground-state wave functions of these nuclei up to at least around 40 have dominant component of good SU(3) symmetry. According to the Bayman-Bohr theorem [
151], the wave function with good SU(3) symmetry has dual character of mean-field-type structure and cluster structure. A typical example which we referred to in previous sections is the duality property of the double closed-shell wave function
of
16O. This wave function with
symmetry of SU(3) can be rewritten in the form of cluster-model wave function of
12C +
α configuration,
Here r is the relative coordinate between 12C and α clusters and stands for the harmonic oscillator function with the oscillator quanta , namely the nodal number , and with the size parameter . is the internal wave function of 12C with SU(3) symmetry and is given by the -shell shell-model wave function as
The relative wave function has SU(3) symmetry . stands for the SU(3) vector-coupling of with symmetry and with symmetry to the SU(3) symmetry by the use of the SU(3) reduced Clebsch-Gordan coefficient ; .
The angular-momentum-coupled wave function of 12C + α configuration can be expressed by the linear combination of the SU(3)-coupled wave function of 12C + α configuration as
The number of total oscillator quanta, Q, of the double-closed-shell wave function is . Since the 16O wave function having is restricted only to this double-closed-shell wave function having SU(3) symmetry, we know that there hold
Therefore Eq.(139) can be written as
Thus the double-closed-shell wave function can be equivalently expressed in three ways by the use of 12C + α cluster-model wave functions, 12C() + α(S-wave), 12C() + α(D-wave), and 12C() + α(G-wave).
Besides the case of the double closed-shell wave function of
16O, many other cases of duality property have been discussed rather in detail in many review papers of nuclear clustering [
6,
152]. They include the ground-state wave functions of
12C,
20Ne,
28Si,
40Ca, and
44Ti. The shell-model wave function with duality property is not restricted to the ground state. A typical example is the
excited state of
32S which is the main component of the superdeformed state and has the configuration with 4 particles raised from
-shell to
-shell. The intrinsic wave function of this state has the duality character that it is equivalent to the
16O-
16O cluster model wave function as shown below:
Here expresses the Cartesian harmonic oscillator single-particle function with nk being the number of oscillator quanta along k direction (, y, or z). The coordinate rO−O stands for the relative coordinate between two 16O clusters and stands for the Cartesian harmonic oscillator function of coordinate r with size parameter γ. This intrinsic wave function belongs to the SU(3) symmetry .
6.1.2 Formation of cluster states from mean-field state with duality property
The duality property of the shell-model wave function of the ground state was already known in 1950’s among people who studied microscopic cluster model. For example, Perring and Skyrme wrote at the beginning of the abstract of their paper [
153] “It is shown that it is possible to write down
α-particle wave functions for the ground states of
8Be,
12C, and
16O, which become, when antisymmetrized, identical with shell-model wave functions.” The microscopic cluster model study by Wildermuth which he started with his paper at 1958 [
154,
155] was on the basis of the duality property of the shell-model wave function of the ground state. The year 1958 was about 10 years after the proposal of the shell-model by Mayer and Jensen at 1949 and therefore the shell-model had already been well established [
156]. Wildermuth aimed to discuss the cluster physics which coexists with mean-field physics. Wildermuth considered that the formation of cluster states is due to the activation of cluster degrees of freedom embedded in the ground state due to its duality property. Namely he considered that, since the ground state wave function possesses the duality property, there occurs two kinds of excitation of the ground state, one is the activation of the mean-field degrees of freedom of the ground state yielding the mea-field-type excited states and the other is the activation of the cluster degrees of freedom of the ground state yielding the cluster states. His philosophy can be shown as in Fig. 52.
Wildermuth assigned cluster structures to many excited states in light nuclei [
157]. For example, in
16O he assigned
12C(0
+) +
α structure to 0
+(6.06 MeV), 2
+(6.91 MeV), 4
+(10.36 MeV), 1
-(9.58 MeV), and 3
-(11.62 MeV). He also assigned
12C(2
+) +
α structure to 2
+(9.84 MeV). Wildermuth’s assignments are mostly supported by the present nuclear structure study.
In Refs. [
6,
152], detailed discussions are given that, in many self-conjugate
nuclei from
12C up to
44Ti, many cluster states are really existent as excited states in accordance with the prediction by the duality property of the ground state. We can say that the existence of cluster states in addition to mean-field-type states is an inevitable consequence of the duality of the ground state.
We discussed the duality property of the excited state by using the example in
32S given in Eq. (142). The AMD study of Ref. [
158] showed that the excitation of the
16O-
16O clustering degree of freedom embedded in the
excited state of
32S expressed in Eq. (142) yields the molecular resonance states of
16O-
16O. This result means that the formation of the
16O-
16O molecular resonance states by the evolution of nuclear structure from the
32S ground state is through two steps, one is the formation of the
excited state from the ground state by activating mean-field degree of freedom and the other is the formation of the
16O-
16O molecular resonance states from the
excited state by activating the cluster degree of freedom.
Mean-field state can have the duality property when the state has a dominant component of good SU(3) symmetry. Therefore the ground state of heavy nucleus is considered to have no duality property since the SU(3) symmetry of Elliott is considered to be nonexistent in such ground state. The cluster states we know until now except for lighter nuclei than 12C can be regarded as being formed by the excitation of cluster degree of freedom from some respective mean-field states, but it does not necessarily mean the nonexistence of cluster state which are not related to any mean-field state through the duality nature. In lighter nuclei than 12C, the formation of the mean field of nucleons in the ground state is not so firm as is seen in the ground state of 8Be which has cluster structure.
6.1.3 Observed large monopole transition strengths between ground state and cluster states
Cluster states have very different structure from the ground state with mean-field-type structure. For example, the Hoyle state of
12C is a 3
α gas-like state and its density is about 1/3 of the ground state. Also, the 6th 0
+ state of
16O is predicted to be a 4
α gas-like state and its density is calculated to be more dilute [
19,
20]. However, the observed strengths of the monopole (
E0) transitions between cluster states and the ground state are large and comparable with the single-nucleon strength in its order of magnitude, in spite of large difference of structures of initial and final states of the transition [
159].
The
E0 single-nucleon strength [
79] is roughly given by
with
R standing for the nuclear radius. It is obtained by calculating
by using uniform-density approximation for
and
,
for
. If we adopt
for
R for light nuclei like
12C and
16O, we have
. The observed
E0 value between the Hoyle state and the ground state in
12C is 5.4±0.2 fm
2, and that between the first-excited 0
+ state and the ground state in
16O is 3.55±0.21 fm
2. The first-excited 0
+ state in
16O is a well-known cluster state which has the structure of
+
with
standing for an
α cluster moving around
with
S-wave
.
It looks very contradictory that we have, on one hand, the large difference of structures between cluster states and the ground state with mean-field-type structure, but, on the other hand, we have the large magnitude of the monopole (E0) transitions between cluster states and the ground state. The reason why it looks contradictory comes from our way of understanding of cluster states as being described by superpositions of many particle-hole excited configurations which are very much complicated compared with the ground state. The matrix element of the monopole transition operator between the ground state and the cluster state expressed by the superposition of many particle-hole configurations looks to be much smaller than the E0 single-nucleon strength. Therefore in order to resolve this seeming contradiction we have to abondon above-described way of understanding of cluster states in the language of mean-field theory.
It was shown in Ref. [
79] that the resolution of this seeming contradiction is very easy if we rewrite the ground-state wave function in the form of cluster-model wave function according to the duality property of the ground state. When we express the ground-state wave function in the form of cluster-model wave function, it is immediately clear that the monopole transition between the ground state and a cluster state is nothing but the monopole transition of the inter-cluster relative motion of the ground state to that of the cluster state. For example, we consider the monopole transition in
16O between the first-excited
(
) state and the ground state (
). Because of the duality property the ground state is expressed as
with
. On the other hand, since the first-excited 0
+ state has the
+
structure, it is expressed as
. Thus the
E0 transition between two states is the transition between
and
. Since the
E0 transition operator
can be written as
we see that only the relative-motion part,
, contributes to the
E0 transition. We can easily show that both the operators
and
do not contribute [
79]. The fact that the
E0 transition reduces to the
E0 transition of the inter-cluster relative motion means that the
E0 transition is due to only single degree of freedom, inter-cluster relative motion. Therefore it is very natural that the observed strength of
E0 transition is comparable to the single-nucleon strength. One may think that the effect of the antisymmetrization operator
A invalidates the above argument. But in Ref. [
79], it is shown that the effect of the antisymmetrization operation is small. Actually the result of the analytical calculation of the
E0 matrix element is given as
The quantity represents the effect of the antisymmetrization and actually is fairly smaller than unity in general for non-large N. However, in the above analytical formula, quantities appear in the form of ratio, , and the magnitude of the ratio is close to unity, which implies that the effect of antisymmetrization has only little influence on the value. The quantity is the coefficient of the 2-jump component contained in , and its magnitude is around 0.4. Note that is not percentage quantity, . The E0 matrix element of the relative motion, is larger than the corresponding matrix elements of the single-nucleon motion, and , by about 50%. We thus see that the order of magnitude of is the same as the single-nucleon strength.
In
16O, the
E0 transition from the second-excited
(
) state to the ground state has been also observed and its observed strength is large as 4.03 fm
2. The dominant configuration of this excited state is the cluster structure of
+
and the double-closed-shell wave function of the ground state can be expressed, because of its duality property, by the cluster-model wave function of the configuration of
+
as we explained in Section 6.1.1. Therefore the
E0 transition in this case can also be regarded as being the
E0 transition of the relative motion between
cluster and
α cluster. Here also the
E0 transition is due to only single degree of freedom, inter-cluster relative motion and hence it is very natural that the observed strength of
E0 transition is comparable to the single-nucleon strength. In Ref. [
79], there is shown the analytical formula of the
E0 matrix element between
and
states which is quite similar to that between
and
states.
The numerical values of
and
calculated by the analytical formulae mentioned above with suitable parameter values are 1.97 fm
2 and 3.89 fm
2, respectively. These values are somewhat smaller than the observed values although the order of magnitudes are well reproduced. In Ref. [
79] it was shown that the inclusion of the ground-state correlation into the ground-state wave function makes the reproduction of the observed values by theory fairly satisfactory. The ground-state correlation adopted is that of the inter-cluster relative wave function; namely the
12C-
α relative wave function of the ground state is not simply
but has contributions from higher
N components,
with
. The important point of this result is that the cluster degree-of-freedom due to the duality property induces the ground-state correlation.
There have been many cluster-model calculations of
E0 transition strengths by many authors. As far as the microscopic wave functions like those of RGM, Brink-GCM, AMD, and THSR are used, and also so far as the OCM is adopted, the duality property of the ground state is automatically taken into account, and hence there have been no difficulty in obtaining the
E0 transition strengths between the ground state and cluster states with reasonably good agreement with the observed data. What was done in Ref. [
79] was the clarification of the reason why the
E0 transition strengths between the ground state and cluster states have the magnitude comparable with the single-nucleon strength. The most important point is the fact that the
E0 transition is due to few degrees of freedom, namely the inter-cluster relative motion. Especially, for cluster states composed of two clusters, the number of degrees of freedom is only one. Another important point is the fact that the antisymmetrization operation by
A has only small effect. This point was proved by obtaining the explicit analytical formulas of the
E0 matrix elements between the ground state and cluster states. An example in the two-cluster system is given in Eq. (146).
In the above we explained that large E0 transition strengths comparable with the single-nucleon strength between the ground state and cluster states can be considered as being quite natural because of the duality property of the ground state. Conversely if some excited state has a large E0 transition strength to the ground state which is comparable to the single-nucleon strength, we can say that it is quite probable that this state is a cluster state if this state can be judged not to be a mean-field state having one-particle one-hole excited configuration from the ground state. Therefore large E0 transition strength is a good probe for identifying cluster states just like large E2 transition strength which is a good probe for identifying states with large quadrupole deformation.
6.2 Pauli-allowed states of many-cluster systems
6.2.1 Two-cluster Pauli-forbidden states and many-cluster Pauli-allowed states
In the container model, clusters make nonlocalized motion in containers under inter-cluster Pauli repulsion. The inter-cluster Pauli repulsion between clusters Cp and Cq originates from the Pauli-forbidden states defined by
where
is the internal wave functions of cluster
Cj and
is the relative coordinate between
Cp and
Cq. In the case of
system, the Pauli-forbidden states are harmonic oscillator wave functions
with
[
40]. The wave function of the relative motion between clusters
Cp and
Cq should be orthogonal to the Pauli-forbidden states. This orthogonality makes the relative wave function to avoid the spatial region where the Pauli-forbidden wave functions are existent. This is the effect of the repulsive force by the Pauli-forbidden states. A little more detailed explanation of the repulsive force effect of the Pauli-forbidden states is as follows. The orthogonality to the Pauli-forbidden states creates the oscillatory behavior of the relative wave function in the inner spatial region where the Pauli-forbidden wave functions are existent. The oscillatory behavior of the wave function means the large value of the kinetic energy. Therefore in order to avoid the large kinetic energy, the amplitude of the wave function becomes small in the inner spatial region. The small amplitude of the relative wave function in the inner spatial region is nothing but the repulsion between two clusters avoiding to come close to each other. The effective repulsive force coming from the orthogonality to the Pauli-forbidden states is called structural repulsive core [
160]. The structural repulsive core does not kill the inner-part wave function and so, if the attractive force is strong, the inner-part wave function can have large amplitude in order to gain the potential energy. Its difference from the rigid repulsive core is in many ways. For example at high scattering energy of clusters, high frequency oscillatory behavior of the relative wave function satisfies automatically the orthogonality to the Pauli-forbidden states, which means the fadeaway of the repulsive force effect. Another example which is very important is the weakening of the repulsive effect in three or many cluster systems where inter-cluster attractive force is strengthened by the Borromean effect. In the system of many clusters composed of
Cj (
–
n), every pair of clusters of
Cp and
Cq (
–
n) should have their relative wave function orthogonal to their Pauli-forbidden states
. It means that the relative wave function
of the many-cluster system should be orthogonal to all the Pauli-forbidden states of all pairs of clusters. Here
(
) are relative coordinates among clusters like Jacobi coordinates. In the case of Jacobi coordinates, we have
, where
and
Aj are center-of-mass coordinate and mass number of cluster
Cj, respectively.
In discussing the effect of the Pauli principle in the many-cluster system, the RGM norm kernel
[
44] is useful which is defined as
The eigen functions and eigen values of
are denoted as
and
, respectively, which satisfy [
44,
161]
The eigen functions constitute an orthonormal complete set and any wave function can be expanded by this complete set. The eigen functions with zero eigen value are Pauli forbidden because for , while the eigen functions with non-zero eigen value are Pauli allowed because for . The functional space spanned by all the Pauli-allowed eigen functions is designated as Sn.
Two-cluster Pauli-forbidden state is also a many-cluster Pauli-forbidden state because there holds . Then we have
From this relation we know that every Pauli-allowed eigen function
with
0 is orthogonal to all the two-cluster Pauli-forbidden state
. Therefore any Pauli-allowed state which is expressed by a linear combination of Pauli-allowed eigenfunctions
is orthogonal to all the two-cluster Pauli-forbidden states
. Now we consider the functional space
which is spanned by many-cluster wave functions that are orthogonal to all the two-cluster Pauli-forbidden state
. Clearly the space
contains in it the space
Sn. However explicit comparisons of
and
Sn in many kinds of many-cluster systems have shown that
and
Sn are equal to each other (
) [
44,
161,
162]. In the OCM for many-cluster system, the total wave function is required to be orthogonal to all the two-cluster Pauli-forbidden state
. Thus the functional space of many-cluster OCM is
. The fact that
ensures that the functional space of many-cluster OCM is the correct Pauli-allowed space
Sn spanned by all the Pauli-allowed eigenstates of the norm Kernel
.
6.2.2 Many-cluster Pauli-allowed state and evolution of cluster structure
In Table 16 we show the Pauli-allowed eigen states of the norm kernel of 4α system. Since the eigen states of the norm kernel have, as their quantum numbers, the total number of harmonic oscillator quanta N and SU(3) label , the allowed states are indicated by using N and . Their eigen values are given below their . The Pauli-allowed eigen state with lowest harmonic oscillator quanta has . This eigen state which we denote as has the following property,
This is of course the duality property of the double closed-shell wave function. When we denote the Pauli-allowed eigen state with as , we have the equality that is the same as the - shell-model wave function . This is the duality property of the 1p-1h shell-model wave function. The formation of this state by the excitation of inter-cluster wave function from the ground state is equivalent to the formation of this state by the excitation of mean-field degree of freedom from the ground state.
As we discussed in Section 4, in two-cluster system, the inter-cluster Pauli repulsion works so as to separate two clusters, which means that the inter-cluster Pauli repulsion works so as to separate spatial localization of clusters. The spatial localization-cluster Pauli repulsion does not necessarily generate spatial localization of clusters. A typical example is the Hoyle state in 12C which has the condensate-like gas structure of three α clusters with no deformation. We also know that it is highly probable that the sixth 0+ state at 15.1 MeV in 16O has condensate-like 4α gas structure. Moreover in general there can be variety of deformations in many-cluster systems. We can see this point in Table 16. Except the compact shell-model states for and 13, there are variety of Pauli-allowed states with different deformations ranging from prolate to oblate. For example, we see many Pauli-allowed states with SU(3) label corresponding to prolate deformation and also many Pauli-allowed states with SU(3) label corresponding to oblate deformation.
The most compact Pauli-allowed state has the lowest number of oscillator quanta
and expresses the shell-model wave function which has the duality property and is a dominant component of the ground state. The Pauli-allowed state with
often expresses a dominant component of a shell-model state with
–
character. According to the Ikeda diagram the cluster state which appears first as the excitation energy rises is two-cluster state. In general the cluster state composed of two or more clusters has the wave function which is a linear combination of Pauli-allowed state with many different
N. The two-cluster state has spatially localized clusters because of inter-cluster Pauli repulsion and its deformation is prolate. As the excitation energy goes up further higher, there appears three-cluster state. Three-cluster states do not necessarily have spatial localization of clusters as is typically seen in the case of the Hoyle state which has 3
α condensate-like structure. When a three-cluster state has a spatial localization of clusters, its shape or deformation can be both prolate and oblate. The shape or deformation of a three-cluster state is largely influenced by the orthogonality requirement to the states lying in the energy region lower than the three-cluster state. A recent study of this point is given in Ref. [
17] in the case of the three-
α states above the Hoyle state in
12C.
6.3 Container model and evolution of cluster structure
The so-called Ikeda diagram [
31] discusses the necessary condition of the excitation energy for the formation of the cluster state. It says that a cluster state appears around or above the threshold energy of the breakup into constituent clusters. The argument is that if the excitation energy is fairly below the threshold energy, attractive interaction between clusters is no more weak and clusters will lose their identity and dissolve into mean-field-like structure. Fig. 53 displays the Ikeda diagram given in Ref. [
31].
This diagram shows what kind of cluster structures are formed along the increase of excitation energy. This diagram may be understood such that they express the evolution path of cluster structure along the increase of excitation energy. In Fig. 54 a new diagram of the cluster formation is also shown, which is along the container picture discussed throughout this review paper. The path of cluster formation along the increase of excitation energy can be described by dilatation of the container, from spatially compact ground states to more dilute cluster structures. This point of view of cluster evolution, which was referred to as the container evolution, was demonstrated in the case of 16O, in Section 5.4. The concept of the container evolution was here generalized in other cluster formations appearing in the Ikeda diagram.
The evolution path of cluster structure along the increase of excitation energy, which can be identified with the container evolution process, is closely associated by the duality property discussed in the previous section. For example, in 12C, because of the duality property of the ground state, the activation of cluster degree-of-freedom leads to dilatation of container for 3α clusters and generates the 3α cluster states around or above the 3α breakup threshold energy including the Hoyle state. In 16O, because of the duality property of the ground state, the activation of 12C + α clustering degree-of-freedom leads to dilatation of container for clusters and generates the 12C + α cluster states around or above the 12C + α breakup threshold energy. Then the 12C cluster of the 12C + α cluster state evolves into 3α cluster structure due to the 12C duality property forming the 4α cluster states around or above the 4α breakup threshold energy. We see that each step of the evolution process of the cluster structure, firstly from the ground state to the 12C + α cluster states and next from the 12C + α cluster states to the 4α cluster states, can be understood according to the duality property of 16O and 12C. We here say that, in these nuclei, the evolution of cluster structure is governed by the duality path along the evolution path.
However, there exist evolution paths of cluster structure which contain non-duality path. A typical example is the formation of the 16O + 16O molecular resonance states. The dominant component of the ground state of 32S has the SU(3) symmetry and this component does not have the clustering degree of freedom of 16O + 16O configuration whose SU(3) symmetry is (N, 0) with . Thus the formation of the 16O + 16O molecular resonance states is not made only by the activation of the cluster degree of freedom from the ground state. As we already discussed in the previous Section 6.1.1, the 16O + 16O molecular resonance states are formed by the activation of the 16O + 16O cluster degree of freedom from the superdeformed mean-field state which are formed from the ground state by the activation of mean-field degree of freedom. Therefore the Ikeda diagram, as well as the container evolution diagram, does not describe the evolution path of cluster structure which involves non-duality path and does not show the mean-field-type excited state like the superdeformed mean-field state in 32S.
Special attention needs to be paid for the ground states of
8Be,
20Ne, and
44Ti which are the nuclei having two protons and two neutrons outside the double-closed shell nuclei
4He,
16O, and
40Ca, respectively. In the case of
8Be, as we already pointed out, the ground state is not a mean-field-type state but a cluster state. We can conjecture that the shell-model wave function
is spread into ground state and also into many excited states, and there is no excited state to which this shell-model wave function concentrates as a dominant component.
8Be in the Ikeda diagram indicates the existence of a very important change of nuclear structure, namely the formation of the nucleon mean field in the ground state is not energetically favourable in lighter nuclei than
12C. In the case of
20Ne, the ground state rotational band constitutes an inversion-doublet together with the negative-parity rotational band with
built upon the
at 5.80 MeV. The intrinsic state of the inversion-doublet bands is the
16O +
α cluster state. Thus it looks similar to the
8Be case. However in Refs. [
21,
31] the authors argue that the clustering feature is much weaker in the ground band than in the
band. The AMD calculation of
20Ne of Ref. [
163] supports the argument of Refs. [
21,
31]. For example the intrinsic wave function of the ground state obtained by AMD does not have density distribution of
16O +
α clustering while that of the
state of the
band shows clearly density distribution of
16O +
α clustering. Although the density distribution of the ground state does not show
16O +
α clustering, the component of
16O +
α cluster wave function contained in the ground state wave function has the squared amplitude of 70%, which is due to the duality property of the shell-model wave function,
. In Table 17 we show the properties of the levels of the inversion-doublet bands of
20Ne obtained by the AMD study of Ref. [
163]. We see that the squared amplitude of the
band component
WJ in the ground band decreases rather rapidly as the spin value
J increases from 70% of the 0
+ state to 28% of the 8
+ state. This decrease of the cluster character is due to the increase of the nucleon-spin alignment for larger
J [
164]. It can be seen in the increase of the expectation value of the two-body spin–orbit interaction
for larger
J. This feature of the spin-alignment well reflects the mean-field-type character of the ground bend of
20Ne. On the other hand, in the
band, the
WJ value of the 1
- state is as large as 95% and the smallest value of
WJ for the 9
- state is still as large as 70%. This result means that the
band has good clustering character of
16O +
α. It is supported by the small expectation values
for all the states of the
band.
In the case of
44Ti, just like the
20Ne case, there exists inversion-doublet bands composed of the ground rotational band and the
band built upon the
at 7.0 MeV. The negative-parity band levels were found by the
α-transfer reaction [
100]. In Table 18 we show the squared amplitude of the
40Ca +
α component
WJ of each level of the inversion-doublet bands of
44Ti obtained by the AMD study of Ref. [
165]. We see that the
40Ca +
α clustering component is very small in the ground band. Even the largest clustering percentage for the ground 0
+ state is 40%, which means that non-clustering component is 60% for the ground state. The 10
+ and 12
+ states have almost no clustering component. As for the
band, the largest clustering percentage is 56% for the 1
- state, which is fairly smaller than 95% of the 1
- state of
20Ne. The small clustering component in the inversion doublet bands is considered to be due to the strong spin–orbit potential and also to the strong attraction from the
40Ca core nucleus. Although the
40Ca +
α percentage of the ground band is not larger than 40%, the AMD calculation does not show the existence of
excited state whose
40Ca +
α component is larger than 40% and has relative-motion nodal number
. Here the relative-motion nodal number
nr means the nodal number of the reduced width amplitude
defined as
The relative-motion nodal number is the smallest nodal number of which is allowed by the inter-cluster Pauli principle. The value of nr of the of the band is also . The fact that the percentage of 40Ca + α component of the ground state is 40% means that the percentage of the SU(3) shell-model component Ca-core, is less than 40% because has the duality property that it is equivalent the cluster-model wave function of 40Ca + α. The decrease of the shell-model component with good spatial and SU(3) symmetry seen in the 44Ti ground state is considered to be more pronounced in the ground states of heavier nuclei.
7 Summary
In this report, we mainly focus on the new concept of nonlocalized clustering and the container picture for the evolution of cluster structure in nuclei, which can be another important fruition from the creative THSR idea in addition to the first important achievement of α condensation. In the recent five years, various cluster systems from nα nuclei, neutron-rich nuclei, to Λ-hypernuclei, are investigated with different THSR-type wave functions, which all provide the concept of nonlocalized clustering and the container picture. It was surprising that, even for the very old and historical subjects, like the α+16O clustering, α-linear-chain states, etc., there is still a large room for deepening the basic understanding of them.
We start with the famous Hoyle state of 12C. This novel 3α gas-like state triggers great interests in a wide field in physics and it has been studied for half a century by numerous models in nuclear physics. It came as a surprise that the THSR wave function with a simple structure gives a quite accurate description of this state and this Hoyle state is finally proved to be the α condensation state. Now, the search for the α condensation in 16O and 20Ne nuclei is becoming a very important subject both from theory and experiment. The great successful description of the gas-like states by the THSR wave function encourages us to extend this idea to the general cluster systems. By the study of the inversion doublet bands of 20Ne one significance underlying the THSR wave function is clarified, i.e., the clusters make a nonlocalized motion occupying the lowest orbit of the cluster mean-field potential instead of the localized motion from the traditional understanding. Subsequent calculations show that, this kind of nonlocalized motion, characterized by the size parameter in the THSR wave function, is not only the essential feature of the Hoyle state and typical α+ 16O cluster states but also of the alpha-linear-chain structures, molecular structures in neutron-rich nuclei like 9,10Be, and even cluster structures in hypernuclei.
The container picture was developed further for expressing this essential concept of nonlocalized clustering. In the container picture, the clusters are making almost independent motion in a nonlocalized way from a dynamical aspect, while the Pauli principle plays a key role for the formation of clusters from a kinematical aspect. The evolution of cluster structure along the increase of excitation energy can be naturally described by solving the Hill–Wheeler equation with the size parameter being the generator coordinate. This description can be identified with the container evolution, in which the size parameter characterizing the container expands along the path of the cluster evolution. We discussed the diverse clusters 0+ states in 12C and 16O. These evolution of cluster structures are well described and explained in our container picture in a unified way.
To investigate various complex nuclear systems in the container model, some new generations of the THSR wave functions and also new techniques are developed. For instance, the different widths of clusters can be treated exactly in an extended THSR wave function. The Monte Carlo arithmetic has been applied to deal with the complicated integration problem for multi-cluster systems. These extensions show the original THSR wave function has a very flexible character and it is becoming a powerful tool in nuclear cluster physics. In the future, we are extending this container idea to study the cluster structure of nuclei in an ab initio way.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2020