Nonlocalized clustering and evolution of cluster structure in nuclei

Bo Zhou , Yasuro Funaki , Hisashi Horiuchi , Akihiro Tohsaki

Front. Phys. ›› 2020, Vol. 15 ›› Issue (1) : 14401

PDF (57452KB)
Front. Phys. ›› 2020, Vol. 15 ›› Issue (1) :14401 DOI: 10.1007/s11467-019-0917-0
REVIEW ARTICLE
Nonlocalized clustering and evolution of cluster structure in nuclei
Author information +
History +
PDF (57452KB)

Abstract

We explain various facets of the THSR (Tohsaki–Horiuchi–Schuck-Röpke) wave function. We first discuss the THSR wave function as a wave function of cluster-gas state, since the THSR wave function was originally introduced to elucidate the 3α-condensate-like character of the Hoyle state (02+ state) of 12C. We briefly review the cluster-model studies of the Hoyle state in 1970’s in order to explain how there emerged the idea to assign the α condensate character to the Hoyle state. We then explain that the THSR wave function can describe very well also non-gaslike ordinary cluster states with spatial localization of clusters. This fact means that the dynamical motion of clusters is of nonlocalized nature just as in gas-like states of clusters and the localization of clusters is due to the inter-cluster Pauli principle which is against the close approach of two clusters. The nonlocalized cluster dynamics is formulated by the container model of cluster dynamics. The container model describes 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. When we notice that fact that the THSR wave function with the smallest value of size parameter is equivalent to the shell-model wave function, we see that the container model describes the evolution of cluster structure from the ground state with shell-model structure up to the gas-like cluster state via ordinary non-gaslike cluster states. For the description of various cluster structure, more generation of THSR wave function have been introduced and we review some typical examples with their actual applications.

Graphical abstract

Keywords

Hoyle state / alpha condensation / THSR wave function / nonlocalized clustering / container model / evolution of cluster structure

Cite this article

Download citation ▾
Bo Zhou, Yasuro Funaki, Hisashi Horiuchi, Akihiro Tohsaki. Nonlocalized clustering and evolution of cluster structure in nuclei. Front. Phys., 2020, 15 (1) : 14401 DOI:10.1007/s11467-019-0917-0

登录浏览全文

4963

注册一个新账户 忘记密码

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 (N=Z=18) 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 P=0.05 MeV/fm3 line. When the energy is low (E*/A=8 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 (E*/A=12–16 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 E*/A=12–16 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 0+ state which is known as the Hoyle state is the 3α-cluster gas-like state having the 3α-condensate-like structure [5, 6]. The excitation energy 7.65 MeV of this 02+ state is slightly above the 3α and 8Be(01+)+α 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 3α system has the form

(1)Φ3αTHSR(B)=N(B)A{exp[−2B2∑k=13(  X k−  X CM)2]∏i=13ϕ(αi)}.

Here Xk and XCM stand for the center-of-mass coordinates of the αk cluster and the total system, respectively. The intrinsic wave function of the kth cluster is represented by the ϕ(αk). 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 3α cluster model wave functions of the Hoyle state obtained by solving 3α 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]:

(2)|〈Φ3αTHSR|Φ3αRGM/GCM〉|2≈100%.

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];

(3)|〈Φ16O+αTHSR|Φ16O+αRGM/GCM〉|2≈100%.

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 ΦTHSR(B) which is just the size parameter of the container,

(4)∑B〈ΦTHSR(B′)|(H−Ek)|ΦTHSR(B)〉fk(B)=0.

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 3α-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 3α 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 3α 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 3α 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 3α THSR wave function namely that of the old 3α 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 03+ state of 12C. Next we review the research of the 4α condensate-like state in 16O which is assigned to the observed 6th 0+ state at 15.1 MeV. We explain how this assignment of the 4α condensate-like state is supported by theoretical calculations and their comparison with experiments. Theoretical calculations are those by 4α OCM (orthogonality condition model) [19] and by the GCM with 4α 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 3α and 4α systems. The solutions of the Hill–Wheeler equation especially in the case of the 4α 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 3α 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 nα cluster states in self-conjugate 4n nuclei are expected to appear around the nα 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 12C, whose excitation energy is 7.65 MeV, located at 0.38 MeV above the 3α 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 3α 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 3α breakup threshold and is well known as one of the mysterious 0+ 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 7.65 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 3α 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 8Be 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 8Be(0+)+α channel was successfully reproduced by semi-microscopic three-body calculation [12], called the Orthogonality Condition Model (OCM) [38-40]. This 3α OCM calculation gave the S-wave dominant structure between 8Be(01+) and α. Since 8Be(01+) is composed of loosely coupled two α clusters in relative S wave, the Hoyle state was concluded to have a loosely coupled 3α structure in relative S waves with large spatial extension, namely a gas-like structure of the 3α clusters, which is quite different from the 3α linear-chain structure. This understanding of the structure of the Hoyle state was confirmed, a few years later, by full microscopic 3α calculations, i.e., the 3α RGM calculation [14] and the 3α Brink-GCM calculations [13], where the Brink basis function is used [41]. The 3α Brink-GCM calculation with resonance boundary condition was also performed by Brussels group [42].

Both methods, the 3α RGM and 3α 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 3α RGM wave function is obtained by solving the following equation of motion,

(5)〈ϕ3(α)|(H−E)|A{χ3α( ξ 1, ξ 2)ϕ3(α)}〉=0,

where ϕ3(α)≡ϕ(α1)ϕ(α2)ϕ(α3) and ϕ(αi) is the internal wave function of the α cluster with the four nucleons inside it being constrained to (0s)4 harmonic oscillator configuration, and χ3α( ξ 1, ξ 2) is the relative wave function between the 3α 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 E0 and E2 transition properties. In Table 1 we show that the experimental data available in 12C are nicely reproduced by the 3α 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 8Be+α in 12C calculated with the 3α GCM, which is defined below:

(6)Yc=[L,l]J(r)=12!8!4!〈δ(r−ξ)r2[Φ8BeL,Yl(ξ^)]Jϕ(α)|Φ12CJ〉.

For the Hoyle state (02+ state), the S-wave component, which is denoted as [0,0], 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 4n nuclei in 2001 [8], and later it was extended so as to include deformation [45], in which the following operator form is assumed:

(7)|ΨnαTHSR〉=(Cα†)n|vac〉,

(8)Cα†=∫d  R exp[−∑k=x,y,zRk2βk2]∫d r1⋯d r4×​φ0s( r1​−​  R ,ν)aσ1τ1†( r1)⋯φ0s( r4​−​  R ,ν)aσ4τ4†(r4),

where aστ†(  r ) is a creation operator of one nucleon with spin σ, isospin τ at a spatial point r, and φ0s(  r −  R ,ν) is a 0s harmonic oscillator (H.O.) wave function centered at R with the form of φ0s(  r ,ν)=(2ν/π)3/4exp[−ν  r 2]. The coordinate representation of this THSR state |ΨnαTHSR〉 is expressed as

(9)​〈  r 1σ1τ1,⋯  r 4nσ4nτ4n|ΨnαTHSR〉​=​​∫d  R 1⋯d  R nexp​(​​−​∑i=1n∑kx,y,zRik2βk2​)​ΨnαB(R 1,⋯,  R n),

where ΨnαB(  R 1,⋯,  R n) is the Brink wave function in the nα cluster system [41]. The nα Brink wave function can be written in a form of Slater determinant and another equivalent form as shown below:

(10)ΨnαB(  R 1,⋯,  R n)=​det[φ0s(  r 1​​−​​  R 1,ν)σ1τ1⋯φ0s(  r 4n​​−​​  R n,ν)σ4nτ4n]

(11)=​A[φ0s(  X 1​​−​​  R 1,4ν) ​​⋯φ0s(  X n​​−​​  R n,4ν)ϕ(α1) ​​⋯ϕ(αn)],

where Xi and ϕ(αi)(i=1,⋯,n) are the center-of-mass (c.o.m.) coordinates and the internal wave function of the i-th α cluster, respectively, and ϕ(αi) is described by the (0s)4 H.O. shell-model wave function as follows:

(12)ϕ(αi)​∝​exp​[−​ν∑j=14(  r 4(i−1)+j​−​  X i)2]​σ4i−3τ4i−3⋯σ4iτ4i,

and A 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

(13)ΨnαTHSR=〈  r 1σ1τ1,⋯  r 4nσ4nτ4n|ΨnαTHSR〉 =​A{∏i=1n∫​d  R iexp(−​​∑k=x,y,zRik2βk2)​φ0s(  X i​​−​​  R i,4ν)ϕ(αi)} ∝​A{∏i=1nexp(−​∑k=x,y,z2Xik2Bk2)​ϕ(αi)} =​g(XG)ΦnαTHSR,Bk2=b2+2βk2, b=1/2ν

with

(14)ΦnαTHSR∝A{∏i=1nexp[−∑k=x,y,z2(Xik−XGk)2Bk2]ϕ(αi)}=A{∏i=1n−1exp(−μi∑k=x,y,z2ξik2Bk2)​∏i=1nϕ(αi)},

and g( XG)∝exp(−2n∑k=x,y,zXGk2/Bk2), where   ξ i=  X i−∑j=i+1n  X j/(n−i) and μi=(n−i)/(n−i+1) for i=1,⋯,n−1, and XG=∑j=1n Xj/n. g( XG) and ΦnαTHSR 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, exp{−2∑k=x,y,z(Xk−XGk)2/Bk2} with the size parameter Bk (k=x,y,z), under the condition that the effect of the antisymmetrizer A is negligible. This condition is satisfied when the size parameter Bk is much larger than the size parameter of the α cluster, b, or equivalently βk≫0. 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 βk 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 A, so that there only remains the α-cluster-like correlation. Furthermore in the limit of Bk→b(k=x,y,z), or βk→0(k=x,y,z), 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 ΨnαB( R1,⋯,  Rn) coincides with the shell-model Slater determinant in the limit of R i→0(i=1,⋯,n). By using an alternative expression of the delta function, limβ→0(1/πβ2)exp(−R2/β2)=δ(R), this limiting operation βk→0(k=x,y,z) in Eq. (9) gives the Brink wave function ΨnαB( R1,⋯,  Rn) in the limit of R i→0(i=1,⋯,n), 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 α+16O inversion doublet bands of 20Ne 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 4n-body microscopic wave function so that the 4n nucleons should be explicitly antisymmetrized by the operator A. 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),

(15)ΨnαTHSR=∫dRGexp[−n∑kx,y,z(RGk)2βk2]×φ0s( XG−RG,4nν)∫d R˜1⋯dR˜n−1×exp[−∑i=1n−1μi∑k=x,y,zR˜ik2βk2]ΦnαB( R˜ 1,⋯,  R˜n−1),

where { R˜ } is the set of Jacobi coordinates transformed from the position-parameter set { R }, with the relation, R˜ i=  Ri−∑j=i+1n Rj/(n−i) for i=1,⋯,n−1 and RG=∑j=1n Rj/n, and ΦnαB( R˜ 1,⋯,  R˜ n−1) is the internal part of the nα Brink wave function defined by

(16)ΦnαB(R˜1,⋯​,  R˜n−1)​​=​​φ0s−1( XG​​−​​  R˜n,4nν)ΨnαB( R1,⋯​,  Rn).

While in Eq. (15) the first half which is in the first line is nothing but the spurious total c.o.m. wave function, g( XG) 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.,

(17)ΦnαTHSR=∫d R˜1⋯d R˜n−1exp[−∑i=1n−1μi∑kx,y,zR˜ik2βk2]×ΦnαB( R˜1,⋯,  R˜n−1).

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 ΦnαB( R˜1,⋯, R˜n−1) in the above equation deviates from the pure Slater determinant, ΨnαB( R1,⋯, Rn), 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 O, with a pure Slater determinant, like 〈ΨB|O|ΨB′〉, the matrix element to be desired, 〈ΦB|O|ΦB′〉, is also easily derived, since the transformation matrix U:{ R }→{ R˜ } 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 ΦnαTHSR 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:

(18)ΦnαJM(β )=NPMKJΦnαTHSR=N∫dΩDMKJ∗(Ω)A{∏i=1n−1exp[−μi∑kx,y,z2(R^(Ω)  ξ i)k2Bk2]∏i=1nϕ(αi)}=N∫dΩDMKJ∗(Ω)∫d R˜1⋯d R˜n−1exp(−∑i=1n−1μi∑kx,y,zR˜ik2βk2)ΦnαB(R^−1(Ω)  R˜1,⋯,R^−1(Ω)  R˜n−1),=N∫dΩDMKJ∗(Ω)∫d R˜1⋯d R˜n−1exp[−∑i=1n−1μi∑kx,y,z(R^(Ω)  R˜ i)k2βk2]ΦnαB( R˜1,⋯,  R˜n−1),

with β ≡(βx,βy,βz), Ω, R^(Ω), DMKJ(Ω), and N, 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 O for the THSR wave function ΦnαTHSR,(JM), i.e., 〈ΦnαJM|O|ΦnαJM〉, one should first calculate the matrix element for the intrinsic Brink wave function ΦB, i.e., 〈ΦB|O|ΦB〉. This is basically represented as a form of Gaussians. Thus integrating over the set of variables { R˜ } 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: β⊥​≡​βx​=​βy​≠​βz. The above relation of Eq. (18) is simplified and ΦnαJM can be redefined as

(19)ΦnαJM(β⊥,βz)=N∫dcosθdM0J(θ)∫d R˜1⋯d R˜n−1 ×exp[−∑i=1n−1μi∑kx,y,z(R^y(θ) R˜i)k2βk2]ΦnαB(R˜1,⋯, R˜n−1),

with a redefined normalization constant N.

The energy eigenstates can be represented by superposing the above wave functions, Eq. (19), as follows:

(20)Φnα,λJM=∑β⊥,βzfλJ(β⊥,βz)ΦnαJM(β⊥,βz).

The coefficients fλJ can then be determined by solving the following Hill–Wheeler equation,

(21)∑β⊥′,βz′〈ΦnαJM(β⊥,βz)|(H​−​E)|ΦnαJM(β⊥′,βz′)〉fλJ(β⊥′,βz′)​=​0,

where the microscopic Hamiltonian H=T−TG+∑i<jAVij(NN)+∑i<jAVij(C) is composed of the kinetic energy with the c.o.m. kinetic energy subtracted, T-TG, nucleon-nucleon effective interaction Vij(NN), and Coulomb interaction Vij(C). 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 3α RGM and 3α GCM were performed for 12C 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 8Be(0+)+α 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 3α-cluster structure, which are loosely coupled in relative S-waves, rather than the 3α-linear-chain structure. In the 3α RGM, or equivalently in the 3α 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 3α RGM wave function can be expressed in the following form:

(22)Ψ3αRGM=A[χ3α(  ξ1,  ξ2)ϕ(α1)ϕ(α2)ϕ(α3)],

where the Jacobi coordinates   ξ i= Xi−∑j=i+12 Xj/(n−i) for i=1,2. On the other hand, as we can easily see from the r-space representation of the 3α THSR wave function in Eq. (13) for n=3, the relative wave function between the 3α clusters in the 3α THSR wave function is written like

(23)χ3αTHSR( ξ1,  ξ2)∝PM=K=0J=0×exp[−2∑i=13{(Xix−XGx)2+(Xiy−XGy)2b2+2β⊥2+(Xiz−XGz)2b2+2βz2}]=PM=K=0J=0exp[−2∑i=12μi{ξix2+ξiy2b2+2β⊥2+ξiz2b2+2βz2}].

