计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1274|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[其它量化程序] 求助:Turbomole程序中基函数球谐部分的存放顺序

[复制链接 Copy URL]

7

帖子

0

威望

295

eV
积分
302

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 wangyongchaocc 于 2022-12-8 10:58 编辑

最近需要读Turbomole程序输出的分子轨道(MO系数存在mos或者alpha/beta文件中),调用外部积分库计算一些量。为此需要知道Turbomole程序输出的MO系数中基函数球谐部分的存放顺序,请问大家有知道这一信息的吗,谢谢!

5万

帖子

99

威望

5万

eV
积分
112384

管理员

公社社长

2#
发表于 Post on 2022-12-8 11:53:31 | 只看该作者 Only view this author
如果是输出molden文件,格式定义了确切的顺序:https://www.theochem.ru.nl/molden/molden_format.html
对于其它的输出文件,也可以用简单体系,和输出顺序已知的程序得到的系数对照一下来判断

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
wangyongchaocc + 4 赞!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

7

帖子

0

威望

295

eV
积分
302

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-12-8 16:08:45 | 只看该作者 Only view this author
sobereva 发表于 2022-12-8 11:53
如果是输出molden文件,格式定义了确切的顺序:https://www.theochem.ru.nl/molden/molden_format.html
对 ...

好的,非常感谢您!

3808

帖子

4

威望

8000

eV
积分
11888

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2024-2-20 17:40:52 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2024-2-20 19:17 编辑

提醒一下,molden文件强烈依赖于产生其的量化程序,经常不遵循molden标准,因此molden文件现在虽然看似通用,但其实有点“一团糟”,普通用户难以获知传的东西是否严格正确。

Turbomole的轨道可以放在好几处,最常见的是mos文件,其5个球谐函数与Gaussian的除了顺序不对应外,有的还相差一个负号,详见MOKIT代码
https://gitlab.com/jxzou/mokit/- ... ref_type=heads#L982
如果你的AO单/双电子积分与Turbomole的角动量顺序不对应,或者相差一个(些)负号,那么运算结果是不对的。

当然,也可以不用了解这些,直接使用MOKIT的fch2tm和molden2fch小程序传轨道,前者从fch文件产生control和mos文件,后者将Turbomole的molden文件转化为fch文件,例子见http://bbs.keinsci.com/forum.php ... 281430&fromuid=2632。fch文件格式大家是很标准统一的。这样也便于你使用fch文件进行波函数分析,或者进一步传轨道给其他程序。

MOKIT的molden2fch小程序特地区分了不同的量化程序,例如
  1. molden2fch a.molden -orca
  2. molden2fch a.molden -tm
复制代码
自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 14:27 , Processed in 0.164571 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list