기술 요약: 행렬 곱 상태(Matrix Product States)를 이용한 파인만 다이어그램 계산

1. 문제 정의

본 논문은 양자 나노 전자 시스템, 특히 단일 불순물 안데르센 모델(Single Impurity Anderson Model, SIAM)과 같은 비평형 상호작용 모델에서 물리적 관측량을 계산하는 과제를 다룹니다. 목표는 상호작용 강도 $U$에 대한 체계적인 섭동 전개로서 관측량 $Q(U, t)$ (전하 또는 전류와 같은)를 계산하는 것입니다:
$$Q(U, t) = \sum_{n=0}^{\infty} Q_n(t) U^n$$
계수 $Q_n(t)$를 계산하는 것은 비상호작용 그린 함수(Green's functions)의 곱에 대한 고차원 적분을 평가하는 과정을 포함합니다. 저자들은 이 접근 방식에서 세 가지 근본적인 어려움을 식별했습니다:

• 문제 A (조합 폭발): 섭동 차수 $n$이 증가함에 따라 파인만 다이어그램의 수가 계수적으로($O(n!)$) 증가하여, $n \approx 6-7$을 넘어서는 직접적인 합산이 불가능해집니다.
• 문제 B (부호 문제와 적분): 이러한 다이어그램의 피적분 함수는 종종 급격하게 진동하여 몬테카를로 방법의 "부호 문제(sign problem)"를 야기합니다. 또한, 필요한 함수 호출 횟수가 매우 많아(흔히 $\sim 10^9$번) 표준 샘플링은 비효율적이며, 특히 $n$에 따라 함수 호출당 비용이 지수적으로 증가합니다.
• 문제 C (수렴 및 재합산): 계수를 계산할 수 있다 하더라도, 섭동 급수는 유한한 수렴 반경을 가질 수 있으며(특히 정상 상태 극한 $t \to \infty$에서), 강한 상호작용에서의 물리적 결과를 추출하기 위해서는 정교한 재합산(resummation) 기술이 필요합니다.

2. 방법론

저자들은 전통적인 다이어그램적 몬테카를로(Diagrammatic Monte-Carlo)에서 벗어나 "비다이어그램적 비몬테카를로(Non-diagrammatic Non-Monte-Carlo)" 접근법으로 나아가는 3단계 프레임워크를 제안합니다.

1단계: 상호작용 전개 (문제 A 해결)

개별 파인만 다이어그램을 합산하는 대신, 저자들은 **켈디시 형식론(Keldysh formalism)**과 **윅 정리(Wick's theorem)**를 활용합니다.

• 그들은 전개 계수를 시간 변수와 켈디시 인덱스에 대한 다차원 적분으로 표현합니다.
• 결정적으로, 페르미온 연산자의 모든 순열에 대한 합(이 $n!$개의 다이어그램을 생성함)이 비상호작환 그린 함수로 구성된 **행렬의 행렬식(determinant)**과 수학적으로 동일하다는 점을 인식합니다.
• 이는 복잡도를 $O(n!)$에서 $O(2^n)$으로 감소시킵니다 (구체적으로는 $(2n+1) \times (2n+1)$ 행렬의 행렬식을 평가하는 것으로 감소).
• 핵심 통찰은 시간 적분 전에 켈디시 인덱스를 먼저 합산하면 연결되지 않은 다이어그램(disconnected diagrams)이 상쇄되어, 측정 시간 주변에 집중되는 잘 정의된 양의 정부호(positive-definite) 피적분 함수가 된다는 점입니다.

2단계: 텐서 교차 보간 (문제 B 해결)

몬테카를로 샘플링 없이 결과로 나온 고차원 적분을 평가하기 위해, 저자들은 **텐서 교차 보간(Tensor Cross Interpolation, TCI)**을 사용합니다.

• 개념: TCI는 고차원 텐서(피적분 함수)를 **행렬 곱 상태(Matrix Product State, MPS)**로 근사하는 학습 알고리즘입니다. 신경망과 달리, TCI는 경사 하강법을 사용하지 않고, 서브 행렬의 행렬식(maxvol 원리)을 최대화하는 "피벗(pivots, 특정 격자점)"을 선택하는 능동 학습 전략을 사용합니다.
• 메커니즘: 알고리즘은 텐서 차원을 반복적으로 훑으며, 텐서를 슬라이싱하여 형성된 행렬에 대해 **교차 보간(Cross Interpolation, CI)**을 수행합니다. CI는 적은 수의 행과 열(피벗)만을 사용하여 저계수(low-rank) 행렬 $A$를 근사합니다: $A \approx A_{I,J} A_{I,I}^{-1} A_{I,J}$.
• 통합: 피적분 함수가 MPS 표현으로 압축되면, 다차원 적분은 행렬-벡터 곱을 통해 정확하고 효율적으로 계산될 수 있으며, 이를 통해 차원의 저주와 부호 문제를 완전히 우회합니다.
• 성능: 이 방법은 $n \approx 20-25$까지의 계산을 가능하게 하여 이전의 한계를 훨씬 뛰어넘습니다.

3단계: 교차 외삽 (문제 C 해결)

유한한 수의 계산된 계수 $Q_n(t)$로부터 전체 함수 $Q(U, t)$를 재구성하기 위해, 저자들은 **교차 외삽(Cross Extrapolation)**을 사용합니다.

• 전략: 함수 $Q(U, t)$를 행은 상호작용 $U$, 열은 시간 $t$에 대응하는 행렬로 취급합니다. 저자들은 이 행렬이 대략적으로 저계수(low-rank) 구조를 가진다는 점에 주목합니다.
• 실행: TCI와 동일한 수학적 메커니즘을 사용하여, 알려진 "코너" 데이터(작은 $U$ 또는 작은 $t$)를 바탕으로 미지의 영역(큰 $U$ 및 큰 $t$)의 행렬을 외삽합니다.
• 타당성: 이는 $Q(U, t)$가 분리 가능한 함수들의 합 $\sum g_k(U)h_k(t)$로 근사될 수 있다는 가정에 기초합니다. 오차는 외삽의 랭크 $\chi$를 변화시킴으로써 제어됩니다.

3. 주요 기여

• 행렬식 재정의: 외평형 관측량의 섭동 전개가 행렬식의 적분으로 환원될 수 있음을 입증하여, 다이어그램 개수의 계수적 증가를 제거했습니다.
• 적분을 위한 TCI: 양자 다체 물리에서 발생하는 진동하는 고차원 함수를 적분하기 위한 몬테카를로의 우수한 대안으로서 텐서 교차 보간을 도입했습니다. 이는 확률적 노이즈 없이 정확한 적분을 제공합니다.
• 교차 외삽: $(U, t)$ 평면에서 관측량의 저계수 구조를 활용하여 섭동 급수를 재합산하는 강력한 방법을 개발함으로써, 급수가 본래 발산할 수 있는 영역에 접근할 수 있게 했습니다.
• 오픈 소스 구현: 이 연구는 tensor4all 라이브러리를 통해 지원되며, 텐서 네트워크 기술을 사용할 수 있는 도구를 제공합니다.

4. 결과

이 방법론은 단일 불순물 안데르센 모델(양자 점이 두 개의 리드에 결합된 형태)에 적용되어 전하와 전류를 계산하는 데 사용되었습니다.

• 정밀도: 저자들은 최대 8자리의 정밀도로 $n=18$까지의 전개 계수를 계산하였으며, 이는 정상 상태에서의 정확한 베테 이론(Bethe Ansatz) 결과와 일치합니다.
• 물리적 관측량: 이 방법은 다음을 포함한 SIAM의 전체 상도표(phase diagram)를 성공적으로 재현했습니다:
  • 쿨롱 다이아몬드(Coulomb Diamonds): 충전 에너지로 인해 전도도가 억제되는 영역.
  • 콘도 리지(Kondo Ridge): 콘도 효과에서 비롯되는 미분 전도도의 제로 바이어스 이상 현상.
• 전류-전압 특성: 저자들은 강한 상호작용($U=12$) 상태에서 콘도 리지를 가로지르는 전류 $I(V_b)$를 계산했습니다. 결과는 다음과 같은 예상되는 물리적 영역들을 포착했습니다:
  • 높은 바이어스($V_b > U$)에서의 포화.
  • 낮은 바이어스($V_b < T_K$)에서 완벽한 투과($I \approx V_b$)가 기대되는 콘도 온도 $T_K$와 관련된 플래토(plateau).
  • 이 영역 사이의 전이는 여러 에너지 스케일($U$, $\sqrt{\Gamma U}$, $\Gamma$, $T_K$)을 가로지릅니다.

5. 의의 및 주장

본 논문은 이 작업이 "컴퓨터에게 파인만 다이어그램을 계산하는 법을 가르치는 것"의 중요한 단계라고 위치시킵니다.

• 제어된 정확도: 평균장 근사(mean-field approximation)와 달리, 이 기술은 계산 자원(시간 및 랭크)을 늘림으로써 체계적으로 개선할 수 있는 알려진 오차 범위를 제공하는 "제어된(controlled)" 방식입니다.
• 패러다임의 전환: 저자들은 다이어그램 규칙을 수동으로 유도하는 것에서 벗어나, 문제의 기저에 깔린 저계수 구조를 자동으로 발견하는 알고리즘(TCI)을 사용하는 방식으로의 전환을 강조합니다.
• 확장성: 행렬식 감소와 TCI 적분의 결합은 다이어그램적 방법으로는 접근할 수 없었던 영역까지 계산을 가능하게 하여, 이론적 이해와 실험적 관측량 사이의 간극을 메워줍니다.
• 겸손한 태도: 저자들은 이 기술이 강력하지만, 모든 다체 문제를 해결하는 보편적인 솔버는 아님을 인정합니다 (예: 2D 허바드 모델은 향후 목표로 언급됨). 그들은 이 성공이 나노 전자 시스템의 피적분 함수가 갖는 특정한 저계수 구조에 의존한다는 점을 강조합니다.