[논문 리뷰] The {\it victory} project v1.0: an efficient parquet equations solver

Gang Li, Anna Kauch외 2명|arXiv (Cornell University)|2017. 08. 24.

Physics of Superconductivity and Magnetism참고 문헌 11인용 수 4

한 줄 요약

이 논문은 허버드 모델에 대한 파라케트 방정식을 위한 고성능 포트란 구현인 Victory v1.0을 제시한다. 이는 단일 및 이차 입자 상관 함수의 자기일관적 해를 가능하게 하며, 두 수준의 커널 근사와 마츠부라 주파수 공간에서의 고주파수 정규화를 통해 수렴성과 정확성을 확보한다. 고유값 분석을 통해 상전이 신호를 확인한다.

ABSTRACT

{\it Victory}, i.e. \underline{vi}enna \underline{c}omputational \underline{to}ol deposito\underline{ry}, is a collection of numerical tools for solving the parquet equations for the Hubbard model and similar many body problems. The parquet formalism is a self-consistent theory at both the single- and two-particle levels, and can thus describe individual fermions as well as their collective behavior on equal footing. This is essential for the understanding of various emergent phases and their transitions in many-body systems, in particular for cases in which a single-particle description fails. Our implementation of {\it victory} is in modern Fortran and it fully respects the structure of various vertex functions in both momentum and Matsubara frequency space. We found the latter to be crucial for the convergence of the parquet equations, as well as for the correct determination of various physical observables. In this release, we thoroughly explain the program structure and the controlled approximations to efficiently solve the parquet equations, i.e. the two-level kernel approximation and the high-frequency regulation.

연구 동기 및 목표

• 강한 상관 전자 시스템 맥락에서 파라케트 방정식을 위한 수치적으로 효율적이고 확장 가능한 솔버를 개발하기 위해.
• 마츠부라 공간에서 세 개의 운동량과 세 개의 주파수에 의존하는 큰 이차 입자 꼬리함수로 인한 계산 블로킹 문제를 해결하기 위해.
• 유한한 주파수 상자 내에서 주파수 공간의 열린 경계 조건을 적절히 처리하여 수렴성과 물리적 정확성을 확보하기 위해.
• 다양한 채널에서 베티-살파터 방정식의 고유값을 계산하여 상전이 연구를 가능하게 하기 위해.
• 허버드 모델에서 나타나는 잠재적 현상(예: 반자성 및 초전도성)을 연구하기 위한 병렬 처리, 오픈소스 도구를 제공하기 위해.

제안 방법

• 프로그램은 전체 파라케트 체계를 구현하여, 상호작용하는 베티-살파터, 슈윙거-다이슨, 그리고 파라케트 방정식을 자기일관적으로 풀어낸다.
• 이중 수준의 커널 근사를 사용하여 두 입자 꼬리함수의 복잡성을 줄이며, 더 적은 변수의 함수로 근사한다.
• 고주파수 정규화를 적용하여 꼬리함수의 점 渐진 주파수 의존성을 모델링하고, 유한한 주파수 상자 내에서 수렴성을 보장한다.
• MPI를 사용한 병렬 처리를 통해 점 대 점 메시지 교환을 최소화하고 확장성을 향상시킨다.
• 수치적 구현에서 이동 대칭성, SU(2) 스핀 대칭성, 시간 역행 대칭성을 명시적으로 적용한다.
• 자기일관적 해를 통해 다양한 채널(자기성, 쌍성, 전하)에서 베티-살파터 커널의 최대 고유값을 계산하여 임박한 상전이를 탐지한다.

실험 결과

연구 질문

• RQ1 어떻게 하면 운동량 및 마츠부라 주파수 공간에서 큰 이차 입자 꼬리함수를 가진 허버드 모델에 대해 파라케트 방정식을 효율적으로 풀 수 있는가?
• RQ2 유한한 주파수 및 운동량 격자에서 파라케트 방정식을 풀 때 수렴성과 물리적 정확성을 확보하기 위해 필요한 수치 근사는 무엇인가?
• RQ3 베티-살파터 방정식의 고유값은 어떻게 사용하여 강한 상관 시스템의 상전이를 탐지하고 특성화할 수 있는가?
• RQ4 고주파수 정규화는 마츠부라 주파수 공간에서 파라케트 방정식의 해를 안정화시키는 데 어떤 역할을 하는가?
• RQ5 6×6 크기의 클러스터에서 프로그램의 병렬 처리 전략은 다수의 코어에서 어떻게 스케일업되는가?

주요 결과

• 두 수준의 커널 근사와 고주파수 정규화는 유한한 주파수 상자 내에서 파라케트 방정식의 수렴성을 확보하는 데 필수적이다.
• 반도체 격자에서 반필링 상태에서 온도가 감소함에 따라 자기 채널의 고유값 λ는 단조적으로 증가하고 1에 수렴하며, 이는 반자성 불안정성을 시사한다.
• 쌍성 채널에서는 고온에서 s-파동 쌍성이 지배적이지만, 저온에서는 d-파동 쌍성($d_{x^2-y^2}$)이 지배적이 된다.
• 자기 채널에서 최대 고유값은 ${\bf q} = (\pi, \pi)$에서 나타나며, 이는 정사각형 격자에서의 반자성 질서와 일치한다.
• 프로그램은 베티-살파터 방정식의 고유함수를 성공적으로 계산하여, 운동량 공간에서의 파동함수 패턴을 통해 쌍성 대칭성을 식별할 수 있다.
• Victory v1.0은 추가 최적화 없이도 6×6 클러스터에서 임의의 전자 수를 가진 허버드 모델을 연구하는 데 이미 가능하다.