[논문 리뷰] An Electronically Non-Adiabatic Generalization of Ring Polymer Molecular Dynamics

Timothy J. H. Hele|arXiv (Cornell University)|2013. 08. 19.

Advanced Physical and Chemical Molecular InteractionsChemistry참고 문헌 7인용 수 17

한 줄 요약

이 논문은 전자적으로 비단일성인 시스템에서 양자 효과와 재교차 동역학을 정확하게 모델링할 수 있도록 레이어드 고리 고분자 분자역학(RPMD) 속도 이론의 비단일성 일반화를 제안한다. 이 방법은 고전적 RPMD를 전자 상태 결합을 반영한 수정된 잠재 에너지 표면을 통해 확장하여, 깊은 터널링 영역에서도 정확한 양자 속도와 요인 2 이내로 일치시키며, 계의 크기에 따라 선형적으로 증가하는 계산 비용을 보인다.

ABSTRACT

In this thesis I generalize Ring Polymer Molecular Dynamics (RPMD) rate theory to electronically non-adiabatic systems, followed by application to two one-dimensional curve crossing models and a multidimensional spin-boson model.

연구 동기 및 목표

• 핵의 운동이 여러 전자 잠재 에너지 표면 위에서 일어나는 전자적으로 비단일성인 시스템에 대해 레이어드 고리 고분자 분자역학(RPMD) 속도 이론을 확장하기.
• 양자 효과인 터널링과 재교차 동역학을 비단일성 과정에서 포괄하는 계산적으로 효율적인 방법 개발하기.
• 전자 상태가 단일 상태일 경우 기존의 단일성 RPMD로 수렴하고 고온 근처에서 물리적 일관성을 유지하는 방법 확보하기.
• 정확한 양자 기준과 함께 일차원 곡선 교차 모델과 다차원 스피너-보스온 모델에서 접근법 검증하기.
• 계산 비용이 계의 크기에 따라 선형적으로 증가하면서도 다른 근사 방법보다 정확도가 뛰어난 것을 보여주기.

제안 방법

• 고전적 이sovomerism을 일반화하여 양자 시스템을 조화 스프링으로 연결된 고전적 비드들의 고리 고분자로 표현하고, 전자 자유도를 포함한 확장된 위상공간을 도입한다.
• 해밀토니안에 수정된 결합 항을 포함시켜 레이어드 고분자에 대한 새로운 효과적 잠재 에너지 표면을 유도한다.
• 벤넷-캐놀러 인수분해를 적용하여 RPMD 속도를 분리 표면에서의 중심점 확률 밀도와 생성 영역에서의 양운동량 중심점의 유량의 곱으로 분해한다.
• 비대칭 곡선 교차 모델에서 최적의 분리 표면을 결정하기 위해 변분 방법을 사용한다.
• 확장된 위상공간에서 고전적 분자역학을 구현하고, 레이어드 고분자 중심점의 시간 상관 함수로부터 반응 속도를 계산한다.
• 스핀-보스온 모델에서 양자 전이 상태 이론(QTST) 속도가 발산하는 문제를 해결하기 위해, 두 개의 물리적으로 의미 있는 유한한 기여로 수렴 가능한 인수분해를 사용한다.

실험 결과

연구 질문

• RQ1 RPMD 속도 이론은 정확도와 계산 효율성을 유지하면서 전자적으로 비단일성 시스템으로 일반화될 수 있는가?
• RQ2 대칭 및 비대칭 결합이 있는 일차원 곡선 교차 모델에서 비단일성 RPMD 속도는 정확한 양자 속도와 어떻게 비교되는가?
• RQ3 반드시 반응 좌표 질량이 0에 수렴할 때, 디바이 스펙트럼 밀도를 가진 스피너-보스온 모델을 정확히 기술할 수 있는가?
• RQ4 전자 상태가 분리되었을 경우 비단일성 RPMD 방법은 정상적인 단일성 RPMD로 올바르게 복원되는가?
• RQ5 표준 QTST가 질량 요인의 발산 또는 정의되지 않은 속도로 실패하는 시스템에서 RPMD 속도를 신뢰성 있게 계산할 수 있는가?

주요 결과

• 대칭 일차원 곡선 교차 모델에서, 고전적 속도가 수개의 주기만큼 실패하는 100 K에서도 RPMD 속도는 정확한 양자 속도와 요인 2 이내로 일치한다.
• 비대칭 곡선 교차 모델에서는 대칭 케이스보다도 정확한 양자 속도를 더 잘 근사한다. 이는 단일성 근사의 경우와 유사한 편차를 보인다.
• 디바이 스펙트럼 밀도를 가진 스피너-보스온 모델에서, QTST는 질량 요인이 발산하여 정의되지 않지만, RPMD 속도는 수렴 가능한 인수분해를 통해 계산 가능하다.
• 비단일성 RPMD 속도는 거의 모든 테스트된 근사 방법(대표적인 평균장 및 양자역학적 반고전적 접근 포함)보다 정확한 양자 속도에 더 가깝게 나타난다.
• 전자 결합을 끄면 기존의 단일성 RPMD로 정확히 수렴함을 확인하여 기존 이론과의 물리적 일관성을 입증한다.
• 계산 비용이 계의 크기에 따라 선형적으로 증가하여 전자적으로 비단일성 동역학의 대규모 시뮬레이션에 적합하다.