[交流] 用高斯计算垂直电离能和电子亲和能的方法

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2013-05-31 00:51:13

回帖支持 ( 显示支持度最高的前 50 名 )

beefly

专家顾问 (职业作家)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
luckone2013: 金币+2 2013-05-31 12:17:31
gmy1990: 金币+3 2013-06-02 22:07:39

第一种方法：通过给分子加、减电荷做多个独立的计算。
第二种方法：SAC-CI。不过，虽然SAC-CI是一种加上近似的EOM-CCSD方法，但是计算效率太低了。不知道是Gaussian的原因还是SAC-CI的原因。建议改用其它程序比如Q-Chem、CFour，GAMESS做EOM-IP-CCSD或EOM-EA-CCSD。

赞一下(1人)

3楼2013-05-31 04:27:27

普通回帖

gzl9901

铁杆木虫 (文坛精英)

应助: 3 (幼儿园)
金币: 6279.5
帖子: 11263
在线: 1769.2小时
虫号: 1294464

祝福，努力发论文！

2楼2013-05-31 01:19:13

luckone2013

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 495.5
帖子: 171
在线: 45.2小时
虫号: 2379784

内容已删除

4楼2013-05-31 12:07:37

kris4691

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 522.8
帖子: 57
在线: 24小时
虫号: 1348141

★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2013-06-02 22:08:01

所谓垂直就是不考虑结构变形。
1垂直电离能计算：优化分子构型，（同时也就得到分子能量了）使用这个构型减少电子计算离子的能量，离子减分子能量就是电离能了。
2垂直电子亲和能：和电离相反的过程，先优化离子的构型，再以此构型直接计算分子的能量，再求能量差。

我没做过，但是听老师讲过，可能有不正确的地方，楼主可以再查一下。

5楼2013-05-31 16:35:47

luckone2013

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 495.5
帖子: 171
在线: 45.2小时
虫号: 2379784

引用回帖:

5楼: Originally posted by kris4691 at 2013-05-31 16:35:47
所谓垂直就是不考虑结构变形。
1垂直电离能计算：优化分子构型，（同时也就得到分子能量了）使用这个构型减少电子计算离子的能量，离子减分子能量就是电离能了。
2垂直电子亲和能：和电离相反的过程，先优化离子的 ...

那意思就是垂直的要保持构型不变，求分子和离子的能量差值吧？
还有最关键的是优化好构型，减少电子计算离子的能量时，怎么让分子的构型保持不变？

[ Last edited by luckone2013 on 2013-5-31 at 16:49 ]

6楼2013-05-31 16:44:13

kris4691

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 522.8
帖子: 57
在线: 24小时
虫号: 1348141

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
luckone2013: 金币+5, 谢谢！ 2013-05-31 23:17:38

引用回帖:

6楼: Originally posted by luckone2013 at 2013-05-31 16:44:13
那意思就是垂直的要保持构型不变，求分子和离子的能量差值吧？
还有最关键的是优化好构型，减少电子计算离子的能量时，怎么让分子的构型保持不变？
...

区别在于，电离能是用分子的构型计算能力差；亲和能是用离子构型。
我很久没用高斯了，大概讲一下，可能有错，你最好去看看高斯手册。
把优化结构的构型取出来，就是用gaussview打开优化的log文件再另存为gjf，在calculate》calculation setup里面设置电荷，回车以后他会自动修改自旋。general里面沟选ignore symmetry。就可以了。
方法和基组按老师要求选吧。
例如
%chk=c6.chk
%mem=1800MB
%nproc=12
# b3lyp/gen nosymm

Title Card Required

1 2

7楼2013-05-31 18:07:37

luckone2013

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 495.5
帖子: 171
在线: 45.2小时
虫号: 2379784

引用回帖:

7楼: Originally posted by kris4691 at 2013-05-31 18:07:37
区别在于，电离能是用分子的构型计算能力差；亲和能是用离子构型。
我很久没用高斯了，大概讲一下，可能有错，你最好去看看高斯手册。
把优化结构的构型取出来，就是用gaussview打开优化的log文件再另存为gjf，在 ...

是不是在优化好的构型里，在选用的基组后面加上nosymm，然后将电荷和多重度改掉就可以啦？

9楼2013-06-01 17:27:49

kris4691

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 522.8
帖子: 57
在线: 24小时
虫号: 1348141

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
luckone2013: 金币+1 2013-06-02 23:14:38

引用回帖:

9楼: Originally posted by luckone2013 at 2013-06-01 17:27:49
是不是在优化好的构型里，在选用的基组后面加上nosymm，然后将电荷和多重度改掉就可以啦？...

是的，注意自旋多重度千万别弄错了。

10楼2013-06-02 15:28:26

luckone2013

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 495.5
帖子: 171
在线: 45.2小时
虫号: 2379784

引用回帖:

10楼: Originally posted by kris4691 at 2013-06-02 15:28:26
是的，注意自旋多重度千万别弄错了。...

但是，命令行加了nosymm后构型还是从阴离子型变成了中性分子的构型。感觉nosymm不能保持分子的构型，而只是和分子的对称性有关的项。是不是应该把命令行的 opt 删掉？

11楼2013-06-02 23:26:53

kris4691

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 522.8
帖子: 57
在线: 24小时
虫号: 1348141

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
luckone2013: 金币+1, 这个真的没人教，谢谢 2013-06-03 12:27:36

引用回帖:

11楼: Originally posted by luckone2013 at 2013-06-02 23:26:53
但是，命令行加了nosymm后构型还是从阴离子型变成了中性分子的构型。感觉nosymm不能保持分子的构型，而只是和分子的对称性有关的项。是不是应该把命令行的 opt 删掉？...

我上面不是给你写了个例子你还不明白，看来你是对高斯软件没有基础。多看看高斯手册和教程，网上很好找的。
opt命是对结构进行优化，肯定要改变构型，而你说会优化成中性分子的构型，可能是你设置哪儿又出什么问题了，也可能离子和分子的构型差不多，你仔细看看原子坐标。这些都是很基础的东西，没有人教你吗？先把基础操作和理论知识弄扎实了再去做计算，你现在这样现学现用很容易出错。

12楼2013-06-03 10:43:30