[논문 리뷰] Analytic influence functionals for numerical Feynman integrals in most open quantum systems

Nikesh S. Dattani, Felix A. Pollock외 1명|arXiv (Cornell University)|2012. 03. 20.

Spectroscopy and Quantum Chemical Studies참고 문헌 34인용 수 7

한 줄 요약

이 논문은 개방 양자 시스템에서 이산화된 영향 함수를 위한 완전한 해석적 공식을 제시하며, 스펙트럼 밀도의 전통적인 수치적으로 불안정한 비정적 적분을 대체로, 환경 반응 함수 α(t)의 유한하고 해석적으로 다룰 수 있는 적분으로 전환한다. 높은 물리적 의미를 지닌 스펙트럼 분포들—예를 들어 로렌츠-드루드 및 거듭제곱 법칙 형태—에 대해, α(t)가 지수 함수의 합인 경우, 이러한 적분은 닫힌 형태의 표현식을 제공하며, 이는 더 빠르고 정확하며 더 엄밀한 수치적 파인만 경로 적분 계산을 가능하게 한다.

ABSTRACT

Fully analytic formulas, which do not involve any numerical integration, are derived for the discretized influence functionals of a very extensive assortment of spectral distributions. For Feynman integrals derived using the Trotter splitting and Strang splitting, we present general formulas for the discretized influence functionals in terms of proper integrals of the bath response function. When an analytic expression exists for the bath response function, these integrals can almost always be evaluated analytically. In cases where these proper integrals cannot be integrated analytically, numerically computing them is much faster and less error-prone than calculating the discretized influence functionals in the traditional way, which involves numerically calculating integrals whose bounds are both infinite. As an example, we present the analytic discretized influence functional for a bath response function of the form α(t)=\sum_{j}^{K}p_{j}e^{Ω_{j}t}, which is a natural form for many spectral distribution functions (including the very popular Lorentz-Drude/Debye function), and for other spectral distribution functions it is a form that is easily obtainable by a least-squares fit . Evaluating our analytic formulas for this example case is much faster and easier to implement than numerically calculating the discretized influence funcitonal in the traditional way. In the appendix we provide analytic expressions for p_{j} and Ω_{j} for a variety of spectral distribution forms, and as a second example we provide the analytic bath response function and analytic influence functionals for spectral distributions of the form J(ω)\proptoω^{s}e^{-(ω/ω_{c})^{q}}. The value of the analytic expression for this bath response function extends beyond its use for calculating Feynman integrals.

연구 동기 및 목표

• 파인만 경로 적분 계산에서 전통적인 비정적 적분의 계산 부담과 수치적 불안정성을 제거하기 위해.
• 스펙트럼 밀도 J(ω)가 아닌 환경 반응 함수 α(t)를 기반으로 한 이산화된 영향 함수(DIF)의 일반적 해석적 표현식을 유도하기 위해.
• 다양한 물리적으로 관련된 스펙트럼 분포에 대해 정확하고, 신속하며 강건한 수치적 DIF 평가를 가능하게 하기 위해.
• 즉각적인 구현을 위한 오픈소스 MATLAB 및 Mathematica 코드를 제공하기 위해.
• HEOM 및 나카지마-츠반징 방정식과 같은 근사 양자 동역학 방법의 벤치마킹을 지원하기 위해 일관되고 해석적으로 정확한 환경 결합 모델을 확보하기 위해.

제안 방법

• 환경 반응 함수 α(t)의 적절한 적분을 통해 이산화된 영향 함수의 일반적 표현식을 유도하여, J(ω)의 무한 구간 적분을 대체한다.
• α(t) = ∑ⱼ pⱼ e^{Ωⱼ t} 형태로 표현되는 경우에 집중하며, 이는 로렌츠-드루드 및 디바 함수와 같은 스펙트럼 밀도에서 자연스럽게 유도된다.
• α(t)가 해석적일 경우, DIF에 대한 해당 적분은 유한한 적분 구간 덕분에 정확히 또는 높은 수치 효율성으로 평가될 수 있음을 보여준다.
• 스펙트럼 분포 매개변수(예: J(ω) = ω^s e^{-(ω/ω_c)^q})와 역온도 β에 따라 pⱼ 및 Ωⱼ의 명시적 해석적 표현식을 제공한다.
• 페데 근사(예: [N-1/N])를 사용하여 페르미-디랙 및 보즈-아인슈타인 분포를 정확하게 표현함으로써 α(t) 분해를 수행한다.
• MATLAB 및 Mathematica에서 구현하여 오픈소스 코드를 제공하여 스펙트럼 매개변수에서 직접 해석적 DIF를 계산할 수 있도록 한다.

실험 결과

연구 질문

• RQ1 파인만 경로 적분에서 이산화된 영향 함수는 환경 반응 함수 α(t)의 유한하고 해석적으로 풀 수 있는 적분으로 표현될 수 있는가?
• RQ2 어느 정도의 물리적으로 관련된 스펙트럼 분포에서 α(t)가 ∑ⱼ pⱼ e^{Ωⱼ t} 형태를 띠며 해석적 DIF 평가가 가능한가?
• RQ3 J(ω)의 무한 구간 적분을 α(t)의 유한 적분으로 대체함으로써 수치적 파인만 적분 계산의 속도, 정확도, 안정성이 어떻게 향상되는가?
• RQ4 해석적 DIF는 HEOM 및 나카지마-츠반징 방정식과 같은 근사 개방 양자 시스템 동역학 방법의 벤치마킹을 어느 정도 지원할 수 있는가?
• RQ5 표준 수치 소프트웨어를 사용하여 해석적 DIF를 실제로 효율적으로 계산하고 구현할 수 있는가?

주요 결과

• 논문은 스펙트럼 밀도 J(ω)의 수치적으로 불안정한 비정적 적분이 필요 없도록, 환경 반응 함수 α(t)의 적절한 적분을 기반으로 한 이산화된 영향 함수의 일반 공식을 유도한다.
• 로렌츠-드루드/디바 및 거듭제곱 법칙 형태(J(ω) ∝ ω^s e^{-(ω/ω_c)^q})와 같은 스펙트럼 분포에 대해, 해당 α(t)는 ∑ⱼ pⱼ e^{Ωⱼ t} 형태를 띠며, 이는 완전한 해석적 DIF를 가능하게 한다.
• 해석적 DIF는 계산적으로 열등하다: 기존의 무한 구간 수치 적분에 의존하는 방법보다 평가 속도가 더 빠르고 수치 오차의 위험이 적다.
• 페르미-디랙 및 보즈-아인슈타인 형태를 페데 근사를 통해 표현함으로써, 세 가지 주요 스펙트럼 분포 유형에 대해 pⱼ 및 Ωⱼ의 명시적 해석적 표현식을 도출하였다.
• 복잡한 환경 스펙트럼 밀도를 가진 시스템, 특히 거듭제곱 법칙 형태 J(ω) ∝ ω^s e^{-(ω/ω_c)^q}에 대해서도 정확한 DIF 계산이 가능하며, 상부 비정적 감마 함수 및 삼각함수 항을 포함한 α(t)에 대한 해석적 표현식을 유도하였다.
• 저자들은 스펙트럼 매개변수에서 직접 해석적 DIF를 계산할 수 있도록 오픈소스 MATLAB 및 Mathematica 코드를 제공하여 광범위한 재현 가능성과 사용성을 보장한다.