MPS와 TTNS의 비아다믹 비교: 결합 차원 확대가 가져온 정확도 회복

초록

본 연구는 Renormalizer 통합 프레임워크를 이용해 P3HT:PCBM 이종접합 모델에서 MPS와 TTNS의 비아다믹 시뮬레이션을 재검증한다. 이전 보고서에서 나타난 60% 수준의 차이는 충분히 큰 결합 차원(M)을 사용하지 않아 발생한 것이며, M을 확대하면 약한 전자 상태 점유율의 상대 오차가 10% 이하, 절대 오차가 0.005 이하로 감소한다. 최적화된 TTNS 트리 구조는 계산 효율성을 높이며, 두 방법 모두 결합 차원을 충분히 키우면 수치적으로 정확한 해에 수렴함을 확인하였다.

상세 분석

본 논문은 텐서 네트워크 기반 양자 동역학 시뮬레이션에서 가장 널리 사용되는 두 가지 ansatz, 즉 matrix product states (MPS)와 tree tensor network states (TTNS)의 수렴 특성을 체계적으로 비교한다. 핵심 모델은 P3HT:PCBM 이종접합에서의 exciton dissociation이며, 전자 26개 상태와 113개의 진동 모드(로컬 OT 진동 8×13, 풀러렌 진동 8, 인터분자 진동 R)로 구성된 고차원 해밀토니안을 포함한다. 기존 연구에서는 장시간 시뮬레이션 후 MPS와 TTNS 사이에 최대 60%에 달하는 점유율 차이가 보고되었는데, 이는 두 방법이 동일한 물리적 현상을 기술함에도 불구하고 구현상의 차이(소프트웨어, 트리 구조, 결합 차원 설정) 때문이라고 추정되었다.

연구팀은 Renormalizer 패키지를 활용해 MPS와 TTNS를 동일한 알고리즘(시간 의존 변분 원리 기반 프로젝트 스플리팅, TDVP‑PS)으로 구현함으로써 소프트웨어 차이를 원천적으로 제거하였다. 이후 결합 차원 M을 단계적으로 증가시켜 (M=64, 128, 256 등) 시뮬레이션을 수행했으며, 각 단계에서 (1) 전자 상태 점유율, (2) 전체 에너지 보존, (3) 엔트로피 기반 결합 차원 요구량을 정량적으로 평가했다. 결과는 다음과 같다.

1. 오차 감소: M을 64에서 256으로 확대함에 따라 약한 전자 상태(예: CS 2, CS 3)의 점유율 상대 오차가 60% → 9% 수준으로 급격히 감소했고, 절대 오차는 0.05 → 0.0048으로 수렴했다. 이는 결합 차원이 충분히 크면 양자 얽힘을 정확히 포착한다는 이론적 기대와 일치한다.

2. 엔트로피와 결합 차원 관계: TTNS 트리 구조에서 각 결합을 기준으로 von Neumann 엔트로피 S를 계산한 결과, 최대 S≈3.2 nats(≈ln 25) 수준이었다. 따라서 M≥e^S≈25이면 충분히 정확한 근사가 가능함을 확인했으며, 실제 시뮬레이션에서는 M≈64가 안전한 여유를 제공한다는 결론을 얻었다.

3. 구조 최적화: 기존 연구에서 사용된 “선형” MPS와 “균등 3‑branch” TTNS보다, 전자-진동 상호작용이 집중되는 부위(LE 1↔CS 1 전이와 R 모드) 주변에 높은 차원을 배치하고, 먼쪽(멀리 떨어진 OT 진동)에는 낮은 차원을 할당하는 비균등 트리 구조를 설계했다. 이 최적화 TTNS는 동일한 M에서도 연산량이 30% 감소하면서 오차는 기존 구조보다 2배 이하로 억제되었다.

4. 계산 효율성: MPS는 선형 연결 구조 덕분에 메모리 복잡도가 O(N M²)이며, TTNS는 트리 깊이에 따라 O(N M³) 정도로 약간 더 무겁다. 그러나 최적화된 TTNS는 결합 차원을 비균등하게 배분함으로써 실제 실행 시간은 MPS와 비슷하거나 오히려 빠르게 진행되었다. 이는 TTNS가 얽힘이 집중된 부분을 효율적으로 표현한다는 점을 실증한다.

5. 수치적 정확도 검증: 최종적으로 M=256, 최적화 TTNS 구조를 적용한 결과, 전체 파동함수의 정규화 오차는 10⁻⁸ 이하이며, 에너지 보존율도 0.999999 수준을 유지했다. 이는 TD‑DMRG와 ML‑MCTDH가 이론적으로 동일한 정확도 한계에 도달했음을 의미한다.

결론적으로, 두 방법 사이의 차이는 “결합 차원 부족”이라는 전형적인 텐서 네트워크 수렴 문제에서 비롯된 것이며, 충분히 큰 M과 물리적 상호작용 구조를 반영한 트리 설계가 이루어지면 MPS와 TTNS는 동일한 수치적 정확도를 제공한다. 이는 향후 복잡한 유기·무기 이종계, 광전지, 초고속 전자 전달 현상 등에 두 방법을 자유롭게 선택해 적용할 수 있는 기반을 마련한다.