Høyer Nicolai Machholdt, Christiansen Ove
Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark.
J Chem Theory Comput. 2024 Jan 23;20(2):558-579. doi: 10.1021/acs.jctc.3c00962. Epub 2024 Jan 6.
We present a new quasi-direct quantum molecular dynamics computational method which offers a compromise between quantum dynamics using a precomputed potential energy surface (PES) and fully direct quantum dynamics. This method is termed the time-dependent adaptive density-guided approach (TD-ADGA) and is a method for constructing a PES on the fly during a dynamics simulation. This is achieved by acquisition of new single-point (SP) calculations and refitting of the PES, depending on the need of the dynamics. The TD-ADGA is a further development of the adaptive density-guided approach (ADGA) for PES construction where the placement of SPs is guided by the density of the nuclear wave function. In TD-ADGA, the ADGA framework has been integrated into the time propagation of the time-dependent nuclear wave function and we use the reduced one-mode density of this wave function to guide when and where new SPs are placed. The PES is thus extended or updated if the wave function moves into new areas or if a certain area becomes more important. Here, we derive equations for the reduced one-mode density for the time-dependent Hartree (TDH) method and for multiconfiguration time-dependent Hartree (MCTDH) methods, but the TD-ADGA can be used with any time-dependent wave function method as long as a density is available. The TD-ADGA method has been investigated on molecular systems containing single- and double-minimum potentials and on single-mode and multi-mode systems. We explore different approaches to handle the fact that the TD-ADGA involves a PES that changes during the computation and show how results can be obtained that are in very good agreement with results obtained by using an accurate reference PES. Dynamics with TD-ADGA is essentially a black box procedure, where only the initialization of the system and how to compute SPs must be provided. The TD-ADGA thus makes it easier to carry out quantum molecular dynamics and the quasi-direct framework opens up the possibility to compute quantum dynamics accurately for larger molecular systems.
我们提出了一种新的准直接量子分子动力学计算方法,该方法在使用预先计算的势能面(PES)的量子动力学和完全直接量子动力学之间提供了一种折衷方案。这种方法被称为时间相关自适应密度引导方法(TD-ADGA),是一种在动力学模拟过程中即时构建PES的方法。这是通过根据动力学的需要获取新的单点(SP)计算并重新拟合PES来实现的。TD-ADGA是用于PES构建的自适应密度引导方法(ADGA)的进一步发展,其中SP的放置由核波函数的密度引导。在TD-ADGA中,ADGA框架已被集成到时间相关核波函数的时间传播中,并且我们使用该波函数的约化单模密度来指导新SP的放置时间和位置。因此,如果波函数移动到新区域或某个区域变得更加重要,则PES会被扩展或更新。在这里,我们推导了时间相关哈特里(TDH)方法和多组态时间相关哈特里(MCTDH)方法的约化单模密度的方程,但只要有密度可用,TD-ADGA就可以与任何时间相关波函数方法一起使用。TD-ADGA方法已在包含单势阱和双势阱的分子系统以及单模和多模系统上进行了研究。我们探索了不同的方法来处理TD-ADGA涉及在计算过程中变化的PES这一事实,并展示了如何获得与使用精确参考PES获得的结果非常吻合的结果。使用TD-ADGA的动力学本质上是一个黑箱过程,其中只需要提供系统的初始化以及如何计算SP。因此,TD-ADGA使得进行量子分子动力学更容易,并且准直接框架为准确计算更大分子系统的量子动力学开辟了可能性。