| 查看: 1713 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
dftb中默认温度是0K，计算单点能怎么设置温度呢？
|
||
dftb中默认温度是0K，要计算300K下的单点能，怎么设置温度呢？输入文件见附件，输出结果显示，Driver 里的Thermostat有问题，这该怎么解决呢？ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : dftb_in.hsd
2013-08-22 21:04:43, 2.83 K
» 猜你喜欢
国家级人才课题组招收2026年入学博士
已经有2人回复
-
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有140人回复
-
散金
已经有87人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有34人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有1人回复
-
KAUST（阿卜杜拉国王科技大学）MXene 器件方向博士后招聘
已经有0人回复
-
沙特阿拉伯阿卜杜拉国王科技大学（KAUST）电池方向博士后招聘
已经有0人回复
-
福州大学新能源材料与工程研究院招收2026年入学博士
已经有0人回复
-
南通大学电气与自动化学院2026年优秀博士招聘公告
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】场及量子化
已经有13人回复
送红花一朵 
3楼2013-08-27 13:08:38
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