Obviously the special form is assumed for the relative wave function between the 3α 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 3α RGM/GCM wave function that is the full solution of the microscopic 3α-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 12C 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 3α 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 M=0.575 is adopted and in the latter a slightly modified version of Volkov No. 2 with M=0.59 [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 12C. 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 0.7 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 3α 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 3α threshold is expected to have, than the GCM wave function, in which the finite number of spatially localized 3α-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 2α+α-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 12C 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, q≥2.5 fm−1, 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 β⊥,βz, with the Force 2,

(24)E(β⊥,βz)=〈Φ3αJ=0(β⊥,βz)|H|Φ3αJ=0(β⊥,βz)〉,

In Fig. 6(a), we show the energy contour map in the two-parameter space β⊥ and βz for the above quantity. A minimum appears at a spherical position, (β⊥,βz)=(1.5, 1.5 fm). 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, Φ3α,λ=1J=0 is calculated to be |〈Φ3αJ=0(β⊥=1.5, βz=1.5 fm)|Φ3α,λ=1J=0〉|2=0.93. This means that the single THSR wave function with a rather small value of the parameter (β⊥,βz)=(1.5, 1.5 fm) very well approximates the ground state wave function Φ3α,λ=1J=0, though the minimum energy −87.68 MeV is about 1.8 MeV higher than the energy −89.52 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:

(25)E⊥(β⊥,βz)=〈P^⊥g.s.Φ3αJ=0(β⊥,βz)|H|P^⊥g.s.Φ3αJ=0(β⊥,βz)〉〈P^⊥g.s.Φ3αJ=0(β⊥,βz)|P^⊥g.s.Φ3αJ=0(β⊥,βz)〉,

P^⊥g.s.=1−|Φ3α,λ=1J=0〉〈Φ3α,λ=1J=0|, 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 β⊥,βz. We see an energy minimum at (β⊥, βz)=(5.2, 1.5 fm) in the oblate region of the map and a second energy minimum at (β⊥,βz)=(2.6, 7.5 fm) in the prolate region. The minimum energy value is −81.75 MeV. This value is almost the same as the energy of −81.79 MeV given by solving the Hill–Wheeler equation. The minimum energy −81.75 MeV is close to the second minimum energy of −81.67 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 0.9. 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

(26)Φ3α,⊥J=0(β⊥,βz)​​=​P^⊥g.s.Φ3αJ=0(β⊥,βz)​〈P^⊥g.s.Φ3αJ=0(β⊥,βz)|P^⊥g.s.Φ3αJ=0(β⊥,βz)〉.

We then show in Fig. 6(c) the contour map of the squared overlap surface defined below:

(27)O(β⊥,βz)=|〈Φ3α,⊥J=0(β⊥,βz)|Φ3α,λ=2J=0〉|2,

where Φ3α,λ=2J=0 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 0.993 at (β⊥,βz)=(5.3, 1.5 fm). This means that the Hoyle state wave function obtained by solving the Hill–Wheeler equation is almost identical to the single THSR wave function Φ3α,⊥J=0(β⊥=5.3, βz=1.5 fm). Since we keep in mind that the wave function Φ3α,λ=2J=0 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 3α 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 ab−initio 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 12C 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 3α 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 22+ state in 12C, and discuss the physical observables of the states. The calculated r.m.s. radius of the Hoyle state is 2.93 fm, which is sizably larger than that of the ground state (2.4 fm). The monopole transition strength is also calculated, which is 2.1 efm2, compared with the experimental value, 5.3 efm2 (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 8Be 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 v18 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 12C, 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 E0 transition form factor is compared to the experimental data, which leads to M(E0)=5.29±0.14 efm2 at k2=0, 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 5.36 efm2 and 6.45 efm2 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, 10.4 versus 7.65 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 3.0−3.5 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 0+ state above the Hoyle state, at 10.3 MeV with a width of 2.7 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 Jπ=2+ state around 10 MeV almost 40 years ago by using the microscopic cluster models [13, 14, 42]. Amazingly the Jπ=2+ state has been observed recently, as a different state from the broad 10.3 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 0+ peak is decomposed into two peaks, giving the 03+ and 04+ states at 1.77 MeV and 3.29 MeV above the 3α threshold, with the widths of 1.45 MeV and 1.42 MeV, respectively [69]. Theoretically it is very difficult to obtain both 0+ 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 0+ states [70, 71]. Though the interpretation of these states is not yet converged, the result of this study suggests that the 03+ state has a higher nodal structure between 8Be+α in relative S-wave [70]. Also the AMD and FMD calculations suggest that the 04+ state has a bent-armed shape of the 3α clusters, resembling the 3α linear-chain configuration, though in these calculations the 03+ state is missing [51, 77].

The THSR ansatz could also reproduce the 03+ and 04+ 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 2α+α correlation [50], as follows:

(28)Φ3αJ=0( β1,  β2)​​=​​P^J=0A{∏i=12exp​​[−​μi​∑k=x,y,z​2ξix2Bik2]​∏i=13​ϕ(αi)},

where μi=i/(i+1) and Bik2=b2+2βik2 with i=1,2 and k=x,y,z, and P^J=0 is the angular-momentum projection operator onto J=0. The degree of freedom of the cluster motions characterized by β  is divided into β 1 and β 2, corresponding to 2α (8Be) and α, respectively. Thus the case β1=  β 2 comes back to the original THSR wave function, and this natural extension allows for including the 2α+α 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, β1=(β1x=β1y,β1z) and β2=(β2x=β2y,β2z) 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 03+ and 04+ states. We can see that the 03+ and 04+ states are reasonably reproduced. Amazing feature of the 03+ state is the large monopole transition to the Hoyle state, whose strength is around 47 efm2. The r.m.s. radius of this state is calculated to be much larger than the Hoyle state. These results suggest that the 03+ state is excited from the Hoyle state by monopole-like dilatation.

We should also emphasize that for the formation of the 03+ 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 P^⊥=1−|01+〉〈01+|−|02+〉〈02+|. Second we construct the THSR wave function in the orthogonal space, like Ψ⊥(  β 1,  β 2)=P^⊥Φ3αJ=0(  β 1,  β 2). We then calculate the energy of this THSR state, which is defined as

(29)E⊥=〈Ψ⊥(  β 1,  β 2)|H|Ψ⊥(  β 1,  β 2)〉〈Ψ⊥(  β 1,  β 2)|Ψ⊥(  β 1,  β 2)〉.

In Fig. 10, the contour map of the energy surface in two-parameter space β1x=β1y=β1z and β2x=β2y=β2z is shown. In construction of the projection operator P^⊥, optimal single configurations in the THSR states are taken as the 01+ and 02+ 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 03+ 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 03+ 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 03+ state is excited from the Hoyle state with the strong monopole strength, the following quantities are calculated:

(30)Y(a)​=​12!4!8!〈δ(ξ2​−​a)ξ22[Φ2α,Y0( ξ^2)]0ϕ(α)|ΦGCM(0λ+)〉,

(31)χ(a)​=​12!4!4!4!〈δ(ξ1​−​a)ξ12Φ2α−αY0( ξ^1)|ΦGCM(0λ+)〉,

where Φ2α=N0P^J=0A[exp{−∑k=x,y,zξ1k2/Bk2}ϕ2(α)] and Φ2α−α=N˜0P^J=0A[exp{−∑k=x,y,zξ2k2/B2k2}ϕ3(α)], with N0 and N˜0 normalization factors, and Bk2=b2+βk2 and B2k2=3/4b2+β2k2. ξ1 and ξ2 are the Jacobi coordinates. The upper equation is nothing but the RWA of 8Be+α channel, where Bx, By, and Bz in Φ2α is best parameterized to describe 8Be nucleus. The lower equation is defined as representing the correlation of the 2α clusters associated with the coordinate   ξ 1.

In Fig. 11(a), the RWAs in the 8Be(0+)+α(S) channel for the 01+–04+ states are shown. The RWAs for the ground state and the 03+ 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 8Be 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 8Be core, and hence formation of the loosely coupled 3α cluster state.

On the contrary, the 03+ state recovers the distinct nodal behavior and, with one additional node, forms a higher nodal 8Be(0+)+α structure. Obviously the relative motion between the 8Be 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 2α part associated with   ξ 1 is shown, where the relative motion between the 2α and α is chosen to best simulate the 03+ state. In this calculation, this part is taken as Φ2α−α with β2x=β2y=13.8 fm and β2z=13.7 fm, since the 03+ state obtained by solving the Hill–Wheeler equation has the largest squared overlap (around 30%) with the single THSR wave function for this β2 value. While for the 03+ state the two nodes in inner region remain, which come from the repulsive core between the 2α clusters, an oscillatory behavior is also clearly seen in outer region more than 4 fm. This indicates that the 2α relative motion associated with ξ1 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 O(E0,12C)=∑i=112(ri− XG)2=∑k=13∑i∈αk(  r i−  Xk)2+2ξ12+8ξ22/3 [79]. The strong monopole transition strength from the Hoyle state can be from the two relative motions associated with the coordinates ξ1 and ξ2, giving the picture that the 03+ state is featured as monopole dilatation from the Hoyle state with a vibrational mode.

2.2 4α condensate-like state in 16O

2.2.1 4α OCM calculation

We discuss in this section about the Hoyle-analog state in 16O calculated by the 4α OCM, in which the 4α 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 4α-cluster system are written in the following way, respectively,

(32)〈ϕ4(α)|(H−E)|A[χ4α( ξ1,  ξ2,  ξ3)ϕ4(α)]〉=0.ϕ4(α)≡ϕ(α1)ϕ(α2)ϕ(α3)ϕ(α4).

and

(33)Ψ4αRGM( r1,⋯,  r16)=A[χ4α( ξ1,  ξ2,  ξ3)ϕ4(α)],

where ξ1, ξ2, and ξ3 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,

(34)∫dξ ′[H(ξ,  ξ ′)−EN(ξ, ξ ′)]χ4α(ξ′)=0,

with

(35){H(ξ,  ξ′)N(ξ,  ξ′)}​=​〈δ(ξ​−η )ϕ4(α)|{H1}|A[δ( ξ ′−η )ϕ4(α)]〉,

where ξ , ξ′, and  η  denote the sets  ξi,  ξi′, and ηi, (i=1,2,3), respectively, which are all the Jacobi coordinates between the 4α 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 4α clusters and Jacobi coordinates ηi(i=1,2,3) between the 4α 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 H and N. This Eq. (35) can further be transformed into the following form:

(36)(N−1/2HN−1/2−E)ΦJ(B)=0,

where Φ(B) can be defined as a bosonic wave function of the 4α clusters with only the Jacobi coordinates ξ left, as follows:

(37)ΦJ(B)(ξ1, ξ2, ξ3)=N1/2χ4α=∫​dξ′N1/2(ξ, ξ′)χ4α(ξ′).

In the OCM, the following approximation with respect to the treatment of the antisymmetrization in these kernels is taken,

(38)(N−1/2HN−1/2−E)ΦJ(B)=0⇔Λ(N−1/2HN−1/2−E)ΛΦJ(B)=0⇒Λ(H(OCM)−E)ΛΦJ(OCM)=0,

where Λ is a projection operator onto Pauli allowed space [40, 44] and H(OCM) is the effective local Hamiltonian of the OCM written below:

(39)ΛH(OCM)Λ=∑i4Ti−Tcm+∑i<j4[V2α(N)(i,j)+V2α(C)(i,j) +V2α(P)(i,j)]+∑i<j<k4V3α(i,j,k)+V4α(1,2,3,4),

where Ti, V2α(N)(i,j), V2α(C)(i,j), V3α(i,j,k) and V4α(1,2,3,4) 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 Tcm is subtracted from the Hamiltonian. V2α(P)(i,j) is the Pauli exclusion operator, by which the Pauli forbidden states between the 2α clusters in 0S, 0D, and 1S 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 12C and 16O 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 Jπ=0+ spectrum with the calculated ones by using the 4α 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 0+ states calculated by the 4α OCM and the four 0+ states calculated by the 4α THSR ansatz as (01+)OCM−(06+)OCM and the (01+)THSR−(04+)THSR, respectively. We can see that the 4α OCM calculation gives a satisfactory one-to-one correspondence with the experimental spectrum. The low-lying states (01+–03+) of this spectrum are actually in good agreement with earlier 12C+α OCM calculations [85, 86] (see also Ref. [28]). We should emphasize that the orthogonality to the low-lying 12C+α states and the ground state is indispensable for the formation of the 4α-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 Rrms, monopole matrix element to the ground state M(E0), and α-decay width Γ between the experimental data and the results of the 4α OCM and THSR ansatz is shown. The α-decay width is calculated based on the R-matrix theory [28, 87]. The (06+)OCM state has the largest r.m.s. radii among the calculated states and this large value 5.4 fm indicates that this state is the 4α gas-like state and considered to be the Hoyle-analog state. We can see that the observed widths of the experimental 04+, 05+, and 06+ states are nicely reproduced by the 4α OCM calculations. In particular, the observed width of the 06+ state at 15.1 MeV is very well reproduced by (06+)OCM state, and hence the Hoyle-analog state corresponds to the 06+ state at 15.1 MeV.

We then investigate the Hoyle-analog state by calculating the RWAs for the 12C+α channels, which are defined by

(40)Y[L,L]0(r)=4!3!1!×〈[δ(r′−r)r′2[ΦL(OCM)(12C),YL(r^′)]0ϕ(α)]|ΦJ=0(OCM)(16O)〉,

with ΦL(OCM)(12C) and ΦJ=0(OCM)(16O) being the OCM states of 12C and 16O, respectively. We show in Fig. 13 the RWAs of the (06+)OCM 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 4α condensate state must have a large overlap with the 12C(02+)+α structure.

The α-condensate character of the (06+)OCM is also clearly seen in the momentum distribution of the c.o.m. motion of α-clusters, which is defined by

(41)ρ(k)=1(2π)3∫drdr′ρ˜(B)(r, r′)e−i k ⋅(r−r′),

where ρ˜(B)( r , r′) 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 (01+)OCM−(06+)OCM states. Only the one of the (06+)OCM state has a markedly sharp peak near zero momentum compared with the other 0+ states. In Fig. 14(b), we show for comparison the momentum distribution of the Hoyle state by 3α OCM (dotted curve), which is also sharply peaked near zero momentum.

In addition to the α condensate state, we also discussed the other 0+ states. The structures of the 02+ and 03+ states are well established as having the 12C(01+)+alpha(S) and 12C(21+)+alpha(D) cluster structures, respectively. These structures of the 02+ and the 03+ states are confirmed in the present 4α OCM calculation. The present calculations also show that the 04+ and 05+ states mainly have 12C(01+)+α(S) structure with higher nodal behavior and 12C(1−)+α(P) structure, respectively [28].

2.2.2 4α THSR calculation

The Hoyle-analog state in 16O is also investigated with the 4α THSR wave function. While in the 4α OCM discussed in the previous subsection the 4α problem is fully solved without any model assumption with respect to the α+α relative motions within the semi-microscopic treatment, in the 4α THSR the THSR-type constraint is imposed on the α+α relative motions. In this calculation, we only consider the spherical case, i.e., β0≡β⊥=βz, 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 (01+)THSR–(04+)THSR [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 (03+)OCM and (05+)OCM states are shown to have a large fraction of S-factor from the 12C(21+)+α(D) and 12C(11−)+α(P) 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 12C+α cluster states, it is necessary to further extend the THSR ansatz in Eq. (14) so as to separate the degree of freedom between the 12C core and the remaining α cluster, by assigning different parameter sets of β  to the first 3α’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 03+ and 05+ states to be completely missing and the (02+)THSR and (03+)THSR states to be averaged representations of the 02+ and 04+ states, respectively. We can therefore say that the (04+)THSR state corresponds to the 06+ state with the 4α condensate character. In fact, in Table 3 we can see that the (04+)THSR 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 4α OCM.

We also calculate the RWAs for the 12C+α channels in the THSR ansatz, which are defined by

(42)Y[0,0]0(r)=16!12!4!〈[δ(r′−r)r′2[Φ3α,μL=0,Y0( r^′)]J=0ϕ(α)]|Φ4α,λJ=0〉,

with Φ3α,μJ=0 and Φ4α,λJ=0 being the wave functions of the 12C(0μ+) and 16O(0λ+) states, respectively, calculated with the (spherical) THSR ansatz in Eqs. (20) and (21). We show in Figs. 15(a)–(d) the amplitudes for the (01+)THSR−(04+)THSR, respectively, in the 12C(01+)+α and 12C(02+)+α channels. The ones of the (04+)THSR well correspond to the amplitudes for the (06+)OCM shown in Fig. 13. This again indicates that the (04+)THSR state is the counterpart of the (06+)OCM state and hence the Hoyle-analog state. We can also see that the (02+)THSR and (03+)THSR states have large amplitudes from the 12C(01+)+α channel and their peak positions are in outer region, beyond 5 fm, which are also consistent with the results of the 4α OCM calculation (see Ref. [28]).

In the same way as the 4α OCM, we calculate the momentum distribution of the α clusters in the THSR ansatz. The 4α THSR wave function is, however, the fermionic wave function of the 16 nucleons. We therefore have to extract a bosonic degree of freedom from Φ4αTHSR 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],

(43)Φ(B)=N1/2χ4αTHSR≈Nχ4αTHSR〈Nχ4αTHSR|Nχ4αTHSR〉,

with

(44)χ4αTHSR(ξ1,ξ2,ξ3)​​=​​∑βfλJ=0(β)exp[​​−​2​∑i=13μi{ξi2b2​+​2β2}]

where fλJ=0(β) 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 0λ+ states with λ=1,⋯,4 are shown. We can see that the 04+ state has a sharp peak near zero momentum and its shape is in fairly good agreement with the one of the (06+)OCM state in Fig. 14, where both curves have a peak height around 2.0 fm and width around 0.5 fm−1.

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 nα 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 B→b (or β →0). 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 3α THSR wave function as an example,

(45)limB→bΨ3α,intTHSR​=limB→bC0(B)exp(−12νXG2)×A{exp[−(1B2ξ12+43B2ξ22)]ϕ3(α)}

(46)​=​c0A{R2,0(ξ1,1b2)R2,0(ξ2,43b2)ϕ3(α)}

(47)​=​Φshell((0s)4(0p)8;[444](λ,μ)​=​(0,4),J​=​0),

where ξ1=X2-X1 and ξ2=X3-(X1+X2)/2. C0(B) and c0 are normalization constants. RN,L(r,γ0) is the radial harmonic oscillator function with size parameter γ0 with and N stands for the number of harmonic oscillator quanta, N=2n+L. The obtained Eq. (47) shows that the 3α THSR wave function at the limit of B→b 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 12C [9, 11].

Furthermore, if the composed clusters of a nuclear system are closed-shell nuclei like α and 16O, the THSR wave function for this cluster system also has the simple shell-model limit. In the 16O+α system, after angular-momentum projection with even L, the THSR wave function can be written, exp(−20νXG2)A{rLexp(−γr2)YLM(r^)ϕ(16O)ϕ(α)}. The shell-model limit of this THSR wave function can be expressed as follows [11]:

(48)limB→bΨ20Ne,LTHSR=limB→bCL(B)exp(−20νXG2)×A{rLexp(−85B2r2)YLM(r^)ϕ(16O)ϕ(α)}

(49)=cLA{R8,L(r,85b2)YLM(r^)ϕ(16O)ϕ(α)}

(50)=Φshell((0s)4(0p)12(0d1s)4;[4](λ,μ)=(8,0),LM),

where CL(B) and cL are normalization constants. The derived Eq. (50) shows that at the limit of B→b, the THSR wave function becomes the most important sd shell-model wave function with spatial symmetry [4] and SU(3) symmetry (λ,μ)=(8,0) [92]. This point is very important for the description of shell-model-like states in 20Ne 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 B→+∞, 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 0+ THSR wave function with 0+ THSR wave functions with various deformations in 20Ne [11]. The deformed THSR wave function ΦNe of 20Ne has the form A[χ(  r )ϕ(α)φ(16O)], where χ(r) is exp[−∑k=x,y,z(8/5Bk2)rk2] and Bk2=b2+2βk2. We see in this figure that the prolate THSR wave function with βx = βy = 0.9 fm, βz = 2.5 fm is almost 100% equivalent to oblate THSR wave functions with βx=βy≈2.1 fm, βz≈0∼1.2 fm after angular-momentum projections onto the 0+ state. The equivalence of prolate and oblate THSR wave functions after angular-momentum projection is also true for other the spin-parity states of 20Ne.

In spite of the equivalence of prolate and oblate THSR wave functions after angular-momentum projection, we can say that the 20Ne 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 20Ne states expressed by THSR wave functions have negative signs. From the well-known formula Q(J)=−[J/(2J+3)]Q (intrinsic), we know that when the expectation value Q(J) 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 ΦNeJ=A[χJ( r)ϕ(α)ϕ(16O)] and we can prove that this type of wave function ΦNeJ gives us the following formula for Q(J),

(51)Q(J)=−J2J+3165〈r2〉,

(52)165〈r2〉=〈ΦNeJ|∑j=120( rj− XG)2|ΦNeJ〉​−​R2(16O)​−​R2(α),

(53)R2(Ck)=〈ϕ(Ck)|∑j∈Ck( rj− XG(Ck))2|ϕ(Ck)〉.

This formula shows that Q(J) has negative value and explains why the calculated values of Q(J) by THSR wave functions have all negative signs. Of course the negative sign of Q(J) 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 Φave(βx=βy;βz) generated from a prolate THSR wave function Φpro(βx=βy;βz) is defined as follows:

(54)Φave(βx=βy;βz)​=​(12π∫02πdθeiθJx)Φpro(βx​=​βy;βz).

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 0+ state, when we construct the rotation-average wave function Φave(βx=βy=0.9;βz=2.5 fm) from the prolate THSR wave function with (βx,βy,βz)=(0.9,0.9,2.5 fm) which gives the minimum energy for 0+, it is almost 100% equivalent to the oblate THSR wave function with (βx,βy,βz)=(0.9,2.1,2.1 fm),

(55)Φobl(βx=0.9,βy=βz=2.1fm)  ​=Φave(βx=βy=0.9,βz=2.5 fm).

When we construct the good-spin wave function from Φave(βx=βy;βz), it is clear that the resulting wave function is the same as the good-spin wave function constructed from Φpro(βx=βy;βz). Therefore in the above example, the oblate THSR wave function Φobl(βx=0.9;βy=βz=2.1 fm) is almost 100% equivalent to the prolate THSR wave function Φpro(βx=βy=0.9;βz=2.5 fm) after angular momentum projection to the 0+ 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 20Ne case. We proposed a new type of microscopic cluster wave function [10], which we call hybrid-Brink–THSR wave function,

(56)Φcluster( βi,  Si)=∫d R 1…d Rnexp[−∑i=1n∑k=x,y,zRik2βik2]ΦclusterB(R1+S1,…,  Rn+ Sn)

(57)∝A{∏i=1nexp[−Ai∑k=x,y,z(Xik−Sik)22Bik2]ϕ(Ci)},

(58)ΦclusterB(S1,…,  Sn)∝A{∏i=1nexp[−Ai( Xi− Si)22b2]ϕ(Ci)}.

Here Bik2=βik2+2b2 and Xi and ϕ(Ci) 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 βi. The oscillator parameter of the cluster Ci is called b, which is taken as the same size value for different clusters. ΦclusterB 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 βi or Bi becomes the size parameter. When βik=0, i.e., Bik=b(k=x,y,z), 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]

(59)Φn(rn)​=∫dRnexp(​−Rn,x2βn,xy2−Rn,y2βn,xy2−Rn,z2βn,z2)eimϕRn×(πb2)−3/4e−(rn−Rn)22b2.

When we change the integration variables (Rn,x,Rn,y,Rn,z) to (−Rn,x,−Rn,y,−Rn,z), in the integral representation of Φn(−rn) by Eq. (59), we obtain Φn(−rn)=−Φn(rn). It is because the azimuthal angle of (−Rn,x,−Rn,y,−Rn,z) is (π+ϕRn) and exp[i(π+ϕRn)]=−exp(iϕRn). When m=0, 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,

(60)π=πα(1)×πα(2)×πn={+,m=0,−,m=±1.

From Eq. (59) we can also see that the z component of orbital angular momentum lz,n=m is a good quantum number for the extra neutron. Because of the rotational symmetry of the two-α-cluster subsystem about z-axis, we have lz,α=0 for α clusters. Thus the z-component of the orbital angular momentum lz of total system is

(61)lz=lz,α(1)+lz,α(2)+lz,n=m.

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

(62)Φl(λ)(R)=∑BCB(λ)(R)Φl(B,R)

with

(63)Φl(B,R)=A[Φ3αTHSR(B)φl( rn−X3α−R )],

where Φ3αTHSR(B) is the spherical THSR wave function for 12C with B=Bx=By=Bz in Eq. (14) and φl( r )=(2ν/π)3/4exp(−νr2)Ylm( r^ ), with rn and X3α being the coordinates of the valence neutron and the c.o.m. coordinate of the 3α system, respectively. The coefficient of the superposition in Eq. (62), CB(λ)(R), and the corresponding energy eigenvalue El(λ)(R) can be obtained as a function of the 3α-n relative distance parameter R by solving the following Hill–Wheeler-type equation,

(64)∑B′〈Φl(B,R)|H−El(λ)(R)|Φl(B′,R)〉CB′(λ)(R)=0.

The adiabatic energy curves El(λ)(R) with l=0,1 for the lowest (λ=1) and the second lowest (λ=2) eigenstates are shown in Fig. 19. It is noted that the cases of λ=1 and λ=2 correspond to the ground state and the Hoyle state for the core 3α system, respectively, and hence demonstrate, respectively, that the valence neutron couples with the ground state and the Hoyle state in l=0 or l=1. 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 l=0 is energetically favored for both λ=1,2 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 El=0(λ=1)(R) becomes higher than El=1(λ=1)(R) inside the surface region of the ground state of 12C, i.e., for R less than around 3 fm. At the center of the core, the neutron sits in the 0p1/2 orbit rather than in the 1s1/2 orbit.

However, a more interesting result is for the λ=2 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 El=0(λ=2) goes higher than that of the P-wave neutron, El=1(λ=2), at a larger value of R, according to the larger size for the Hoyle state. Again at the center of the core, the 1s1/2 orbit for the valence neutron becomes higher than the 0p1/2 orbit, though the energy spacing between them is much smaller than in the λ=1 case. This result means that even for the Hoyle state, which is a dilute object with about 1/3 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 α+ 16O cluster system of 20Ne has been extensively studied by many nuclear models [95-99]. The observed ground-state band and Kπ=01− band in 20Ne were regarded as being an inversion doublet band arising from the hetero-polar di-nucleus configuration [21] of α+ 16O cluster structure. See the schematic diagram in Fig. 20 This kind of inversion doublet band in 20Ne can be considered as the direct manifestation of cluster structures in 20Ne. 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 20Ne, a hybrid-Brink–THSR wave function wave function was proposed as follows [10]:

(65)Φ^Ne(β , S )=A{exp[−∑kx,y,z8( ξ −S )k25Bk2]ϕ(α)ϕ(16O)},

where Bk2=b2+2βk2, ξ = X2−X1,  XG=(4 X1+16X2)/20. X1 and X2 represent the center-of-mass coordinates of the α cluster and the 16O cluster, respectively. All calculations are performed with restriction to axially symmetric deformation, that is, βx=βy≠βz and S ≡(0,0,Sz). 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 M=0.59, and the oscillator parameter b=1.46 fm.

Figure 21 shows the contour map of the energy surface of the intrinsic wave function of 20Ne in the two-parameter space, Sz and βx=βy=βz. We can see the minimum energy −159.66 MeV appears at Sz=0 and βx=βy=βz=1.8 fm. The value Sz=0 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 20Ne. 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 20Ne than the Brink wave function.

After angular-momentum projection on the intrinsic hybrid-Brink–THSR wave function with the spherical case βx=βy=βz, the optimal projected wave function can be obtained by variational calculations. Figures 22 and 23 show the contour maps for the 0+, 2+ positive-parity states and 1−, 3− negative-parity states of 20Ne in the two-parameter space, Sz and βx=βy=βz, respectively. It is surprising to find that the obtained minimum energies with respect to the projected states all appear at Sz=0. For instance, the obtained minimum energy −159.66 MeV for Jπ=0+ state appears at Sz=0 and βx=βy=βz=1.8 fm. For the Jπ=1− state, the minimum energy −155.33 MeV appears at Sz≈0 and βx=βy=βz=2.4 fm in the contour map. For the other states of inversion doublet bands of 20Ne, they have the same character after variational calculations. The inter-cluster distance parameter Sz→0 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 20Ne. 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 Φ^Ne(  β ,  S ) 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 20Ne appears at Sz=3.0 fm. For the Jπ=1− state, the optimum position appears at Sz=3.9 fm. The non-zero values of Sz seem to indicate that the α+ 16O structure of 20Ne 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, βx=βy=βz=1.8 fm, 1.5 fm, 2.4 fm, and 1.9 fm for Jπ=0+, 2+, 1−, and 3− states, respectively, according to their minimum positions in the contour maps, then we find that the minima appear at Sz=0 in Fig. 24. This indicates that the separation distance parameter Sz does not play any physical role in describing the α+ 16O 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 20Ne.

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 20Ne, we can make further variational calculations in the two-parameter space, βx=βy and βz using the projected hybrid-Brink–THSR wave function with the parameter Sz=0 (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 βx=βy and βz using the THSR-type wave functions are listed in Table 4.

On the other hand, the exact solution of the α+ 16O cluster system can be obtained by superposing the Brink wave functions, that is the Brink-GCM wave function,

(66)∑j〈ΦBrinkJπ(Ri)|H^−E|ΦBrinkJπ(Rj)〉f(Rj)=0.

Here, ΦBrinkJπ(Ri) can be obtained directly from the projected Brink wave function ΦBrinkJπ(45 R ,−15 R ) with  R =(0,0,Ri). By solving the Hill–Wheeler equation Eq. (66), we can obtain the following Brink-GCM wave function,

(67)ΦGCMJπ=∑if(Ri)ΦBrinkJπ(Ri).

Thus, we can compare the single THSR-type wave function with the exact Brink-GCM wave function by calculating their squared overlap |〈Φ^minTHSR|Φ^GCMBrink〉|2. 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 Jπ=0+, 2+, 4+, 1−, and 3− states of 20Ne, 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 20Ne. 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 α+ 16O cluster structure in 20Ne.

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 8Be and 20Ne 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 20Ne. 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 20Ne 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 (βx,βy,βz)=(1.78 fm,1.78 fm,7.85 fm). 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 (βx,βy,βz)=(0.01 fm,0.01 fm,5.1 fm) 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 βx=βy, 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 (βx,βy,βz)=(0.01 fm,0.01 fm,5.1 fm).

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 12C but also non-gaslike localized cluster structure states like the inversion doublet band states in 20Ne. 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 4α (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 (02+)THSR, (03+)THSR, and (04+)THSR states in the THSR ansatz correspond to the (02+)OCM, (04+)OCM, and (06+)OCM states in the OCM calculation, which have the structures of 12C(0+) and α in an S-wave, of 12C(0+) and α in a higher nodal S-wave, and of the 4α condensate, respectively. The THSR wave functions of these states, Φ4α,λJ=0 with λ=1,⋯,4 of Eq. (20), are obtained by solving the Hill–Wheeler equation Eq. (21), with a constraint of β0≡β⊥=βz, which are represented as follows: Φ4α,λJ=0=∑β0fλJ=0(β0)Φ4αJ=0(β0). 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 β0.

We first construct orthonormal wave functions from the THSR wave functions as follows:

(68)Φ˜λ(β0)=NλP^λ−1Φ4αJ=0(β0), λ=1,⋯,4,

where Nλ are normalization constants and P^λ−1 are projection operators defined by

(69)P^λ−1=1−∑k=1λ−1|Φ4α,kJ=0〉〈Φ4α,kJ=0|, λ=1,⋯,4.

Here the case λ=1 gives P^0=1. The wave function Φ˜λ(β0) depends on the parameter β0 and is orthogonal to the wave functions Φ4α,λ′J=0 of the (0λ′+)THSR states satisfying λ′≤λ−1. We then calculate the following squared overlap amplitudes between the wave functions Φ˜λ(β0) and Φ4α,λJ=0,

(70)Θλ(β0)=|〈Φ4α,λJ=0|Φ˜λ(β0)〉|2, λ=1,⋯,4.

In Fig. 28 we show Θλ(β0) of Eq. (70) as a function of β0. For any wave functions Φ4α,λJ=0 (λ=1,⋯,4), Θλ(β0) (λ=1,⋯,4) take values close to 100% when β0 is optimally chosen. This result means that all the (0λ+)THSR states are almost represented by single THSR wave functions parametrized by single β0 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 3α condensate state. In particular, in this case for the (04+)THSR state the maximum of Θ4(β0) amounts to 96% at β0=6.5 fm, for which this large value of β0 again supports the picture that this state is considered to be the 4α condensate state, in which the 4α clusters are trapped into a very large (spherical) container characterized by the nuclear size parameter β0.

However, this figure shows that not only the 4α dilute gas-like state but also the lower states (0λ+)THSR with λ=1,2,3, in which 12C+α cluster dynamics prevails, are expressed as single THSR wave functions in the proper orthogonal spaces. This corresponds to the case of 20Ne discussed in the previous section, Section 4, in which the 16O+α 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 β0 values for the four states. We can see that the optimal β0 values increase from the (01+)THSR to the (04+)THSR states. The r.m.s. radii are also proportional to the β0 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 β0 value, i.e., dilatation of the container, from the most compact (01+)THSR state to the most dilute gas-like (04+)THSR state with the 4α condensate structure, through the intermediate density (02+)THSR and more dilute (03+)THSR states with the 12C+α cluster structures. The evolution of the cluster excitations can be described by an evolution of the container, from a smaller value of β0 to a larger one. This also supports the correctness of the container picture.

We should also note that the projection operator P^λ, 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 8Be, where the antisymmetrization operator A removes the Pauli forbidden states, to construct a structural repulsive core between the two α clusters.

The more sophisticated treatment of 12C+α 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 20Ne in Section 4 and even of lower-lying states of 16O in Section 5.1. The ideal situation of the nα 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:

(71)ΦJ(B−GCM)=∑R˜1z,⋯,R˜nzfJ(R˜1z,⋯,R˜nz)P^MKJΦnαB(R˜1z,⋯,R˜nz),

where ΦnαB(R˜1z,⋯,R˜nz) is the Brink wave function in Eq. (16) with R˜ i=(0,0,R˜iz), and the coefficients fJ(R˜1z,⋯,R˜nz) are determined by solving the Hill–Wheeler equation like Eq. (21) taking the variables, R˜1z,⋯,R˜nz, as the generator coordinates instead of β parameter in this equation. In Ref. [104], we calculated the 3α and 4α 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 3α and 4α 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: |〈ΦnαJ(  β )|ΨJ(B−GCM)〉|2. 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 Jπ=0+, 2+, and 4+ states in 12C and 16O nuclei, respectively. The optimal β parameter values in this situation all indicate extremely-prolate deformation, namely, β⊥∼0 and βz≫0.

In Table 6, we show the maximum squared overlap values and the corresponding β values (denoted as 3α and 4α LCS). These large overlap values are to be compared with the maximum square overlap values with 3α and 4α 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 16O+α 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 3α and 4α 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 A, 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 8Be [45, 48], 12C [9, 48, 50], 20Ne [[10, 11, 22], and Λ9Be [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 03+ state in 12C 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 8Be+α-type cluster structure. On the other hand, in Section 4, we showed that the single THSR wave functions excellently describe the 16O+α cluster states in 20Ne, 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 8Be+α structure in 12C, or even 12C+α structure in 16O.

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 β1 and β2 are associated to the two Jacobi coordinates   ξ 1 and   ξ 2, corresponding to α+α and 2α-α motions, respectively. We show it again below,

(72)Φ3αJM(β1,  β2)​=P^JMA[∏i=12exp(−μi∑k=x,y,z2ξix2b2+2βik2)∏i=13ϕ(αi)].

This generalization is a natural extension of the original THSR wave function, since taking   β 1=  β 2 results in it. This also follows the container concept that is discussed in the previous sections, since the case of |  β 1|≤|  β 2| corresponds to the 8Be+α cluster structure, in which the 8Be core and the remaining α cluster are confined in a small-size container, while in the core 8Be the 2α clusters are confined in a large-size container. The shapes and sizes of the two containers are then characterized by the corresponding two parameters,   β 1 and   β 2. We should note that the fact that the THSR wave function includes the 8Be+α asymptotic feature for |  β 1|≪|  β 2| 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 3α 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:

(73)∑  β 1′,  β 2′[〈Φ3αJM(β1,  β2)|112∑i=112(ri−XG)2|Φ3αJM(β 1′,  β 2′)〉−{R(γ)}2〈Φ3αJM(β1,  β2)|Φ3αJM(β 1′,  β 2′)〉]g(γ)(β 1′,  β 2′)​=​0.

We then remove out of the present model space the eigenstates belonging to unphysically large eigenvalues. By taking the following bases,

(74)ΦJM(γ)=∑β1,β2g(γ)( β1,  β2)Φ3αJM( β1,  β2),

with γ satisfying R(γ)≤Rcut, we diagonalize Hamiltonian as follows:

(75)∑γ′〈ΦJM(γ)|H|ΦJM(γ′)〉fλ(γ′)=Eλfλ(γ).

Since the present extended THSR wave function can include 8Be+α asymptotic form by taking the large values of the two width parameters   β 1 and   β 2, the 8Be+α continuum components, as well as the 3α continuum components, can be successfully removed by imposing the cutoff for the mean square radius R(γ)≤Rcut. Although we could not obtain the excited states except for the 02+ and 22+ states by using the original THSR wave function [47], we can now obtain the other observed 03+, 04+, and 42+ 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 03+ 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 12C. 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 12C. In particular, the nature of the series of states (02+, 22+, 42+) and the 03+ and 04+ states have recently been much discussed from the experimental side. As was mentioned in Section 2.1.5, the 22+ state that theoretically has been predicted at a few MeV above the Hoyle state already in the early works of 3α Brink-GCM [13, 42] and 3α RGM [14] was recently confirmed by several experiments [64-68]. A strong candidate for a member of the Hoyle family of states with Jπ=4+ was also reported by Freer et al. [108]. Itoh et al. recently pointed out that the broad 0+ resonance at 10.3 MeV should be decomposed into two states: 03+ and 04+ [69]. This finding is consistent with theoretical predictions where the 03+ state is considered as a breathing excitation of the Hoyle state [16-18, 70] and the 04+ state as the bent arm or linear chain configuration [17, 51, 77].

In Fig. 30(a), the E2 transition strengths between J and J±2 states and monopole transitions between 0+ states are also shown with corresponding arrows. We can note the very strong E2 transitions inside the Hoyle band, B(E2;42+→22+)=591 e2fm4 and B(E2;22+→02+)=295 e2fm4. The transition between the 22+ and 03+ states is also very large, B(E2;22+→03+)=104 e2fm4. In Fig. 30(b), the calculated energy levels are plotted as a function of J(J+1), 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 22+→02+ and 22+→03+ 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 22+,42+ form a rotational band not with the 02+ but with the 03+ 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 02+ and 03+ states. The similar effect is also argued in the study of the 4α condensate and 12C(02+)+α rotational band in 16O [111, 112]. Furthermore, we show in Fig. 31 the S2 factors of the 8Be+α components, which are defined below:

(76)S[I,l]2(Jλ+)=∫​dr[rY[I,l]J(r)]2,

where Y[I,l]J is the RWAs of the 12C(Jλ+) states in the 8Be(I)+α(l) channel. Except for the 04+ state, all the states have dominant components from the channels [0,J]J, which is consistent with 8Be(0+) and α rotation. However, In the 22+ and 42+ states, the components from the channels [0, 2]2 and [0, 4]4, respectively, are dominant. However, in the 22+ and 42+ 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 04+ state in 12C. In the AMD and FMD calculations [51, 77], the dominant intrinsic configuration of the 04+ state is a bent-armed structure of the 3α clusters, resembling the linear-chain structure. However, any firm conclusions cannot still be drawn, since in the AMD, FMD, or even the 3α Brink-GCM calculation, the four 0+ states are not consistently reproduced. As we mentioned above, in these calculations, the 03+ state is missing and only the 04+ 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 04+ 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 8Be(2+) 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 04+ state.

In Fig. 32, we show the contour map of the squared overlap between the 04+ state and the single THSR configuration in two-parameter space,   β 1=  β 2=(βx=βy,βz). 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 3α linear-chain structure. The parameter values giving the maximal squared overlap is calculated to be   β 1=  β 2=(0.6, 6.7 fm), which is close to   β 1=  β 2=(0.1, 5.1 fm) 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 3α clusters, which is fairly different from the ordinary picture of the linear-chain state with rigid-body 3α-cluster configuration arranged in a line in a spatially localized way. Thus we can say that the present 04+ state has the one-dimensional α condensate structure by around 50%, where the 3α 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 12C also have many well-developed cluster states [13, 14, 114]. Quite recently, the 1− 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 12C 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 12C was assumed to have D3h 3α cluster structure and many cluster states 0+,2+,3−,4±,5− were predicted. Particularly, the 5− state at 22.4(2) MeV in 12C 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 12C were not favored, especially the Hoyle state and other developed positive-parity cluster states. The positive-parity states of 12C 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 12C, we focus on the 3- and 4- states in 12C 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 20Ne in Section 4. We start from the general two-β THSR wave function,

(77)Φ( β ,  S)=∫ R1 d R2 exp[−( R1− S1)2β12−( R2− S2)2β22]ΦB( R1,  R2)∝ϕGA{exp[−(ξ1− S1)2B12−(ξ2− S2)2B22ϕ(α1)ϕ(α2)ϕ(α3)]},

(78)ΦB( R1,  R2)∝ϕGA{exp(−( ξ1− R1)2b2−(ξ2−R2)234b2)ϕ(α1)ϕ(α2)ϕ(α3)},

where B12=b2+  β 12, B22=34b2+  β 22. β1≡(β1x=β1y,β1z),   β2 =(β2x=β2y,β2z). ξ1=X2−X1 and  ξ2=X3−(X1+X2)/2. The  X1, X2, and  X3 are the center-of-mass coordinates of the three clusters in 12C. In the above equation, ΦB(R1,  R2) is the 3α Brink wave function with the corresponding generator coordinates  R1 and  R2. The R1 represents the specified distance parameter of two clusters and R2 these two clusters and the third α cluster. The  S1 and S2 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  β1=2β  and  β2=3/2 β  in the above equations. And the introduced shift parameter  S1≡(S1x,S1y=S1z=0) and  S2≡(S2x,S2y,S2z=0) 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),   S 2=(0,0,3/2).

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:

(79)     Φ4αJ=0( β1,  β2)=P^J=0×A{exp[−2∑kx,y,z(μ1ξ1k2+μ2ξ2k2b2+2β1k2−μ3ξ3k2b2+2β2k2)]∏i=14ϕ(αi)},

where ξi is again the Jacobi coordinates between the α particles, and μi=i/(i+1), 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 βi⊥≡βix=βiy, so as to deal with the four parameters, β1⊥,β1z,β2⊥,β2z, 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 Jπ=0+ 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 0V+ 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 06+ state has the 4α condensate character and the 02+-05+ states all have α+12C cluster structures, i.e., α(S)+12C(01+), α(D)+12C(21+), α(S)+12C(01+), and α(P)+12C(1−) cluster structures, respectively. The difference between the 02+ and 04+ 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 α(P)+12C(1−) cluster structure, like the 05+ 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 0I+ to the 0V+ 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:

(80)Oλ( β1,  β2)=|〈Φ˜λJ=0( β1,  β2)|Ψλ〉|2,

with Ψλ the λth eigenfunction and Φ˜λJ=0(β1,  β2) normalized single eTHSR wave function in a space orthogonal to the lower eigenstates, i.e.,

(81)Φ˜λJ=0(β1,  β2)=NλP^λΦ4αJ=0( β1,  β2),

where P^λ=1−∑i=1λ−1|Ψi〉〈Ψi| with λ=I,⋯,IV, and P^V=1−∑i=16|Ψi〉〈Ψi|, and Nλ are the corresponding normalization constants.

This quantity indicates how these five states Ψλ(λ=I,⋯,V) 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 ΨI–ΨV with single configurations in the β2⊥ and β2z parameter space in Eq. (80) are shown. Here the β1 parameter, i.e., (β1⊥,β1z) is fixed at the position denoted by ⊗ in these figures, so that the maximum value of the squared overlap in the four parameter space (β1⊥,β1z,β2⊥,β2z) 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 01+, 21+ and 02+ 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 β1⊥≫β1z, while the remaining α cluster is put into a prolate-deformed and very compact container with β2⊥≪β2z. 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 0I+ state by 98%.

In the 0II+ state, shown in Fig. 41(b), the 3α clusters are in a spherical container with β1⊥∼β1z. The fourth α cluster is put into a larger size container with spherical shape, i.e., β2⊥∼β2z>β1⊥∼β1z. In particular, the parameter set (β1⊥,β1z)=(1.8,1.8 fm) is almost the same as that for 12C in Table 9, i.e., (β⊥,βz)=(1,9,1.8 fm). This means that the first 3α clusters are confined in a compact container to form the ground state of 12C, since the 12C(01+) 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 (β2⊥,β2z)=(3.5,3.6 fm), 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 12C(01+) state.

The 0III+ state, shown in Fig. 41(c), is similar to the 0II+ state but both containers are not spherical but deformed. The β1 parameter takes almost the same value as that of the isolated 12C(2+) state, as shown in Table 9, which means that the first 3α clusters form the 12C(2+) state, since the state is described by the single parameter value of β . The configuration of the remaining α cluster (β2⊥,β2z)=(5.1,0.1 fm), giving the largest value 76%, means that the α cluster moves in a deformed and larger container. This is present understanding of the 03+ state, which is conventionally considered to have the α(D)+12C(2+) structure.

In the 0IV+ state, shown in Fig. 41(d), one can see that the 3α clusters are put in slightly larger container than that for the 12C(01+) 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 0IV+ state dominantly by 84%. This means that the second container characterized by  β2 is further evolved from that in the 0II+ state. We can say that this state corresponds to the 04+ state in the former 4α OCM calculation, which predicts the α+12C(01+) higher nodal structure for the state.

The 0V+ state, shown in Fig. 41(e), is the most striking. All the α clusters occupy an identical orbit, with (β1⊥,β1z,β2⊥,β2z)=(5.3,1.9,5.3,1.8 fm). 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 (β1⊥,β1z)=(5.6,1.4 fm) in Table 9. This means that the 0V+ 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 06+ 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 Φ˜kJ=0 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 nα 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 02+ 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 02+ 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]

(82)|Φ9Be〉=(Cα†)2cn,σ†|vac〉,

(83)cn,σ†​=​∫d Rnexp(imϕRn)exp(−Rn,x2+Rn,y2βn,xy2−Rn,z2βn,z2)   ×∫drn(πb2)−3/4e−(rn−Rn)22b2an,σ,†(rn),

where Cα† and cn,σ are creation operators of α-particle and neutron, respectively. rn represents the coordinate of the extra neutron and its corresponding generate coordinate is Rn. aσ,τ†(rn) is the creation operator of the extra neutron with spin σ. Here, spin is assumed to be up. The ϕRn is the azimuthal angle of Rn, 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

(84)|Φ10Be〉=(Cα†)2cn,↑†cn,↓†|vac〉.

Here, Cα† is the creation operator of the α cluster. cn,↑† and cn,↓† are the creation operators of the two valence neutrons with spin up and down, respectively. See details in Ref. [125].

As we know, 02+ state in 10Be was reported to be a developed molecular cluster state from AMD calculations [124]. The calculated energy spectrum of the 02+ 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 02+ 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-nα 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.

· 2α+Λ cluster structure in Λ9Be. 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 nα THSR wave function.

The H-THSR wave function describing the 2α+Λ hypernucleus with good angular momentum is then introduced as follows:

(85)ΦJH−THSR(β)=P^MKJΦ2αTHSR(  β )∑κfΛ( β ,κ)φΛ(κ),

where P^MKJ is the angular-momentum projection operator and the Λ particle simply couples to the 2α core nucleus in an S wave. Its radial part is expanded in terms of Gaussian basis functions, φΛ(κ)=(π/2κ)−3/4exp(−κr2α−Λ2) with respect to the width parameter κ, where  r2α−Λ=XΛ−(X1+X2)/2 with XΛ a position vector of the Λ particle. In the practical calculations, we assume an axially-symmetric deformation for Φ2αTHSR(β ), βx=βy≡β⊥≠βz.

The coefficients of the expansion fΛ(β⊥,βz,κ) are then determined by solving the following Hill–Wheeler-type equation of motion,

(86)∑κ′〈P^M0JΦ2αTHSR(β⊥,βz)φΛ(κ)|H^−E(β⊥,βz)|P^M0JΦ2αTHSR(β⊥,βz)φΛ(κ′)〉fΛ(β⊥,βz,κ′)=0.

This H-THSR wave function is characterized only by the parameters β⊥ and βz, 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  Λ13C.

In order to compare the single H-THSR wave function with the Brink-GCM wave function, we calculate the following squared overlap:

(87)OJ(β⊥,βz)=|〈ΦJH − THSR(β⊥,βz)|ΦJB − GCM〉|2〈ΦJH − THSR(β⊥,βz)|ΦJH − THSR(β⊥,βz)〉〈ΦJB − GCM|ΦJB − GCM〉.

The above wave function ΦJB−GCM is the Brink-GCM wave function projected onto the model space with the angular-momentum channel (L,λ)=(J,0), in which L, λ are the angular momentum between the 2α clusters and the one between the Λ particle and 8Be, respectively.

In Fig. 45, we show the contour maps of this quantity in the two-parameter space β⊥ and βz, calculated for the Jπ=0+ (a), 2+ (b), and 4+ (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 Jπ=0+,2+ states are 99.5% and 99.4%, respectively. For the Jπ=4+ state, the maximum value is slightly down to 97.7%. This reduction may originate from a mixture of spurious scattering-state components. These practically 100% squared overlap values of course mean that, at least on the subspace with L=J and λ=0, the Brink-GCM wave function obtained by solving the 2α+Λ 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 (J,λ)=(L,0) 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 Λ9Be the 2α clusters are trapped into a container, which is specified by the optimal value of the size parameter β, under the influence of the antisymmetrizer A acting on the nucleons, such as realized in the form of the wave function Eq. (85).

We here discuss the intrinsic structure of Λ9Be, together with that of 8Be. 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,

(88)ρN(  r )=〈Φ2α−ΛH − THSR(β⊥,βz)|∑i=18δ( r i− XG−  r )|Φ2α−ΛH − THSR(β⊥,βz)〉〈Φ2α−ΛH − THSR(β⊥,βz)|Φ2α−ΛH − THSR(β⊥,βz)〉.

This density is normalized as usual to the total number of nucleons, ∫​drρN( r)=8.

In Fig. 46, we show the intrinsic density profiles of Λ9Be defined by Eq. (88) and of 8Be. The single optimal β values giving the minimum energies are adopted, i.e., (β⊥,βz)=(1.5,2.8 fm) for Λ9Be and (1.8,7.8 fm) for 8Be. While both ones clearly show the 2α-cluster structure with the prolate-deformed shape, 8Be has a gas-like tail of the 2α clusters and Λ9Be not. In Λ9Be the Λ particle gives rise to strong shrinkage and the gas-like tail in 8Be disappears. The r.m.s. radius of the 8Be core is accordingly changed from Rrms=2.87 fm for 8Be to Rrms(c)=2.31 fm for Λ9Be. 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 A 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 ΛN interaction artificially, with overall factor δ multiplied with V(ΛN). The β values to give the energy minima after projection onto J=0 space, which are listed in Table 11, are adopted. As the ΛN interaction is strengthened, from δ=0.4 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 ΛN potential, rms radii of the 8Be core, are also shown in Table 11. δ=0 and 1 correspond to the cases of 8Be and Λ9Be, respectively. Starting from the case of 8Be with very large value of βz and small value of β⊥, i.e., β⊥=1.8 fm and βz=7.8 fm, only the βz value drastically gets smaller as the increase of ΛN interaction, and eventually for Λ9Be, βz becomes much smaller while β⊥ is almost unchanged, i.e., β⊥=1.5 fm and βz=2.8 fm.

We should note that at each step of the change of ΛN 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 8Be a single optimal Brink wave function can only approximate the Brink-GCM wave function by 0.722 [48], while for Λ9Be the former can approximate the latter by 0.940 [106].

· Multi-cluster dynamics in Λ13C. We discuss  Λ13C in this subsection by using the H-THSR ansatz.  Λ13C is theoretically investigated and the second 1/2+ state is found as a bound state, which is considered to correspond to the 12C(02+)⊗Λ state [133, 134], although experimentally the second 1/2+ state in  Λ13C 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 Λ9Be, as follows:

(89)Φ3αH−THSR( β1,  β2,  βΛ)=A[exp{−∑i=12μi∑k=x,y,z2ξik2b2+2βik2}ϕ3(α)]φΛ(βΛ),

with

(90)φΛ(βΛ)=exp(−μΛ∑k=x,y,z2ξΛk2b2+2βΛk2),

where the φΛ(  β Λ) is the wave function of the Λ particle with its width parameter,   β Λ, which can be out of the antisymmetrization, μΛ=3mΛ/(12mN+mΛ) with mN and mΛ being nucleon and Λ-particle masses, respectively, and ξΛ=XΛ−( X1+X2+X3)/3. 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., βix=βiy≠βiz and βΛx=βΛy≠βΛz. 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 Jπ=0+ spectra of 12C and  Λ13C and the observed Jπ=0+ 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 03+ state is concluded to be the excited state from the Hoyle state by monopole vibration and the 04+ state dominantly has a 3α linear-chain configuration. While the 03+ state is considered to be a family of the Hoyle state, the 04+ state has a quite different structure from the 02+ and 03+ states. In  Λ13C, first we should mention that we simply denote the spectrum of  Λ13C as Jπ=0+ states, neglecting the intrinsic spin 1/2 of the Λ particle, for simplicity. From this figure, one can see that  Λ13C gives a much larger number of decay channels than 12C, i.e., Λ+12C(01+), α+Λ9Be(01+), Λ+12C(21+), α+Λ9Be(21+), 8Be(0+)+Λ5He, α+Λ9Be(02+), and 8Be(2+)+Λ5He. Up to the 3α+Λ threshold, we have the four Jπ=0+ states, corresponding to the four Jπ=0+ states in 12C. We can see that the 02+ in  Λ13C is obtained as a bound state which is located below the lowest threshold, Λ+12C(01+), by 0.8 MeV. The 03+ and 04+ states are connected to the 04+ and 03+ states of 12C, respectively, as you can also see from the figure. This assignment is because the 03+ and 04+ states in  Λ13C largely include 12C(04+)+Λ and 12C(03+)+Λ components, respectively [131].

In fact, the 03+ state in  Λ13C is found to keep essentially the same features as the 04+ state in 12C, which are shown in Section 5.3. For example, we obtain very similar result as in Fig. 32, also for the 04+ state in  Λ13C, 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 04+ state in 12C, but the β parameter values to give the maximum is shifted from (βx=βy, βz)=(0.6,6.7 fm) in the 04+ state of 12C, to (βx=βy,βz)=(0.4,4.6 fm) in the 03+ state of  Λ13C. The linear-chain configuration in  Λ13C then shrinks along z-direction, compared with the case of 12C, due to the additional Λ particle.

In Table 12, we show Λ binding energy BΛ, the root mean square (r.m.s.) radius of the core Rrms(c), the r.m.s. distance between the core and Λ particle Rrms(c−Λ), and the r.m.s. radius of  Λ13C, Rrms, together with the r.m.s. radii of isolated 12C calculated in Ref. [16]. Here the Λ binding energy is defined as the binding energy measured from the corresponding 12C(0+)+Λ threshold. From the calculated r.m.s. radii for these states, we can expect that the 04+ state is qualified to be a gas-like cluster state, since only this state has a larger r.m.s. radius 4.3 fm than the one of the Hoyle state 3.7 fm. From the RWA analysis, we find that the 04+ state has the largest component of α+Λ9Be(02+), in which Λ9Be(02+) has a dilute α+Λ5He structure [131]. Therefore this state is a loosely coupled 2α+Λ5He state, i.e., a Hoyle-analog state. In this table, we can also see that the 04+ state suffers from much less shrinkage effect from the additional Λ particle than the other 02+ and 03+ states, except for the ground state that has saturation density and gives almost no shrinkage effect. The BΛ value for this state is also smaller than the one of the other states. This is because for the 04+ state the Λ particle couples with only one α cluster, to form Λ5He. The little overlap with nucleons reduces the shrinkage as well as the gain of the Λ binding energy. In fact, this state appears around the 2α+Λ5He threshold energy, for the threshold rule. The slightly higher energy position, 1.6 MeV above the 2α+Λ5He threshold, allows for the diluteness of this state.

On the other hand, the dominant component included in the 02+ state is α+Λ9Be(01+) configuration. The Λ particle then overlaps with the 2α clusters and therefore gives much stronger shrinkage effect and larger Λ binding energy. The binding energy of the Λ9Be is so large that it is energetically favored that the injected Λ particle shrinks the core and forms α+Λ9Be(01+) cluster configuration. In fact, this state is located at around the α+Λ9Be(01+) threshold as well as around the Λ+12C(01+) threshold (see Fig. 47), which follows the so-called threshold rule.

We further show in Fig. 48 the RWAs in the 8Be(0+)+Λ5He channel of the 02+ (solid curve in red) and the 04+ (solid curve in blue) states in  Λ13C. In comparison, the corresponding RWAs of the Hoyle state (dotted curve in red) and the 03+ state (dotted curve in blue) in 12C for the 8Be(01+)+α 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 8Be core, while in the 03+ state in 12C (dotted curve in blue) there exist four nodes. On the other hand, in  Λ13C, the most prominent feature is the disappearance of nodes for the 04+ state, though the innermost node only remains. This means that the 8Be core is dissolved in the 04+ state, completely unlike the case of the 03+ state in 12C. We should also be aware that for the 04+ state a long tail is still very much developed, which is similar to the Hoyle state. This is also completely different from the 02+ state, in which the whole amplitude is pushed inside with no long tail any more. All these results give a strong support that the 04+ state has a gas-like structure of 2α+Λ5He clusters, as a Hoyle analog state. It is important to mention that this Hoyle analog state uniquely appears in  Λ13C and we cannot expect in 13C 2α+5He gas, since 5He is not bound [136], whereas Λ5He 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

(91)Φ0( βα,  βp)​=∫d R1d R2exp​[−​∑i=12​(Ri,x2+Ri,y2βα,xy2​+​Ri,z2βα,z2)​]∫d Rpexp[−Rp,x2−Rp,y2βp,xy2−Rp,z2βp,z2]×eimϕRpΦ2α+pB(R1,R2,Rp),

where βα is the size parameter for the nonlocalized motion of 2α clusters and βp is the size parameter for the valence proton around the 2α clusters. The eimϕRp is the phase factor for the treatment of the negative-parity states. Φ2α+pB(R1,R2,Rp) is the Brink wave function, in which the ( R1,  R2), and Rp are corresponding generator coordinates for the α clusters and the valence proton, respectively.

In the 2α+n 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],

(92)Φj( βα,  βp)​=​∫​d R1dR2exp​[​−​∑i=12​(Ri,x2+Ri,y2βα,xy2​+​Ri,z2βα,z2)]  ×​∫dRpexp[−Rp,x2−Rp,y2βp,xy2−Rp,z2βp,z2]  ×​eimϕRp+RjΦ2α+pB(R1,R2,Rp​+​Rj), j​=​1,2.

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 Φj(  β α,  β p) 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 Φj(  β α,  β p) and large spreading configuration Φ0(  β α,  β p)

(93)Ψ=c(Φ1(βα,  βp′)+Φ2( βα,  βp′))+dΦ0( βα,  βp).

Here Φ1 and Φ2 can include cluster-correlated configurations of 9B describing the correlated motion of valence proton around each α-cluster. The wave function Φ0 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 Φ0, 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],

(94)G^n(β0)=∫d R1⋯d Rnexp[−∑i=1nAiRi2β02−2bi2]D^(R ).

Here, the D^ is a shift operator. With the auxiliary generator coordinate  R (R1… Rn), this G^n(β0) 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:

(95)Ψc(β0,bc1,bc2)=∫d R1dR2exp(−A1R12β02−2bc12−A2R22β02−2bc22)ΦcB(R1,R2)

(96)∝exp(−Aβ02Xcm2)A[exp(−A1A2Aβ02Xrel2)ϕc1int(bc1)ϕc2int(bc2)].

(97)ΦcB(R1,R2)∝A{exp[−A12bc12(X1−R1)2−A22bc22(X2−R2)2]ϕc1int(bc1)ϕc2int(bc2)}.

Here the clusters C1 and C2 have the atomic mass number A1 and A2, respectively. X1 and X2 are the corresponding center-of-mass coordinate, respectively. Xcm=(A1X1+A2X2)/A and Xrel=X2−X1. bc1 and bc2 are the width parameters or the size parameters of cluster C1 and C2. ΦcB 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 20Ne, 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 20Ne wave function can be written directly,

(98)ΨNe(β0,bα,b16O)∝exp(−20β02Xcm2)A[exp(−165β02Xrel2)ϕαint(bα)ϕ16Oint(b16O)].

Here Xcm=(4X1+16X2)/20 and Xrel=X2−X1. X1 and X2 are the center-of-mass coordinate of α cluster and 16O cluster, respectively. bα and b16O are the width parameters or the size parameters of α cluster and 16O cluster, respectively. The ΦNeB(R1,R2) is a generalized Brink wave function [41, 144] for 20Ne 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., bα=b16O in Eq. (97), to avoid the center-of-mass problem.

By making variational calculations for widths bα, b16O, and β0 variables, the energy curve of 20Ne and the variation of width variables can be obtained. Figure 50 shows the energy curve of 20Ne along with the obtained optimum widths bα and b16O using the intrinsic wave function in Eq. (98). Clearly, it can be seen that with the increase of the size parameter β0 indicating the cluster system expands, the gap between bα and b16O is becoming larger and larger. Another saying, imagine a situation where the stable double-closed-shell nuclei α and 16O are approaching each other to form a stable nuclear state (the ground state of 20Ne), in the process Fig. 50 tells us the 16O is slightly shrinking while the α nucleus is becoming more loosely bound. It has been known that the ground state of 20Ne 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 β0=4.2 fm and slightly different values of widths bα=1.45 fm and b16O=1.46 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 E(bα=b16O=1.47 fm, R=3.0 fm) =−158.43 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

(99)〈Ψnα(β )|H^|Ψnα(  β ′)〉​=∫d R1…d Rn′exp(−∑i=1n(Ri2β2+R′i2β′2))

(100)×〈ΦnαB(R1,⋯,  Rn)|H^|ΦnαB( R1′,⋯,  Rn′)〉.

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 20Ne (α+ 16O) 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 6×3×2 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

(101)Ψnα( β )=∫d R1…d Rnexp[−∑i=1nRi2β2]A[Φα1( R1)⋯Φαn(Rn)]

(102)=A{∏i=1n[∫d Riexp[−Ri2β2]Φαi( Ri)]}.

Here Φαi is the wave function for the ith α cluster centered at Ri, which has a Gaussian form Exp[−ν(r − Ri)2] with the width ν. To simplify the calculation of this THSR wave function, the integral over the generator coordinate parameters Ri in the original THSR wave function can be replaced by the summation of many Slater determinants as follows:

(103)Ψnα(β)≈∑k=1mA[Φα1(R1)⋯Φαn( Rn)].

The generator coordinate parameter {Ri} can be randomly produced, however, the random numbers are generated by the weight function W(Ri) with the same Gaussian shape:

(104)W(Ri)∝exp(−Ri2β2).

With the increasing ensemble number, the distribution of {Ri} 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 12C [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 16O core in 24Mg. Figure 51 shows the energy convergence of 24Mg measured from the 16O +2α threshold, where two α-clusters occupy the α-condensed state around the 16O core. The solid and dotted lines (β=2, 3 fm) converge to the energies below the threshold, however, the dashed- and dash-dotted-lines (β=4, 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 24Mg around the 16O +α +α 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 α+ 16O 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:

(105) Φ(R1,R2,R3)​=​∏p​|ϕp(1)(R1)ϕp(1)(R2)ϕp(1)(R3)ϕp(2)(R1)ϕp(2)(R2)ϕp(2)(R3)ϕp(3)(R1)ϕp(3)(R2)ϕp(3)(R3)|p,

where suffix p means the spin and the isospin part. and Ri represents the coordinate of the ith α-cluster in the parameter space R. The spatial part of a single nucleon wave function is written by

(106)ϕp(i)(Rj)=1(πb2)3/4exp{−12b2{rp(i)−Rj)2},

where i,j=1,2,3. The variable rp(i) stands for the real coordinate for ith nucleon with pth 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 3α clusters, is given by

(107)C^(β)=1(πβ2)3×3/2∫∏i=13dRiexp(−Ri2β2),

where the integral is carried out over the infinite parameter space. The THSR wave function is represented by

(108)Ψ(β)=C^(β)Φ(R1,R2,R3).

The size parameter of the container is denoted as β.

Let us obtain the norm kernel by the analytical expression:

(109)​〈Ψ(β)|Ψ(β′)〉=C^(β)C^(β′)×〈Φ(R1,R2,R3)|Φ(R1′,R2′,R3′)〉.

The overlap of the wave function (105) is given after carrying the integral with respect to real variable rp(i):

(110)〈Φ(R1,R2,R3)|Φ(R1′,R2′,R3′)〉=|G11G12G13G21G22G23G31G32G33|4,

where the overlap between single nucleon wave functions is written by

(111)Gij=〈ϕp(k)(Ri)|ϕp(k)(Rj′)〉=exp{−14b2(Ri−Rj′)2},

where the suffix p is for all the spin-isospin states, where k=1,2,3. 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 (G11G22G33)4 has the exponential form:

(112)exp{−(R˜)t((A˜)(B˜)(B˜)t(A˜′))(R˜)},

where matrices are given by

(113)(R˜)t=(R1,R2,R3,R1′,R2′,R3′),(A˜)=(1β2+1b20001β2+1b20001β2+1b2),(A˜′)=(1β′2+1b20001β′2+1b20001β′2+1b2),(B˜)=−(1b20001b20001b2),

where only the matrix (B˜) is varied depending on the exchange property, but the matrices (A˜),(A˜′) are unchanged during all the permutations.

The multiple integral with respect to the parameters, R1,R2,R3,R1′,R2′,R3′ 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:

(114)N1(β,β′)=1(β2β′2)3×3/2|(A˜)(B˜)(B˜)t(A˜′)|−3/2={(1​+​β2b2)​3(1​+​β′2b2)​3|(A˜˜)(B˜˜′)(B˜˜)t(A˜˜)|}−3/2,

where

(114)(A˜˜)=(100010001),(B˜˜)=−(x000x000x) with x=β2b21+β2b2,(B˜˜′)=−(x′000x′000x′) with x′=β′2b21+β′2b2.

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):

(115)|(A˜˜)(B˜˜′)(B˜˜)t(A˜˜)|=|100−x′000100−x′000100−x′−x001000−x001000−x001|.

At first, we have to divide all the columns into two parts; that is, for instance, the first column is written by

(116)(100−x00)=(100000)−x(000100).

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 x′. Almost all of the separated determinants are exactly 0 except for the power of xx′ 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:

(117)Nall(x,x′)=f(x,x′)∑k=19Wkfk−3/2(xx′),

where

(117)f(x,x′)={(2−x1−x)(2−x′1−x′)}−3×3/2,

(117)fk(xx′)=∑l=03ckl(−xx′16)l.

The weight Wk, which comes from the same exchange character of

(118)(G1i1G2j1G2k1)↑↑(G1i1G2j1G2k1)↑↓,(G1i1G2j1G2k1)↓↑(G1i1G2j1G2k1)↓↓,

and the coefficients ckl are listed in Table 13. All the 9 terms with respect to the polynomial have a common factor (1−xx′) 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 x,x′ relating to the α condensation strength are defined in the range of 0<x,x′<1. In the case of x,x′→0, the system goes to shell-model state, on the other hand, for x,x′→1 the system shows a free α gas aggregate.

· The integrand with distance parameter. Let us assume 2α–α system. The distance parameter d between 2α and α clusters is introduced; namely, In Eq. (110) the matrix elements of overlapping are explicitly written by

(119)Gij=exp{−14b2(Ri−Rj′−13d+13d′)2}            for i,j≤2,Gi3=exp{−14b2(Ri−R3′−13d−23d′)2}            for i≤2,G3j=exp{−14b2(R3−Rj′+23d+13d′)2}            for j≤2,G33=exp{−14b2(R3−R3′+23d−23d′)2}.

Concerning Eq. (113) related to (G11G22G33)4, we write the power part in exponential as

(120){−(R˜d)t((A˜)(B˜)(D˜11)(D˜12)(B˜)t(A˜′)(D˜21)(D˜22)(D˜11)t(D˜21)td11d12(D˜12)t(D˜22)td21d22)(R˜d)},

where matrices are given by

(120)(R˜d)t=(R1,R2,R3,R1′,R2′,R3′,d,d′),

(120)(D˜11)t=1b2(−13,−13,23),

(120)(D˜12)t=1b2(13,13,−23),

(120)(D˜21)t=1b2(13,13,−23),

(120)(D˜22)t=1b2(−13,−13,23),

and the other scalars are

(121)d11=−d12=−d21=d22=69b2.

The D matrices and the scalars are independent of the permutation property. The matrices (A˜), (A˜′) and (B˜) are common to those in Eq. (113).

The multiple integral with respect to the parameters, R1,R2,R3,R1′,R2′,R3′ 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:

(122)N1(β,β′)exp{−(d,d′)(Δ˜11Δ˜12Δ˜21Δ˜22)(dd′)},

where

(122)Δ˜ij=|(A˜)(B˜)(D˜1j)(B˜)t(A˜′)(D˜2j)(D˜1i)t(D˜2i)tdij||(A˜)(B˜)(B˜)t(A˜′)|.

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 Δ11 as

(123)Δ˜11=19b2|(A˜˜)(B˜˜)(D˜˜11)(B˜˜)t(A˜˜′)(D˜˜21)(D′~~11)t(D′~~21)t6||(A˜˜)(B˜˜)(B˜˜)t(A˜˜′)|,

where

(123)(D˜11)t=(−1,−1,2)

(123)(D˜21)t=(1,1,−2)

(123)(D˜11′)t=x(−1,−1,2)

(123)(D˜21′)t=x′(1,1,−2).

The explicit form of the numerator in Eq. (123) is

(124)|100−x′00−10100−x′0−100100−x′2−x0010010−x0010100−x001−2−x−x2xx′x′−2x′6|.

The same procedure as that of the original THSR derivation is carried out. As a consequence, the analytical form is written by

(125)Nall(x,x′:d,d′)=f(x,x′)∑k=19fk−3/2(xx′)∑nWkn(1)gkn(1),

where

(125)gkn(1)=exp{−1fk(xx′)(d,d′)(Δ11(kn)Δ12(kn)Δ21(kn)Δ22(kn))(dd′)}.

The coefficients Wkn(1) are listed in Table 14, where Wk=∑nWkn(1) is satisfied for each type of k in the original THSR ansatz. The property of the function Δij(kn) with respect to the polynomial of x,x′ is as follows:

(125)type  s:∑lcl(xx′)l+(x+x′)∑lcl′(xx′)l,

(125)type  s0:∑lcl(xx′)l,

(125)type  ns:∑lcl(xx′)l+x∑lcl′(xx′)l+x′∑lcl′′(xx′)l,

(125)type  nsx:∑lcl(xx′)l+x∑lcl′(xx′)l,

(125)type  nsx′:∑lcl(xx′)l+x′∑lcl′(xx′)l.

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 (d⋅d′).

· 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

(126)C^(β1,β2)=1(πβ12)2×3/21(πβ22)3/2×∫∏i=12dRiexp{−(Ri−13d)2β12}×∫dR3exp{−(R3+23d)2β22}.

In the same way as the previous discussions, the Eq. (123) is

(127){−(R˜d)t((A˜)(B˜)(D˜11)(D˜12)(B˜)t(A˜′)(D˜21)(D˜22)(D˜11)t(D˜21)td11d12(D˜12)t(D˜22)td21d22)(R˜d)},

where matrices are

(127)(R˜d)t=(R1,R2,R3,R1′,R2′,R3′,d,d′),

(127)(D˜11)t=(−13β12,−13β12,23β22),

(127)(D˜12)t=(0,0,0),

(127)(D˜21)t=(0,0,0),

(127)(D˜22)t=(−13β′12,−13β′12,23β′22).

The scalars dij are as follows:

(128)d11=29β12+49β22, d12=d21=0, d22=29β′12+49β′22.

First of all, we should scrutinize two-strength THSR ansatz as the case of d=0. Namely, for the first term of (G11G22G33)4 the deteminant is given by

(129)|100−x1′000100−x1′000100−x2′−x1001000−x1001000−x2001|,

where

(129)xi=βi2b21+βi2b2, xi′=β′i2b21+β′i2b2  with  i=1,2.

The function of f(x,x′) should be also exchanged by f(x1,x2,x1′,x2′) as

(129)f(x1,x2,x1′,x2′)={(2−x11−x1)(2−x1′1−x1′)}−3×3/2×{(2−x21−x2)(2−x2′1−x2′)}−3/2.

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 d≠0. However, the function of fk(xx′) is written by the combination of 4 kinds of x parameter. As for the case of d≠0, the function for each term is analytically written by the same form as that of Eq. (122). The numerator Δ11 for the term of (G11G22G33)4 is given by dividing the final column as

(130)(−1β12−1β122β220002β12+4β22)=1β12(−1−100002)+1β22(0020004).

Thus, we can easily write the numerator such as

(131)19β12|100−x1′00−10100−x1′0−100100−x2′0−x10010000−x10010000−x20010−1+x1−1+x12+2x20002|+19β22|100−x1′0000100x1′0000100−x2′2−x10010000−x10010000−x20010−1+x1−1+x12+2x20004|,

where we use the relation of

(131)11+βi2b2=1−xi,  with  i=1,2

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 R3,R3′ to the origin (=0), we carry out the integral with respect to R1,R2,R1′,R2′. For the first term of (G11G22G33)4, the determinant is written by

(132)|10−x′0010−x′−x0100−x01|.

The resultant norm kernel is as follows:

(133)Nall(p)(x,x′)=f(p)(x,x′)∑k=120Wk(p){fk(p)(xx′)}−3/2,

where

(133)f(p)(x,x′)={(2−x1−x)(2−x′1−x′)}−2×3/2,

(133)fk(p)(xx′)=∑l=02ckl(p)(−xx′16)l.

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 x,x′ 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 (10!)4≈1.73×1026 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, r1,⋯,rn 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 4n nuclei (“N=Z=even” 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 det|(0s)4(0p)12| of 16O. This wave function with (λ,μ)=(0,0) symmetry of SU(3) can be rewritten in the form of cluster-model wave function of 12C + α configuration,

(134)116!det|(0s)4(0p)12|​=Ng1C164A{[R4( r ,3νN)ϕ(0,4)(12C)](0,0)ϕ(α)},

(135)g(r (16O),16νN)[R4( r ,3νN)ϕ(0,4)(12C)](0,0)​=∑L=0,2,4CL[R4,L(r ,3νN)ϕ(0,4)L(12C)]J=0.

(136)g( X ,γ)=(2γπ)3/4exp(−γ X 2),  r (16O)=116∑i=116  ri,C164=16!12!4!, CL=〈(4,0)L,(0,4)L||(0,0)0〉.

Here r is the relative coordinate between 12C and α clusters and R4,L( r ,3νN) stands for the harmonic oscillator function R4,LM( r ,3νN)=R4,L(r,3νN)YLM(r^) with the oscillator quanta N=4, namely the nodal number n=(4−L)/2, and with the size parameter 3νN. ϕ(0,4)L(12C) is the internal wave function of 12C with SU(3) symmetry (0,4) and is given by the 0p-shell shell-model wave function as

(137)ϕ(0,4)L(12C)g( r (12C),12νN)​=​|(0s)4(0p)8;(0,4)L〉,

(138)r (12C)=112∑i=112  r i.

The relative wave function R4,LM( r,3νN) has SU(3) symmetry (4,0). [R4(r ,3νN)ϕ(0,4)(12C)](0,0) stands for the SU(3) vector-coupling of R4,L(r ,3νN) with (4,0) symmetry and ϕ(0,4)L(12C) with (0,4) symmetry to the SU(3) symmetry (0,0) by the use of the SU(3) reduced Clebsch-Gordan coefficient 〈(4,0)L,(0,4)L||(0,0)0〉; (4,0)​×​(0,4)​→(0,0).

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

(139)A{[R4,L( r ,3νN)ϕ(0,4)L(12C)]J=0ϕ(α)}=∑(λ,μ)〈(4,0)L,(0,4)L||(λ,μ)0〉A{[R4( r ,3νN)ϕ(0,4)(12C)](λ,μ)ϕ(α)}.

The number of total oscillator quanta, Q, of the double-closed-shell wave function det|(0s)4(0p)12| is Q=8. Since the 16O wave function having Q=8 is restricted only to this double-closed-shell wave function having (0,0) SU(3) symmetry, we know that there hold

(140)A{[R4(r ,3νN)ϕ(0,4)(12C)](λ,μ)ϕ(α)}=0  for (λ,μ)≠(0,0).

Therefore Eq.(139) can be written as

(141)A{[R4,L( r ,3νN)ϕ(0,4)L(12C)]J=0ϕ(α)}=〈(4,0)L,(0,4)L||(0,0)0〉A{(0,0)ϕ(α)[R4(r ,3νN)ϕ(0,4)(12C)]}.

Thus the double-closed-shell wave function can be equivalently expressed in three ways by the use of 12C + α cluster-model wave functions, 12C(L=0) + α(S-wave), 12C(L=2) + α(D-wave), and 12C(L=4) + α(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 4ℏω excited state of 32S which is the main component of the superdeformed state and has the configuration with 4 particles raised from 1s0d-shell to 1p0f-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:

(142)(0,0,0)4(1,0,0)4(0,1,0)4(0,0,1)4(1,0,1)4​×(0,1,1)4(0,0,2)4(0,0,3)4​=n0A{X(0,0,24)(rO − O,8νN)ϕ2(16O)}g(r (32S),32νN),

(143)  r (32S)=132∑i=132  ri.

Here (nx,ny,nz) expresses the Cartesian harmonic oscillator single-particle function with nk being the number of oscillator quanta along k direction (k=x, y, or z). The coordinate rO−O stands for the relative coordinate between two 16O clusters and X(nx,ny,nz)(r,γ) stands for the Cartesian harmonic oscillator function of coordinate r with size parameter γ. This intrinsic wave function belongs to the SU(3) symmetry (λ,μ)=(24,0).

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 4n 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 4ℏω 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 4ℏω 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 4ℏω 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 (3/5)R2 with R standing for the nuclear radius. It is obtained by calculating 〈uf|r2|ui〉 by using uniform-density approximation for uf(r) and ui(r), u(r)≈(3/R3)1/2 for 0≤r≤R. If we adopt R≈3 fm for R for light nuclei like 12C and 16O, we have (3/5)R2≈5.4 fm2. The observed E0 value between the Hoyle state and the ground state in 12C is 5.4±0.2 fm2, and that between the first-excited 0+ state and the ground state in 16O is 3.55±0.21 fm2. The first-excited 0+ state in 16O is a well-known cluster state which has the structure of 12C(01+) + α(S) with α(S) standing for an α cluster moving around 12C(01+) with S-wave (L=0).

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 0+ (02+) state and the ground state (01+). Because of the duality property the ground state is expressed as A{R4,0(r,3νN)h0} with h0=Y0,0(r^)ϕ(0,4),L=0(12C)ϕ(α). On the other hand, since the first-excited 0+ state has the 12C(01+) + α(S) structure, it is expressed as A{χ0(r)h0}. Thus the E0 transition between two states is the transition between R4,0(r,3νN) and χ0(r). Since the E0 transition operator O(E0,16O) can be written as

(144)O(E0,16O)=12∑i=116(ri−r(16O))2=O(E0,12C)+O(E0,α)+1212×416r2,

we see that only the relative-motion part, (1/2)((12×4)/16)r2, contributes to the E0 transition. We can easily show that both the operators O(E0,12C) and O(E0,α) 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

(145)M(E0,02+↔01+)=〈(0s)4(0p)12,J=0|O(E0,16O)|​×A{χ0(r)[Y0(r^)ϕ0(12C)]J=0ϕ(α)}g(r(16O),16ν)〉

(146)​=12τ0,4τ0,6η6〈R40(r,ν)|r2|R60(r,ν)〉,

(147)τL,N=〈ΨL,N|A{ΨL,N}〉,

(148)ΨL,N=RN,L(rC−α,3ν)[YL(r^C−α)ϕL(12C)]J=0ϕ(α),

(149)A{χ0(r)[Y0(r^)ϕ0(12C)]J=0ϕ(α)}​​=​​∑N=6∞ηN(CNA{Ψ0,N}),

(150)||CNA{Ψ0,N}||=1.

The quantity τL,N 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 τ0,N appear in the form of ratio, τ0,4/τ0,6, and the magnitude of the ratio is close to unity, which implies that the effect of antisymmetrization has only little influence on the M(E0) value. The quantity η6 is the coefficient of the 2ℏω-jump component contained in A{χ0(r)[Y0(r^)ϕ0(12C)]J=0 ϕ(α)}, and its magnitude is around 0.4. Note that η6 is not percentage quantity, (η6)2. The E0 matrix element of the relative motion, 〈R40(r,ν)|r2|R60(r,ν)〉 is larger than the corresponding E0 matrix elements of the single-nucleon motion, 〈R00(r,ν)|r2|R20(r,ν)〉 and 〈R11(r,ν)|r2|R31(r,ν)〉, by about 50%. We thus see that the order of magnitude of M(E0,02+→01+) is the same as the single-nucleon strength.

In 16O, the E0 transition from the second-excited 0+ (03+) state to the ground state has been also observed and its observed strength is large as 4.03 fm2. The dominant configuration of this excited state is the cluster structure of 12C(21+) + α(D) 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 12C(21+) + α(D) 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 12C(21+) 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 03+ and 01+ states which is quite similar to that between 02+ and 01+ states.

The numerical values of M(E0,02+→01+) and M(E0,03+→01+) calculated by the analytical formulae mentioned above with suitable parameter values are 1.97 fm2 and 3.89 fm2, 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 R4,L(r,3ν) but has contributions from higher N components, RN,L(r,3ν) with N>4. 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 χuF,pq(  r pq) defined by

(151)A{χuF,pq( rpq)ϕ(Cp)ϕ(Cq)}=0,

where ϕ(Cj) is the internal wave functions of cluster Cj (j=p,q) and  rpq is the relative coordinate between Cp and Cq. In the case of α+α system, the Pauli-forbidden states are harmonic oscillator wave functions RN,l(rαα,νN/2)Yl,m(r^αα) with N=2n+l<4 [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 (j=1–n), every pair of clusters of Cp and Cq (p,q=1–n) should have their relative wave function orthogonal to their Pauli-forbidden states χuF,pq(  r pq). It means that the relative wave function χ(  ξ 1,⋯,  ξ n−1) of the many-cluster system should be orthogonal to all the Pauli-forbidden states of all pairs of clusters. Here   ξ k (k=1−(n−1)) are relative coordinates among clusters like Jacobi coordinates. In the case of Jacobi coordinates, we have   ξ k= Xk+1−(1/∑j=1kAj)∑j=1kAj Xj, where Xj 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 N [44] is useful which is defined as

(152)N( a1,⋯,  an−1;  a1′,⋯,  an−1′)​=〈δ( ξ1− a1)⋯δ( ξn−1− an−1)ϕ(C1)⋯ϕ(Cn)|A{δ( ξ1− a1′)⋯δ( ξn−1−  an−1′)ϕ(C1)⋯ϕ(Cn)}〉

The eigen functions and eigen values of N are denoted as χi( ξ1,⋯,  ξn−1) and μi, respectively, which satisfy [44, 161]

(153)∫d  a1′⋯d  an−1′N( a1,⋯,  an−1;  a1′,⋯,  an−1′)χi(  a1′,⋯,  an−1′)=μiχi( a1,⋯,  an−1),

(154)〈ϕ(C1)⋯ϕ(Cn)|A{χi(ξ1,⋯,  ξn−1)ϕ(C1)⋯ϕ(Cn)}〉​=μiχi( ξ1,⋯,  ξn−1).

The eigen functions {χi( ξ1,⋯,  ξn−1)} constitute an orthonormal complete set and any wave function χ( ξ1,⋯,  ξn−1) can be expanded by this complete set. The eigen functions with zero eigen value μi=0 are Pauli forbidden because A{χi(ξ1,⋯,  ξn−1)ϕ(C1)⋯ϕ(Cn)}=0 for μi=0, while the eigen functions with non-zero eigen value μi≠0 are Pauli allowed because A{χi( ξ1,⋯,  ξn−1)ϕ(C1)⋯ϕ(Cn)}≠0 for μi≠0. The functional space spanned by all the Pauli-allowed eigen functions is designated as Sn.

Two-cluster Pauli-forbidden state χuF,pq( rpq) is also a many-cluster Pauli-forbidden state because there holds A{χuF,pq( rpq)ϕ(C1)⋯ϕ(Cn)}=0. Then we have

(155)0=〈χuF,pq(rpq)ϕ(C1)⋯ϕ(Cn)}|A{χi( ξ1,⋯,  ξn−1)〉,

(156)ϕ(C1)⋯ϕ(Cn)}〉=μi〈χuF,pq( rpq)|χi( ξ1,⋯,  ξn−1)〉.

From this relation we know that every Pauli-allowed eigen function χi( ξ1,⋯,  ξn−1) with μi≠ 0 is orthogonal to all the two-cluster Pauli-forbidden state {χuF,pq( rpq)}. Therefore any Pauli-allowed state which is expressed by a linear combination of Pauli-allowed eigenfunctions {χi} is orthogonal to all the two-cluster Pauli-forbidden states {χuF,pq( rpq)}. Now we consider the functional space Sn′ which is spanned by many-cluster wave functions that are orthogonal to all the two-cluster Pauli-forbidden state {χuF,pq( rpq)}. Clearly the space Sn′ contains in it the space Sn. However explicit comparisons of Sn′ and Sn in many kinds of many-cluster systems have shown that Sn′ and Sn are equal to each other (Sn′=Sn) [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 {χuF,pq( rpq)}. Thus the functional space of many-cluster OCM is Sn′. The fact that Sn′=Sn 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 N.

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 μi are given below their (λ,μ). The Pauli-allowed eigen state with lowest harmonic oscillator quanta has N=12. This eigen state which we denote as χ12,(0,0)(ξ1,  ξ2,  ξ3,) has the following property,

(157)A{χ12,(0,0)( ξ1,  ξ2,  ξ3)ϕ(α1)⋯ϕ(α4)}​∝​det|(0s)4(0p)12|.

This is of course the duality property of the double closed-shell wave function. When we denote the Pauli-allowed eigen state with N=13 as χ13,(2,1)( ξ1,  ξ2,  ξ3,), we have the equality that A{χ13,(2,1)(  ξ 1,  ξ 2,  ξ 3)ϕ(α1)⋯ϕ(α4)} is the same as the 1p-1h shell-model wave function C|(0p)−1(1s0d)1;(λ,μ)=(2,1)〉. 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 N=12 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 (λ≠0,μ=0) corresponding to prolate deformation and also many Pauli-allowed states with SU(3) label (λ=0,μ≠0) corresponding to oblate deformation.

The most compact Pauli-allowed state has the lowest number of oscillator quanta N=Nlowest 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 N=Nlowest+1 often expresses a dominant component of a shell-model state with 1p–1h 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 12C+α 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 (λ,μ)=(0,12) and this component does not have the clustering degree of freedom of 16O + 16O configuration whose SU(3) symmetry is (N, 0) with N≥24. 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 |(0s)4(0p)4,[44]J=0〉 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 Kπ=0− built upon the Jπ=1− 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 Kπ=0− 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 Jπ=1− state of the Kπ=0− 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, |(0s)4(0p)12(1s0d)4;[44444](λ,μ)=(8,0),J=0〉. 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 Kπ=0− 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 〈V^ls〉 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 Kπ=0− 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 Kπ=0− band has good clustering character of 16O + α. It is supported by the small expectation values 〈V^ls〉 for all the states of the Kπ=0− band.

In the case of 44Ti, just like the 20Ne case, there exists inversion-doublet bands composed of the ground rotational band and the Kπ=0− band built upon the Jπ=1− 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 Kπ=0− 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 Jπ=0+ excited state whose 40Ca + α component is larger than 40% and has relative-motion nodal number nr=6. Here the relative-motion nodal number nr means the nodal number of the reduced width amplitude Y(r) defined as

(158)Y(r)​=​r〈δ(rrel​−​r)r2Y0.0(r^rel)ϕ(40Ca)ϕ(α)|ΦAMD〉.

The relative-motion nodal number nr=6 is the smallest nodal number of Y(r) which is allowed by the inter-cluster Pauli principle. The value of nr of the Jπ=1− of the Kπ=0− band is also nr=6. 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 Φ(12,0)=|40Ca-core, (1p0f)4;[410](λ,μ)=(12,0)〉 is less than 40% because Φ(12,0) 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.

References

[1]

K. Wildermuth and Y. C. Tang, A Unified Theory of the Nucleus, Vieweg, 1977

[2]

Y. Fujiwara, H. Horiuchi, K. Ikeda, M. Kamimura, K. Katō, Y. Suzuki, and E. Uegaki, Comprehensive Study of Alpha-Nuclei, Prog. Theor. Phys. Suppl. 68, 29 (1980)

[3]

M. Freer, H. Horiuchi, Y. Kanada-En’yo, D. Lee, and U.-G. Meißner, Microscopic clustering in light nuclei, Rev. Mod. Phys. 90, 035004 (2018)

[4]

T. Furuta and A. Ono, Nuclear liquid-gas phase transition studied with antisymmetrized molecular dynamics, Phys. Rev. C 74, 014612 (2006)

[5]

T. Yamada, Y. Funaki, H. Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, Nuclear Alpha-Particle Condensates, in: Clusters in Nuclei, Vol.2, Lecture Notes in Physics, edited by C. Beck, Springer Berlin Heidelberg, Berlin, Heidelberg, 2012, pp 229–298

[6]

H. Horiuchi, K. Ikeda, and K. Katō, Recent developments in nuclear cluster physics, Prog. Theor. Phys. Suppl. 192, 1 (2012)

[7]

Y. Funaki, H. Horiuchi, and A. Tohsaki, Cluster models from RGM to alpha condensation and beyond, Prog. Part. Nucl. Phys. 82, 78 (2015)

[8]

A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, Alpha cluster condensation in 12C and 16O, Phys. Rev. Lett. 87, 192501 (2001)

[9]

Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, Analysis of previous microscopic calculations for the second 0+ state in 12C in terms of 3α particle Bose-condensed state, Phys. Rev. C 67, 051306 (2003)

[10]

B. Zhou, Y. Funaki, H. Horiuchi, Z. Ren, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, and T. Yamada, Nonlocalized clustering: A new concept in nuclear cluster structure physics, Phys. Rev. Lett. 110, 262501 (2013)

[11]

B. Zhou, Y. Funaki, H. Horiuchi, Z. Ren, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, and T. Yamada, Nonlocalized cluster dynamics and nuclear molecular structure, Phys. Rev. C 89, 034319 (2014a)

[12]

H. Horiuchi, Three-alpha model of 12C orthogonality condition model as an approximation of resonating group treatment, Prog. Theor. Phys. 51, 1266 (1974)

[13]

E. Uegaki, S. Okabe, Y. Abe, and H. Tanaka, Structure of the Excited States in 12C (I), Prog. Theor. Phys. 57, 1262 (1977)

[14]

M. Kamimura, Transition densities between the 01+, 21+, 41+, 02+, 22+, 11- and 31- states in 12C derived from the three-alpha resonating-group wave functions, Nucl. Phys. A 351, 456 (1981)

[15]

J. Carlson, S. Gandolfi, F. Pederiva, S. C. Pieper, R. Schiavilla, K. E. Schmidt, and R. B. Wiringa, Quantum Monte Carlo methods for nuclear physics, Rev. Mod. Phys. 87, 1067 (2015)

[16]

Y. Funaki, Hoyle band and alpha condensation in 12C, Phys. Rev. C 92, 021302 (2015)

[17]

Y. Funaki, Monopole excitation of the Hoyle state and linear-chain state in 12C, Phys. Rev. C 94, 024344 (2016)

[18]

B. Zhou, A. Tohsaki, H. Horiuchi, and Z. Ren, Breathing-like excited state of the Hoyle state in 12C, Phys. Rev. C 94, 044319 (2016)

[19]

Y. Funaki, T. Yamada, H. Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, α-particle condensation in 16O studied with a full four-body orthogonality Condition model calculation, Phys. Rev. Lett. 101, 082502 (2008)

[20]

Y. Funaki, T. Yamada, A. Tohsaki, H. Horiuchi, G. Röpke, and P. Schuck, Microscopic study of 4α-particle condensation with inclusion of resonances, Phys. Rev. C 82, 024312 (2010)

[21]

H. Horiuchi and K. Ikeda, A molecule-like structure in atomic nuclei of 16O* and 20Ne, Prog. Theor. Phys. 40, 277 (1968)

[22]

B. Zhou, Z. Ren, C. Xu, Y. Funaki, T. Yamada, A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, New concept for the ground-state band in 20Ne within a microscopic cluster model, Phys. Rev. C 86, 014301 (2012)

[23]

P. Navrátil, J. P. Vary, and B. R. Barrett, Large-basis ab initio no-core shell model and its application to 12C, Phys. Rev. C 62, 054311 (2000)

[24]

S. Quaglioni, C. Romero-Redondo, and P. Navrátil, Three-cluster dynamics within an ab initio framework, Phys. Rev. C 88, 034320 (2013)

[25]

J.-P. Ebran, E. Khan, T. Nikšić, and D. Vretenar, How atomic nuclei cluster, Nature 487, 341 (2012)

[26]

P. W. Zhao, N. Itagaki, and J. Meng, Rod-shaped nuclei at extreme spin and isospin, Phys. Rev. Lett. 115, 022501 (2015)

[27]

Y. Iwata, T. Ichikawa, N. Itagaki, J. A. Maruhn, and T. Otsuka, Examination of the stability of a rod-shaped structure in 24Mg, Phys. Rev. C 92, 011303 (2015)

[28]

T. Yamada, Y. Funaki, T. Myo, H. Horiuchi, K. Ikeda, G. Röpke, P. Schuck, and A. Tohsaki, Isoscalar monopole excitations in 16O: Alpha-cluster states at low energy and mean-field-type states at higher energy, Phys. Rev. C 85, 034315 (2012b)

[29]

W. B. He, Y. G. Ma, X. G. Cao, X. Z. Cai, and G. Q. Zhang, Giant dipole resonance as a fingerprint of α clustering configurations in 12C and 16O, Phys. Rev. Lett. 113, 032506 (2014)

[30]

Y. Chiba, M. Kimura, and Y. Taniguchi, Isoscalar dipole transition as a probe for asymmetric clustering, Phys. Rev. C 93, 034319 (2016)

[31]

K. Ikeda, N. Takigawa, and H. Horiuchi, The systematic structure-change into the molecule-like structures in the self-conjugate 4n nuclei, Prog. Theor. Phys. Suppl. E68, 464 (1968)

[32]

F. Hoyle, On nuclear reactions occuring in very hot STARS (I): The synthesis of elements from carbon to nickel, Astrophys. J., Suppl. Ser. 1, 121 (1954)

[33]

C. W. Cook, W. A. Fowler, C. C. Lauritsen, and T. Lauritsen, 12B, 12C, and the red giants, Phys. Rev. 107, 508 (1957)

[34]

J. P. Vary, P. Maris, H. Potter, M. A. Caprio, R. Smith, S. Binder, A. Calci, S. Fischer, J. Langhammer, R. Roth, H. M. Aktulga, E. Ng, C. Yang, D. Oryspayev, M. Sosonkina, E. Saule, and U. Çatalyürek, Ab initio no core shell model – recent results and further prospects, arXiv: 1507.04693[nucl-th] (2015)

[35]

H. Morinaga, Interpretation of some of the excited states of 4n self-conjugate nuclei, Phys. Rev. 101, 254 (1956)

[36]

H. Morinaga, On the spin of a broad state around 10 MeV in 12C, Phys. Lett. 21, 78 (1966)

[37]

Y. Suzuki, H. Horiuchi, and K. Ikeda, Study of α chain states through their decay widths, Prog. Theor. Phys. 47, 1517 (1972)

[38]

S. Saito, Effect of Pauli principle in scattering of two clusters, Prog. Theor. Phys. 40, 893 (1968)

[39]

S. Saito, Interaction between clusters and Pauli principle, Prog. Theor. Phys. 41, 705 (1969)

[40]

S. Saito, Theory of resonating group method and generator coordinate method, and orthogonality condition model, Prog. Theor. Phys. Suppl. 62, 11 (1977)

[41]

D. Brink, In: Proceedings of the International School of Physics ”Enrico Fermi”, Course 36, Academic Press, New York, London, 1966

[42]

P. Descouvemont and D. Baye, Microscopic theory of the 8Be(α,γ)12C reaction in a three-cluster model, Phys. Rev. C 36, 54 (1987)

[43]

H. Horiuchi, Generator coordinate treatment of composite particle reaction and molecule-like structures, Prog. Theor. Phys. 43, 375 (1970)

[44]

H. Horiuchi, Kernels of GCM, RGM and OCM and their calculation methods, Prog. Theor. Phys. Suppl. 62, 90 (1977a)

[45]

Y. Funaki, H. Horiuchi, A. Tohsaki, P. Schuck, and G. Röpke, Description of 8Be as deformed gas-like two-alpha-particle states, Prog. Theor. Phys. 108, 297 (2002)

[46]

A. Tohsaki-Suzuki, M. Kamimura, and K. Ikeda, Microscopic study of the interaction between complex nuclei, Prog. Theor. Phys. Suppl. 68, 359 (1980)

[47]

Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, Resonance states in 12C and α-particle condensation, Eur. Phys. J. A 24, 321 (2005)

[48]

Y. Funaki, H. Horiuchi, W. von Oertzen, G. Röpke, P. Schuck, A. Tohsaki, and T. Yamada, Concepts of nuclear α-particle condensation, Phys. Rev. C 80, 064326 (2009)

[49]

A. Volkov, Equilibrium deformation calculations of the ground state energies of 1p shell nuclei, Nucl. Phys. 74, 33 (1965)

[50]

B. Zhou, Y. Funaki, A. Tohsaki, H. Horiuchi, and Z. Ren, The container picture with two-alpha correlation for the ground state of 12C, Prog. Theor. Exp. Phys. 2014, 101D01 (2014b)

[51]

M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann-Cosel, and A. Richter, Structure of the Hoyle state in 12C, Phys. Rev. Lett. 98, 032501 (2007)

[52]

M. Chernykh, H. Feldmeier, T. Neff, P. von Neumann-Cosel, and A. Richter, Pair decay width of the Hoyle state and its role for stellar carbon production, Phys. Rev. Lett. 105, 022501 (2010)

[53]

Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, Inelastic form factors to alpha-particle condensate states in 12C and 16O: What can we learn? Eur. Phys. J. A 28, 259 (2006a)

[54]

S. Elhatisari, D. Lee, G. Rupak, E. Epelbaum, H. Krebs, T. A. Lähde, T. Luu, and U.-G. Meißner, Ab Initio alpha–alpha scattering, Nature 528, 111 (2015)

[55]

S. Elhatisari, N. Li, A. Rokash, J. M. Alarcón, D. Du, N. Klein, B.-n. Lu, U.-G. Meißner, E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, and G. Rupak, Nuclear binding near a quantum phase transition, Phys. Rev. Lett. 117, 132501 (2016)

[56]

E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, and U.-G. Meißner, Structure and rotations of the Hoyle state, Phys. Rev. Lett. 109, 252501 (2012)

[57]

E. Epelbaum, H. Krebs, T. A. Lähde, D. Lee, and U.-G. Meißner, Viability of carbon-based life as a function of the light quark mass, Phys. Rev. Lett. 110, 112502 (2013)

[58]

A. C. Dreyfuss, K. D. Launey, T. Dytrych, J. P. Draayer, and C. Bahri, Hoyle state and rotational features in carbon-12 within a No-Core shell-model framework, Phys. Lett. B 727, 511 (2013)

[59]

T. Dytrych, K. D. Launey, J. P. Draayer, P. Maris, J. P. Vary, E. Saule, U. Catalyurek, M. Sosonkina, D. Langr, and M. A. Caprio, Collective Modes in Light Nuclei from First Principles, Phys. Rev. Lett. 111, 252501 (2013)

[60]

G. K. Tobin, M. C. Ferriss, K. D. Launey, T. Dytrych, J. P. Draayer, A. C. Dreyfuss, and C. Bahri, Symplectic No-Core Shell-Model approach to intermediate-mass nuclei, Phys. Rev. C 89, 034312 (2014)

[61]

R. B. Wiringa, S. C. Pieper, J. Carlson, and V. R. Pandharipande, Quantum Monte Carlo calculations of A=8 nuclei, Phys. Rev. C 62, 014001 (2000)

[62]

P. Schuck, Y. Funaki, H. Horiuchi, G. Röpke, A. Tohsaki, and T. Yamada, Alpha particle clusters and their condensation in nuclear systems, Phys. Scr. 91, 123001 (2016)

[63]

M. Itoh, H. Akimune, M. Fujiwara, U. Garg, H. Hashimoto, T. Kawabata, K. Kawase, S. Kishi, T. Murakami, K. Nakanishi, Y. Nakatsugawa, B. K. Nayak, S. Okumura, H. Sakaguchi, H. Takeda, S. Terashima, M. Uchida, Y. Yasuda, M. Yosoi, and J. Zenihiro, Study of the cluster state at Ex=10.3 MeV in 12C, Nucl. Phys. A 738, 268 (2004)

[64]

M. Freer, H. Fujita, Z. Buthelezi, J. Carter, R. W. Fearick, S. V. Förtsch, R. Neveling, S. M. Perez, P. Papka, F. D. Smit, J. A. Swartz, and I. Usman, 2+ excitation of the 12C Hoyle state, Phys. Rev. C 80, 041303 (2009)

[65]

M. Itoh, H. Akimune, M. Fujiwara, U. Garg, N. Hashimoto, T. Kawabata, K. Kawase, S. Kishi, T. Murakami, K. Nakanishi, Y. Nakatsugawa, B. K. Nayak, S. Okumura, H. Sakaguchi, H. Takeda, S. Terashima, M. Uchida, Y. Yasuda, M. Yosoi, and J. Zenihiro, Candidate for the 2+ excited Hoyle state at Ex ~10 MeV in 12C, Phys. Rev. C 84, 054308 (2011)

[66]

H. O. U. Fynbo and M. Freer, Rotations of the Hoyle State in 12C, Physics 4, 94 (2011)

[67]

W. R. Zimmerman, N. E. Destefano, M. Freer, M. Gai, and F. D. Smit, Further evidence for the broad 22+ state at 9.6 MeV in 12C, Phys. Rev. C 84, 027304 (2011)

[68]

W. R. Zimmerman, M. W. Ahmed, B. Bromberger, S. C. Stave, A. Breskin, V. Dangendorf, T. Delbar, M. Gai, S. S. Henshaw, J. M. Mueller, C. Sun, K. Tittelmeier, H. R. Weller, and Y. K. Wu, Unambiguous identification of the second 2+ state in 12C and the structure of the Hoyle state, Phys. Rev. Lett. 110, 152502 (2013)

[69]

M. Itoh, H. Akimune, M. Fujiwara, U. Garg, T. Kawabata, K. Kawase, Tetsuya Murakami, K. Nakanishi, Y. Nakatsugawa, H. Sakaguchi, S. Terashima, Makoto Uchida, Y. Yasuda, M. Yosoi, and J. Zenihiro, Nature of 10 MeV state in 12C, J. Phys.: Conf. Ser. 436, 012006 (2013)

[70]

C. Kurokawa and K. Katō, New broad 0+ state in 12C, Phys. Rev. C 71, 021301 (2005)

[71]

S.-I. Ohtsubo, Y. Fukushima, M. Kamimura, and E. Hiyama, Complex-scaling calculation of three-body resonances using complex-range Gaussian basis functions: Application to 3α resonances in 12C, Prog. Theor. Exp. Phys. 2013, 073D02 (2013)

[72]

J. Aguilar and J. M. Combes, A class of analytic perturbations for one-body Schrödinger Hamiltonians, Commun. Math. Phys. 22, 269 (1971)

[73]

E. Balslev and J. M. Combes, Spectral properties of many-body Schrödinger operators with dilatation-analytic interactions, Comm. Math. Phys. 22, 280 (1971)

[74]

A. T. Kruppa, R. G. Lovas, and B. Gyarmati, Complex scaling in the cluster model: Resonances in 8Be, Phys. Rev. C 37, 383 (1988)

[75]

A. T. Kruppa and K. Katō, Resonances in complex-scaled orthogonality condition model of nuclear cluster system, Prog. Theor. Phys. 84, 1145 (1990)

[76]

S. Aoyama, T. Myo, K. Katō, and K. Ikeda, The complex scaling method for many-body resonances and its applications to three-body resonances, Prog. Theor. Phys. 116, 1 (2006)

[77]

Y. Kanada-En’yo, The structure of ground and excited states of 12C, Prog. Theor. Phys. 117, 655 (2007)

[78]

Y. Funaki, H. Horiuchi, and A. Tohsaki, New treatment of resonances with a bound state approximation using a pseudo-potential, Prog. Theor. Phys. 115, 115 (2006b)

[79]

T. Yamada, Y. Funaki, H. Horiuchi, K. Ikeda, and A. Tohsaki, Monopole excitation to cluster states, Prog. Theor. Phys. 120, 1139 (2008a)

[80]

T. Yamada and P. Schuck, Single α-particle orbits and Bose-Einstein condensation in 12C, Eur. Phys. J. A 26, 185 (2005)

[81]

A. Hasegawa and S. Nagata, Ground state of 6Li, Prog. Theor. Phys. 45, 1786 (1971)

[82]

F. Tanabe, A. Tohsaki, and R. Tamagaki, α-α scattering at intermediate energies applicability of orthogonality condition model and upper limit of isoscalar meson-nucleon coupling constants inferred from potential tail, Prog. Theor. Phys. 53, 677 (1975)

[83]

F. Ajzenberg-Selove, Energy levels of light nuclei A = 16–17, Nucl. Phys. A 460, 1 (1986)

[84]

T. Wakasa, E. Ihara, K. Fujita, Y. Funaki, K. Hatanaka, H. Horiuchi, M. Itoh, J. Kamiya, G. Röpke, H. Sakaguchi, N. Sakamoto, Y. Sakemi, P. Schuck, Y. Shimizu, M. Takashina, S. Terashima, A. Tohsaki, M. Uchida, H. P. Yoshida, and M. Yosoi, New candidate for an alpha cluster condensed state in 16O(α,A′) at 400 MeV, Phys. Lett. B 653, 173 (2007)

[85]

Y. Suzuki, Structure Study of 16O by 12C + alpha cluster-coupling model (I), Prog. Theor. Phys. 55, 1751 (1976a)

[86]

Y. Suzuki, Structure study of 16O by 12C + α cluster-coupling model (II), Prog. Theor. Phys. 56, 111 (1976b)

[87]

A. M. Lane and R. G. Thomas, R-matrix theory of nuclear reactions, Rev. Mod. Phys. 30, 257 (1958)

[88]

T. Yamada, Y. Funaki, H. Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, Criterion for Bose–Einstein condensation in traps and self-bound systems, Phys. Rev. A 78, 035603 (2008b)

[89]

T. Yamada, Y. Funaki, H. Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, Internal one-particle density matrix for Bose–Einstein condensates with finite number of particles in a harmonic potential, Phys. Rev. C 79, 054314 (2009)

[90]

A. Tohsaki, New effective internucleon forces in microscopic α-cluster model, Phys. Rev. C 49, 1814 (1994)

[91]

H. Matsumura and Y. Suzuki, A microscopic analysis of the amount of α-condensation in 12C, Nucl. Phys. A 739, 238 (2004)

[92]

A. Arima, Elliott’s SU(3) model and its developments in nuclear physics, J. Phys. G Nucl. Part. Phys. 25, 581 (1999)

[93]

M. Lyu, Z. Ren, B. Zhou, Y. Funaki, H. Horiuchi, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, and T. Yamada, Investigation of 9Be from a nonlocalized clustering concept, Phys. Rev. C 91, 014313 (2015)

[94]

A. Tohsaki, The interaction of a nucleon and alpha condensate, Int. J. Mod. Phys. E 17, 2106 (2008)

[95]

J. Hiura, Y. Abe, S. Saitō, and O. Endō, Alpha-cluster plus 16O-core model for 20Ne and neighboring nuclei, Prog. Theor. Phys. 42, 555 (1969)

[96]

T. Matsuse, M. Kamimura, and Y. Fukushima, Study of the alpha-clustering structure of 20Ne based on the resonating group method for 16O+α, Prog. Theor. Phys. 53, 706 (1975)

[97]

M. LeMere, Y. C. Tang, and D. R. Thompson, Study of the A+16O system with the resonating-group method, Phys. Rev. C 14, 23 (1976)

[98]

S. Hara, K. T. Hecht, and Y. Suzuki, Alpha cluster formation in 20Ne in the cluster-orbital shell model (I), Prog. Theor. Phys. 84, 254 (1990)

[99]

S. Hara, K. Ogawa, and Y. Suzuki, Alpha Cluster Formation in 20Ne in the cluster-orbital shell model (II), Prog. Theor. Phys. 88, 329 (1992)

[100]

T. Yamaya, K. Katori, M. Fujiwara, S. Katō, and S. Ohkubo, Alpha-cluster study of 40Ca and 44Ti by the (6Li,d) reaction, Prog. Theor. Phys. 132, 73 (1998)

[101]

A. Astier, P. Petkov, M.-G. Porquet, D. S. Delion, and P. Schuck, Novel manifestation of α-clustering structures: New α+208Pb states in 212Po revealed by their enhanced E1 decays, Phys. Rev. Lett. 104, 042701 (2010)

[102]

Z. Ren and B. Zhou, Alpha-clustering effects in heavy nuclei, Front. Phys. 13, 132110 (2018)

[103]

Y. Kanada-En’yo, Description of an α-cluster tail in 8Be and 20Ne: Delocalization of the α cluster by quantum penetration, Prog. Theor. Exp. Phys. 2014, 103D03 (2014)

[104]

T. Suhara, Y. Funaki, B. Zhou, H. Horiuchi, and A. Tohsaki, One-dimensional α condensation of α-linear-chain states, J. Phys.: Conf. Ser. 569, 012008 (2014)

[105]

T. Ichikawa, J. A. Maruhn, N. Itagaki, and S. Ohkubo, Linear chain structure of four-α clusters in 16O, Phys. Rev. Lett. 107, 112501 (2011a)

[106]

Y. Funaki, T. Yamada, E. Hiyama, B. Zhou, and K. Ikeda, Container structure of alpha-alpha-Lambda clusters in 9-lambda–Beryrium, Prog. Theor. Exp. Phys. 2014, 113D01 (2014)

[107]

V. I. Kukulin and V. M. Krasnopol’sky, Description of few-body systems via analytical continuation in coupling constant, J. Phys. A: Math. Gen. 10, L33 (1977)

[108]

M. Freer, S. Almaraz-Calderon, A. Aprahamian, N. I. Ashwood, M. Barr, B. Bucher, P. Copp, M. Couder, N. Curtis, X. Fang, F. Jung, S. Lesher, W. Lu, J. D. Malcolm, A. Roberts, W. P. Tan, C. Wheldon, and V. A. Ziman, Evidence for a new 12C state at 13.3 MeV, Phys. Rev. C 83, 034314 (2011)

[109]

R. Bijker and F. Iachello, Cluster states in nuclei as representations of a U(N+1) group, Phys. Rev. C 61, 067305 (2000)

[110]

D. J. Marín-Lámbarri, R. Bijker, M. Freer, M. Gai, T. Kokalova, D. J. Parker, and C. Wheldon, Evidence for triangular D3h symmetry in 12C, Phys. Rev. Lett. 113, 012502 (2014)

[111]

S. Ohkubo and Y. Hirabayashi, Alpha-particle condensate states in 16O, Phys. Lett. B 684, 127 (2010)

[112]

Y. Funaki, T. Yamada, H. Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, Alpha cluster states and condensation in 16O, Prog. Theor. Phys. 196, 439 (2012)

[113]

T. Suhara and Y. Kanada-En’yo, Effects of alpha-cluster breaking on 3α-cluster structures in 12C, Phys. Rev. C 91, 024315 (2015)

[114]

E. Uegaki, Y. Abe, S. Okabe, and H. Tanaka, Structure of the excited states in 12C (II), Prog. Theor. Phys. 62, 1621 (1979)

[115]

R. Imai, T. Tada, and M. Kimura, Real-time evolution method and its application to the 3α cluster system, Phys. Rev. C 99, 064327 (2019)

[116]

R. Bijker and F. Iachello, The algebraic cluster model: Three-body clusters, Ann. Phys. 298, 334 (2002)

[117]

B. Zhou, Y. Funaki, H. Horiuchi, M. Kimura, Z. Ren, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, and T. Yamada, Nonlocalized motion in a two-dimensional container of α particles in 3- and 4- states of 12C, Phys. Rev. C 99, 051303 (2019)

[118]

R. Bijker and F. Iachello, Evidence for tetrahedral symmetry in 16O, Phys. Rev. Lett. 112, 152501 (2014)

[119]

Y. Kanada-En’yo, Tetrahedral 4α and 12C+α cluster structures in 16O, Phys. Rev. C 96, 034306 (2017)

[120]

Y. Kanada-En’yo and H. Horiuchi, Structure of light unstable nuclei studied with antisymmetrized molecular dynamics, Prog. Theor. Phys. Suppl. 142, 205 (2001)

[121]

M. Freer, The clustered nucleus — cluster structures in stable and unstable nuclei, Rep. Prog. Phys. 70, 2149 (2007)

[122]

Z. H. Yang, Y. L. Ye, Z. H. Li, J. L. Lou, J. S. Wang, D. X. Jiang, Y. C. Ge, Q. T. Li, H. Hua, X. Q. Li, F. R. Xu, J. C. Pei, R. Qiao, H. B. You, H. Wang, Z. Y. Tian, K. A. Li, Y. L. Sun, H. N. Liu, J. Chen, J. Wu, J. Li, W. Jiang, C. Wen, B. Yang, Y. Y. Yang, P. Ma, J. B. Ma, S. L. Jin, J. L. Han, and J. Lee, Observation of enhanced monopole strength and clustering in 12Be, Phys. Rev. Lett. 112, 162501 (2014)

[123]

T. Baba and M. Kimura, Structure and decay pattern of the linear-chain state in 14C, Phys. Rev. C 94, 044303 (2016)

[124]

Y. Kanada-En’yo, M. Kimura, and A. Ono, Antisymmetrized molecular dynamics and its applications to cluster phenomena, Prog. Theor. Exp. Phys. 2012, 01A202 (2012)

[125]

M. Lyu, Z. Ren, B. Zhou, Y. Funaki, H. Horiuchi, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, and T. Yamada, Investigation of 10Be and its cluster dynamics with the nonlocalized clustering approach, Phys. Rev. C 93, 054308 (2016)

[126]

R. Tamagaki, Potential models of nuclear forces at small distances, Prog. Theor. Phys. 39, 91 (1968)

[127]

N. Yamaguchi, T. Kasahara, S. Nagata, and Y. Akaishi, Effective interaction with three-body effects, Prog. Theor. Phys. 62, 1018 (1979)

[128]

F. Kobayashi and Y. Kanada-En’yo, Novel cluster states in 10Be, Phys. Rev. C 86, 064303 (2012)

[129]

T. Suhara and Y. Kanada-En’yo, Quadrupole deformation β and γ constraint in a framework of antisymmetrized molecular dynamics, Prog. Theor. Phys. 123, 303 (2010)

[130]

E. Hiyama and T. Yamada, Structure of light hypernuclei, Prog. Part. Nucl. Phys. 63, 339 (2009)

[131]

Y. Funaki, M. Isaka, E. Hiyama, T. Yamada, and K. Ikeda, Multi-cluster dynamics in Λ13C and analogy to clustering in 12C, Phys. Lett. B 773, 336 (2017)

[132]

Y. Yamamoto, T. Motoba, T. Fukuda, M. Takahashi, and K. Ikeda, Formation and transition of strangeness =−2 nuclear systems, Prog. Theor. Phys. 117, 281 (1994)

[133]

T. Yamada, T. Motoba, K. Ikeda, and H. Bandō, Structure study of typical light hypernuclei, Prog. Theor. Phys. 81, 104 (1985)

[134]

E. Hiyama, M. Kamimura, T. Motoba, T. Yamada, and Y. Yamamoto, Three- and four-body cluster models of Hypernuclei Using the G-Matrix ΛN Interaction Λ9Be, Λ13C, ΛΛ6He and ΛΛ10Be, Prog. Theor. Phys. 97, 881 (1997)

[135]

O. Hashimoto and H. Tamura, Spectroscopy of Λ hypernuclei, Prog. Part. Nucl. Phys. 57, 564 (2006)

[136]

T. Yamada and Y. Funaki, Alpha-cluster structures and monopole excitations in 13C, Phys. Rev. C 92, 034326 (2015)

[137]

Q. Zhao, Z. Ren, M. Lyu, H. Horiuchi, Y. Funaki, G. Röpke, P. Schuck, A. Tohsaki, C. Xu, T. Yamada, and B. Zhou, Investigation of the 9B nucleus and its cluster-nucleon correlations, Phys. Rev. C 97, 054323 (2018)

[138]

W. Sünkel, A generator coordinate technique for arbitrary cluster widths, Phys. Lett. B 65, 419 (1976)

[139]

M. Colonna and P. Chomaz, Spinodal decomposition in nuclear molecular dynamics, Phys. Lett. B 436, 1 (1998)

[140]

T. Maruyama, K. Niita, and A. Iwamoto, Extension of quantum molecular dynamics and its application to heavy-ion collisions, Phys. Rev. C 53, 297 (1996)

[141]

Y. Tang, M. LeMere, and D. Thompsom, Resonating-group method for nuclear many-body problems, Phys. Rep. 47, 167 (1978)

[142]

M. Dufour and P. Descouvemont, Molecular states in a microscopic four α model, Nucl. Phys. A 927, 134 (2014)

[143]

B. Zhou, New trial wave function for the nuclear cluster structure of nuclei, Prog. Theor. Exp. Phys. 2018, 041D01 (2018)

[144]

A. Tohsaki-Suzuki, New computational methods of analytical derivation of exchange kernels, Prog. Theor. Phys. Suppl. 62, 191 (1977)

[145]

Y. Tang, Microscopic cluster theory for nuclear systems, Nucl. Phys. A 463, 377 (1987)

[146]

K. Toshitaka, M. Takehiro, and A. Akito, Effect of breathing excitations of the triton nucleus on the alpha-triton cluster structure of 7Li, Nucl. Phys. A 414, 185 (1984)

[147]

N. Itagaki, M. Kimura, C. Kurokawa, M. Ito, and W. von Oertzen, Alpha condensed state with a core nucleus, Phys. Rev. C 75, 037303 (2007)

[148]

N. Itagaki, T. Kokalova, and W. Von Oertzen, Three alpha state around 40Ca, Mod. Phys. Lett. A 25, 1947 (2010)

[149]

T. Ichikawa, N. Itagaki, T. Kawabata, T. Kokalova, and W. von Oertzen, Gas-like state of α clusters around a 16O core in 24Mg, Phys. Rev. C 83, 061301 (2011b)

[150]

A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, Colloquium status of α-particle condensate structure of the Hoyle state, Rev. Mod. Phys. 89, 011002 (2017)

[151]

B. Bayman and A. Bohr, On the connection between the cluster model and the SU3 coupling scheme for particles in a harmonic oscillator potential, Nucl. Phys. 9, 596 (1958)

[152]

H. Horiuchi, Coexistence of cluster states and mean-field-type states, J. Phys. G Nucl. Part. Phys. 37, 064021 (2010)

[153]

J. K. Perring and T. H. R. Skyrme, The alpha-particle and shell models of the nucleus, Proc. Phys. Soc. A 69, 600 (1956)

[154]

K. Wildermuth and T. Kanellopoulos, The “cluster model” of the atomic nuclei, Nucl. Phys. 7, 150 (1958a)

[155]

K. Wildermuth and T. Kanellopoulos, On the structure of rotational levels of atomic nuclei, Nucl. Phys. 9, 449 (1958b)

[156]

M. G. Mayer and J. H. D. Jensen, Elementary Theory of Nuclear Shell Structure, Wiley, 1960

[157]

R. K. Shelline and K. Wildermuth, Experimental evidence for cluster structures in light and medium weight nuclei, Nucl. Phys. 21, 196 (1960)

[158]

M. Kimura and H. Horiuchi, 16O +  16O molecular band and their relation to the superdeformed band of 32S, Nucl. Phys. A 738, 236 (2004)

[159]

T. Kawabata, H. Akimune, H. Fujita, Y. Fujita, M. Fujiwara, K. Hara, K. Hatanaka, M. Itoh, Y. Kanada-En’yo, S. Kishi, K. Nakanishi, H. Sakaguchi, Y. Shimbara, A. Tamii, S. Terashima, M. Uchida, T. Wakasa, Y. Yasuda, H. Yoshida, and M. Yosoi, Cluster structure in 11B, Phys. Lett. B 646, 6 (2007)

[160]

R. Tamagaki and H. Tanaka, Repulsive core of effective α-α potential, Prog. Theor. Phys. 34, 191 (1965)

[161]

H. Horiuchi, Multi-cluster allowed states and spectroscopic amplitude of cluster transfer, Prog. Theor. Phys. 58, 204 (1977b)

[162]

Y. Fujiwara and H. Horiuchi, Generator coordinate theory of normalization kernels of cluster systems (II), Prog. Theor. Phys. 65, 1632 (1981)

[163]

M. Kimura, Deformed-basis antisymmetrized molecular dynamics and its application to 20Ne, Phys. Rev. C 69, 044319 (2004)

[164]

Y. Kanada-En’yo and H. Horiuchi, Clustering in Yrast States of 20Ne studied with antisymmetrized molecular dynamics, Prog. Theor. Phys. 93, 115 (1995)

[165]

M. Kimura and H. Horiuchi, Coexistence of cluster structure and superdeformation in 44Ti, Nucl. Phys. A 767, 58 (2006)

Rights & permissions

Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2020

PDF (57452KB)

2426

Accesses

0

Citation

Detail

Sections
Recommended

/

〈 〉