時間反演不對稱投影后變分方法

圖 雅, 何 艷, 高早春, 陳永壽, 田永威

(1. 沈陽師范大學 物理科學與技術學院, 沈陽 110034; 2. 中國原子能科學研究院, 北京 102413)

?

時間反演不對稱投影后變分方法

圖 雅1, 何 艷1, 高早春2, 陳永壽2, 田永威1

(1. 沈陽師范大學 物理科學與技術學院, 沈陽 110034; 2. 中國原子能科學研究院, 北京 102413)

在HF(Hartree-Fock)平均場的基礎上,建立了新的投影后變分方法(Variation after projection即VAP方法)。由于算法中考慮了HF平均場的時間反演對稱性破缺,使得新的VAP方法能夠很自然地推廣到偶偶核的奇自旋態,以及奇奇核和奇-A核的所有暈態計算中。采用這種新的VAP方法,初步計算了sd-殼原子核24Mg,25Mg,26Mg和26Al的暈帶,并將其與殼模型(Shell Model,SM)給出的準確值進行比較,結果發現二者非常接近。這表明VAP方法是一種非常有效的殼模型近似方法。同時,根據計算結果,分析認為角動量投影在獲得良好的殼模型近似中起到了關鍵的作用。

超越平均場; 投影后變分; 原子核結構

0 引 言

在原子核理論中,原則上,量子多體體系的波函數可以通過求解薛定諤方程而獲得。傳統殼模型[1]就是基于這種思想,但其在對重核的計算中遇到了很大困難。對于重核,其組態空間異常巨大,即使當今最快的超級計算機也難以實現重核的殼模型計算,這直接阻礙了人們對于重核結構性質的深入研究。為了回避這一困難,人們引入比較簡單的平均場近似方法[2],以揭示原子核的粗塊性質。其中比較成功且應用廣泛的平均場方法有HF (Hartree-Fock)[3]和HFB(Hartree-Fock-Bogoliubov)[3]方法。但是平均場方法破壞了體系的對稱性,同時很多物理效應被忽略,因此無法給出好的波函數。發展超越平均場方法,不僅可以大大改善平均場波函數,而且為重核的殼模型近似計算提供了一條非常有效的途徑。

超越平均場方法[4-7]已經成為研究量子多體體系的一類重要理論方法。其基本思想是,通過投影方法將平均場的對稱性進行恢復,得到的波函數具有好的量子數,進而能夠更好地描述量子多體體系的結構性質。為了盡可能好地優化波函數,理想的辦法是將投影波函數進行變分,即投影后變分(Variation After Projection,VAP)。到目前為止,在核理論中,VAMPIR[8-10]是比較成功的一種VAP方法。在VAMPIR計算中,同時考慮了粒子數投影和角動量投影,因此這需要很大的計算量。為了減少計算量,在通常的VAMPIR計算中,采用的平均場具有時間反演對稱性和軸對稱性。

在前人工作的基礎上,我們發展了自己的VAP[11]方法。在文獻[12]中,首次在VAP方法中引入三軸形變,并且首次實現了角動量(J),同位旋(T),以及質量數(A)的同時投影。計算結果表明,角動量投影是獲得很好殼模型近似的關鍵,同時也引起了三軸形變的普遍出現。在該工作中,仍然保留了時間反演對稱性。這就使得VAP的計算僅限于偶偶核的偶自旋態。為了將VAP的計算推廣到偶偶核的奇自旋態,奇-A 核和奇奇核,必須破壞平均場的時間反演稱性。

目前,在VAP中進一步考慮時間反演對稱性破缺。同時為簡化計算,將HFB平均場變為HF平均場。因為HF平均場中粒子數是確定的,因此不必進行粒子數投影。這不僅大大地優化了計算,同時也能更清楚的揭示純粹角動量投影帶來的新現象。采用新的VAP方法,計算了sd-殼一些核的暈態能量。通過比較發現,VAP能量值與殼模型給出的準確值十分接近。

1 VAP理論框架

(1)

式中:N是歸一化常數;而d是一個反對稱矩陣,其矩陣元dμν是決定|Φ〉的變分參量。可以通過改變dμν矩陣元來改變波函數|Φ〉。在這里設定dμν是復數,形式如下(其中,xμν和 yμν均是實數。)

(2)

(3)

(4)

為簡單起見,記

(5)

(6)

其相對應的投影能量EJ則可以展開為如下的形式

(7)

事實上,投影能量EJ及其相關的系數fk,可以通過求解Hill-Wheeler(HW)[13]方程同時獲得,

(8)

同時在這里,fk滿足歸一化條件〈ΨJ,M|ΨJ,M〉=1,

(9)

由于dμν是復數,我們要對實部xμν和虛部yμν分別求一階偏導。對于實部xμν,EJ的一階偏導數為

(10)

(11)

對于虛部yμν處理方法與上式相同,只需將xμν換成yμν。

(12)

求二階偏導數的優點在于:

