[論文レビュー] Identification of the invariant manifolds of the LiCN molecule using Lagrangian descriptors

F. Revuelta, R. M. Benito+1 more|arXiv (Cornell University)|Oct 19, 2021

Advanced Chemical Physics Studies参考文献 57被引用数 8

ひとこと要約

本研究では、Lagrangian記述子（LDs）を用いてLiCN⇌LiNCの異性化反応における不変多様体を同定し、LDにおける特異点として多様体を解像するには長い積分時間が必要であることを示した。2次元のMorse型ポテンシャルエネルギー表面（PES）を構築し、ab initio結果を高精度で再現するとともに、断熱近似を用いて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.

研究の動機と目的

• Lagrangian記述子を用いて、LiCN⇌LiNC異性化反応における力学的起源の不変多様体を同定すること。
• ab initio結果を高精度で再現できる2自由度のMorse振動子に基づく簡略化された2次元PESを構築すること。
• 高周波数のCN伸張振動モードを分離する断熱近似を適用し、反応障壁付近の主要な力学的特徴を保持したまま、2次元系を1次元有効ポテンシャルに還元すること。
• 分岐が発生する直前における不変多様体の出現をLDがどのように明らかにするかを調査し、近似によって生じる誤った構造についても検討すること。
• 簡略化された1次元モデルが、元の2次元系で観察された位相空間構造を定性的に捉えられるかを検証すること。

提案手法

• 有限の積分時間でLagrangian記述子（LDs）を用い、位相空間における特異点として不変多様体を検出する。
• LiCN異性化をモデル化する2自由度（Rおよびϑ）のMorse型PESを構築し、許容誤差範囲内でab initio PESと一致させる。
• 高周波数のCN伸張振動を分離する断熱近似を適用し、1自由度の有効ポテンシャルに系を簡略化する。
• 1次元モデルを用いて、分岐がLD構造に与える影響を調査し、特に架空の不安定点に起因する誤ったサドル点および関連する多様体の出現を分析する。
• 2次元および1次元モデル間のLDパターンを比較し、簡略化モデルが力学的障壁および分離子を的確に再現できるかを検証する。
• PES上での古典的軌道の数値積分と、以下の式を用いたLDの計算を実施：$ M_{\tau}^{(q)}(x_0) = \int_{-\tau}^{\tau} \| \dot{\mathbf{x}}(t; x_0) \| dt $、ここで$ \dot{\mathbf{x}} $は軌道の時間微分を表す。

実験結果

リサーチクエスチョン

• RQ1 Lagrangian記述子は、LiCN異性化反応における障壁上部の周期的軌道に関連する不変多様体を信頼性高く同定できるか？
• RQ2 積分時間の長さが、この系におけるLDにおける不変多様体の解像度にどのように影響するか？
• RQ3 Morse型2次元PESは、元のab initio PESの力学的挙動をどの程度正確に再現できるか？
• RQ4 断熱近似により得られる1次元モデルは、不変多様体および力学的障壁の定性的な構造を保持できるか？
• RQ5 簡略化された1次元モデルにおけるLDは、実際に分岐が発生する前段階で、分岐の動的予兆をどの程度明らかにするか？

主な発見

• Lagrangian記述子は、周期的軌道の安定性指数に対して十分に長い積分時間をとれば、位相空間における特異点として不変多様体を的確に同定できる。
• 2次元のMorse型PESは、定数の慣性モーメントを用いても、ab initio結果を高い精度で再現できる。
• 断熱近似により、反応障壁付近の定性的な力学的挙動を保持した有効な1次元系への還元が可能である。
• 1次元モデルにおいて、nR ≥ 7のとき、ϑ = π radに架空の不安定点に起因する誤った不変多様体がLDに現れるが、これは近似による人工的構造である。
• これらの誤った構造は、実際の分岐が発生する前からLDに現れるため、本手法が近い将来の力学的遷移に敏感であることを示している。
• 簡略化された1次元モデルは、2次元系の本質的現象を捉えており、力学的障壁や分離子の形成を含む。これにより、簡略化された動的解析に有効であることが検証された。