Improving Perturbation Theory with the Sum-of-squares II: Large Density-Density Terms

기술 요약: 합의 제곱(Sum-of-squares) II를 통한 섭동 이론의 개선: 거대 밀도-밀도 항(Large Density-Density Terms)

문제 정의
본 논문은 양자 화학 해밀토니안에 자기 일관적(self-consistent) 합의 제곱(SOS) 분해를 적용할 때 발생하는 중대한 한계를 다룬다. 이전 연구(참조 [1])는 6차 SOS 방법이 섭동적으로 선택된 3차 연산자($\tau_i$)를 사용할 때, 약한 무작위 4차 상호작용을 가진 모델 해밀토니안에서 표준 4차 SOS(2RDM)보다 속도와 정확도 면에서 더 뛰어난 성능을 보임을 입증했다. 그러나 이 방법은 상호작용 항이 작다($\epsilon$ 차수)는 가정에 의존한다. 실제 양자 화학 해밀토니안은 강력한 "밀도-밀도" 상호작용 항($U_{ij} n_i n_j$)뿐만 아니라, 밀도 의존적 호핑(hopping) 및 스핀-스핀 상호작용을 포함하는 경우가 많다. 이러한 항들이 클 때(1 차수), 섭동적 방식의 3차 연산자 선택은 이들을 정확하게 처리하지 못하며, 결과적으로 이전 방법론을 이러한 응용 분야에 효과적이지 않게 만든다. 핵심 과제는 계산 효율성을 유지하고 기저 상태 에너지에 대한 엄격한 하한(lower bound)을 제공하는 능력을 유지하면서, 이러한 강력한 밀도-밀도 항을 처리할 수 있도록 자기 일관적 SOS 프레임워크를 확장하는 것이다.

방법론
저자들은 강력한 밀도-밀도 항을 수용하기 위해 3차 연산자($\tau_i$)의 정의와 SOS 분해의 구조를 수정할 것을 제안한다.

1. 수정된 해밀토니안 구조: 해밀토니안을 $H = H_0 + \epsilon V$로 재정의한다. 여기서 $H_0$는 이제 이차 부분($H_{quad}$)과 강력한 밀도-밀도 상호작용 항($\sum U_{ij} n_i n_j$)을 명시적으로 포함한다. $\epsilon V$ 항은 나머지 더 약한 4차 상호작용들을 포함한다.
2. 재규격화된 에너지 분모: 이전의 섭동적 접근 방식에서는 3차 연산자의 계수가 $1/(e_i + e_j + e_k + e_l)$ 형태의 에너지 분모에 의해 결정되었다. 강력한 밀도-밀도 항을 처리하기 위해, 저자들은 이 분모를 상호작용 에너지를 포함하도록 재정의한다. 구체적으로, 분모는 다음과 같이 된다:
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{pairs}} U_{pair}$$
   이는 연산자 $\psi_i + \tau_i$가 오차 $O(\epsilon^2)$까지 기저 상태를 소멸(annihilate)하게 함으로써, 방법론이 $O(\epsilon^4)$ 정확도의 에너지 하한을 제공하는 데 필요한 "오차 속성(error property)"을 보존하도록 보장한다.
3. 폐쇄성(Closure) 및 밀도 의존 항 처리: "폐쇄성"(유효 해밀토니안이 동일한 4차 해밀토니안 클래스 내에 머물도록 보장함)을 유지하고 밀도 의-존 호핑 항($n_i \psi^\dagger_j \psi_k$)을 처리하기 위해, 저자들은 더 일반적인 분해 구조를 도입한다. 이들은 시도 해밀토니안(trial Hamiltonian)을 $\sum \epsilon n_i V_i$ 형태의 항을 포함하도록 확장한다.
   새로운 분해는 $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$ 형태를 취하며, 여기서:
   $$F_a = \psi_i + n_j \sum_{k \neq j} d_k \psi_k + \theta_a$$
   $$\theta_a = \tau_i + \mu + \epsilon n_j \nu$$
   이다. 여기서 $\mu$와 $\nu$는 오차 속성을 교정하고 필요한 밀도 의존 항을 생성하기 위해 선택된 추가적인 3차 연산자들이다. 이 구조를 통해 반복적인 자기 일관적 절차를 통해 대상 해밀토니안 $H$를 재현할 수 있다.

