Role of unitary correlation operator on high-momentum antisymmetrized molecular dynamics using bare NN interaction for 3H and 4He

We extend the high-momentum antisymmetrized molecular dynamics (HMAMD) by incorporating the short-range part of the unitary correlation operator method (UCOM) as the variational method of finite nuclei. In this HMAMD+UCOM calculation of light nuclei, the HMAMD is mainly in charge of the tensor correlation with up to the four-body correlation, while the short-range correlation is further improved by using the UCOM. The binding energies of the 3H and 4He nuclei are calculated with this HMAMD+UCOM using the AV8' bare nucleon-nucleon (NN) interaction. The different roles of the short-range and tensor correlations from the HMAMD and UCOM are analyzed in the numerical results. Compared with the previous calculations based on the different variational methods, this newly extended HMAMD+UCOM method can almost provide the consistent results with the ab initio results.


I. INTRODUCTION
In the recent decade, the existence of the strong short-range repulsion and tensor interaction in the bare nucleon-nucleon (NN) interaction is confirmed to provide strong short-range and tensor correlations by the ab initio studies of many-body physics in nuclei [1][2][3]. The short-range correlation contributes dominantly to the very high-momentum region around 4 fm −1 , while the tensor correlation is strong in the midium-momentum region around 2 fm −1 [4]. Therefore, the description of the NN correlations including high-momentum component becomes more and more essential for a developed framework to perform the ab initio calculation with the bare NN interaction.
When reproducing the nuclear structures and the binding energies of light mass region, one of the most successful calculations is the Green's function Monte-Carlo method (GFMC), where various correlation functions with many variational parameters are used for the nuclear wave function without any renormalization [5,6]. Besides, in recent years, the frameworks of tensor-optimized shell model (TOSM) [7][8][9] and tensor-optimized antisymmetrized molecular dynamics (TOAMD) [10][11][12][13][14][15] are developed successfully to reproduce the binding energies of s-shell nuclei using the AV8' bare interaction. In TOSM, two-particle two-hole excitations induced by the NN interaction are fully included to get the converging energy. In TOAMD, the correlation functions are introduced for the short-range correlation and tensor correlation and their functional forms are determined to minimize the total energy of the nucleus. In both methods, the high-momentum motion of nucleons in nuclei are explicitly treated. The description of the short-range correlation and tensor correlation basing on the microscopic framework remains an essential topic for the understanding of the inter-nucleon correlations.
In addition to the expansion scheme with the correlation functions for nuclei, another framework named the high-momentum antisymmetrized molecular dynamics (HMAMD) [16][17][18] has been developed for the description of the NN correlation and this method is an extension of the pioneering idea given in Ref. [19]. In HMAMD, the correlated nucleon-pair with high-momentum component is introduced in the nuclear many-body states and this treatment is extendable by increasing the number of correlated pairs at the same time. It is found that the short-range correlation and the tensor correlation are well described within a rather simple form of high-momentum NN pair excitation. Another unique treatment for the correlations between nucleons is the unitary transformation of the uncorrelated state to the correlated many-body states as done by the unitary correlation operator method (UCOM) [9,20,21]. In UCOM, the short-range correlation and the tensor correlation are described with two kinds of unitary correlation operators separately, which have a good property for porting them to the various kinds of nuclear many-body wave functions. In particular, short-range part of UCOM can work nicely to describe the short-range correlation in finite nuclei [9,20,22] and also in nuclear matter [23,24] Thus, it is rather straightforward to combine the HMAMD with UCOM to further improve the description of the inter-nucleon correlations. Under the unitary transformation, the HMAMD still maitains its analytical simplicity while providing a better description of manybody correlations. In our previous work [9], we have incorporated the UCOM as well as the S-wave UCOM (S-UCOM) into the TOSM for the calculation of the 4 He nucleus and obtained better numerical results which are very close to the rigorous calculation. In this work, we will combine the HMAMD with UCOM, namely HMAMD+UCOM, and investigate the applicability of this new variational method by performing the calculations of 3 H and 4 He nuclei using AV8' bare interaction. Since the tensor correlation has already been well described within the framework of HMAMD [4,16], we will adopt only the short-range part of UCOM in the current work. The detailed role of UCOM will be discussed from the numerical results of this new extension. We will also compare these new results with those from previous TOSM+UCOM/S-UCOM calculations as well as the TOAMD calculations. This paper is organized as follows: Section II explains the framework of the HMAMD wave function and the UCOM transformation, which we are aiming to propose. The numerical results of the HMAMD+UCOM calculations for the 3 H and 4 He nuclei are presented and discussed in Sec. III, where the comparison of the current HMAMD+UCOM between the TOSM+UCOM/S-UCOM or TOAMD will also be included. The conclusion is given in Sec.

