计算化学公社

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

[辅助/分析程序] 一个Fortran写的电荷转移积分程序

  [复制链接 Copy URL]

286

帖子

3

威望

5617

eV
积分
5963

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Warm_Cloud 于 2020-7-19 13:31 编辑

attach://20792.x电荷转移积分的计算是很多小伙伴的需求,这两天基于高斯的输出文件做了一个,细节可以参考http://bbs.keinsci.com/thread-11369-1-1.html
使用方法很简单,需要提供dimer和monomers的高斯输出文件和fchk文件,以此按照提示输入就可以。输入文件要加 nosymm iop(3/33=1)。
顺便说一下,以后更强大和好用的电荷转移功能会加入到我自己写的量化软件中。
CT.x (8.96 MB, 下载次数 Times of downloads: 197)
代码地址:https://github.com/zhangyingfeng1993/Charge-Transport-Integral
有问题的小伙伴可以留言讨论。


评分 Rate

参与人数
Participants 11
eV +51 收起 理由
Reason
sailing + 4 GJ!
B1te_theDust + 5 膜拜大佬
corei70715 + 3 谢谢分享
yjmaxpayne + 5
柒月小鱼 + 5 big law good
ggdh + 5 大师终于出手了
978142355 + 5 大师大法呱呱叫
wbn + 5 赞!
linqiaosong + 4 大师好,大师妙
ene + 5
冰释之川 + 5 好物!

查看全部评分 View all ratings

1102

帖子

18

威望

6643

eV
积分
8105

Level 6 (一方通行)

計算化学の社畜

2#
发表于 Post on 2019-8-23 22:33:33 | 只看该作者 Only view this author
大师大法好!先收藏之,待日后使用
Stand on the shoulders of giants

125

帖子

0

威望

3509

eV
积分
3634

Level 5 (御坂)

アルトリア・ペンドラゴン

3#
发表于 Post on 2019-8-24 11:03:37 | 只看该作者 Only view this author
开源大师

490

帖子

2

威望

4881

eV
积分
5411

Level 6 (一方通行)

4#
发表于 Post on 2019-8-24 11:45:33 | 只看该作者 Only view this author
大师大法好

245

帖子

0

威望

2576

eV
积分
2821

Level 5 (御坂)

5#
发表于 Post on 2019-8-24 18:05:42 来自手机 | 只看该作者 Only view this author
大师大法好!

144

帖子

0

威望

2818

eV
积分
2962

Level 5 (御坂)

6#
发表于 Post on 2019-9-16 17:15:38 | 只看该作者 Only view this author
大师也真够懒的,简单写个小命令行的编译方法噻。。。不过,我自己随便捣鼓了一下:

ifort -O3 CT.f90 -L/xxxxx/mkl/intel64/ -lmkl_intel_lp64 -lmkl_core -lmkl_intel_thread -lpthread -lm -qopenmp -o CT.x

对同一系统在同一理论级别下测试结果:

VOTCA CTP (git master): -12.4 meV
大师CT: -18.2meV

差别很小。

有两个小建议:
1) 建议暖大师提供 HOMO,HOMO-1, LUMO,LUMO+1搜索并单独输出的选项,绝大多数时候,这几个轨道的耦合是最主要的
2)耦合部分输出加入每个项的对应说明。

286

帖子

3

威望

5617

eV
积分
5963

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2019-9-17 08:09:09 | 只看该作者 Only view this author
yjmaxpayne 发表于 2019-9-16 17:15
大师也真够懒的,简单写个小命令行的编译方法噻。。。不过,我自己随便捣鼓了一下:

ifort -O3 CT.f90 - ...

啊,吼吼吼

7

帖子

0

威望

545

eV
积分
552

Level 4 (黑子)

8#
发表于 Post on 2019-10-26 13:17:34 | 只看该作者 Only view this author
请问,这四项哪一个是我想要的电荷转移积分啊
J(meV)    S
Je(meV)   Je(kcal/mol)

286

帖子

3

威望

5617

eV
积分
5963

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2019-10-27 13:42:49 | 只看该作者 Only view this author
小白的暖暖 发表于 2019-10-26 13:17
请问,这四项哪一个是我想要的电荷转移积分啊
J(meV)    S
Je(meV)   Je(kcal/mol)

Je

7

帖子

0

威望

545

eV
积分
552

Level 4 (黑子)

10#
发表于 Post on 2019-10-28 10:02:51 | 只看该作者 Only view this author

55

帖子

1

威望

1551

eV
积分
1626

Level 5 (御坂)

11#
发表于 Post on 2019-11-30 21:50:47 | 只看该作者 Only view this author
太牛批了!

286

帖子

3

威望

5617

eV
积分
5963

Level 6 (一方通行)

12#
 楼主 Author| 发表于 Post on 2020-7-19 13:32:27 | 只看该作者 Only view this author
更新了一下输出。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
ene + 5

查看全部评分 View all ratings

1

帖子

0

威望

399

eV
积分
400

Level 3 能力者

13#
发表于 Post on 2021-12-3 14:32:00 | 只看该作者 Only view this author
请问楼主,您计算转移积分所采用的理论方法是基于这篇 JACS,128(30):9882-9886吗?是属于能级劈裂法还是在位能修正法啊?我个人的理解是作者是在能级劈裂法的基础上进行修正?

9

帖子

0

威望

363

eV
积分
372

Level 3 能力者

14#
发表于 Post on 2022-4-24 12:15:05 | 只看该作者 Only view this author
楼主您好,最近我用您写的程序测试电荷转移积分遇到一点问题,请教下面报错会是什么问题呢?谢谢!(在m062x/ma-SVP下计算的dimer,monomer单点。dimer.out载入时出错,但是把monomer.out当作dimer载入又没问题。)

nput for Dimer:
  Please input the Gaussian log file :<e.g. out.log>
pPDM_Dimer1_SP.out

  Spin Multiplicity:             1
  Number of basis functions:  2560
  Readind Orbit...
forrtl: severe (59): list-directed I/O syntax error, unit 1101, file /root/Documents/Study/Charge-Transport-Integral-master/pPDM_Dimer1_SP.fchk
Image              PC                Routine            Line        Source            
CT.x               000000000094B016  Unknown               Unknown  Unknown
CT.x               0000000000967727  Unknown               Unknown  Unknown
CT.x               0000000000965985  Unknown               Unknown  Unknown
CT.x               0000000000404E8B  Unknown               Unknown  Unknown
CT.x               0000000000401AE0  Unknown               Unknown  Unknown
CT.x               0000000000400D29  Unknown               Unknown  Unknown
CT.x               0000000000400CDE  Unknown               Unknown  Unknown
CT.x               00000000009DE5E9  Unknown               Unknown  Unknown
CT.x               0000000000400BBA  Unknown               Unknown  Unknown

286

帖子

3

威望

5617

eV
积分
5963

Level 6 (一方通行)

15#
 楼主 Author| 发表于 Post on 2022-4-24 20:28:19 | 只看该作者 Only view this author
kaiden 发表于 2022-4-24 12:15
楼主您好,最近我用您写的程序测试电荷转移积分遇到一点问题,请教下面报错会是什么问题呢?谢谢!(在m062 ...

文件上传一下我试试。

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 17:37 , Processed in 0.276073 second(s), 31 queries , Gzip On.

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