| 查看: 1748 | 回复: 2 | ||
[求助]
dftb中默认温度是0K，计算单点能怎么设置温度呢？
|
dftb中默认温度是0K，要计算300K下的单点能，怎么设置温度呢？输入文件见附件，输出结果显示，Driver 里的Thermostat有问题，这该怎么解决呢？ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : dftb_in.hsd
2013-08-22 21:04:43, 2.83 K
» 猜你喜欢
1
已经有21人回复
-
0703 调剂
已经有6人回复
-
物理化学论文润色/翻译怎么收费?
已经有56人回复
-
西南某211招27届博士
已经有33人回复
-
大连工业杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有0人回复
-
深圳大学高等研究院赵伟课题组招收化学调剂生
已经有3人回复
-
青岛大学固态电池研究团队招收材料、化学、物理等方向博士研究生
已经有31人回复
-
化学工程085602 305分求调剂
已经有17人回复
-
武汉纺织大学全国重点实验室杨应奎教授团队拟招17-19位硕博生！
已经有0人回复
-
材料与化工专硕接受跨专业调剂（重庆市属高校）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】场及量子化
已经有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那篇文献中报道的。 |
毛毛虫_2楼2013-08-23 01:01:38
3楼2013-08-27 13:08:38
