arxiv: v1 [nucl-th] 15 Jan 2015

Transcription

1 Accuracy of the new pairing theory and its improvement L. Y. Jia, Department of Physics, University of Shanghai for Science and Technology, Shanghai 93, P. R. China Department of Physics, Hebei Normal University, Shijiazhuang, Hebei 54, P. R. China (Dated: September 6, 8) Recently I proposed[] a new method for solving the pairing Hamiltonian with the pair-condensate wavefunction ansatz based on the Heisenberg equations of motion for the density matrix operators. In this work an improved version is given by deriving the relevant equations more carefully. I evaluate both versions in a large ensemble with random interactions, and the accuracy of the methods is given statistically in terms of root-mean-square derivations from the results. The widely used variational calculation is also done and the results and computing-time costs are compared. arxiv:5.3684v [nucl-th] 5 Jan 5 PACS numbers:.6.ev,..re, Keywords: I. INTRODUCTION Pairing correlation has long been recognized in nuclei [] and influences all the properties of the latter, such as mass, gap of excitation energy, and moment of inertia [3]. In general, any mean-field treatment of the nuclei needs to account it somehow to get reasonable results. Among the methods [4] the most popular one may be the Bardeen-Cooper-Schrieffer (BCS) theory [5], or its advanced version the Hartree-Fock-Bogoliubov(HBF) theory [3], where pairing correlation is considered by introducing quasi-particles and writing the ground state as a vacuum of the latter. But there are disadvantages of breaking the particle number and need for a unphysical minimum pairing strength [6, 7]. A common improvement is to use the pair-condensate [Eq. (), the BCS wavefunction projected onto good particle number] as the ground state wavefunction [8]. Usually the criteria to determine the variational parameters is minimizing the energy in the variation principle (variation after projection BCS) [8 3]. In Ref. [] I proposed a new criteria based on the Heisenberg equations of motion (EOM) for density matrix operators, as the lowest-order (mean-field) results of the generalized density matrix (GDM) formalism [4 ]. The method was applied to the calcium isotopes with the well known FPD6 interaction []. In this work I give an improved version of the method. In Ref. [] the relevant equations are derived assuming that neighboring even-even nuclei have the same Hartree- Fock (HF) single-particle energies and occupation numbers; here I do it more carefully allowing the latters to be different. To see the validity of the approaches, both versions are applied to a large ensemble with examples and random interactions; the results are good in almost all the examples. I also compare the results and necessary computing-time costs with those of the variation principle. In Sec. II I derive the improved version of the GDM pairing theory. Then in Sec. III the approaches are evaluated in the large random ensemble. Finally Sec. IV summarizes the work and discusses further directions. II. FORMALISM The pairing theory in Ref. [] was derived as the lowest-order (mean-field) results of the complete GDM formalism. In fact the derivation could be simpler and more clear if we focus on the mean fields and do not introduce collective (quadrupole) phonons. Also, in Ref. [] I made the assumption that neighboring even-even nuclei have the same HF single-particle energies and occupation numbers. Below I derive an improved version of the pairing theory abandoning the above assumption. As before, the ground state of the N-particle system is assumed to be an N-pair condensate, φ N = χn (P ) N, () where χ N is the normalization factor, and P is the pair creation operator P = v a. () a In Eq. () the summation runs over the entire singleparticle space. is the time-reversed level of the singleparticle level. The pair structure v are parameters to be determined by the theory. With the antisymmetrized fermionic Hamiltonian H = ǫ a a + V 34 a 4 a a 3a 4, (3) 34 I calculate the Heisenberg equations of motion for the density matrix operators R a a and K a a, [R,H] = [ǫ,r] V 543 a 5 a 4 a 3a + V 345 a a 3 a 4a 5, (4) [K,H] = (ǫk) (ǫk) + V 34 K V 543 a 5 a 4a 3 a + V 543 a 5 a 4a 3 a. (5)

