24小时热门版块排行榜

返回列表

当前主题已经存档。

zhangwei_suda

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
贵宾: 0.25
金币: 8369.6
帖子: 728
在线: 25.9小时
虫号: 99273

[资源] 血的教训_关于pseudo=read

不知道咱们这里是怎么做过渡金属的，我是刚入手，所以犯了很多低级错误，在这里说一下，大家随便BS，关键是防着后面有人犯同样的错误，hoho

我现在做的是Pt的体系，cis-Cl2-Pt-R2，两个Cl原子处于cis-位置，但是我刚开始的时候，画分子时很粗心的把cis画成了trans结构，开始尝试做计算。。。。然后，当天晚些时候，重新检查分子结构，发现图画错了，于是改正，同时同组里面的同学告诉我做过渡金属，一般需要在路径部分输入“pseudo=read", 问他们是什么作用，他们也说的不是很清楚，反正有用，我就加了。

  但是今天无聊的时候用"SCF D",检查两个log文件发现，能量相差了

2000 hartree

  !!!!

后来请教老师，发现输出部分有如下不同：
未加入”pseudo=read"时， 178 alpha electrons    178 beta electrons
加入 "pseudo=read" 时， 128 alpha electrons    128 beta electrons

也就是在用了"pseudo=read"关键字以后，电子数少了整整100个。后来老师帮我分析了这个原因：pseudo=read, 是一个对赝势场的定义（具体定义什么的详见gaussian手册），我只说直观的，加入这个关键字以后，饱和壳层的原子不再纳入考虑，把它们当做是核的一部分处理。而不加这个关键字，计算的是全电子结构，所有的电子波函数都参与一次计算。简单的说，对于我做的体系，含有一个Pt原子，两个P原子，两个Cl原子。对于Pt而言，不加pseudo=read时，有78个电子，加时是18个电子（只算外层的，价电子)，Cl，不加时17电子，加7电子；P分别是15和5，刚好(60+10*2+10*2)=100个电子。

总而言之，在加 "pseudo=read" 之后，整个分子的能量大大降低，但这种降低并非真实情况。所以我们在做有关于过渡金属的时候，一定要加这个关键字。

现学现卖，说的比较粗俗，希望大家能够看懂。

回复此楼