计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Markmahao
打印 Print 上一主题 Last thread 下一主题 Next thread

[算法与编程] 请教波函数与空间的对应关系

[复制链接 Copy URL]

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

16#
发表于 Post on 2022-7-7 18:49:11 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2022-7-7 19:13 编辑
Markmahao 发表于 2022-7-7 18:36
最后还是旋转了系数矩阵,s和p的部分很快完成但又遇到一个问题,就是高角动量的笛卡尔型(比如d轨道), ...

5个展开成6个,展开系数是已知的,写程序时可以直接写进去,比较方便。但6个变回5个的时候,需要解线性方程组(示例),没法简单通过已知的一些系数线性组合得到,不知道你是不是这么做的。
自动做多参考态计算的程序MOKIT

123

帖子

0

威望

2617

eV
积分
2740

Level 5 (御坂)

17#
 楼主 Author| 发表于 Post on 2022-7-7 22:11:20 | 只看该作者 Only view this author
zjxitcc 发表于 2022-7-7 18:49
5个展开成6个,展开系数是已知的,写程序时可以直接写进去,比较方便。但6个变回5个的时候,需要解线性方 ...

完全不是这么做的,我学习一下,谢谢老师指点

123

帖子

0

威望

2617

eV
积分
2740

Level 5 (御坂)

18#
 楼主 Author| 发表于 Post on 2022-9-13 17:18:46 | 只看该作者 Only view this author
zjxitcc 发表于 2022-7-7 18:49
5个展开成6个,展开系数是已知的,写程序时可以直接写进去,比较方便。但6个变回5个的时候,需要解线性方 ...

游客,本帖隐藏的内容需要积分高于 150 才可浏览,您当前积分为 0

HCl_x1.fch

34.49 KB, 下载次数 Times of downloads: 0

旋转后计算结果

HCl_z.fch

34.49 KB, 下载次数 Times of downloads: 0

旋转前计算结果

rotated_MO.txt

13.85 KB, 下载次数 Times of downloads: 0

我旋转的结果

4289

帖子

4

威望

9550

eV
积分
13919

Level 6 (一方通行)

MOKIT开发者

19#
发表于 Post on 2023-6-26 00:11:00 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2023-6-26 00:12 编辑

MOKIT已支持轨道系数随分子朝向旋转,此功能无需计算AO基重叠积分,没用到Wigner-D矩阵。对球谐型和Cartesian型基函数都支持,以下是一个例子(在python中使用)
  1. from mokit.lib.mirror_wfn import rotate_atoms_wfn
  2. rotate_atoms_wfn('ZnMe2F.fch', 'ZnMe2F_new.gjf')
复制代码
第一个文件含一套已知的轨道系数,在旧的朝向下;第二个文件含新朝向下的坐标,支持gjf和xyz两种格式。运行后会产生ZnMe2F_r.fch文件,内含新朝向的坐标和轨道系数,SCF可1圈收敛。

评分 Rate

参与人数
Participants 2
eV +10 收起 理由
Reason
Markmahao + 5 好物!
wzkchem5 + 5

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 01:57 , Processed in 0.271195 second(s), 23 queries , Gzip On.

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