[논문 리뷰] Symmetry-projected variational calculations with the numerical suite TAURUS I. Variation after particle-number projection

B. Bally, Adrián Sánchez-Fernández외 1명|arXiv (Cornell University)|2020. 10. 27.

Nuclear physics research studies참고 문헌 6인용 수 6

한 줄 요약

이 논문은 핵물리학에서 대칭성 프로젝션을 통한 변분 계산을 위한 공개된 Fortran 2008 코드인 TAURUS_vap을 제시한다. 이 코드는 구형 조화진동자 기저에서 일반적인 보골리우브 준입자 상태에 대해 입자수 프로젝션 이후 변분 계산을 가능하게 하며, 일반적인 하모닉 오실레이터 기저에서의 입자수 프로젝션을 통한 변분 계산을 지원한다. 이는 밸런스 스페이스 및 노 코어 계산을 모두 지원하며, 모든 핵 종류(짝-짝, 홀-짝, 홀-홀)를 처리할 수 있고, 탈형, 쌍화, 각운동량에 대한 동시 제약 조건을 허용하여 현대 하미르톤을 사용한 아비니티 및 셸 모델 연구에 매우 융통성 있는 도구를 제공한다.

ABSTRACT

We present the numerical code TAURUS_vap that solves the variation after particle-number projection equations for symmetry-unrestricted real Bogoliubov quasiparticle states represented in a spherical harmonic oscillator basis. The model space considered is invariant under spatial and isospin rotations but no specific set of orbits is assumed such that the code can carry out both valence-space and no-core calculations. In addition, no number parity is assumed for the Bogoliubov quasiparticle states such that the code can be used to describe even-even, odd-even and odd-odd nuclei. The variational procedure can be performed under several simultaneous constraints on the expectation values of a variety of operators such as the multipole deformations, the pairing field or the components of the angular momentum. To demonstrate the potential and versatility of the code, we perform several example calculations using an empirical shell-model interaction as well as a chiral interaction. The ability to perform advanced variational Bogoliubov calculations offered by this code will, we hope, be beneficial to the shell model and ab initio communities.

연구 동기 및 목표

• 핵many-body 시스템에서 입자수 프로젝션 이후 변분 계산을 위한 공개적이고 수치적으로 효율적인 코드 개발
• 기존 코드에서 널리 지원되지 않는 효과적 장 이론에서 유래한 하미르톤을 포함한 일반 하미르톤을 사용한 계산 가능화
• 고정된 활성 공간이나 수의 기현성 가정 없이 밸런스 스페이스 및 노 코어 구성 모두 지원
• 다중 관측량(다중극 탈형, 쌍화장, 각운동량 성분 등)에 대한 동시 제약 조건을 통합된 프레임워크에서 허용
• BLAS 및 LAPACK를 활용한 확장성 있고 병렬화된 수치적 프레임워크를 통한 고성능 계산을 위한 설계

제안 방법

• 코드는 구형 조화진동자 기저에서 일반적인 보골리우브 준입자 상태에 대해 실수값의 입자수 프로젝션 이후 변분 계산을 수행한다.
• 에너지 기대값을 최소화하기 위한 변분 절차를 사용하며, 연산자 기대값에 대한 제약 조건을 수반한다.
• 기저 대칭성(회전 대칭, 우도, 시간 역전 대칭, 입자 교환에 대한 반대칭성)을 활용하여 행렬 요소를 효율적으로 저장함으로써 메모리와 계산 자원을 절약한다.
• 편미분 인덱스와 위상 인자들을 사용하여 누락된 행렬 요소를 실시간으로 재구성함으로써 정확도를 유지하면서 최소한의 저장소를 사용한다.
• MPI를 통한 분산 메모리 병렬화를 지원하여 단일 노드의 메모리 한계를 초월한 큰 모델 공간을 처리할 수 있다.
• 수치적 안정성과 성능 향상을 위해 최적화된 BLAS 및 LAPACK 루틴을 활용한다.

실험 결과

연구 질문

• RQ1 일반적인 목적의 공개 코드가 구형 조화진동자 기저에서 임의의 보골리우브 상태에 대해 입자수 프로젝션 이후 변분 계산을 수행할 수 있는가?
• RQ2 이러한 코드가 현대 효과적 장 이론에서 유래한 일반 하미르톤을 얼마나 효율적으로 처리할 수 있는가?
• RQ3 대칭성 감소 전략과 행렬 요소 재구성 기법이 계산 속도를 희생시키지 않으면서 얼마나 메모리 사용을 줄일 수 있는가?
• RQ4 이 코드는 탈형, 쌍화, 각운동량과 같은 다중 관측량에 대한 동시 제약 조건을 통합된 프레임워크에서 지원할 수 있는가?
• RQ5 대규모 핵구조 계산에서 여러 프로세서에 걸쳐 코드의 성능가 어떻게 스케일링되는가?

주요 결과

• TAURUS_vap 코드는 모든 핵 종류(짝-짝, 홀-짝, 홀-홀)를 포함한 일반적인 보골리우브 상태에 대해 구형 조화진동자 기저에서 입자수 프로젝션 이후 변분 계산을 성공적으로 수행한다.
• 코드는 활성 공간이나 수의 기현성에 대한 사전 가정 없이도 일반 하미르톤(예: 경험적 셸 모델, 카이랄 상호작용 포함)을 처리할 수 있다.
• 회전 대칭, 우도, 시간 역전 대칭, 입자 교환에 대한 반대칭성 등의 대칭성을 활용함으로써 메모리 사용이 크게 감소하였으며, 행렬 요소의 일부만 저장된다.
• 편미분 인덱스의 사용으로 누락된 행렬 요소를 실시간으로 신속히 재구성할 수 있어, 큰 모델 공간에서도 계산 오버헤드를 최소화한다.
• MPI 기반 병렬화를 통해 효율적인 확장이 가능하여 단일 노드의 메모리 한계를 초월한 계산이 가능하다.
• 기존 솔버와의 벤치마킹을 통해 코드의 정확성과 성능이 확인되었으며, 고급 아비니티 및 셸 모델 연구에의 적용 가능성을 입증하였다.