PIMD轨迹分析分析质心轨迹还是珠子轨迹需要逐个分析？

Assiduitas

我参考文献《Quantum fluctuations and isotope effects in ab initio descriptions of water》DOI：10.1063/1.4894287进行AIMD模拟。想复现一下文中所提及的“质子漂移”现象，即对于氢键，质子氢原子可能会存在距离施体氧、受体氧距离一样的情况，甚至是距离受体氧的距离小于施体氧。
对于这个问题，想请教一下各位做过PIMD的大佬，分析动力学轨迹时，是仅仅分析质心轨迹就能看到感兴趣的量子效应，还是需要分子每个珠子的轨迹，才应该能看到上述所阐述的量子效应？
就是图中会存在横坐标>0的情况，即V>0

sobereva

仅仅分析质心轨迹

Assiduitas

sobereva 发表于 2025-6-26 05:29
仅仅分析质心轨迹

sob老师，我把质心轨迹和珠子轨迹都分析了一遍。发现：对于H而言，珠子轨迹存在比较多的“质子漂移”现象。而D的珠子轨迹基本没有“质子漂移”现象。这符合常识——H的量子效应要强于D。
而分析H与D的质心轨迹，两者都不存在“质子漂移”现象。我能不能用珠子的轨迹来解释这个现象呢？
模型都是32个水分子，不同点在于一个是H2O，一个是D2O。

Assiduitas

sobereva 发表于 2025-6-26 05:29
仅仅分析质心轨迹

sob老师，如果我研究的这个“质子漂移”是量子效应主导，我是不是应该研究珠子轨迹。而对于质心运动轨迹，我应该把它当作经典动力学轨迹看待。

sobereva

Assiduitas 发表于 2025-6-26 09:51
sob老师，如果我研究的这个“质子漂移”是量子效应主导，我是不是应该研究珠子轨迹。而对于质心运动轨迹 ...

是的。珠子轨迹能体现核波包分布，与核量子效应关系最直接。
但质心运动轨迹也不能说就是不考虑量子效应的经典动力学轨迹，而应该说成是考虑了核量子效应情况下的核波包中心的轨迹。并不总是波包质心轨迹都和经典动力学轨迹相仿佛。