[논문 리뷰] Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry

Mi-Song Dupuy, Gero Friesecke|arXiv (Cornell University)|2020. 02. 18.

Tensor decomposition and applications참고 문헌 28인용 수 6

한 줄 요약

이 논문은 양자화학에서 텐서 트레이스(TT) 근사에 대한 새로운 궤도 순서 정렬 방법을 소개한다. 이 방법은 비상호작용 해밀토니안의 기저 상태에서 나타나는 예상치 못한 역대칭성에 기반하여, 단일값(singular values)의 이론적 대칭성에서 유도된 숨겨진 불변량을 최적화함으로써, 단일값 尾(꼬리) 감쇠를 4–5개 주기 정도 감소시킨다. 이는 기존의 Fiedler 순서 정렬과 비교해 극적으로 뛰어난 성능을 발휘한다.

ABSTRACT

The tensor train approximation of electronic wave functions lies at the core of the QC-DMRG (Quantum Chemistry Density Matrix Renormalization Group) method, a recent state-of-the-art method for numerically solving the $N$-electron Schr\\"odinger equation. It is well known that the accuracy of TT approximations is governed by the tail of the associated singular values, which in turn strongly depends on the ordering of the one-body basis. Here we find that the singular values $s_1\\ge s_2\\ge ... \\ge s_d$ of tensors representing ground states of noninteracting Hamiltonians possess a surprising inversion symmetry, $s_1s_d=s_2s_{d-1}=s_3s_{d-2}=...$, thus reducing the tail behaviour to a single hidden invariant, which moreover depends explicitly on the ordering of the basis. For correlated wavefunctions, we find that the tail is upper bounded by a suitable superposition of the invariants. Optimizing the invariants or their superposition thus provides a new ordering scheme for QC-DMRG. Numerical tests on simple examples, i.e. linear combinations of a few Slater determinants, show that the new scheme reduces the tail of the singular values by several orders of magnitudes over existing methods, including the widely used Fiedler order.

연구 동기 및 목표

• 양자화학에서 텐서 트레이스(TT) 근사의 정확도를 좌우하는 핵심 과제인 기저 순서 정렬 문제를 해결한다.
• 비상호작용 해밀토니안에 대한 TT 근사에서 단일값의 행동, 특히 그들의 꼬리 감쇠가 궤도 순서에 어떻게 의존하는지 조사한다.
• 상관관계 있는 웨이브패킷에 대해 단일값 꼬리 감쇠를 최소화하는 이론적으로 탄탄한 최적화된 궤도 순서 정렬 기법을 개발한다.
• 기존 방법들(예: Fiedler 순서 정렬)보다 뛰어난 단일값 감쇠와 TT 근사 품질을 보이는 실용적이고 수치적으로 검증된 방법을 제공한다.

제안 방법

• 비상호작용 기저 상태를 나타내는 텐서의 단일값에 존재하는 새로운 역대칭성을 식별한다: $ s_1 s_d = s_2 s_{d-1} = \text{const} $, 여기서 $ s_1 \to s_d $ 는 감소하는 순서로 정렬된 단일값이다.
• 이 대칭성에서 유도된 숨겨진 불변량 $ I = s_1 s_d $ 를 정의하며, 이는 주어진 궤도 순서에 따른 단일값의 꼬리 행동을 완전히 특징짓는다.
• 이 불변량을 최소화함으로써 단일값 감쇠를 억제하는 '최적의 계수 순서 정렬' 전략을 제안한다.
• 성분으로서의 슬래터 결정자들로부터 유도된 불변량의 초과합을 사용하여, 이 방법을 상관관계 있는 웨이브패킷으로 확장한다.
• 일부 슬래터 결정자의 선형 조합을 시뮬레이션하여, 일반적인 양자화학 상태를 모의 실험한다.
• 수치 실험을 통해 새로운 방법이 단일값 감쇠 측면에서 Fiedler 순서 정렬 및 기타 표준 기법들과 비교하여 어떻게 성능을 냉각하는지 평가한다.

실험 결과

연구 질문

• RQ1 비상호작용 시스템에서 텐서 트레이스 근사의 단일값 꼬리 감쇠는 일체의 오비탈 궤도 순서에 어떻게 의존하는가?
• RQ2 비상호작용 기저 상태의 단일값 스펙트럼은 개선된 근사에 활용할 수 있는 숨겨진 대칭성을 보이는가?
• RQ3 슬래터 결정자 구성 요소에서 유도된 불변량을 사용하여 상관관계 있는 웨이브패킷의 단일값 꼬리 행동을 상한으로 제한할 수 있는가?
• RQ4 단일 불변량을 최소화하는 궤도 순서 정렬 기법이 Fiedler 순서 정렬과 같은 기존 방법보다 텐서 트레이스 근사 품질을 크게 향상시킬 수 있는가?

주요 결과

• 비상호작용 기저 상태(슬래터 결정자)의 텐서 트레이스 표현에서 단일값은 정확한 역대칭성을 보인다: $ s_1 s_d = s_2 s_{d-1} = \text{const} $, 이는 궤도 순서에 명시적으로 의존한다.
• 숨겨진 불변량 $ I = s_1 s_d $ 는 단일값의 꼬리 행동을 완전히 결정하며, 꼬리 제어의 복잡성을 단일 스칼라 양으로 줄인다.
• '최적의 계수 순서 정렬'을 통해 이 불변량을 최적화함으로써, 기존 표준 방법 대비 단일값 꼬리 감쇠가 4–5개 주기 감소한다.
• 새로운 방법은 단일값 감쇠에서 1개 주기의 향상만을 달성하는 Fiedler 순서 정렬보다 뛰어나다.
• 상관관계 있는 웨이브패킷의 경우, 개별 슬래터 결정자로부터 유도된 불변량의 초과합으로 꼬리 감쇠를 상한으로 제한할 수 있으며, 이는 일반화된 순서 정렬 기법을 가능하게 한다.
• 소수의 슬래터 결정자의 선형 조합에 대한 수치 실험 결과, 새로운 방법은 테스트 케이스 전반에서 일관되게 뛰어난 단일값 감쇠 성능을 보였다.