如何通过molden生成molcas所需要的.Orb文件

铭律凉 · 发表于 Post on 2020-4-16 17:42:26

我在使用molcas得到包含轨道信息的molden文件以后，用Multiwfn打开了molden文件对轨道做了一些修改，然后保存了新的molden文件，却不知道该怎么想办法让molcas通过新的molden文件生成后面所需要的.Orb文件，请问有什么办法或者其他软件可以生成.Orb吗？

zjxitcc · 发表于 Post on 2020-5-5 21:38:33

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

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

铭律凉 · 发表于 Post on 2020-4-25 18:48:32

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

M062X正则轨道：Dimertest.INPORB

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

zjxitcc · 发表于 Post on 2020-4-25 18:33:34

铭律凉发表于 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

铭律凉 · 发表于 Post on 2020-4-25 13:13:07

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

zjxitcc · 发表于 Post on 2020-4-25 12:26:57

本帖最后由 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。

铭律凉 · 发表于 Post on 2020-4-25 11:52:50

本帖最后由铭律凉于 2020-4-25 13:15 编辑

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

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

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

zjxitcc · 发表于 Post on 2020-4-18 11:28:02

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

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

铭律凉 · 发表于 Post on 2020-4-18 08:24:16

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

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

zjxitcc · 发表于 Post on 2020-4-16 21:14:24

本帖最后由 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文件发给我。并告知做局域化的轨道序号。

铭律凉 · 发表于 Post on 2020-4-16 20:00:58

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

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

zjxitcc · 发表于 Post on 2020-4-16 18:36:10

如果你只是调换轨道、做做局域化啥的简单操作，直接在MOLCAS里面也有关键词可行的咧，无需通过molden。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Molcas] 如何通过molden生成molcas所需要的.Orb文件

浏览过的版块