| 查看: 2056 | 回复: 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 ] |
回复此楼» 猜你喜欢
理论计算合作
已经有4人回复
-
小木虫的论文辅导靠谱吗？有没有用过的同学？
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有221人回复
-
QE计算声子ph.out的Γ点出现虚频（-30cm-1）vasp计算没有
已经有0人回复
-
40-50万/年，中国散裂中子源诚聘计算模拟方向博士后
已经有81人回复
-
比利时鲁汶大学与国家留学基金委共同资助博士研究生CSC-KU Leuven PhD Scholarship
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
书籍求助：汽车市场营销理论与实务（电子版）——章小平
已经有0人回复
-
东方理-中科大联合博士生招聘
已经有0人回复
-
伦敦大学学院Benjamin Moss博士招收人工智能、原位光谱、催化相关26年博士生
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - QC强帖: 12
- 应助: 86 (初中生)
- 金币: 12784.9
- 散金: 27307
- 红花: 300
- 帖子: 4664
- 在线: 925.5小时
- 虫号: 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
lifei_dut
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1254.7
- 散金: 3
- 帖子: 180
- 在线: 48.5小时
- 虫号: 527008
- 注册: 2008-03-17
- 专业: 化学环境污染与健康
3楼2009-11-20 16:04:39
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