2 Terms on the right-hand sides of Eqs. (4) and (5) are read as matrix multiplications; for example, [ǫ,r] = (ǫr) (Rǫ) = 3 ǫ 3R 3 3 R 3ǫ 3. On the pair condensate (), the density matrices are diagonal : ρ N φ N a a φ N = δ n N, (6) κ N φ N a a φ N = δ sn. (7) In practical shell-model calculations usually each singleparticle level has distinct spin and parity, thus the mean fields are also diagonal : f N ǫ + V 43 ρ N 34 = δ e N, 34 (8) δ N V 34 κ N 43 gn. (9) 34 The pairing mean field δ N should not be mixed with the Kronecker delta δ. Now I take matrix elements of the equations of motion between the pairing ground states (): equation (4) between those of the same nuclei ( φ N and φ N ); equation (5) between those of neighboring even-even nuclei ( φ N and φ N ). On the left-hand side of Eq. (5) we have φ N [K,H] φ N = φ N K H φ N φ N HK φ N (E N E N )κ N, where E N and E N are the ground state energies, H φ N E N φ N and H φ N E N φ N. Similarly for Eq. (4) we have φ N [R,H] φ N (E N E N )ρ N =. On the right-hand sides the two-body density matrices are approximated in the following way: φ N a 4 a 3 a a φ N φ N a 4 a φ N φ N a 3 a φ N φ N a 4 a φ N φ N a 3 a φ N + φ N a 4 a 3 φ N φ N a a φ N, = ρ N 4ρ N 3 ρ N 4ρ N 3 +κ N 34 κn, () φ N a 4 a 3a a φ N φ N a 4 a φ N φ N a 3 a φ N φ N a 4 a φ N φ N a 3 a φ N + φ N a 4 a 3 φ N φ N a a φ N = ρ N 4 κ N 3 ρ N 4 κ N 3 +ρ N 34 κ N. () Equations () and () would be if the ground states were taken as single-particle Slater determinants: the right-hand sides were just the fully contracted terms in Wick s theorem. Here they are approximations because the ground states are taken as pair condensates (). Finally Eqs. (4) and (5) become = [f N,ρ N ] κ N δ N +δ Nκ N, () (E N E N )κ N = f N κ N +κ N f T N +δ N δ N ρ T N ρ N δ N. (3) Under the diagonal properties(6), (7), (8), and(9), Eq. () is satisfied automatically, and Eq. (3) becomes E N E N = e N +g N n N. (4) Equation (4) is the main equation of the improved theory. It implies that the right-hand side is independent of the single-particle label, by which the parameters v in Eq. () are fixed. The main equation of the old theory, Eq. (9) in Ref. [], corresponds to replacing e N and n N by e N and n N. Equation (4) includes the wellknown particle-particle random phase approximation [7] as its special case of N = (n =, e = ǫ). The normalization factor χ N (), occupation numbers n N (6), pair-transfer amplitudes s N (7), mean fields en (8) and g N (9) are functions of the pair structures v (); their functional forms, as the kinematics of the system, have already been given in Eqs. (3) and (4) of Ref. [] and are not repeated here. In the next section I take the pairing Hamiltonian: ǫ = δ ǫ, V 34 = δ δ 3 4 G 3 (5) in Eq. (3). Consequently the mean fields (8) and (9) become e N = ǫ G n N, gn = G s N. (6) I apply both versions of the theory to a large random ensemble; for convenience I call Eq. (4) the GDM pairing theory, and Eq. (9) in Ref. [] the GDM pairing theory. s N III. RANDOM ENSEMBLE The GDM theorywasapplied in Ref. [] to calcium isotopes with the FPD6 interaction and the results are good. Here I would like to consider both the GDM and GDM theories in a large ensemble with random interactions; consequently we can speak statistically the accuracy of the theories in terms of the root-mean-square derivations from the results. The variational calculation with the trial wavefunction () is also performed and the accuracy and computing-time cost are compared. The random ensemble has examples with different parameters determined in the following way. For each example, I first pick up the single-particle levels randomly from the pool j =,3,5,7,9,,3. Each angular momentum j has a 4% probability of being selected; the selected ones (at least two) constitute the model space. Second, the single-particle energies ǫ (5) are determined randomly following the uniform distribution from MeV to. Third, I pick up the pairing strength G max as a random number following the uniform distribution from to MeV, then the pairing matrix elements G (5) are distributed uniformly from

