[논문 리뷰] Sublinear scaling in non-Markovian open quantum systems simulations

Moritz Cygorek, Jonathan Keeling외 2명|arXiv (Cornell University)|2023. 04. 11.

Tensor decomposition and applications참고 문헌 99인용 수 4

한 줄 요약

이 논문은 가우시안 환경을 나타내는 텐서 네트워크의 자기유사성 특성을 활용하여 비마르코프성 열린 양자 시스템을 수치적으로 정확하게 시뮬레이션하는 알고리즘을 제안한다. 반복 가능한 블록을 기반으로 한 분할 정복 전략을 통해 무한 메모리의 경우 O(n log n) 스케일링, 절단된 메모리의 경우 O(n_c log n_c) 스케일링을 달성하여 수백만 단위의 시간 스텝에 걸친 시뮬레이션을 분 단위로 수행할 수 있다. 이는 양자점 광방출, 초방사성, 강한 결합 영역에서의 성능을 입증하였다.

ABSTRACT

While several numerical techniques are available for predicting the dynamics of non-Markovian open quantum systems, most struggle with simulations for very long memory and propagation times, e.g., due to superlinear scaling with the number of time steps $n$. Here, we introduce a numerically exact algorithm to calculate process tensors -- compact representations of environmental influences -- which provides a scaling advantage over previous algorithms by leveraging self-similarity of the tensor networks that represent Gaussian environments. Based on a divide-and-conquer strategy, our approach requires only $\mathcal{O}(n\log n)$ singular value decompositions for environments with infinite memory. Where the memory can be truncated after $n_c$ time steps, a scaling $\mathcal{O}(n_c\log n_c)$ is found, which is independent of $n$. This improved scaling is enabled by identifying process tensors with repeatable blocks. To demonstrate the power and utility of our approach we provide three examples. (1) We calculate the fluorescence spectra of a quantum dot under both strong driving and strong dot-phonon couplings, a task requiring simulations over millions of time steps, which we are able to perform in minutes. (2) We efficiently find process tensors describing superradiance of multiple emitters. (3) We explore the limits of our algorithm by considering coherence decay with a very strongly coupled environment. The algorithm we present here not only significantly extends the scope of numerically exact techniques to open quantum systems with long memory times, but also has fundamental implications for simulation complexity.

연구 동기 및 목표

• 장기 메모리 시간과 장기간의 전파 기간을 가진 비마르코프성 열린 양자 시스템을 시뮬레이션할 때 발생하는 계산적 병목 현상을 해결한다.
• 기존의 프로세스 텐서 알고리즘에서 발생하는 초선형 스케일링(O(n²) 또는 그 이상)이 시뮬레이션의 가능성을 제한하는 문제를 극복한다.
• 가우시안 환경에서 정확도를 유지하면서도 계산 비용을 극적으로 줄이는 수치적으로 정확한 방법을 개발한다.
• 강한 결합, 초방사성, 장시간 양자 얽힘 붕괴와 같은 복잡한 양자 동역학의 실용적 시뮬레이션을 가능하게 한다.
• 기본적인 시뮬레이션 복잡성에 대한 영향을 미치는 스케일러블한 텐서 네트워크 기반 열린 양자 시스템 시뮬레이션의 기반을 마련한다.

제안 방법

• 환경 영향의 텐서 네트워크 표현에서 나타나는 자기유사성을 활용하여 프로세스 텐서 내 반복 가능한 블록을 식별한다.
• 단일값 분해(SVD)를 사용하여 프로세스 텐서를 재귀적으로 분해하는 분할 정복 알고리즘을 구현하여 필요한 연산 수를 감소시킨다.
• 스핀-보존 모델과 같이 환경 상관 함수가 지수적으로 감쇠하는 가우시안 통계를 가지는 환경에 적용한다.
• 무한 메모리의 경우 알고리즘이 시간 스텝 수 n에 대해 O(n log n) 스케일링을 보이며, 절단된 메모리의 경우 n_c 단계 이후 O(n_c log n_c) 스케일링을 보여, 총 전파 시간과는 독립적이다.
• 프로세스 텐서의 구조를 이용하여 동일하거나 유사한 블록을 재계산하지 않고, 텐서 네트워크 내 재귀적 패턴을 활용한다.
• 기존의 MPO(Matrix Product Operators) 및 프로세스 텐서 방법과의 통합을 통해 효율적인 동역학 전파를 실현한다.

실험 결과

연구 질문

• RQ1 비마르코프성 열린 양자 시스템을 위한 수치적으로 정확한 알고리즘이 시간 스텝 수에 대해 초선형 스케일링을 달성할 수 있는가?
• RQ2 환경 메모리의 텐서 네트워크 표현에서 나타나는 자기유사성을 어떻게 활용하여 필요한 SVD 연산 수를 줄일 수 있는가?
• RQ3 기존의 프로세스 텐서 기법으로는 해결이 어려운 장기 메모리 및 장기간 전파 동역학을 얼마나 잘 시뮬레이션할 수 있는가?
• RQ4 이 알고리즘은 강한 도트-포논 결합을 가진 양자점이나 초방사성 발광체 어레이와 같은 물리적으로 관련된 모델에 적용 가능한가?
• RQ5 수백만 단위의 시간 스텝을 포함한 시뮬레이션에서 계산 시간 측면에서 실용적인 성능 향상은 어느 정도인가?

주요 결과

• 무한 메모리 환경의 경우 O(n log n) 스케일링, 절단된 메모리의 경우 n_c(메모리 절단점)에 대해 O(n_c log n_c) 스케일링을 달성하여 기존의 O(n²) 또는 그 이상의 스케일링에 비해 근본적인 개선을 이룬다.
• 강한 외부 드라이브와 강한 도트-포논 결합 조건에서의 양자점 광방출을 수백만 시간 스텝에 걸쳐 분 단위로 성공적으로 시뮬레이션하여 실용적 타당성을 입증하였다.
• 초방사성 발광체 시스템의 프로세스 텐서를 효율적으로 계산하여 집단적 방출 동역학을 고정밀도로 연구할 수 있도록 하였다.
• 강한 결합 환경에서의 얽힘 붕괴를 성공적으로 시뮬레이션하여, 어려운 매개변수 영역에서도 알고리즘의 강건성을 입증하였다.
• 시간 스텝 수에 따라 행렬 차원이 증가하더라도, 다양한 매개변수 범위에서 관측된 계산 시간은 실질적으로 선형 또는 초선형으로 증가하였다.
• 기존에 스케일링 제약으로 인해 계산적으로 불가능했던 비마르코프성 동역학의 수치적 정확한 시뮬레이션을 가능하게 하였다.