[논문 리뷰] Exact rate calculations by trajectory parallelization and twisting

Eric Vanden‐Eijnden, Maddalena Venturoli|ArXiv.org|2009. 04. 23.

Quantum chaos and dynamical systems인용 수 9

한 줄 요약

이 논문은 전이 경로 이론을 사용하여 평형 및 비평형 시스템에서 정확한 반응 속도를 계산하기 위해 궤적의 병렬 처리 및 비틀림 기법을 제시한다. 베로이 테세레이션과 플럭스 매칭 경계 조건을 활용함으로써, 이 방법은 반응 궤적의 병렬 시뮬레이션을 가능하게 하고 정확하게 정방향 및 역방향 속도 상수를 계산한다. 2D 멜러 포텐셜 테스트 케이스에서 브루트 포스 시뮬레이션과 우수한 일치를 보이며 검증되었다.

ABSTRACT

A sampling procedure to compute exactly the rate of activated processes arising in systems at equilibrium or nonequilibrium steady state is presented. The procedure is a generalization of the method in [A. Warmflash, P. Bhimalapuram, and A. R. Dinner, J. Chem. Phys., 127, 154112 (2007); A. Dickson, A. Warmflash, and A. R. Dinner, J. Chem. Phys., 130, 074104 (2009)] in which one performs simulations restricted into cells by using a reinjection rule at the boundaries of the cells which is consistent with the exact probability fluxes through these boundaries. Our generalization uses results from transition path theory which indicate how to twist the dynamics to calculate reaction rates.

연구 동기 및 목표

• 통계적 평형 상태에 있는 임의의 평형 또는 비평형 시스템에서 정확한 반응 속도를 계산할 수 있는 일반적인 방법을 개발하기 위해.
• 이전의 궤적 병렬 처리 기법을 확장하여 플럭스 매칭 경계 조건을 통합함으로써 비평형 동역학을 다룰 수 있도록 하기 위해.
• 정확한 통계적 성질을 유지하면서 비틀림된 베로이 세포를 통해 반응 궤적의 병렬 계산을 가능하게 하기 위해.
• 마일스토닝, TIS, FFS를 일반화하기 위해 비단조화적 인터페이스와 임의의 세포 구성 조건을 허용하기 위해.
• 직접 브루트 포스 시뮬레이션과 정확한 비교를 통해 2D 멜러 포텐셜 시스템에서 방법을 검증하기 위해.

제안 방법

• 상태 공간을 생성점의 집합을 사용하여 베로이 세포로 분할하여, 상태 공간의 각 점이 가장 가까운 세포에 속하도록 보장한다.
• 각 세포에서 독립적으로 반응 궤적을 시뮬레이션하고, 출구 점을 저장하여 플럭스 매칭 재진입 규칙을 통해 인접한 세포로 궤적을 재초기화한다.
• 전이 경로 이론을 사용하여 동역학을 비틀어 반응물과 생성물 상태를 연결하는 궤적을 샘플링하도록 편향함으로써 정확한 속도 계산이 가능해진다.
• 비편향 및 비틀린 동역학에서 각 세포당 10⁸ 단계를 시뮬레이션하기 위해 Δt = 10⁻⁴인 2차 랑주완 통합기를 사용한다.
• 궤적 집합에서 평형 확률 πα와 παᴬ를 계산하며, A가 소스이고 B가 싱크일 때 παᴬ는 비평형 분포를 나타낸다.
• 시간 역행성이나 단조로운 인터페이스 순서에 대한 가정을 피하기 때문에, 복잡하고 비단조화적인 반응 경로에 적용 가능하다.

실험 결과

연구 질문

• RQ1 비평형 시스템에서 병렬 처리 및 세포 기반 시뮬레이션 프레임워크를 사용하여 정확한 반응 속도를 계산할 수 있는가?
• RQ2 비단조화적 인터페이스 순서를 가정하지 않고, 임의의 베로이 세포에 대해 궤적 병렬 처리를 어떻게 일반화할 수 있는가?
• RQ3 전이 경로 이론을 통한 동역학 비틀림이 속도 계산에서 정확한 통계적 성질을 얼마나 잘 유지하는가?
• RQ4 세포가 최적화되지 않거나 반응 좌표에 정렬되지 않은 경우에도 방법이 정확성을 유지할 수 있는가?
• RQ5 계산된 속도가 직접 장시간 비편향 시뮬레이션으로부터 구한 속도와 어떻게 비교되는가?

주요 결과

• 비틀린 동역학을 사용하여 계산된 반응 플럭스 νR은 5.9×10⁻³이었으며, 브루트 포스 값인 5.8×10⁻³과 매우 유사했다.
• 제안된 방법을 통해 도출된 정방향 속도 상수 kA,B는 7.4×10⁻³이었고, 직접 시뮬레이션 결과인 7.1×10⁻³과 비교하여 유사했다.
• 본 방법을 통해 도출된 역방향 속도 상수 kB,A는 3.1×10⁻³이었으며, 직접 시뮬레이션 결과인 3.2×10⁻³과 뛰어난 일치를 보였다.
• 본 방법을 통해 계산된 평형 확률 πα는 장시간 비편향 궤적을 통한 결과와 매우 우수한 일치를 보였으며, 그림 5에서 확인되었다.
• 비평형 조건에서 비틀린 동역학 하에 παᴬ는 B 쪽으로 단조롭게 감소하였으며, 이는 예상과 일치한다. 특히 π₃₆ᴬ는 구성상 0이었다.
• 세포 배치가 최적화되지 않은 경우에도 방법이 강인함을 입증하였으며, 36개의 무작위로 분포된 베로이 세포가 정확도를 저하시키지 않았다.