3 3 to G max. Finally, the number of pairs N is determined following the uniform distribution from to Ω, where Ω = j (j+) is the maximal particle number allowed by the model space. I perform four sets of calculations for the ensemble: two GDM calculations GDM and GDM, variational calculation VAR, and the calculation. (In this work the calculation is done by diagonalization in spaces with fixed seniority [, 3]. It can also be achieved by the Monte Carlo algorithm [4 6]; or the Richardson s method in some special cases [7 9].) For the variational calculation of the pair-transfer amplitudes s = φ N a a φ N (7), I show two sets of results. In VAR the pair structure v () in φ N and φ N are the same, given by minimizing φ N H φ N ; while in VAR v in φ N and φ N are different, given by minimizing φ N H φ N and φ N H φ N, respectively. In Figs. and I show the complete spectroscopic results for the ensemble. For example, in panel of Fig. there are 354 points (crosses), corresponding to the 354 single-particle levels in the examples of the ensemble. The horizontal coordinate of each point is the value of n of the corresponding single-particle level, while the vertical coordinate is the GDM value. Thus a perfect calculation would have all the points lying onthe y = xstraightline. Similarly, in panels and of Fig. the vertical coordinates are the GDM and VAR values of n, respectively. Figure is plotted in the same way for the pair-transfer amplitudes s. From Figs. and we see that the variational calculation for n and the VAR version of s are generally better than the GDM ones. However, the less-careful VAR calculation of s is worse than the GDM ones. The rootmean-square (σ) derivations from the results are σ GDM n on Fig. and =., σ GDM n σs GDM σs VAR =.88, σ VAR n =.45,(7) =.7, σs GDM =.9, =.373, σs VAR =.9, (8) on Fig.. We note that on Fig. there is a point much worse than others for all the calculations, so let us look at the particular example it belongs to. This example has N = particles on two single-particle levels with angular momenta j = 5, 9 and energies ǫ5 = 8.988, ǫ9 = 9.39 MeV. The pairing two-body matrix elements are G5 =.45,,5 G9 =.456, and,9 G5 =,9 G9 =.6 MeV. The failed point corresponds to the j = 5,5 level. This example is particular in that ǫ5 ǫ9 G5,5 G9,9 G5, thus,9 there is little correlation between the two levels. I have looked at the wavefunction of the daughter nucleus and it is mainly P 5 (P 9 ) 4, which is not representable by Eq. (). Without this example, the root-meansquare derivations for the pair-transfer amplitudes s are σs GDM σ VAR =.77, σ GDM s =.63, σ VAR s =.35, and s =.79. The necessary formula for the GDM and variational calculations were given in Ref. []. In general, in large model spaces the GDM calculation costs less time than the variational one, by a factor of the number of nondegenerate single-particle levels (in order of magnitude), because the former needs to calculate only φ N P φ N while the latter needs φ N P P φ N (P a ). However, as we can see from Figs. and and Eqs. (7) and a (8), the accuracy of the GDM method is close to that of the variational one. The reduction of time cost may be a big advantage when doing ab-initio mean-field calculations for medium and heavy nuclei, especially if we were fitting parameters of the interaction (for example the effort in developing density functionals with spectroscopic accuracy). Next I would like to see the accuracy of the methods depending on different quantities. For this purpose I plot the root-mean-square derivations of n and s as functions of the single-particle angular momentum j, particle number N, and pairing strength G max in Figs. 3, 4, and 5, respectively. We have the following observations:. From Figs. 3 and 4 we see that as a trend the GDM results improve with increasing j and N, which should be expected because the GDM formalism is a collective theory.. There seems to be no obvious trend with the pairing strength shown in Fig On Fig. 4, the GDM calculation for N =, which is the particle-particle random phase approximation, is much better than the GDM calculation. And for small N (until N = 6), the GDM s seems to be slightly better than the GDM s. 4. On Fig. 3, the GDM s is slightly better than the GDM s at the smallest angular momenta j = and 3. For these levels the difference n N n N, which is inversely proportional to j around the Fermi surface, is largest; and the GDM theory with n N seems to be slightly better than the GDM one with n N. 5. The GDM n seems to be consistently better than the GDM n and the reason is still unclear. 6. In panel of Figs. 3, 4, and 5 the worst points are at j = 5, N = 6, and.5 MeV < G max <.6 MeV, respectively; because these three sub-groups contain the worst example mentioned below Eq. (8). At last in Fig. 6 I show the results for the ground state energy by different calculations. In panel there are points (crosses), corresponding to the examples of the ensemble. The horizontal coordinate of each point is the pairing correlation energy of the corresponding example, E pair ǫ n F E shell, where E shell is the ground state energy of the shell model calculation, and n F = or is the occupation number of the naive Fermi distribution. The vertical coordinate shows the variational ground state energy measured from the one, E var = φ var N H φvar N E shell, where φ var N is the pair condensate () with its pair structure v () determined by the variation principle. Similarly for panels and, but the vertical co-

