[논문 리뷰] Stochastic Multi Configuration Time-Dependent Hartree for Dissipative Quantum Dynamics with Strong Intramolecular Coupling

Souvik Mandal, Fabien Gatti외 3명|arXiv (Cornell University)|2022. 06. 16.

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

한 줄 요약

이 논문은 강한 결합된 비가역 양자 시스템에서 정확한 열화를 달성하기 위해 수정된 Lindblad 비가역 연산자를 갖춘 확률적 다구성시간의존하트리(sMCTDH) 방법을 제안한다. 일반 모드 좌표로의 변환과 일반화된 상승/하강 연산자의 사용을 통해 정확한 세부 균형과 버틀름만 통계를 강제함으로써, 조화계에서의 과열 행동을 성공적으로 제거하고 비조화계 및 페르미 공진의 경우 열화를 향상시킨다.

ABSTRACT

In this article, we explore the dissipation dynamics of a strongly coupled multidimensional system in contact with a Markovian bath following a system-bath approach. We use in this endeavour the recently developed stochastic Multi-Configuration Time-Dependent Hartree approach within the Monte Carlo wave packet formalism [J.Chem.Phys.156, 094109 (2022)]. The method proved to yield thermalized ensembles of wave packets when intramolecular coupling is weak. To treat strongly coupled systems, new Lindblad dissipative operators are constructed as linear combinations of the system coordinates and associated momenta. These are obtained by an unitary transformation to a normal mode representation, which reduces intermode coupling up to second order. Additionally, we use combinations of generalized raising/lowering operators to enforce the Boltzmann distribution in the dissipation operators, which yield perfect thermalization in the harmonic limit. The two ansatz are tested using a model two-dimensional hamiltonian parameterized to disentangle the effects of intramolecular potential coupling, of strong mode mixing observed in Fermi resonances, and of anharmonicity.

연구 동기 및 목표

• 강한 다차원 결합 시스템의 비가역 양자 역학에서 비물리적인 과열 행동을 해결하기 위해.
• 체계-환경 접근법을 사용하여 강한 분자 내 결합과 비조화성 존재 시 열화를 향상시키기 위해.
• 조화 근사에서 정확한 세부 균형과 버틀름 분포를 강제하는 비가역 연산자를 개발하기 위해.
• 다양한 결합 강도와 비조화성으로 구성된 모델 시스템에서 방법을 벤치마킹하기 위해.
• 열화 역학에서 모드 혼합과 다중모드 결합의 역할을 조사하기 위해.

제안 방법

• 열린 양자 시스템을 시뮬레이션하기 위해 몬테카를로 웨이브패킷 전파를 사용하는 확률적 MCTDH(sMCTDH) 프레임워크를 채택한다.
• 이중 모드 간 상호작용을 감소시키기 위해 국소 좌표에서 일반 모드로의 선형 좌표 변환을 도입한다.
• 일반 모드 좌표와 운동량의 선형 조합으로 새로운 Lindblad 비가역 연산자를 구성하며, 온도에 의존하는 계수를 사용한다.
• 일반 모드 표현에서 일반화된 상승/하강 연산자를 사용하여 정확한 세부 균형과 정확한 열 분포를 강제한다.
• MCTDH 파동함수 앤티츠에 호환되는 곱의 합 형태로 비가역 연산자를 구성한다.
• 비조화 강결합, 페르미 공진, 이차형 결합된 조화 포텐셜의 세 가지 모델 시스템에 방법을 적용한다.

실험 결과

연구 질문

• RQ1 일반 모드 변환은 강한 결합된 비가역 양자 시스템에서 과열 행동을 제거할 수 있는가?
• RQ2 일반화된 상승/하강 연산자는 기존의 좌표 기반 비가역 연산자에 비해 열화를 얼마나 향상시키는가?
• RQ3 공통의 열역학적 버스를 공유함에도 불구하고, 이차형으로 결합된 조화 모드는 왜 서로 다른 온도에서 열화하는가?
• RQ4 비조화성은 향상된 비가역 연산자 조차도 정확한 열화를 방해하는가?
• RQ5 강한 모드 혼합(예: 페르미 공진에서)은 약한 결합 조건에서도 단일 온도 평형에 도달하게 하는가?

주요 결과

• 일반 모드 변환은 2차까지의 다중모드 상호작용을 크게 감소시켜 조화계 및 비조화계에서 열화를 향상시킨다.
• 이차형으로 결합된 조화 포텐셜의 경우, 시스템은 200 K에서 기대되는 열 분포에 도달하지만, 온도에 의존하는 상수로 인해 400 K에서는 과냉각 행동을 보인다.
• 일반화된 상승/하강 연산자는 조화 근사에서 정확한 열화를 강제하여 페르미 공진의 경우 양 모드가 동일한 평형 온도에 도달하게 한다.
• 페르미 공진 시스템에서 강한 모드 혼합은 약한 결합 조건에서도 완전한 단일 온도 평형에 도달하게 한다.
• 이차형으로 결합된 조화 시스템은 차별적 열화를 보이며(Z-모드 ~200 K, R-모드 ~370 K) 비공명 주파수로 인해 별개의 열역학적 집합이 존재함을 시사한다.
• 비조화계는 여전히 작은 과열 이격을 보이며, 주로 에너지에 독립적인 비가역 연산자 때문이므로 향후 연구에서는 에너지에 의존하는 형태의 필요성이 제기된다.