주요 기여

• 강력한 상호작용으로의 확장: 본 논문은 강력한 밀도-밀도 상호작용 항을 포함하는 해밀토니안에 대한 SOS 분해를 생성하는 방법을 제공하며, 이는 기존에 자기 일관적 SOS 방법을 양자 화학 문제에 적용하는 데 장벽이 되었던 요소였다.
• 일반화된 연산자 정의: 저자들은 상호작용 강도를 고려한 재규격화된 에너지 분모를 사용하여 3차 연산자 $\tau_i$의 수정된 정의를 도입함으로써, 상호작용이 섭동적으로 작지 않은 경우에도 "오차 속성"이 유지되도록 하였다.
• 폐쇄성 보전: 수 연산자($n_i$)와 특정 3차 교정 항을 포함하는 보조 항들을 도입함으로써, 유효 해밀토니안이 여전히 반복적 풀이에 적합한 4차 연산자로 남도록 하는 폐쇄성을 어떻게 유지하는지 입증하였다.
• 네트워크 흐름 정식화: 본 논문은 이 방법으로 분해할 수 있는 해밀토니안의 범위를 규정한다. 음의 상호작용 항($U_{ij} < 0$)을 가진 해밀토니안에 대한 유효한 분해를 구성하는 문제를 네트워크 흐름 문제로 환원한다. 이는 주어진 상호작용 파라미터 집합을 제안된 분해가 처리할 수 있는지 판단하는 구체적인 기준을 제공한다.

결과 및 주장
본 논문은 제안된 방법이 다음을 수행한다고 주장한다:

• 정확도 유지: 밀도-밀도 항이 1 차수일 때도 실제 기저 상태 에너지에 대해 $O(\epsilon^4)$의 정확도를 갖는 기저 상태 에너지의 하한을 산출한다.
• 효율성 보전: 모든 항이 작은($O(\epsilon)$) 경우와 비교하여 계산 노력이 크게 높지 않은데, 이는 이 방법이 고차의 완전한 준정부호 계획법(SDP)을 푸는 대신 행렬 대각화와 반복적 업데이트에 의존하기 때문이다.
• 특정 항 처리: 이 방법은 밀도-밀도 상호작용을 성공적으로 처리하며, 부수적으로 밀도 의존 호핑 항($n_i(\psi^\dagger_j \psi_k + h.c.)$)과 스핀-스핀 상호작용(밀도-밀도 항의 합으로 매핑 가능)을 처리할 수 있다.
• 한계점: 저자들은 이 방법이 현재 특정 형태(입자-홀 공액 이후)의 밀도 의존 호핑만을 다룬다는 점을 언급한다. 쌍 생성($n_i \psi^\dagger_j \psi^\dagger_k$)과 관련된 항은 포함되지 않으나, 저자들은 이러한 항이 일반적으로 더 작으며 섭동적으로 처리될 수 있다고 제안한다. 또한, 싱글렛 호핑(singlet hopping) 항은 이차 기저 상태를 소멸시키는 성질을 공유하지 않으므로 이 특정 확장에 의해 다뤄지지 않는다.

의의
본 연구의 의의는 무작위 모델 해밀토니안에 대한 자기 일관적 SOS 방법의 이론적 성공을 실제 양자 화학 문제로 연결하는 데 있다. 강력한 밀도-밀도 상호작용이라는 장벽을 극로함으로써, 본 논문은 완전한 SDP 솔버의 계산 비용 없이 엄격한 에너지 하한을 제공하는 방법의 적용 범위를 확장한다. 저자들은 이 작업이 이전 연구의 필수적인 후속편이며, 이 비섭동적이고 섭동적으로 유도된 접근법을 통해 이전에는 도달할 수 없었던 물리적으로 유의미한 해밀토니안을 다룰 수 있게 한다고 위치시킨다.