4 4 ordinates are the ground state energies of the GDM and GDM calculations, respectively, measured from the one (E GDM = φ GDM N H φ GDM N E shell and E GDM = φ GDM N H φ GDM N E shell ). We see that all the three calculations give good ground state energies: the errors are small relative to the pairing correlation energy. The average values are Ēvar =.45, Ē GDM =.6, Ē GDM =.7, and Ēpair = 3.64 MeV. It is well known that the variational calculation finds the best and lowest ground state energy for a set of restricted wavefunctions of the form (). However, from the above average values we see that the energies of the GDM wavefunctions φ GDM N and φ GDM N are close to the variational minimum φ var N H φvar N. IV. SUMMARY In summary, I derive a physically improved version of the GDM pairing theory proposed in Ref. []. Both versions are checked in a large random ensemble, and the accuracy is given statistically in terms of root-meansquare derivations from the results. Consequently, we could consider the theories to be correct and apply them with confidence to realistic systems. Based on the results, the GDM theories are not as accurate as (although close to) the variation principle. However, the reduction of computing-time cost is huge for large model spaces (by a factor of the number of non-degenerate single-particle levels). This should be interesting for ab-initio mean-field calculations of medium and heavy nuclei (especially if deformed Nilsson singleparticle levels were used), or the effort in fitting parameters of an interaction. Comparingthe twoversionsofthe GDM theory, we see that in general the new one (physically more reasonable) is slightly better in calculating pair-transfer amplitudes s, while theoldoneproducesbetteroccupationnumbers n and slightly better ground state energy. The reason for the latter is still unclear. The key approximation of the current GDM methods is the factorization or linearization of the two-body density matrix on the pair condensate [Eqs. () and ()]. It would be interesting to see its validity in other circumstances, in particular, whether we could use it in the variational formalism when calculating the twobody part of the average energy, which would reduce the time cost to the same level of the current GDM methods. Support is acknowledged from the startup funding for new faculty member in University of Shanghai for Science and Technology. The numerical calculations of this work are done at the High Performance Computing Center of Michigan State University. [] L. Y. Jia, Phys. Rev. C 88, 4433 (3). [] A. Bohr, B. R. Mottelson, and D. Pines, Phys. Rev., 936 (958). [3] S. T. Belyaev, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 3, () (959). [4] Ricardo A Broglia, and Vladimir Zelevinsky, Fifty Years of Nuclear BCS: Pairing in Finite Systems (World Scientific, 3). [5] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 6, 6 (957); Phys. Rev. 8, 75 (957). [6] A. Bohr and B. Mottelson, Nuclear Structure (Benjamin, New York, 975). [7] P. Ring and P. Schuck, The Nuclear Many-Body Problem (Springer-Verlag, Berlin, 98). [8] K. Dietrich, H. J. Mang, and J. H. Pradal, Phys. Rev. 35, B (964). [9] L.M. Robledo, G.F. Bertsch, arxiv: [nucl-th] (). [] H. Flocard and N. Onishi, Annals of Physics 54, 75 (997). [] Fabian Braun and Jan von Delft, Phys. Rev. Lett. 8, 47 (998). [] J. Dukelsky and G. Sierra, Phys. Rev. B 6, 3 (). [3] Jan von Delft and D.C. Ralph, Physics Reports 345, 6 (). [4] A. Kerman and A. Klein, Phys. Rev. 3, 36 (963). [5] S.T. Belyaev and V.G. Zelevinsky, Yad. Fiz., 74 (97) [Sov. J. Nucl. Phys., 46 (97)]; Yad. Fiz. 6, 95 (97) [Sov. J. Nucl. Phys. 6, 657 (973)]; Yad. Fiz. 7, 55 (973) [Sov. J. Nucl. Phys. 7, 69 (973)]. [6] V.G. Zelevinsky, Prog. Theor. Phys. Suppl , 5 (983). [7] M.I. Shtokman, Yad. Fiz., 479 (975) [Sov. J. Nucl. Phys., 47 (976)]. [8] L. Y. Jia, Phys. Rev. C 84, 438 (). [9] L. Y.Jia, and V.G. Zelevinsky,Phys. Rev. C84, 643 (). [] L. Y.Jia, and V.G. Zelevinsky,Phys. Rev. C86, 435 (). [] W.A. Richter, M.G. Van Der Merwe, R.E. Julies, and B.A. Brown, Nucl. Phys. A53, 35 (99). [] H. Molique and J. Dudek, Phys. Rev. C 56, 795 (997). [3] V. Zelevinsky and A. Volya, Nucl. Phys. A75, 35 (5). [4] N. Cerf and O. Martin, Phys. Rev. C 47, 6 (993); [5] N. Cerf, Nucl. Phys. A564, 383 (993). [6] Abhishek Mukherjee, Y. Alhassid, and G. F. Bertsch, Phys. Rev. C 83, 439 (). [7] R. W. Richardson, Phys. Lett. 3, 77 (963); Phys. Rev. 4, 949 (966). [8] J. Dukelsky, S. Pittel, and G. Sierra, Rev. Mod. Phys. 76, 643 (4). [9] J. Dukelsky, S. Pittel, arxiv:4.95 [nucl-th] ().

