[논문 리뷰] iGVPT2 : an interface to computational chemistry packages for anharmonic corrections to vibrational frequencies

Loïc Barnes, Baptiste Schindler외 2명|arXiv (Cornell University)|2017. 04. 07.

Advanced Chemical Physics Studies참고 문헌 9인용 수 3

한 줄 요약

iGVPT2는 두 번째 순서의 진동 양자역학 섭동 이론(VPT2) 및 그 변종(VPT2+K, DCPT2, HDCPT2)을 사용하여 비조화 진동 주파수 보정을 계산하는 무료이고 오픈소스 소프트웨어입니다. 다양한 양자화학 패kg와 연동되며, 하이브리드 QM//MM 계산과 힘 상수의 병렬 계산을 지원하여 실제 계산 시간을 약 f³(여기서 f는 정규 모드의 수)만큼 감소시켜, 생체분자 포함 큰 분자에서의 비조화적 인프라레드 스펙트럼 예측을 효율적으로 가능하게 합니다.

ABSTRACT

iGVPT2 is a program for computing anharmonic corrections to vibration frequencies, based on force field expansion of the potential energy surface in normal mode coordinates. It includes second order vibrational perturbation theory (VPT2) algorithm and its derived methods (VPT2+K, DCPT2, HDCPT2). iGVPT2 is interfaced with several computation chemistry packages to compute the potential energies and dipoles derivatives. The second, third and quartic derivatives can be computed at the same level of theory but they can be also computed using different methods via one or two computational packages. iGPVT2 includes also a very fast hybrid QM//MM approach for biomolecules. It is provided free-of-charge for non-commercial research (see \url{https://sites.google.com/site/allouchear/igvpt2}).

연구 동기 및 목표

• VPT2 및 그 변종을 사용하여 비조화 진동 주파수 보정을 계산하는 유연하고 효율적인 인터페이스 개발
• 고수준 이론에서의 조화 주파수 계산과 저수준에서의 비조화 보정을 결합한 혼합 수준의 양자화학 계산 지원
• ORCA, Gaussian, Molpro, OpenMOPac 등 다양한 계산 화학 패kg를 플러그인 방식의 인터페이스로 연동하여 상호운용성 향상
• 대규모 생체분자에서 계산 비용을 줄이면서도 정확도를 유지하는 빠른 하이브리드 QM//MM 방법 구현
• 세차 및 사차 도함수 계산을 위한 고도로 병렬화된 프레임워크 제공로 실제 계산 시간을 약 f³(=정규 모드 수)만큼 감소

제안 방법

• iGVPT2는 VPT2를 포함한 변종(VPT2+K, DCPT2, HDCPT2)을 사용하여 비조화 보정을 계산합니다.
• ORCA, Gaussian, Molpro, FireFly, GAMESS-US, DFTB+ 등 다양한 양자화학 패kg와 연동하여 위치 에너지, 전기 dipole, 힘 도함수를 계산합니다.
• 소프트웨어는 하이브리드 계산을 지원합니다: 조화 주파수 계산은 고수준 방법(B3LYP/def2-TZVPP 등)에서 수행하고, 삼차 및 사차 힘 상수는 저수준(MMFF94 또는 반경험적 방법 등)에서 계산합니다.
• 조화 주파수의 수치 도함수는 사용자 정의 단계 크기를 가진 유한 차분 방법으로 계산되며, 원자와 정규 모드에 따라 병렬 처리됩니다.
• 하이브리드 QM//MM 방법은 QM 영역에 DFT를, 나머지 부분에 MMFF94를 사용하여 대규모 생체분자에서 빠르고 정확한 비조화 주파수 예측을 가능하게 합니다.
• 코드는 실시간 계산 시간을 약 f³(f = 정규 모드 수)만큼 감소시키기 위해 병렬화되어 있으며, 이는 세차 및 사차 도함수 계산을 크게 가속화합니다.

실험 결과

연구 질문

• RQ1 다양한 양자화학 패kg를 통해 단일 인터페이스로 비등방성 진동 주파수 보정을 효율적으로 계산할 수 있는가?
• RQ2 iGVPT2의 하이브리드 QM//MM 방법은 대규모 생체분자에서 계산 비용을 크게 줄이면서도 정확도를 유지하는가?
• RQ3 병렬 처리를 통해 세차 및 사차 도함수의 실시간 계산 시간을 약 f³만큼 감소시킬 수 있는가?
• RQ4 iGVPT2로 계산된 비조화 주파수는 실험 데이터 및 다른 소프트웨어 구현과 비교해 얼마나 정확한가?
• RQ5 예를 들어 고수준에서 조화 주파수, 저수준에서 비조화 보정을 계산하는 혼합 수준 계산은 전체 도함수를 고수준에서 계산하지 않아도 신뢰할 수 있는 결과를 도출할 수 있는가?

주요 결과

• 프로톤화된 디펩티드 GlyGlyH+의 경우, 전체 DFT(B3LYP/def2-TZVPP) 계산으로 비조화 주파수를 계산하는 데 109,888 CPU 시간이 소요되었지만, 2000개의 코어를 사용하여 실시간 계산 시간은 56시간으로 단축되었습니다.
• 하이브리드 QM//MM 방법(B3LYP/def2-TZVPP로 조화 주파수, MMFF94로 비조화 보정)을 사용할 경우 동일한 계산은 단 0.5시간(세차 및 사차 도함수 계산은 6초)에 완료되었습니다.
• 이전 연구에서 검증된 바와 같이, 하이브리드 QM//MM 방법은 실험적 비조화 주파수를 높은 정확도로 재현했습니다.
• 시험 분자들에 대한 최대 편차(MAX)는 1224 cm⁻¹ 이하였고, 기준 계산 결과와 비교한 평균 절대 편차(MAD)는 2.5 cm⁻¹ 이하였습니다.
• 일관된 단계 크기를 사용할 경우 비조화 주파수에 대한 평균 제곱 편차(RMSD)는 4.8 cm⁻¹ 이하였으며, 이는 높은 일관성과 신뢰성의 가능성을 시사합니다.
• 병렬 처리와 하이브리드 방법을 통해 대규모 시스템, 특히 생체분자에서도 정확한 비조화 주파수 예측이 가능하며, 빠른 계산 속도를 제공합니다.