计算化学公社

标题: 求助：激光模拟的温度，步长等参数设置为多少合适 [打印本页]

作者
Author: laser-z 时间: 2024-5-23 16:57
标题: 求助：激光模拟的温度，步长等参数设置为多少合适
求助：我在做飞秒激光方面的第一性原理分子动力学模拟，现在在复现一篇文章，在此文章中，整个体系是在电子温度T=4ev下进行模拟，4ev转化为开氏温度为46420K，这温度会不会很高？我看到培训班例子中最高的为3500K（"高温下烷烃裂解的模拟"）。另外由于温度过高，电子速度会变得很快，所以步长是否应设置的极小，达到10的-2次方或-3次方？此外热浴时间常数TIMECONE的设置，50fs是否足够，是否要减小到10fs的层次？

感谢各位大佬指点

作者
Author: sobereva 时间: 2024-5-23 17:25
看那篇文章具体用的模拟方法说明
当前情况不属于一般的AIMD模拟范畴，CP2K本身都未必适合

作者
Author: 乐平 时间: 2024-5-23 17:35
飞秒激光作用在体系上，需要用 real-time TDDFT 来模拟，CP2K 没有这个功能

作者
Author: laser-z 时间: 2024-5-23 18:28

乐平发表于 2024-5-23 17:35
飞秒激光作用在体系上，需要用 real-time TDDFT 来模拟，CP2K 没有这个功能

老师，那篇文章是使用CPMD方法进行模拟的，将飞秒激光的作用改变为电子获得假想动能，才有4ev的电子温度。在此电子温度下，使用NVE系综，固定电子占据数，进行CPMD模拟，文章使用的软件是Quantum Espresso软件包。如果我想计算飞秒激光相关，是否就要更换其他量子化学软件才可以？

作者
Author: laser-z 时间: 2024-5-23 18:40

sobereva 发表于 2024-5-23 17:25
看那篇文章具体用的模拟方法说明
当前情况不属于一般的AIMD模拟范畴，CP2K本身都未必适合

老师，那篇文章使用的是CPMD方法模拟，确实跟AIMD是两种不同的模拟方法。文章中将飞秒激光的作用设定为电子在固定原子位置的场中传播50 fs，从而获得假想动能并平衡配置和动能温度，然后在此状态下，使用NVE系综，固定电子占据数，进行CPMD模拟。文章使用软件为Quantum Espresso。
这样的话，是否就无法使用CP2K做激光与物质在长时间尺度下相互作用的模拟？

作者
Author: sobereva 时间: 2024-5-24 04:30

laser-z 发表于 2024-5-23 18:40
老师，那篇文章使用的是CPMD方法模拟，确实跟AIMD是两种不同的模拟方法。文章中将飞秒激光的作用设定为电 ...

CPMD属于AIMD的一种，和BOMD、ADMP并列，只不过CPMD明确考虑了电子波函数的运动，而CP2K做BOMD不考虑，而是每一步都得到完全自恰的电子波函数
CP2K不支持CPMD，不适合干这个。如果要重现，可以用QE或者目前已经不再继续开发的以前比较流行的CPMD程序

欢迎光临计算化学公社 (http://bbs.keinsci.com/)