[논문 리뷰] Spin relaxation dynamics with a continuous spin environment: the dissipaton equation of motion approach

Wenxiang Ying, Su Yu외 3명|arXiv (Cornell University)|2023. 02. 01.

Spectroscopy and Quantum Chemical Studies인용 수 4

한 줄 요약

이 논문은 연속 스핀 배경에서 스핀 리듬화 동역학을 위한 소산 방정식 운동(DEOM) 접근법을 개발하며, 임의의 스핀 양자수 $S$에 대해 변동-소산 정리(FDT)를 유도하고, 시간 영역 프로니 피팅을 사용하여 지수 감쇠 기저 함수를 구성한다. 이 방법은 표준 스핀-보존 모델을 초월한 비가우시안, 비조화 환경에 대한 정확하고 효율적인 시뮬레이션을 가능하게 한다.

ABSTRACT

We present the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the su(2) Lie algebra. This represents an extreme class of anharmonic environment. Based on the conclusion drawn by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we derive the fluctuation-dissipation theorem (FDT) of the spin bath from a microscopic perspective, and generalize the discussion to the case of arbitrary bath spin quantum number S. Next, the time-domain Prony fitting decomposition scheme is applied to the bare-bath time correlation function (TCF) given by FDT to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy has been justified by a variety of numerical results. We envision this work provides new insights to extend the hierarchical equations of motion (HEOM) and DEOM approach to certain types of anharmonic enviroments with arbitrary TCF or spectral density

연구 동기 및 목표

• 두 수준 시스템이 임의의 스핀 양자수 $S$를 가진 연속 스핀 배경에 결합된 경우를 위한 양자 동역학 프레임워크를 개발하기 위해.
• 기존 결과를 비영인 $S$로 일반화하기 위해, 미세한 시각에서 이러한 스핀 배경에 대한 변동-소산 정리(FDT)를 유도하기 위해.
• 스핀 배경의 시간 상관 함수(TCF)의 시간 영역 프로니 피팅(t-PFD)을 통해 지수 감쇠 기저를 구성함으로써 DEOM 형식을 비조화 환경에 효과적으로 확장하기 위해.
• 다양한 배경 매개변수와 온도 영역에서 방법의 정확성과 효율성을 검증하기 위해.

제안 방법

• 임의의 스핀 양자수 $S$를 가진 독립적인 $\mathfrak{su}(2)$-대수적 수준 시스템으로 구성된 스핀 배경에 대해 FDT를 유도하여, 효과적 스펙트럴 밀도 $J_{\text{eff}}(\omega;\beta,S)$ 를 도출한다.
• 시간 영역 프로니 피팅 분해(t-PFD)를 사용하여 순수한 배경 TCF를 지수 감쇠의 합으로 근사함으로써, DEOM 구성의 기초를 마련한다.
• 유도된 가짜 모드 표현을 바탕으로 소산자 계층을 구성하여, 비조화 스핀 배경이 존재하는 비마르코프 동역학을 기술한다.
• 선형 반응 한계와 효과적 스펙트럴 밀도 매핑을 활용하여, 임의의 스펙트럴 밀도와 비가우시안 환경으로의 접근을 일반화한다.
• 기존 기준과의 비교를 통해 수치적으로 방법을 검증하며, 특히 스핀 배경이 가우시안(보존) 배경으로 감소하는 고스핀 근사에서의 결과를 중심으로 검증한다.
• 배경의 커프오프 주파수와 결합 강도가 다양하게 설정된 모델 시스템에 프레임워크를 적용하여 수렴성과 계산 비용을 평가한다.

Figure 1: Population dynamics of the zero-temperature SSB models. The spin baths are parameterized as (a) $\alpha=0.5$ , $\omega_{c}/\Delta=1$ . (b) $\alpha=0.1$ , $\omega_{c}/\Delta=6$ . (c) $\alpha=0.2$ , $\omega_{c}/\Delta=10$ . (d) $\alpha=0.5$ , $\omega_{c}/\Delta=10$ . (e) $\alpha=0.75$ , $\om

실험 결과

연구 질문

• RQ1 임의의 스핀 양자수 $S$를 가진 스핀 배경에 대해 미세한 시각에서 변동-소산 정리를 어떻게 유도할 수 있는가?
• RQ2 시간 영역 프로니 피팅 전략이 다양한 $S$와 온도에서 스핀 배경의 순수한 배경 시간 상관 함수를 얼마나 정확하게 표현하는가?
• RQ3 표준 보존 배경 가정을 초월하여 비가우시안, 비조화 환경인 스핀 배경을 기술하기 위해 DEOM 형식이 효과적으로 확장될 수 있는가?
• RQ4 t-PFD 기반 DEOM의 정확도가 지수 항의 수와 $\alpha$ 및 $\omega_c/\Delta$와 같은 배경 매개변수에 따라 어떻게 달라지는가?
• RQ5 유 end-온도 스핀 배경에서 다양한 t-PFD 기법을 사용할 경우 계산 비용-정확도의 상호 보완 관계는 어떻게 되는가?

주요 결과

• 스핀 배경의 FDT는 $ C(t) = \frac{1}{\pi} \int_{-\infty}^{\infty} d\omega \, e^{-i\omega t} \frac{J_{\text{eff}}(\omega;\beta,S)}{1 - e^{-\beta\omega}} $ 로 유도되며, 여기서 $ J_{\text{eff}}(\omega;\beta,S) = J(\omega) \zeta(\omega;\beta,S) $ 이고, $ \zeta(\omega;\beta,S) $ 는 스핀 통계를 반영한다.
• $ S = 1/2 $ 일 경우, $ \zeta(\omega;\beta,1/2) = \tanh(\beta\omega/2) $ 로 유도되며, 문헌에서 알려진 결과를 회복한다.
• 고스핀 근사 $ S \gg 1 $ 에서 스핀 배경은 홀스타인-프리마코프 변환에 의해 보존 배경과 동형이 되며, 표준 스핀-보존 모델을 회복한다.
• 4+4 항을 사용한 t-PFD 전략은 온도가 유한한 모델에서 TCF의 실수부와 허수부를 모두 정확하게 피팅하며, 시각적 검토를 통해 높은 정확도를 확인한다.
• 적은 수의 지수 항으로도 높은 정확도를 달성하여 복잡한 비조화 환경에 대해 계산 효율성과 확장 가능성을 입증한다.
• 수치적 결과는 다양한 배경 매개변수, 특히 높은 커프오프 주파수와 강한 결합($ \alpha = 10 $)에서도 방법이 강인함을 확인하며, 온도가 유한한 조건에서도 정확도를 유지한다.

Figure 2: Convergence test of the localization model ( $\omega_{c}/\Delta=1,\ \alpha=10$ ) using $t$ -PFD strategies $2+2$ , $3+3$ and $4+4$ . Comparisons are made against ML-MCTDH (digitized from Ref. 33 ).