[논문 리뷰] A parallel multistate framework for atomistic non-equilibrium reaction dynamics of solutes in strongly interacting organic solvents

David R. Glowacki, Andrew J. Orr‐Ewing외 1명|Bristol Research (University of Bristol)|2014. 12. 13.

Spectroscopy and Quantum Chemical Studies인용 수 34

한 줄 요약

이 논문은 강한 상호작용을 보이는 용매에서 용질의 비평형 반응 동역학을 시뮬레이션하기 위해 CHARMM 내에서 병렬 처리 및 선형 스케일링 다중상태 실증 밸런스 궤도(MS-EVB) 프레임워크를 제시한다. CCSD(T)-F12 전자 구조 데이터와 Hellmann-Feynman 힘을 조합하여, CD3CN 용매에서 F + CD3CN → DF + CD2CN 반응을 정확하게 모델링하는 64상태 PES를 구축한다. 이는 실험적 일시적 스펙트럼과 리스크링 동역학을 고도로 정확하게 재현하며, 진동 상태가 자극된 경우의 경쟁적 파란 이동과 빨간 이동, 그리고 마이크로용매화 효과를 포함한다.

ABSTRACT

We describe a parallel linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM. Forces are obtained using the Hellman-Feynmann relationship, giving continuous gradients, and excellent energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to CCSD(T)-F12 electronic structure theory, we built a 64-state MS-EVB model designed to study the F + CD3CN -> DF + CD2CN reaction in CD3CN solvent. This approach allows us to build a reactive potential energy surface (PES) whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We use our PES to run MD simulations, and examine a range of transient observables which follow in the wake of reaction, including transient spectra of the DF vibrational band, time dependent profiles of vibrationally excited DF in CD3CN solvent, and relaxation rates for energy flow from DF into the solvent, all of which agree well with experimental observations. Immediately following deuterium abstraction, the nascent DF is in a non-equilibrium regime in two different respects: (1) it is highly excited, with ~23 kcal mol-1 localized in the stretch; and (2) not yet Hydrogen bonded to the CD3CN solvent, its microsolvation environment is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of its microsolvation environment results in a red shift. These two competing effects result in a post-reaction relaxation profile distinct from that observed when DF vibration excitation occurs within an equilibrium microsolvation environment. The parallel software framework presented in this paper should be more broadly applicable to a range of complex reactive systems.

연구 동기 및 목표

• 강한 상호작용을 보이는 용매에서 비평형 반응 동역학을 시뮬레이션하기 위한 확장 가능하고 병렬 처리 가능한 계산 프레임워크를 개발하기 위해.
• 특히 진동 자극과 변화하는 용매 환경을 수반하는 화학 반응 후의 일시적 비평형 상태를 정확하게 모델링할 수 있도록 하기 위해.
• 복잡한 반응성 시스템에 대해 정확성과 효율성을 균형 잡은 다중상태 실증 밸런스 궤도 에너지 표면(PES)을 구축하기 위해.
• 용액에서 시간에 따라 변화하는 진동 스펙트럼과 리스크링 속도와 같은 실험적 관측 가능성을 재현하기 위해.
• 반응 후 리스크링에서 진동 자극과 마이크로용매화 동역학 간의 경쟁을 포착하기 위해.

제안 방법

• CHARMM 분자 동역학 엔진 내에서 병렬 처리 및 선형 스케일링 다중상태 MS-EVB 프레임워크의 구현.
• 역동학에서 연속적인 기울기와 뛰어난 에너지 보존을 보장하기 위해 Hellmann-Feynman 힘의 사용.
• 정확한 PES 구축을 위해 고수준의 CCSD(T)-F12 전자 구조 데이터에 다차원 가우시안 결합 요소를 피팅.
• CD3CN 용매에서 F + CD3CN → DF + CD2CN 반응을 위한 64상태 MS-EVB 모델의 구축.
• 결과로 얻어진 PES를 사용하여 비평형 동역학을 시뮬레이션하고 시간에 따라 변화하는 관측 가능성을 계산.
• 반응 후 일시적 진동 스펙트럼, 진동 에너지 리스크링 속도, 마이크로용매화 진화의 분석.

실험 결과

연구 질문

• RQ1 강한 상호작용을 보이는 용매에서 비평형 반응 동역학을 모델링하기 위해 확장 가능하고 병렬 처리 가능한 계산 프레임워크를 어떻게 설계할 수 있는가?
• RQ2 반응 후 생성물의 비평형 동역학을 형성하는 데 있어, 경쟁적 리스크링 경로—진동 자극과 마이크로용매화—의 역할은 무엇인가?
• RQ3 다중상태 실증 밸런스 모델이 실험적 일시적 진동 스펙트럼과 리스크링 동역학을 정확하게 재현할 수 있는가?
• RQ4 진동 자극과 변화하는 용매 구조가 함께 작용할 때, CD3CN 용매에서 DF의 스펙트럼 진화에 어떻게 영향을 미치는가?
• RQ5 비평형 마이크로용매 환경은 평형 용매 환경과 비교해 리스크링 동역학을 어느 정도 변화시키는가?

주요 결과

• 64상태 MS-EVB 모델은 CD3CN 용매에서 DF의 실험적 일시적 진동 스펙트럼을 높은 정확도로 재현하였다.
• 반응 직후 H-탈리화 후 생성된 DF는 약 23 kcal mol⁻¹의 진동 에너지를 가지며, 매우 높은 진동 자극 상태에 있고, 용매화가 불량한 것으로 나타났다.
• DF의 진동 리스크링은 스펙트럼의 파란 이동을 유도하지만, 마이크로용매화 진화는 빨간 이동을 유도하여 복잡하고 단조롭지 않은 리스크링 프로파일을 초래한다.
• 이 두 효과 간의 경쟁은 DF가 평형 용매 환경에서 진동 자극 상태에 있을 때 관찰되는 경로와는 다름을 보였다.
• DF에서 용매로의 에너지 전달은 실험적 리스크링 속도와 일치하는 시간 스케일에서 발생하여, 모델의 역학적 정확성을 검증하였다.
• 병렬 MS-EVB 프레임워크는 복잡한 용매에서 대규모 비평형 반응 동역학을 효율적이고 정확하게 시뮬레이션할 수 있도록 하였다.