计算化学公社

标题: 如何通过molden生成molcas所需要的.Orb文件 [打印本页]

作者
Author: 铭律凉 时间: 2020-4-16 17:42
标题: 如何通过molden生成molcas所需要的.Orb文件
我在使用molcas得到包含轨道信息的molden文件以后，用Multiwfn打开了molden文件对轨道做了一些修改，然后保存了新的molden文件，却不知道该怎么想办法让molcas通过新的molden文件生成后面所需要的.Orb文件，请问有什么办法或者其他软件可以生成.Orb吗？

作者
Author: zjxitcc 时间: 2020-4-16 18:36
如果你只是调换轨道、做做局域化啥的简单操作，直接在MOLCAS里面也有关键词可行的咧，无需通过molden。

作者
Author: 铭律凉 时间: 2020-4-16 20:00

zjxitcc 发表于 2020-4-16 18:36
如果你只是调换轨道、做做局域化啥的简单操作，直接在MOLCAS里面也有关键词可行的咧，无需通过molden。

我对所有轨道局域化，multiwfn很快就做完了，molcas跑了半天还告诉我失败了

作者
Author: zjxitcc 时间: 2020-4-16 21:14
本帖最后由 zjxitcc 于 2020-4-16 21:17 编辑

铭律凉发表于 2020-4-16 20:00
我对所有轨道局域化，multiwfn很快就做完了，molcas跑了半天还告诉我失败了

局域化其实用2*2写都很快，不知道为啥你的体系咋收敛不了。
我写过fchk文件转MOLCAS/OpenMolcas的.input和.INPORB文件的程序，转MOLCAS轨道很复杂。你有兴趣可以在高斯里加nosymm int=nobasistransform关键词然后得到的fchk文件发给我。并告知做局域化的轨道序号。

作者
Author: 铭律凉 时间: 2020-4-18 08:24

zjxitcc 发表于 2020-4-16 21:14
局域化其实用2*2写都很快，不知道为啥你的体系咋收敛不了。
我写过fchk文件转MOLCAS/OpenMolcas的.input ...

本身我这个体系就挺怪，因为是我人为构造的，不过后来通过加大圈数收敛了，现在我希望使用scf模块对齐做dft，初猜为localse后的轨道，其中有一些轨道想forzen住，请问该怎么输入？
（为了方便我已经将整个体系点群设为c1）

作者
Author: zjxitcc 时间: 2020-4-18 11:28

铭律凉发表于 2020-4-18 08:24
本身我这个体系就挺怪，因为是我人为构造的，不过后来通过加大圈数收敛了，现在我希望使用scf模块对齐做d ...

用MOLCAS的人都是奔着多参考方法去的。。。用DFT还不如高斯和ORCA。初猜嘛你产生好了就是xxx.Orb文件，重命名成xxx.INPORB，输入文件里(下属&SCF内）写FILEORB = xxx.INPORB就行。冻结轨道嘛，看看手册4.2.47 SCF部分的FROZen关键词

作者
Author: 铭律凉 时间: 2020-4-25 11:52
本帖最后由铭律凉于 2020-4-25 13:15 编辑

zjxitcc 发表于 2020-4-16 21:14
局域化其实用2*2写都很快，不知道为啥你的体系咋收敛不了。
我写过fchk文件转MOLCAS/OpenMolcas的.input ...

按你的要求我已经做好了.fch，请问能不能帮我把它全部转换位.Orb?如果可以的话能不能分享一下你的程序，因为以后我经常会用到

（请管理员删除此楼层和此附件）

作者
Author: zjxitcc 时间: 2020-4-25 12:26
本帖最后由 zjxitcc 于 2020-4-25 12:29 编辑

铭律凉发表于 2020-4-25 11:52
按你的要求我已经做好了.fch，请问能不能帮我把它全部转换位.Orb?如果可以的话能不能分享一下你的程序， ...

上面说过了，在高斯里加nosymm int=nobasistransform关键词然后得到的fch(k)文件发给我，你这个不是。你无需做自己局域化，直接把高斯算完的fch(k)文件发给我即可。如果你只是想要初始轨道，后面到MOLCAS里做多参考计算，那你其实这一步用RHF算个单点就行。当然，如果硬要用B3LYP也行，也请按照要求产生。并且再加上5D 7F关键词，因为高斯6-31G(d,p)默认是反人类的6D 7F。

作者
Author: 铭律凉 时间: 2020-4-25 13:13

zjxitcc 发表于 2020-4-25 12:26
上面说过了，在高斯里加nosymm int=nobasistransform关键词然后得到的fch(k)文件发给我，你这个不是。你 ...

我用的是m062x 6-31g（d，p）这两个关键词确实都加了，然后我是把这个扔到multiwfn里局域化的，应该没问题吧？我把原fch发给你好了，这次加了5D 7F需要分别对占据和非占据局域化，方法为pipek-mezey with Lowdin population 收敛程度位10^-6

作者
Author: zjxitcc 时间: 2020-4-25 18:33

铭律凉发表于 2020-4-25 13:13
我用的是m062x 6-31g（d，p）这两个关键词确实都加了，然后我是把这个扔到multiwfn里局域化的，应该没问 ...

文件描述如下

M062X正则轨道：Dimertest.INPORB
M062X局域轨道：Dimertest_5D7F_RM062X_PM.fchk，Dimertest_5D7F_RM062X_PM.INPORB
输入文件：Dimertest.input （可自行修改关键词满足自己的计算需求。若想用Dimertest.INPORB里的轨道重复M062X计算亦可，SCF会在3圈内收敛，能量第一圈起即不再变化。但能量会与高斯的有所差别，这是不同程序内部泛函写法导致的，并非轨道问题，因为轨道一开始就收敛了）
链接：https://pan.baidu.com/s/1x0Cq6Q0xhreu-u-PbGgKAg
提取码：ov7m

作者
Author: 铭律凉 时间: 2020-4-25 18:48

zjxitcc 发表于 2020-4-25 18:33
文件描述如下

M062X正则轨道：Dimertest.INPORB

谢谢，其实我更想看看程序是怎么写的，毕竟以后我可能还要用到，总不能一直来麻烦你

作者
Author: zjxitcc 时间: 2020-5-5 21:38

铭律凉发表于 2020-4-25 18:48
谢谢，其实我更想看看程序是怎么写的，毕竟以后我可能还要用到，总不能一直来麻烦你

我把程序开源放到gitlab上了，欢迎使用
https://gitlab.com/jxzou/mokit

欢迎光临计算化学公社 (http://bbs.keinsci.com/)