| 查看: 1433 | 回复: 8 | ||||
[求助]
gaussian中那些算法不能应用维里定理
|
|
刚接触gaussian没多久， 用半经验法PM3优化分子时，所有分子优化收敛后， -V/T都是 接近于1。因为之前一直用密度泛函算都是接近2的，老师也说一定要接近2才算合理，但基本都没出现虚频。得出的能量算合理吗？还是我哪里算错了？求各位指点一下！！！！！还有就是分子总能量为应该可以为正值的吧？势能比动能大，总能量就一定是负值这是真的吗？ ****************************************** # opt pm3 geom=connectivity --------------------------- 1/14=-1,18=20,19=15,26=1,38=1,57=2/1,3; 2/9=110,12=2,17=6,18=5,40=1/2; 3/5=2,16=1,25=1,41=800000,71=1/1,2,3; 4/35=1/1; 5/5=2,35=1,38=5/2; 。。。。。。。。。。。 Requested convergence on RMS density matrix=1.00D-08 within 128 cycles. Requested convergence on MAX density matrix=1.00D-06. Requested convergence on energy=1.00D-06. No special actions if energy rises. Overlap will be assumed to be unity. Keep J ints in memory in canonical form, NReq=941156. SCF Done: E(RPM3) = -0.159845171082 A.U. after 10 cycles Convg = 0.4528D-08 -V/T = 0.9993 Calling FoFJK, ICntrl= 2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0. |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : YHPM3.docx
2015-03-16 10:01:50, 635.14 K
» 猜你喜欢
PDF在线压缩软件
已经有5人回复
-
材料/化学相关专业2026级学术/专业型硕士研究生
已经有19人回复
-
物理化学论文润色/翻译怎么收费?
已经有139人回复
-
求催化作用基础第三版 （甄开吉) 课件
已经有3人回复
-
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
304求调剂
已经有3人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2015-03-16 14:38:22
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13000.7
- 散金: 27307
- 红花: 300
- 帖子: 4672
- 在线: 926.9小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
3楼2015-03-16 20:31:40
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3510.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2015-03-17 11:15:53
5楼2015-03-17 20:11:01
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3510.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★ ★
gmy1990: 金币+2, 感谢应助 2015-03-17 23:06:13
gmy1990: 金币+2, 感谢应助 2015-03-17 23:06:13
|
能量算出来应该是没问题的，Gaussian给出的结果普遍应该是可靠的（不考虑方法本身的缺陷）。 “势能比动能大，总能量就一定是负值”从能量守恒的方程就很简单能看出，核对电子吸引项是负的，因为电子间的排斥通常相对动能及吸引势相比要小。 Virial定理2T=V要求势能必须服从1/r，如果势能偏离这个应该就不满足了，引进参数化就会导致这些问题。Xalpha如果alpha参数化之后可能也会偏离一些，其他的如Beefly前辈说的那样：半经验，ecp，相对论效应也不符合。 |
» 本帖已获得的红花（最新10朵）
黑木秋草6楼2015-03-17 20:38:10
7楼2015-03-18 10:25:19
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 应助: 47 (小学生)
- 金币: 3510.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
8楼2015-03-18 10:36:04
9楼2015-03-20 22:10:46
