| 查看: 1834 | 回复: 2 | |||
[求助]
dftb中默认温度是0K，计算单点能怎么设置温度呢？
|
dftb中默认温度是0K，要计算300K下的单点能，怎么设置温度呢？输入文件见附件，输出结果显示，Driver 里的Thermostat有问题，这该怎么解决呢？ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : dftb_in.hsd
2013-08-22 21:04:43, 2.83 K
» 猜你喜欢
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
求助各位大佬，球球了，在这科萨头
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有73人回复
-
七嗪类物质合成求助
已经有0人回复
-
分子动力学模拟合作需求
已经有4人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
湖北文理学院材料与化工专硕还有调剂名额，欢迎大家推荐考生！
已经有1人回复
-
微纳hub，微纳加工的宝库
已经有0人回复
-
微纳加工里的 lift-off 工艺
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】场及量子化
已经有13人回复
ChemiAndy
木虫 (正式写手)
- QC强帖: 3
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★
gmy1990: 金币+5 2013-08-23 03:25:07
毛毛虫_: 回帖置顶 2013-08-27 22:32:41
gmy1990: 金币+5 2013-08-23 03:25:07
毛毛虫_: 回帖置顶 2013-08-27 22:32:41
|
一般所说的单点能，是指0K下的势能，包括了“原子核，电子”体系之间左右的相互作用能，以及电子动能(电子动能通常很小)。所以计算单点能不需要指定温度，你的输出文件中去掉Driver那一项即可得到单点能。 0K下的分子与300K下的分子仅差一个原子核动能，每个原子在每个方向具有平均动能1/2*kT，所以总动能即为所有N个原子在x,y,z三个方向动能之和为：3*N*(1/2*kT)。这是能均分定理。你就是用MD跑300K下的动力学模拟，其结果也是如此。 另外，我注意到你的体系里有H原子，那么你应该在Hamiltonian里添加如下2项： DampXH = Yes DampXHExponent = 4.53 这两项是指定进行H键校正。这对SCC-DFTB2的能量计算非常重要。 请参考manual和下面文献，其中在方法名称中有“h”的，就是加了这一项氢键校正： M. Gaus, Q. Cui, and M. Elstner. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB). J. of Chem. Theory Comput., 7:931–948,2011. 19, 20 此外，建议你做计算使用SCC-DFTB3方法，否则发文章的话，审稿人会问你，明明DFTB3被证明比DFTB2准确很多，你为什么不用DFTB3。 指定DFTB3是在Hamiltonian中加如下关键词： ThirdOrderFull = Yes HubbardDerivs = { C = -0.1492 O = -0.1575 H = -0.1857 } 这些参数都是上面Cui Qiang那篇文献中报道的。 |
» 本帖已获得的红花（最新10朵）
毛毛虫_2楼2013-08-23 01:01:38
3楼2013-08-27 13:08:38
