| 查看: 1598 | 回复: 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 ] |
回复此楼» 猜你喜欢
中国科学院成都山地灾害与环境研究所2026年特别研究助理（博士后）招聘启事
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有195人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
研究生招生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有6人回复
-
招生材料与化工、环境科学调剂研究生
已经有5人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
-
中国电机工程学报的投稿系统
已经有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
回复