1) 加快迭代的速度,提高計算效率。

2 計算結果與討論

在VAP計算中,采用的是USDB相互作用[14-15]。特別值得注意的是,在這里只用了一個HFSlater行列式進行計算,并且在此基礎上只做了角動量(J)投影。這樣不僅能夠優化算法,減少計算量,大大提高了計算的效率。而且,還能夠充分地表明,角動量投影是VAP近似方法取得成功的關鍵。換言之,角動量投影是獲得良好的殼模型近似的關鍵。

圖1 初期VAP計算結果與SM結果對比

圖2 后期VAP計算結果與SM結果對比

3 結 語

[ 1 ]CAURIER E,MARTINEZ P G,NOWACKI F,et al. The shell model as a unified view of nuclear structure[J]. Rev Mod Phys, 2005,77(2):427-431.

[ 2 ]HOVE D,GARRIDO E,JENSEN A S,et al. Combining Few-Body Cluster Structures with Many-Body Mean-Field Methods[J]. Few-Body Systems, 2017,58(2):33-38.

[ 3 ]RING P,SCHUCK P. The Nuclear Many-Body Problem[M]. New York,Heidelberg, Berlin: Springer Verlag, 1980:189-197,244-248,615-618.

[ 4 ]ROBLEDO L M,BERTSCH G F. Pairing in finite systeme: Beyond the HFB theory[J]. Fifty Years of Nuclear BCS: Pairing in Finite Systems, 2013:89-99.

[ 5 ]HAIDER Q,GOGNY D. Microscopic approach to the generator coordinate method with pairing correlations and density-dependent forces[J]. Phys G: Nucl Part Phys, 1992,18(6):993-1022.

[ 6 ]ENAMI K,TANABE K,YOSHINAGA N. Microscopic description of high-spin states: Quantum-number projections of the cranked Hartree-Fock-Bagoliubov self-consistent,solution[J]. Phys Rev C, 1999,59(1):135-153.

[ 7 ]RODRIGUEZ-GUZMAN R R,SCHMID K W. Spherical Hartree-Fock calculations with linear-momentum projection before the variation[J]. Eur Phys Jour A, 2004,19(1):61-75.

[ 8 ]SCHMID K W. On the use of general symmetry-projected Hartree-Fock-Bogoliubov configurations in variational approaches to the nuclear many-body problem[J]. Prog Part Nucl Phys, 2004,52(2):565-633.

[ 9 ]HJELT T,SCHMID K W,HAMMAREN E,et al. Unrestricted symmetry-projected Hartree-Fock-Bogoliubov calculations for some PF-shell nuclei[J]. Eur Phys J, 2000,7(2):201-208.

[10]PETORVICI A,SCHMID K W,FAESSLE A. Variational approach to shape coexistence in 68Se[J]. Nuclear Physics A, 2002,710(3/4):246-260.

[11]高早春,HOROI M,陳永壽. 三軸形變原子核平均場的投影后變分計算[J]. 中國原子能科學研究院年報, 2015:68-69.

[12]GAO Zaochun,HOROI M,CHEN Yongshou. Variation after projection with a triaxially deformed nuclear mean field[J]. Phys Rev C, 2015,92(6):064310.

[13]HILL D L,WHEELER J A. Nuclear Constitution and the Interpretation of Fission Phenomena[J]. Phys Rev, 1953,89(5):1102-1145.

[14]BROWN B A,RICHTER W A. New “USD” Hamiltonians for the sd shell[J]. Phys Rev C, 2006,74(3):34315-34315.

[15]RICHTER W A,BROWN B A.26Mg Observables for The “USDA” and “USDB” Hamiltonians[J]. Phys Rev C, 2009,80(3):178-190.

The VAP calculations with time-odd Hartree-Fock mean field in nuclear many body theory

TU Ya1, HE Yan1, GAO Zaochun2, CHEN Yongshou2, TIAN Yongwei1

(1. College of Physics Science and Technology, Shenyang Normal University, Shenyang 110034, China; 2. China Institute of Atomic Energy, Beijing 102413, China)

We have implemented a new variation after projection (VAP) calculation based on the Hartree-Fock (HF) mean field. Due to the time reversal symmetry breaking of the mean field,the present VAP can be naturally extended to the yrast states in odd-A and odd-odd nuclei as well as the yrast states in the even-even nuclei. With the new VAP code,we calculated the yrast states in24Mg,25Mg,26Mg and26Al. All the calculated VAP energies are very close to the Shell Model(SM) ones. This clearly shows that the spin projection is very important in achieving a good approximation to the full shell model.

Beyond mean field method; Variation after projection; Nuclei structure

1673-5862(2017)02-0156-05

2016-03-07。

國家自然科學基金青年科學基金資助項目(11305108); 國家自然科學基金理論物理專款項目(11047171)。

圖 雅(1979-),女,遼寧沈陽人,沈陽師范大學副教授,博士。

O571.21

A

10.3969/ j.issn.1673-5862.2017.02.006