xiaoxiong123123 发表于 2024-7-8 21:33 如果涉及到能量或者焓做差，那应该都要用一个程序。 |
卡开发发 发表于 2024-6-21 18:12 哥，如果我用dft＋u计算的话就是用castep计算，他这个模块的话计算的能量是不是就不一样了，还需要在家上▲H吗，是不是其他的物质也要用用这个模块计算了 |
卡开发发 发表于 2024-6-25 21:25 好，谢谢 |
本帖最后由 卡开发发 于 2024-6-25 21:30 编辑 xiaoxiong123123 发表于 2024-6-25 18:34 所以说之前一直在强调你的体系原子数目和化学反应式之间的计量关系。4个Ce8个O单位转换到mol，对应的是1mol Ce4O8或是说4mol CeO2，并不是1mol CeO2。额外需要说明的是之所以能这么算是因为这个晶体当中的Ce和Ce还有O和O是等价的。 |
卡开发发 发表于 2024-6-25 16:12 那这个转换成每mol是什么意思，建立4个Ce，8个O,的意思是1molCe4O8吗，代表的就不是CeO2了吗 |
xiaoxiong123123 发表于 2024-6-25 15:22 不是，晶胞翻一倍能量就翻一倍，你要折合需要自己去除以计量数，比如你的反应式当中写的是CeO2，那么计算Ce4O8的晶胞就得把晶胞能量除以4，而Ce8O16你就得除以8，要跟着你实际反应式去写。 |
卡开发发 发表于 2024-6-20 10:20 老师，他这个单位我还是不太明白，因为我如果做CeO2,模型就比如是4个Ce，8个O，优化后的能量比如是400HA,我扩胞一次变成8个Ce，16个O，优化后的能量成为800多HA，但是他转换单位之后不是成为每CeO2的mol了吗，这个单位转换没搞明白 |
卡开发发 发表于 2024-6-24 16:45 好的，谢谢！ |
xiaoxiong123123 发表于 2024-6-24 14:20 一样，只不过这里经典静电能指的是电子-电子平均排斥、核-电子吸引和核-核排斥，不同写法只是在能量分组或是是否写成扣除重复计入这种差别，这里仍然只是自洽场中得到的能量。 |
卡开发发 发表于 2024-6-21 18:12 好的，谢谢 |
xiaoxiong123123 发表于 2024-6-21 16:42 1、如果你研究的是高温，O2处于气态的（这也是我一开始所提醒的），在DMol3当中最简单的处理方法是直接画个O2的分子不要加盒子去处理，并且使用对称性，此时程序能够正确识别其为气态分子并进行平动和转动贡献的相关热化学计算。也可以加一个足够大的盒子，但平动和转动需要手工进行计算。PS：小心O2基态为三重态。 2、综合来说你的体系不适合用DMol3来研究，理由如下： （1）上面已经提示，DMol3处理晶体的声子频率的结果不一定是对的，虽然我有办法正确处理，但要么是计算量太大，要么需要适当写程序对新手不是非常友好。 （2）你的体系需要进行晶胞优化，DMol3那个stress的处理非常神秘，优化晶胞不太可靠。 （3）Ce这个体系大概率需要DFT+U或者杂化泛函才能算的比较好，DMol3也没办法支持。 综上，建议更换程序，如果有商业程序版权就用一些商业程序，没版权可以考虑Quantum Espresso。但不管更换哪个程序仍然需要注意O2的处理。 |
xiaoxiong123123 发表于 2024-6-21 11:31 1、算吸放热不需要做过渡态。恒温恒压情况下理想气体恒压热等于焓，理想晶体体胀不严重的话可以近似为热力学能，这些都是状态量。要在理想体系下通过谐振近似处理，需要做声子频率计算，但DMol3的声子频率只能做Gamma点的，需要特定处理，但DMol3自身不支持。 2、O2常温下是气相的，你可能需要注意你的研究条件，我不能确定你用晶体O2进行研究的原因。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 11:31 , Processed in 0.196794 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites