[논문 리뷰] Identification of the invariant manifolds of the LiCN molecule using Lagrangian descriptors

F. Revuelta, R. M. Benito외 1명|arXiv (Cornell University)|2021. 10. 19.

Advanced Chemical Physics Studies참고 문헌 57인용 수 8

한 줄 요약

이 연구는 리간그리안 디스크립터(LDs)를 사용하여 LiCN⇌LiNC 이성화 반응에서의 불변 다각형을 식별하며, LD에서의 특이점으로서 다각형을 해소하기 위해 긴 통합 시간이 필수적임을 보여준다. 이는 아비시노 결과를 정확하게 재현하는 2차원 모스 기반의 포텐셜 에너지 표면(PES)을 개발하고, 안정성의 약한 조건을 통해 1차원 모델로 축소시켜 반응 장벽과 분기 효과가 발생하기 이르기 전의 역학적 장벽을 명확히 시각화한다.

ABSTRACT

In this paper, we apply Lagrangian descriptors to study the invariant manifolds that emerge from the top of two barriers existing in the LiCN<->LiNC isomerization reaction. We demonstrate that the integration times must be large enough compared with the characteristic stability exponents of the periodic orbit under study. The invariant manifolds manifest as singularities in the Lagrangian descriptors. Furthermore, we develop an equivalent potential energy surface with 2 degrees of freedom, which reproduces with a great accuracy previous results [Phys. Rev. E 99, 032221 (2019)]. This surface allows the use of an adiabatic approximation to develop a more simplified potential energy with solely 1 degree of freedom. The reduced dimensional model is still able to qualitatively describe the results observed with the original 2-degrees-of-freedom potential energy landscape. Likewise, it is also used to study in a more simple manner the influence on the Lagrangian descriptors of a bifurcation, where some of the previous invariant manifolds emerge, even before it takes place.

연구 동기 및 목표

• 리간그리안 디스크립터를 사용하여 LiCN⇌LiNC 이성화 반응에서 동역학적 기원을 가진 불변 다각형을 식별한다.
• 아비시노 결과를 정확하게 재현하는 모스 진동자 기반의 2자유도 포텐셜 에너지 표면(PES)을 구축한다.
• 고주파수 CN 진동을 분리하기 위해 안정성 근사법을 적용하여 2차원 시스템을 1차원 효과적 포텐셜로 축소시키며, 반응 장벽 근처의 주요 역학적 특징을 유지한다.
• 분기 발생 이전에 LD가 불변 다각형의 출현을 어떻게 드러내는지, 특히 근사법으로 인해 발생하는 오류적 구조를 포함하여 연구한다.
• 완전한 2차원 시스템에서 관측된 위상공간의 구조를 충분히 잘 반영하는지 검증하기 위해 축소된 1차원 모델의 능력을 평가한다.

제안 방법

• 유한한 통합 시간을 가진 리간그리안 디스크립터(LDs)를 사용하여 위상공간에서의 특이점으로서 불변 다각형을 탐지한다.
• LiCN 이성화를 모델링하기 위해 두 개의 자유도(R 및 ϑ)를 가진 2차원 모스 기반 PES를 구축하여 아비시노 PES와 수용 가능한 오차 범위 내에서 일치시킨다.
• 고주파수 CN 진동을 분리하기 위해 안정성 근사법을 적용하여 시스템을 1차원 효과적 포텐셜로 축소시킨다.
• 1차원 모델을 사용하여 분기의 영향, 특히 오류적 안장점과 관련된 다각형의 출현을 연구한다.
• 2차원 및 1차원 모델 간의 LD 패atters를 비교하여 축소된 모델이 역학적 장벽과 분리선을 얼마나 잘 반영하는지 검증한다.
• PES 상에서 고전적 궤적을 수치적으로 통합하고, 공식 $ M_{\tau}^{(q)}(x_0) = \int_{-\tau}^{\tau} \| \dot{\mathbf{x}}(t; x_0) \| dt $ 를 사용하여 LD를 계산한다. 여기서 $ \dot{\mathbf{x}} $ 는 궤적의 시간 도함수이다.

실험 결과

연구 질문

• RQ1 리간그리안 디스크립터는 LiCN 이성화 반응에서 장벽 정상에 위치한 주기 궤도와 관련된 불변 다각형을 신뢰성 있게 식별할 수 있는가?
• RQ2 이 시스템에서 통합 시간은 LD에서 불변 다각형의 해상도에 어떤 영향을 미치는가?
• RQ3 모스 기반 2차원 PES는 원래의 아비시노 PES의 역학을 어느 정도 정확하게 재현하는가?
• RQ4 안정성 근사법을 통해 얻은 1차원 모델은 불변 다각형과 역학적 장벽의 정성적 구조를 유지하는가?
• RQ5 축소된 1차원 모델에서 LD는 실제 분기가 발생하기 이르기 전에 분기의 역학적 전조를 어떻게 드러내는가?

주요 결과

• 리간그리안 디스크립터는 주기 궤도의 안정성 지수에 비례하여 충분히 긴 통합 시간이 확보된 경우, 위상공간에서 특이점으로서 불변 다각형을 성공적으로 식별한다.
• 2차원 모스 기반 PES는 일정한 관성 모멘트를 사용하더라도 아비시노 결과를 매우 정확하게 재현한다.
• 안정성 근사법을 통해 시스템의 유효한 1차원 축소가 가능하며, 반응 장벽 근처의 정성적 역학을 유지한다.
• 1차원 모델에서는 nR ≥ 7 인 경우 ϑ = π 라디안에서 허구의 안장점과 관련된 오류적 불변 다각형이 LD에서 나타나며, 이는 근사법의 산물이다.
• 이러한 오류적 구조는 실제 분기 발생 이전에도 LD에 나타나며, 이는 방법의 임박한 역학 전이에 대한 민감도를 보여준다.
• 축소된 1차원 모델은 2차원 시스템의 핵심 역학 현상, 즉 역학적 장벽과 분리선의 형성 등을 잘 반영하여 간소화된 역학 분석에의 활용 가능성을 검증한다.