| 查看: 1386 | 回复: 11 | ||
| 【奖励】 本帖被评价7次，作者beefly增加金币 7 个 | ||
| 当前主题已经存档。 | ||
[资源]
【beefly个人文集】MOLBO——MOLPRO和NBO的接口程序（2.1版）
|
||
|
MOLPRO是当今计算速度最快的从头量子化学程序之一，但是其NBO模块的开发过于缓慢，功能也极其有限，故自己动手写了MOLBO这个程序。MOLBO从MOLPRO输出文件提取有用的信息，构造用于NBO分析的*.47文件。 MOLBO 2.1版能够提取NBO程序所需的所有数据。MOLBO 2.1版能自动进行Mayer键级分析。 1. 用法： 第一步，在MOLPRO计算中打印基组、重叠矩阵、密度矩阵、Fock矩阵等信息（参见输入模板）； 第二步，鼠标双击MOLBO程序，指定要进行转化的MOLPRO输出文件（扩展名.out或*.output可省略）。如果正常，MOLBO产生*.47文件。 第三步，用GENNBO程序打开*.47文件，进行NBO分析。GENNBO 4.0以上版本是商业软件，需要购买。GENNBO 3.0是免费的，可以从CCL网站下载源代码自己编译： http://www.ccl.net/cca/software/SOURCES/FORTRAN/nbo/index.shtml 也可以从这里下载编译好的GENNBO 3.0 for Windows http://cqb.fc.ul.pt/intheochem/nuno/index.html GENNBO 3.0的$CONTRACT数据区不支持g函数。在基组含有g函数的情况下，为了产生正确的.31文件用于绘图，MOLBO可以自动产生.31文件，用来替换GENNBO 3.0产生的.31文件。 2. 限制和可能存在的问题（一定要读完这一部分！！！） 对称性。对于使用了对称性的计算，若存在对称等价原子，MOLPRO会把重叠矩阵做分块对角化处理，导致对角元不为1.0，无法用于NBO分析。对于这种情况，需要降低分子对称性直到所有原子不等价，并重新运行MOLPRO计算。例如，中心对称的线型分子需要用C2v群，而不能用D2h群；水分子要用Cs群而不能用C2v。对于各种偶极矩分析（$DIPOLE），不能使用对称性。 自定义NBO选项。*.47文件默认加入了键级分析选项（BNDIDX NLMO），还可以选择是否产生绘图文件。也可以自己加入其他选项，请仔细阅读NBO程序的使用手册。 重新开始计算。不能使用1号文件，否则在MOLPRO输出文件中找不到基组信息。 3. 测试过的MOLPRO版本： 2006，2008，2009 4. 测试过的计算类型： SCF: RHF，ROHF，UHF，RDFT，RODFT，UDFT Density Fitting: DF-HF, DF-KS (nosym) Post-HF: MP2, MP3, CCSD, QCISD, QCISD(T), EOM-CCSD, Full CI Local Post-HF: LMP2 (nosym) MC/MR: MCSCF, CASVB, MRCI (including SRCI), CASPT2 (including SRMP2) 参见输入模板的例子，2006和2008格式。若用于2009，部分输入文件需要修改坐标部分。产生$DIPOLE的例子见rhf-dipole.inp。 5. 最后升级： 2010.1.26，2.1版。自动计算Mayer键级。 下载（此链接3个月内有效）： https://webspace.utexas.edu/wz838/molbo.zip [ Last edited by yjcmwgk on 2010-6-16 at 21:36 ] |
回复此楼» 猜你喜欢
085404（计算机技术）293分 希望调剂到材料相关专业
已经有2人回复
-
085404（计算机技术）293分 希望调剂到材料相关专业
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有92人回复
-
085404（计算机技术）293分 希望调剂到材料相关专业
已经有0人回复
-
ELISA试验中不容忽视的细节盘点(二)
已经有0人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
七嗪类物质合成求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
4楼2010-01-12 15:25:17
5楼2010-01-21 20:30:13
6楼2010-01-22 00:55:40
7楼2010-01-22 10:33:05
8楼2010-01-26 16:18:34
9楼2010-01-27 01:33:32
10楼2010-01-27 09:22:55
11楼2010-01-27 09:57:52
12楼2010-01-28 14:51:54
简单回复
wgpig2楼
2009-12-31 14:50
回复
chrinide3楼
2009-12-31 15:47
回复