5 5 occupation numbers n GDM GDM VAR FIG. : (Color online) Occupation numbers n of all the 354 single-particle levels in the ensemble by different calculations. In panel, the single-particle levels and the points are in one-to-one correspondence, with the horizontal coordinate of the point being the n and the vertical coordinate being the GDM n. Similarly for panels and, but the vertical coordinates are n of the GDM and the variational calculations, respectively. See text for details.

6 6 pair-transfer amplitudes s GDM.5 GDM (d) VAR.5 VAR FIG. : (Color online) Pair-transfer amplitudes s of all the 354 single-particle levels in the ensemble by different calculations. The points are plotted in the same way as those in Fig. (), but for the transfer amplitudes s. Panels and (d) plot two sets of variational calculations (see text for details).

7 7 single-particle-level groups with different j group size rms n j GDM GDM VAR rms s j j GDM GDM VAR VAR FIG. 3: (Color online) Results of different calculations grouped by the single-particle-level angular momentum j. I divide all the 354 single-particle levels in the ensemble into different groups according to their angular momentum j, and panel plots the group sizes. Panel plots the root-meansquare derivations from the results of the occupation numbers n by different calculations (GDM, GDM, VAR) within each j group. Similarly, panel plots the root-meansquare derivations of the pair-transfer amplitudes s by four sets of calculations (GDM, GDM, VAR, VAR). See text for details.

8 8 sub-ensemble size rms n 5 5 sub-ensembles with different N N GDM.3 GDM VAR.. rms s N GDM GDM VAR VAR N FIG. 4: (Color online) Results of different calculations in subensembles divided by the particle number N. The examples in the ensemble are divided into sub-ensembles according to their particle number N, and panel plots the sizes of the sub-ensembles. Panels and plot the root-meansquare derivations by different calculations of the occupation numbers n and pair-transfer amplitudes s within each N sub-ensemble. See text for details.

9 9 sub-ensembles with different pairing strength sub-ensemble size G max (MeV) rms n.3.. GDM GDM VAR.5.5 G max (MeV) rms s GDM GDM VAR VAR.5.5 G max (MeV) FIG. 5: (Color online) Results of different calculations in subensembles divided by the pairing strength G max. The examples in the ensemble are divided into sub-ensembles according to their pairing strength G max, and panel plots the sizes of the sub-ensembles. Panels and plot the root-mean-square derivations by different calculations of the occupationnumbersn andpair-transferamplitudess within each sub-ensemble. See text for details.

10 ground state energies Evar (MeV) E pair (MeV) EGDM (MeV) EGDM (MeV) E pair (MeV) E pair (MeV) FIG. 6: (Color online) Ground state energies of all the examples by different calculations. In panel, the examples and the points are in one-to-one correspondence, with the horizontal coordinate of the point being the pairing correlation energy E pair, and the vertical coordinate being the variational ground state energy E var measured from the ground state energy. Similarly for panels and, but the vertical coordinates are ground state energies by the GDM and GDM calculations, respectively, measured from the ground state energy. See text for details.