[논문 리뷰] A numerical method of solving time-dependent Hartree-Fock-Bogoliubov equation with Gogny interaction

Yukio Hashimoto, K. Nodeki|ArXiv.org|2007. 07. 20.

Cold Atom Physics and Bose-Einstein Condensates인용 수 6

한 줄 요약

이 논문은 Gogny 상호작용을 사용하여 시간에 따라 변하는 Hartree-Fock-Bogoliubov (TDHFB) 방정식을 수치적으로 해결하기 위한 방법을 제시한다. 시간 진동수 연산자의 멱급수 전개를 사용하여 운동 방정식을 적분하며, 20O에 적용한 결과, 입자 수와 에너지가 매우 정확하게 보존되며, 소규모 진폭의 2중극 진동수 분포를 성공적으로 재현하여, 이국적 핵에서의 집단 진동수를 연구하는 데의 가능성을 입증한다.

ABSTRACT

A numerical method to solve time-dependent Hartree-Fock-Bogoliubov equation is proposed to treat the case with Gogny effective interaction. To check the feasibility of the method, it is applied to oxygen isotope O20 and small amplitude oscillations are calculated. The conservation of the nucleon numbers as well as energy expectation value is demonstrated. The strength distributions of the small amplitude quadrupole oscillations are also shown.

연구 동기 및 목표

• 시간에 따라 변하는 Hartree-Fock-Bogoliubov (TDHFB) 방정식을 Gogny 효과적 상호작용와 함께 실용적이고 자기일관된 수치적 방법으로 해결하는 것.
• 이론적 방법이 이국적 핵에서 저에너지 집단 진동수를 기술하는 데의 가능성을 시험하는 것, 특히 쌍화 상호작용과 연속체 효과가 존재할 경우에 초점한다.
• TDHFB 프레임워크 내에서 시간 진동수 동안 핵 입자 수와 에너지와 같은 기본 양의 보존을 확보하는 것.
• 결과를 알려진 진동수 스펙트럼과 비교하여, 산소 이sovotope 20O에 적용하여 방법을 검증하는 것.
• 특히 최대 주량수 $N_{\text{shell}}$에 따라 기저 공간 截 截의 수렴성과 정확도를 평가하는 것.

제안 방법

• 각 시간 단계에서 시간 이동 연산자의 멱급수 전개를 사용하여 TDHFB 방정식을 적분함으로써 시간 진동수 연산자의 유니타리성을 보장한다.
• 입자-홀 및 입자-입자 채널 모두에 Gogny 상호작용을 사용하여 동일한 수준에서 다루며, 평균장과 쌍화 효과를 자기일관적으로 기술할 수 있도록 한다.
• 일반화된 밀도 행렬 $\mathcal{R}$ 는 $i\hbar\dot{\mathcal{R}} = [\mathcal{H}, \mathcal{R}]$ 식을 통해 시간에 따라 진화하며, 여기서 $\mathcal{H}$ 는 HFB 해밀토니안 행렬이다.
• 기저 상태는 $N_{\text{max}} = 10$ 까지의 3차원 조화 진동자 상태이며, 수치적 실현 가능성으로 인해 $N_{\text{shell}} = 4$ 에서 공간을 잘라낸다.
• 시간 진동수 계산은 $c\Delta t = 0.2$ fm의 시간 간격을 사용하며, 초기 상태는 2중극 모드를 자극하기 위해 충격형 외란을 통해 준비된다.
• 강도 함수는 2중극 연산자 기대값의 시간 시리즈를 푸리에 변환하여 구하며, 실험 및 QRPA 결과와의 비교를 가능하게 한다.

실험 결과

연구 질문

• RQ1 Gogny 상호작용를 사용한 자기일관된 TDHFB 프레임워크를 실용적인 시간 적분 기법으로 수치적으로 구현할 수 있는가?
• RQ2 Gogny 상호작용를 사용한 TDHFB 계산에서 시간 진동수 동안 입자 수와 에너지의 보존 정도는 어떠한가?
• RQ3 이 방법은 20O의 저에너지 2+ 진동수 스펙트럼을 얼마나 정확하게 재현하는가? 특히 QRPA 및 실험 결과와의 비교에서 어떻게 되는가?
• RQ4 유한한 기저 공간, 특히 $N_{\text{shell}} = 4$ 는 집단 모드의 에너지 및 강도 분포에 어떤 영향을 미치는가?
• RQ5 이 방법은 20O와 같이 탈형이 있고 이국적인 핵에서 소규모 진폭의 2중극 진동수의 본질적인 특징을 포착할 수 있는가?

주요 결과

• 이 방법은 총 양성자 수를 $10^{-11}$ 이내로 보존하여, 양성자 자유도의 시간 진동수에서 매우 높은 정확도를 보여준다.
• 통합 시간 동안 중성자 수는 $10^{-5}$ 이내로 보존되어 쌍화 채널에서의 유니타리성 유지가 양호함을 나타낸다.
• 시간 진동수 전반에 걸쳐 진동수 에너지가 0.2 MeV 이내로 보존되어 에너지 기대값의 안정성을 확인한다.
• 2중극 모멘트 시간 시리즈의 푸리에 변환을 통해 얻은 강도 함수는 약 3 MeV, 10 MeV, 24 MeV 에 세 개의 주요 피크를 보이며, 가장 낮은 피크는 약 3 MeV 근처에 있다.
• 강도 함수의 주요 피크 위치는 실험적 2+ 준위인 1.7, 4, 5, 10 MeV 와 유사하지만, 기저 공간의 제한으로 인해 계산된 에너지는 체계적으로 높게 나타난다.
• 강도 함수의 상대적 피크 높이가 연속체 효과를 포함한 QRPA 결과와 크게 다름을 보이며, 이는 기저 공간이 강도 분포를 정확히 기술하기에 너무 작다는 것을 시사한다.