| 查看: 2183 | 回复: 4 | ||||
[交流]
【分享】gaussian中得到波函数文件的方法在AIM进行下一步的计算 已有1人参与
|
|
用AIM计算分子有关性质，首先需要用gaussian得到其波函数文件，一般情况下对构建的3D分子需要进行优化。得到gaussian的输入文件有很多办法，可以手动输入，也可通过相应软件如gaussview 或chem3D等。 以丙稀C3H6为例。首先在Chem3D中得到其3D分子构型。以PM3进行预优化，创建gaussian输入文件如下: %Chk=c3h6.chk # RHF/6-31G(d,p) Opt c3h6 0 1 C 0 -1.107827 -0.000354 -0.596735 C 0 0.220172 -0.000354 -0.596735 C 0 1.034471 -0.000354 0.639204 H 0 -1.693653 0.000548 -1.510653 H 0 -1.703206 -0.000748 0.312312 H 0 0.782056 -0.000536 -1.538424 H 0 1.681046 -0.887459 0.675106 H 0 0.424156 0.002492 1.552111 H 0 1.684919 0.884001 0.672241 对此文件稍加修改，在＃行后加入 output＝wfn，在坐标最后加上存放的名字，如C3H6.wfn。默认的存放路径为gaussian安装路径下的Scratch文件夹。可将其改为所需的另外的文件夹,如F:\GaussOutput\C3H6.wfn。至此gaussian中输入文件变为 %Chk=c3h6.chk # RHF/6-31G(d,p) Opt output=wfn c3h6 0 1 C 0 -1.107827 -0.000354 -0.596735 C 0 0.220172 -0.000354 -0.596735 C 0 1.034471 -0.000354 0.639204 H 0 -1.693653 0.000548 -1.510653 H 0 -1.703206 -0.000748 0.312312 H 0 0.782056 -0.000536 -1.538424 H 0 1.681046 -0.887459 0.675106 H 0 0.424156 0.002492 1.552111 H 0 1.684919 0.884001 0.672241 F:\GaussOutput\C3H6.wfn 计算完毕，正常结束即可导入AIM2000中进行下一步的计算。 [ Last edited by yjcmwgk on 2009-11-26 at 09:30 ] |
回复此楼» 猜你喜欢
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
070300总分298求调剂
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有184人回复
-
一志愿厦门大学化学学硕307求调剂
已经有13人回复
-
理论与计算化学方向调剂——初试388分（数学144）
已经有1人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
调剂招生
已经有7人回复
-
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
硕士研究生调剂
已经有0人回复
-
投稿期刊求助
已经有1人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- SDS-PAGE样品浓缩的方法
已经有12人回复
- 急！ 关于Ms从Windows提交到Linux计算的问题？
已经有3人回复
- 求一种粉体中除铝粉的方法
已经有6人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- Gaussian计算和其他方法得到的热力学数据哪个更准确呢
已经有11人回复
- 求gaussian计算中理论方法的限制性和非限制性选择原则和差别
已经有3人回复
- gaussian 计算得到的荧光光谱选择其最大发射光谱是看f最大对应的吗？
已经有14人回复
- gaussian 激发态优化后的chk文件和荧光计算后的chk文件有何不同
已经有5人回复
- 【sobereva个人文集】自写Link生成Gaussian的IRC任务中每个点的波函数文件
已经有10人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【原创】并行自洽计算后得到和节点数一样多的波函数文件wfc
已经有4人回复
- 【讨论】含有配位键的金属原子在gaussian中如何计算单点能
已经有7人回复
- 【求助成功】supercell中DOS计算如何得到合理的电荷密度和波函数？
已经有3人回复
liyongxiu
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1874.2
- 散金: 335
- 红花: 1
- 帖子: 384
- 在线: 136.6小时
- 虫号: 442131
- 注册: 2007-10-27
- 专业: 化学动力学
2楼2009-04-11 13:58:28
lifei_dut
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1254.7
- 散金: 3
- 帖子: 180
- 在线: 48.5小时
- 虫号: 527008
- 注册: 2008-03-17
- 专业: 化学环境污染与健康
3楼2009-11-20 16:04:39
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +458 - QC强帖: 12
- 应助: 88 (初中生)
- 金币: 13047.7
- 散金: 27307
- 红花: 300
- 帖子: 4673
- 在线: 927.1小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 11-21 12:27
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 11-21 12:27
|
用gaussian产生wfn文件有几个问题需要注意： 1. 不能用ecp，不能用半经验方法。因为目前的aim程序不支持这些计算。 但是你可以直接用gaussian自带的aim做aim分析（例如 aim=all）。g98的手册里有aim关键词和选项，以后的版本不列这个关键词，但还可以用。gaussian自带的aim有两个问题：对有些元素不支持；用于半经验计算容易出错。 2. 基组不能包含g和g以上的基函数，否则出错。要么删除g以上的基函数，要么换程序（gamess支持g基函数）。此外还可以用checkpoint文件进行aim分析，aimall提供这样的工具。 3. 对于ROHF和RO-DFT计算类型，所有单占据轨道的占据数都是错的，应该是1.0却写成了2.0。需要自己修改wfn文件。 [ Last edited by beefly on 2009-11-23 at 06:39 ] |
4楼2009-11-21 00:25:51
xiaowandouer
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2019.8
- 散金: 90
- 红花: 2
- 帖子: 945
- 在线: 246.5小时
- 虫号: 446338
- 注册: 2007-10-29
- 专业: 理论和计算化学
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
先请教一个问题：用高斯算AIM，可信度怎么样，有虫友说高斯得到的结果可能会引起质疑。 另外，补充下楼主说的第一点，现在AIMALL应该是可以用ecp了。这是摘自aimall FAQs里面的有关ecp的回答（链接：http://aim.tkgristmill.com/faqs.html）。 Can wavefunctions derived from calculations using Effective Core Potentials (ECPs) be used? Yes, as long as the wavefunctions are AIM Extended Wavefunction files (.wfx files) that include additional core density function data for the ECP-modeled core electrons. More information can be found here. Note that Gaussian 09 B.01 and later support writing .wfx files and automatically writes [ Last edited by xiaowandouer on 2011-3-22 at 17:00 ] |
5楼2011-03-22 16:59:31
