24小时热门版块排行榜    

查看: 256  |  回复: 1
当前主题已经存档。

[资源] 血的教训_关于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" 之后,整个分子的能量大大降低,但这种降低并非真实情况。所以我们在做有关于过渡金属的时候,一定要加这个关键字。

   现学现卖,说的比较粗俗,希望大家能够看懂。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yhjiangshiyu

木虫 (正式写手)


这个我刚做MS的时候也遇到过,全电子和价电子,呵呵
2楼2007-11-21 16:54:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangwei_suda 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +3 满头大汗的学生 2026-02-28 3/150 2026-02-28 16:22 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 312求调剂 +4 吃宵夜1 2026-02-28 5/250 2026-02-28 16:11 by gjm133
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[考研] 0856材料求调剂 +7 hyf hyf hyf 2026-02-28 8/400 2026-02-28 12:36 by 52hz~~
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 311求调剂 +3 南迦720 2026-02-28 3/150 2026-02-28 12:24 by houyaoxu
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
[考研] 280求调剂 +6 Qq206./ 2026-02-21 6/300 2026-02-28 11:20 by 我!要一战成硕
[考博] 博士推荐 +3 花儿笑? 2026-02-21 4/200 2026-02-27 18:30 by 黑!在干嘛
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 10/500 2026-02-27 13:22 by ichall
[基金申请] 什么是人一生最重要的? +10 瞬息宇宙 2026-02-21 10/500 2026-02-27 08:46 by tfang
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
[教师之家] 版面费该交吗 +14 苹果在哪里 2026-02-22 17/850 2026-02-26 11:55 by AGI智创机器人
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见