II. FRAMEWORK
We shall begin with the many-body Hamiltonian, whereT i andT c.m. are the kinetic energies of each nucleon and the center of mass, respectively. In this work, we take the bare nucleon-nucleon interaction AV8' [6] forV ij aŝ where v C ij denotes the central term, v T ij denotes the tensor term and v LS ij denotes the spinorbit interaction. In AV8', the central term v C ij includes the short-range repulsive interaction, which describes the short-range correlation when two nucleons are at short distance. Additionally, the tensor interaction v T ij is included in the AV8' comparing with the AV4' central interaction. The Coulomb term is not considered in the current work with the same condition as the other works [9][10][11][12][13][14]25] to compare the results with the AV8' interaction.
The wave function we adopt in this work is the HMAMD wave function, which is based on the antisymmetrized molecular dynamics (AMD). For a nuclear system with A nucleons, the AMD basis can be written as where the single-nucleon wave functions φ(r) are expressed in the Gaussian wave packet multiplied by the spin-isospin wave function χ τ,σ as Here, the range parameter ν represents the size of the nucleon and is common for all nucleons in Eq. (3). It will be determined for each nucleus independently to gain the total energy.
The centroid position Z is a complex vector Z = R + iD, where the real part R represents the spatial position of the Gaussian centroid. We put the condition of A i Z i = 0 to keep the center of mass at the origin of the coordinates and the total momentum at zero. For the s-shell nuclei 3 H and 4 He, the spatial positions of all nucleons are optimized to be R = 0 as the result of the energy variation in TOAMD [13,14]. The imaginary part of Z, which is D, represents the mean momentum of the Gaussian as p = 2 νIm(Z) .
In the framework of HMAMD, the imaginary components to the Gaussian centroids of two nucleons Z 1 and Z 2 are assigned with the same magnitude but opposite directions as In this way, the high-momentum excitation of di-nucleon pairs is introduced into AMD basis with a large magnitude of D [19] and also various vector directions [16], both of which are necessary to describe the tensor correlation without the center-of-mass excitation. The corresponding AMD basis is extended to the single HMAMD basis Φ 1HM , since only a single di-nucleon pair is excited. To extend to the double HMAMD basis Φ 2HM [17], the additional basis with double high-momentum pairs should be included as and is required for the precise description of NN correlations in the order of cluster expansion [25].
The spin-isospin parts of the wave function are generally fixed to be |p ↑ n ↑ n ↓ in the 3 H case and |p ↑ p ↓ n ↑ n ↓ in the 4 He case. However, the spin of the nucleon could be changed due to the tensor operator ofŜ 12 = 3(σ 1 · r)(σ 2 · r) − (σ 1 · σ 2 ) in the tensor force in the NN interaction. It indicates that the spin parallel configuration of two nucleons can contribute to the total wave function satisfying the Pauli exclusion principle. Therefore, we will add the extra wave functions having the spin-parallel configurations, which is technically treated by using the free spin expression of χ σ = η |↑ + (1 − η) |↓ . Here the parameter η is a constant between 0 to 1. The value of η will not influence the final result because of the antisymmetrization.
The final wave function of HMAMD is the superposition of all the basis with grid of the momentum shift D: whereP J π M K is the parity and the angular momentum projector. The corresponding coefficients f i and g j,k will be determined by the diagonalization of the Hamiltonian. The HMAMD method is a powerful method for the description of the short-range correlation and the tensor correlation, especially for the latter one. On the other hand, another way to implement the correlations into the system is to perform a suitable unitary transformation of the uncorrelated many-body states to a correlated basis, Ψ = CΨ. Here the tilde notation represents the correlated wave function. By following the works of Feldmeier et al. [20,21], an unitary correlation operator can be written as a product of two unitary operators: where C Ω is the unitary operator for the tensor correlation and C r is the unitary operator for the short-range correlation. Two-body operators g Ω,ij and g r,ij are the Hermite operator where i and j are the particle index. Since the tensor correlation can be well described by the HMAMD basis as studied in Refs. [16], we only adopt the short-range correlation part C = C r in the current work.
In the calculation of UCOM, the Hermite operator g r,ij is defined as where the momentum p r,ij is the radial component of the relative momentum, which is conjugate to the relative coordinate r ij = |r i − r j |. s(r ij ) is the amount of the shift of the relative wave function of each nucleon pair. In the framework of UCOM, the function R + (r) can be introduced as and one can obtain the relation: The function R + (r) denotes the transformed relative distance and can reduce the amplitude of the short-range part of the relative wave function in nuclei, which should be determined for four spin-isospin channels (singlet-even, triplet-even, singlet-odd and triplet-odd) independently. Thus, the various operators will be transformed to the form of the relative motion as: Here we list the transformation of the relative position operator r, the relative momentum operator p r , the relative orbital angular momentum operator l, the intrinsic spin operator s, the tensor operator S 12 , and the arbitrary scalar function v(r).
The form of R + (r) depends on the potential adopted in the calculation. For the AV8' potential, the following two parameterizations for the short-range correlation functions, which are proposed by Neff-Feldmeier and Roth et al. [20,21,26,27], have proven appropriate: Here, α, β, γ are the variational parameters to minimize the energy of the total system, which will be determined for four channels of the spin-isospin pair separately.
With the above unitary correlation operator C, one can transform the HMAMD wave function to the correlated one as Ψ HMAMD = CΨ HMAMD , and the corresponding Schrodinger equation becomesĤ Ψ HMAMD = E Ψ HMAMD or Ĥ Ψ HMAMD = EΨ HMAMD . Hence, the transformation to the effective Hamiltonian is In addition to the correlated nucleon-nucleon interaction term, t ij , t P ij , and t L ij are the twobody kinetic terms, which come from the correlated kinetic operator with the cluster expansion by the C †T [1] C =T [1] +T [2] +T [3] + .... The number in the square brackets denotes the one-body term, two-body term and so on. For the simplicity of the calculation, we take the approximation to use only up to the two-body operators in the correlated Hamiltonian to treat the short-range correlation, which has been shown to work nicely [9,12]. In this calculation, we apply the Gaussian expansion technique for the radial part of the effective Hamiltonian, so that the forms of these three kinetic operators are t ij = n c n e anr 2 ij t P ij = n c P n (P 2 e a P n r 2 ij + e a P n r 2 where t ij is a central operator type. t P ij is denoted by using the relative momentumP and t L ij is denoted by using the relative orbital angular momentumL. The index n distinguishes the Gaussian range. We also consider the spin-isospin dependence of three operators in Eq. (17).
In the procedure of the current work, we transform the Hamiltonian to the effective one by the UCOM and then solve the eigenvalue problem of the effective Hamiltonian. We will make the variational calculation for the parameters of the R + (r) to minimize the total energy results.

III. RESULTS
In this work, we combine the HMAMD with the UCOM to calculate 3 H and 4 He nuclei using the bare interaction AV8'. Firstly, as a comparison, we show the energy results of 3 H and 4 He nuclei only with the HMAMD in Table. I. Since for the HMAMD+UCOM, the optimized size parameters ν in the single-particle wave function are 0.18 fm −2 and 0.25 fm −2 for 3 H and 4 He, respectively, we adopt the same parameters for the HMAMD calculations without UCOM. In this table, we see clearly the efficiency of the single and   These facts show that the UCOM can surely provide a better description of the short-range correlation in the framework of HMAMD. When comparing with the results shown in Table. I, we notice that the single HMAMD result of 4 He provides larger energy gain when adopting the UCOM than in the 3 H case. It indicates that the UCOM works more effective in the calculation of 4 He than in the calculation of 3 H. This difference can be explained from the property of the short-range correlation. Since the short-range correlation describes the strong repulsion when the nucleons are close to each other, it is natural to understand that the compact nuclei such as 4 He will be more affected by the short-range correlation. This should be the reason why the UCOM for the short-range correlation is more effective in a compact nucleus as 4 He than in a less compact nucleus as 3 H.
Next, we intend to show how the UCOM effect differs between the single and double HMAMD basis in describing the short-range correlation. The form of R + (r) for 3 H calculation are plotted for the even channels in Fig. 2. From the triplet-even (TE) panel of Fig. 2, it is shown that the strength of the TE channel is small in the calculation comparing with other channels. A reasonable explanation is that the tensor force can also be transformed in the TE central part of the Hamiltonian in the UCOM. However, since the strong tensor effect is already implemented by the tensor part of the AV8' interaction, the strength of the TE channel in UCOM should be reduced to survive the bare tensor force as much as possible. When we compare the UCOM effect with the single and the double HMAMD calculations in the SE and TE channels, we can find that the UCOM effect will be reduced when adding the double HMAMD basis. These results again indicate the well description of the short-range correlation by the double HMAMD basis, so that the UCOM effect becomes smaller than the case of single HMAMD. We also show the energy contributions from the even and odd channels of the UCOM for 3 H in Fig. 3. In this figure, when introducing only the even channels of the UCOM, the energy will be rapidly reduced while the odd channels contribute a little to the final result. This result shows that the role of the even channel is important rather than the odd channel in the short-range correlation as the discussion reported in Ref. [25].
Next, we add the spin parallel configurations into the calculated model space. This is because the S 12 tensor operator will change the upward spin parallel states to the downward spin parallel states as |p ↑ n ↑ → |p ↓ n ↓ , so that the spin parallel configuration of the two neutrons or protons should also be considered. In Table. V, we show the comparison of the results of 3 H and 4 He with the spin parallel or only the natural spin configuration.
These calculations are performed with the same parameters for each nucleus and include the double HMAMD basis. From this table, we can observe a slight improvement in the total energy when adding the spin parallel configuration. These results show that the spin parallel configuration is surely required by the bare NN interaction. Furthermore, we can see that the improvement is mainly coming from the tensor interaction, which means that  the spin parallel configuration is favored by the tensor correlation as we discussed before.
Even though the contribution from the spin parallel configuration is not so large, it can not be ignored especially for the benchmark calculation.
Finally we compare the present results with the TOSM+UCOM/S-UCOM work [9] and the TOAMD works [13]. Here, the S-UCOM is the extension of UCOM by introducing the partial-wave dependence in UCOM, where S-wave is only transformed in the even channel.
In S-UCOM, the tensor contribution increases from the ordinary UCOM because of the enhancement of the SD coupling by the tensor force. In Table. VI, we list our current best results for 3 H and 4 He as well as the corresponding TOSM+UCOM/S-UCOM and TOAMD results. The NN interaction adopted in these reference works is also the AV8' interaction. The TOAMD method is basing on the AMD basis as is used in our current work, while the shell model wave functions are used in the work with TOSM with the full optimization of the two-particle two-hole excitations. It is also noted that all methods are variational approach using the bare NN interaction directly. From this comparison, we can is explained in Ref. [9]. Thus, it is reasonable to estimate that HMAMD+S-UCOM instead of the UCOM will get lower energy than the present results. The progress may benefit the description of the tensor correlation as it has been done on the TOSM. And when noticing that the TOAMD gains more energy in the tensor part of the Hamiltonian than the HMAMD+UCOM, the S-UCOM could be the solution for the current HMAMD+UCOM work to reproduce or improve the results in comparison with the TOAMD framework. In S-UCOM, we need the partial-wave expansion of the single-particle wave function in AMD [30], which is considered in the future work. The other possible improvement is the addition of the third pair of nucleons having high-momentum in the HMAMD basis states called the triple HMAMD, which is a natural extension of the present double pair case in Eq. (8).

IV. CONCLUSION
We combine the unitary correlation operator method (UCOM) for short-range part with the double high-momentum antisymmetrized molecular dynamics (HMAMD) as a new variational method for finite nuclei using the bare NN interaction. In HMAMD, up to the two high-momentum pairs excitation is included in the basis states. The UCOM is in charge of the short-range correlation while the HMAMD is expected to be in charge of the tensor correlation. By adopting this new HMAMD+UCOM framework, we calculate the binding energies of the 3 H and 4 He nuclei with AV8' interaction, in which the short-range and the tensor correlation are included.
From the results getting from this work, we can conclude that the UCOM has certainly contributed to the energy gain in the final results. We also find that the UCOM has different efficiency in 3 H and 4 He cases. It is because the short-range correlation described by the UCOM becomes large in the more compact nuclei. By comparing the calculation results from the single HMAMD and double HMAMD with the UCOM, we find that the double high-momentum pair is essential for the description of both the short-range correlation and the tensor correlation. However, even with the double HMAMD basis, the UCOM is still necessary to improve the description of the short-range correlation. When considering the effect of the tensor operator in the NN interaction, we notice that the spin parallel configuration should contribute to the final wave function. Comparing with the present results from the TOSM+UCOM/S-UCOM and TOAMD, the results getting from this work still have the possibility to be further improved by the S-UCOM, so that has the potential to perform the ab initio calculation in future works.