前面就说了电脑很傻逼的,所以对于t的选择只能是离散的点 那么他们之间的时间间隔怎么选取,这个就是步长 为了保障节省计算的时间而不失去其精准性 通常的原则为积分步程应小于系统中最快运动周期的1/ 10 由于分子间的运动是由范德华力影响所以较慢 通常振动频率最高的键是与氢相关的键,周期约10 fs (~3000cm-1) *上述的频率与波数转换,3000cm-1为波数,光速为3x108m/s,光速*波数=频率 ∴频率=3x108m/sx3000cm-1=3x10^10cm/sx3000cm-1=9x10^13/s 频率*周期=1 ∴频率=1/周期=1/(9x10^13/s)=1.1x10^14s≈10fs 设定的步长太长了,就没办法判断正确的移动轨迹 步长短了,就会加大计算工作量 最重要的是时间步长越小体系跨过势垒的能力越差(局部能量最小) 你步长调小的模拟的差异基本可以归结于时间平均不足以达到系综平均, 并不是步长调小了,结果就正确了。 |
Powered by Discuz! X3.5
© 2001-2024 Discuz! Team.