Gabas Fabio, Conte Riccardo, Ceotto Michele
Dipartimento di Chimica, Università degli Studi di Milano , via Golgi 19, 20133 Milano, Italy.
J Chem Theory Comput. 2017 Jun 13;13(6):2378-2388. doi: 10.1021/acs.jctc.6b01018. Epub 2017 May 26.
We present an on-the-fly ab initio semiclassical study of vibrational energy levels of glycine, calculated by Fourier transform of the wavepacket correlation function. It is based on a multiple coherent states approach integrated with monodromy matrix regularization for chaotic dynamics. All four lowest-energy glycine conformers are investigated by means of single-trajectory semiclassical spectra obtained upon classical evolution of on-the-fly trajectories with harmonic zero-point energy. For the most stable conformer I, direct dynamics trajectories are also run for each vibrational mode with energy equal to the first harmonic excitation. An analysis of trajectories evolved up to 50 000 atomic time units demonstrates that, in this time span, conformers II and III can be considered as isolated species, while conformers I and IV show a pretty facile interconversion. Therefore, previous perturbative studies based on the assumption of isolated conformers are often reliable but might be not completely appropriate in the case of conformer IV and conformer I for which interconversion occurs promptly.
我们展示了一项对甘氨酸振动能级的即时从头算半经典研究,该研究通过波包相关函数的傅里叶变换进行计算。它基于一种多相干态方法,并结合了用于混沌动力学的单值矩阵正则化。通过对具有谐波零点能量的即时轨迹进行经典演化所获得的单轨迹半经典光谱,对所有四个最低能量的甘氨酸构象异构体进行了研究。对于最稳定的构象异构体I,还针对每个能量等于第一谐波激发的振动模式运行了直接动力学轨迹。对演化至50000个原子时间单位的轨迹进行分析表明,在此时间跨度内,构象异构体II和III可被视为孤立的物种,而构象异构体I和IV则表现出相当容易的相互转化。因此,先前基于孤立构象异构体假设的微扰研究通常是可靠的,但对于构象异构体IV和I,由于它们会迅速发生相互转化,可能并不完全适用。