| 查看: 1712 | 回复: 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级博士生（与中山大学物理学院联合培养）
已经有22人回复
-
PbS量子点紫外吸收
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有203人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有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