[논문 리뷰] Some Further Results for the Stationary Points and Dynamics of Supercooled Liquids

David J. Wales, Jonathan P. K. Doye|arXiv (Cornell University)|2003. 09. 02.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 4인용 수 98

한 줄 요약

이 논문은 초냉각된 액체의 포텐셜 에너지 표면을 모든 정상점(특히 국소 최소점 외에도)에 대해 유도구역으로 분할하기 위한 뉴턴-라프슨과 고유벡터 추적 기반 방법을 제시한다. 약한 상호작용을 보이는 하위계에서 헤시안 지수에 따른 정상점 수에 대한 해석적 표현을 유도하며, 저지수의 안장점(특히 지수-1)이 더 국소화되고 공간적으로 독립적인 변위를 보임을 발견하여, 동역학 모델에서 효과적인 전이 상태로 사용될 수 있음을 뒷받침한다.

ABSTRACT

We present some new theoretical and computational results for the stationary points of bulk systems. First we demonstrate how the potential energy surface can be partitioned into catchment basins associated with every stationary point using a combination of Newton-Raphson and eigenvector-following techniques. Numerical results are presented for a 256-atom supercell representation of a binary Lennard-Jones system. We then derive analytical formulae for the number of stationary points as a function of both system size and the Hessian index, using a framework based upon weakly interacting subsystems. This analysis reveals a simple relation between the total number of stationary points, the number of local minima, and the number of transition states connected on average to each minimum. Finally we calculate two measures of localisation for the displacements corresponding to Hessian eigenvectors in samples of stationary points obtained from the Newton-Raphson-based geometry optimisation scheme. Systematic differences are found between the properties of eigenvectors corresponding to positive and negative Hessian eigenvalues, and localised character is most pronounced for stationary points with low values of the Hessian index.

연구 동기 및 목표

• 모든 정상점(특히 국소 최소점 외에도)과 관련된 유도구역으로 포텐셜 에너지 표면을 체계적이고 안정적으로 분할하는 방법을 개발하는 것.
• 약한 상호작용을 보이는 하위계에서 각 헤시안 지수를 가진 정상점의 수에 대한 해석적 표현을 유도하여 이전의 통계역학 프레임워크를 확장하는 것.
• 헤시안 고유벡터의 공간적 국소화와 독립성, 특히 양의 고유값과 음의 고유값을 가진 고유벡터 간의 차이를 조사하는 것.
• 기존의 최소점과 진정한 전이 상태(지수-1)에 기반한 동역학 이론이 높은 지수를 가진 안장점이 지배하는 초냉각 액체에서의 표본 추출과도 호환되는지 평가하는 것.

제안 방법

• 뉴턴-라프슨 기반 기하 최적화 기법과 고유벡터 추적을 조합하여 포텐셜 에너지 표면 상의 모든 정상점을 탐색하고 분류하는 데 사용한다.
• 각 정상점에서부터 가장 급강하 경로를 따라 유도구역을 추적하여 그 유도구역을 식별하는 방식으로 유도구역 분할을 수행한다.
• 체계적이고 약한 상호작용을 보이는 하위계 기반의 조합론적 모델을 개발하여 시스템 크기와 헤시안 지수에 따라 정상점의 수를 해석적으로 계산한다.
• 공간적 확장(L)과 공간적 독립성(N)이라는 두 가지 척도를 사용하여 헤시안 고유벡터의 국소화 정도를 정량화하며, 양의 고유값과 음의 고유값을 가진 고유벡터를 비교한다.
• 256원자로 구성된 이원계 레너드존 시뮬레이션에서의 정상점 분석을 통해 해석적 모델의 타당성을 검증하고 고유벡터 성질을 평가한다.
• 다양한 헤시안 지수를 가진 고유벡터의 공간적 특성을 비교하며, 특히 저지수 안장점에 초점을 맞춘다.

실험 결과

연구 질문

• RQ1 모든 정상점(높은 지수의 안장점 포함)에 대해 체계적이고 수치적으로 안정적인 방식으로 포텐셜 에너지 표면을 유도구역으로 분할할 수 있는가?
• RQ2 약한 상호작용을 보이는 시스템에서 총 정상점 수, 국소 최소점 수, 최소점당 평균 전이 상태 수 사이의 해석적 관계는 무엇인가?
• RQ3 헤시안 고유벡터의 공간적 국소화와 독립성은 양의 고유값 모드와 음의 고유값 모드 간에 어떻게 다를까?
• RQ4 음의 고유값에 대응하는 고유벡터가 양의 고유값에 대응하는 고유벡터보다 더 국소화되고 공간적으로 독립적인가, 그리고 이 효과는 낮은 헤시안 지수에서 더욱 두드러지는가?
• RQ5 초냉각 액체에서 높은 지수의 안장점이 지배적인 경우에도 최소점과 진정한 전이 상태(지수-1)에 기반한 기존의 동역학 이론이 여전히 유효한가?

주요 결과

• 뉴턴-라프슨과 고유벡터 추적 기반 방법은 모든 정상점에 대해 포텐셜 에너지 표면을 성공적으로 유도구역으로 분할하여 국소 최소점 외의 분석도 가능하게 한다.
• 약한 상호작용을 보이는 하위계를 위한 해석적 모델은 총 정상점 수, 국소 최소점 수, 최소점당 평균 전이 상태 수 사이에 단순한 관계를 예측한다.
• 음의 헤시안 고유값에 대응하는 고유벡터는 양의 고유값에 대응하는 고유벡터보다 공간적으로 훨씬 더 국소화되고 독립적이다.
• 헤시안 고유벡터의 국소화와 공간적 독립성은 헤시안 지수가 낮을수록 특히 뚜렷하며, 이는 I=1인 경우에 가장 두드러진다.
• 결과는 저지수 안장점이 하위계의 진정한 전이 상태의 조합으로 효과적으로 기술될 수 있으며, Shell 등 이전의 연구와 일치함을 뒷받침한다.
• 초냉각 액체에서 높은 지수의 안장점이 널리 존재하더라도, 동역학은 여전히 최소점과 진정한 전이 상태로 잘 기술될 수 있으며, 가장 낮은 장벽은 여전히 지수-1 안장점을 통해 매개되기 때문이다.