计算化学公社

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

[辅助/分析程序] 调用Multiwfn计算分子间电荷转移积分的程序[测试版]

  [复制链接 Copy URL]

9

帖子

0

威望

383

eV
积分
392

Level 3 能力者

31#
发表于 Post on 2022-4-20 22:29:15 | 只看该作者 Only view this author
本帖最后由 kaiden 于 2022-4-20 22:42 编辑

楼主您好,我最近尝试使用您的程序算电荷转移积分,但是过程遇到些问题,请教下面报错会是什么问题呢?谢谢!(window版程序,已经正确设置Multiwfn环境变量,用楼主的sample测试正常。我的体系中不加-d47 时显示系数矩阵不可逆,加了-d47时出错如下)# Using "dimer.fchk" for dimer.
# Using "monomer1.fchk" for monomer1.
# Using "monomer2.fchk" for monomer2.
# Using "MYDIMER.47" for Fock matrix of dimer.
# Calling Multiwfn to obtain amount of orbitals and electrons of monomer1 ...
# Reading amount of orbitals and electrons of monomer1 ...
# Calling Multiwfn to obtain amount of orbitals and electrons of monomer2 ...
# Reading amount of orbitals and electrons of monomer2 ...
# Calling Multiwfn to obtain amount of orbitals and electrons of dimer ...
# Reading amount of orbitals and electrons of dimer ...
# number of electrons of monomer1/monomer2/dimer:
#                             406     406   812
# number of  orbitals of monomer1/monomer2/dimer:
#                            1280    1280  2560
# Calling Multiwfn to obtain amount of orbitals and electrons of monomer1 ...
# Reading coefficient matrix of monomer1 ...
# Calling Multiwfn to obtain amount of orbitals and electrons of monomer2 ...
# Reading coefficient matrix of monomer2 ...
# Calling Multiwfn to obtain amount of orbitals and electrons of dimer ...
# Reading overlap matrix of dimer ...
# Reading Fock matrix of dimer from an NBO .47 file ...
At line 834 of file calc_coupling.f90
Fortran runtime error: Bad real number in item 1 of list input

Error termination. Backtrace:

Could not print backtrace: libbacktrace could not find executable to open
#0  0xffffffff
#1  0xffffffff
#2  0xffffffff
#3  0xffffffff
#4  0xffffffff
#5  0xffffffff
#6  0xffffffff
#7  0xffffffff
#8  0xffffffff
#9  0xffffffff
#10  0xffffffff
#11  0xffffffff
#12  0xffffffff
#13  0xffffffff
#14  0xffffffff
#15  0xffffffff

1187

帖子

5

威望

2859

eV
积分
4146

Level 6 (一方通行)

32#
 楼主 Author| 发表于 Post on 2022-4-21 00:07:41 | 只看该作者 Only view this author
kaiden 发表于 2022-4-20 22:29
楼主您好,我最近尝试使用您的程序算电荷转移积分,但是过程遇到些问题,请教下面报错会是什么问题呢?谢谢 ...

麻烦把文件发给我,我来测一下。

9

帖子

0

威望

383

eV
积分
392

Level 3 能力者

33#
发表于 Post on 2022-4-21 20:59:32 | 只看该作者 Only view this author
snljty 发表于 2022-4-21 00:07
麻烦把文件发给我,我来测一下。

压缩文件有点大,我用wetransfer传给您。下载链接已经私信发给您了,麻烦啦!

1

帖子

0

威望

73

eV
积分
74

Level 2 能力者

34#
发表于 Post on 2022-6-29 15:02:23 | 只看该作者 Only view this author
你好,我在使用您给的测试文件计算时。出现图中的报错, 显示 “Cannot open file 'Multiwfn_in.txt'”。 但是我在multiwfn文件夹下也没有找到这个文件,想请教下需要怎么处理这种情况。

202206291459293734..png (12.01 KB, 下载次数 Times of downloads: 74)

202206291459293734..png

419

帖子

1

威望

2549

eV
积分
2988

Level 5 (御坂)

35#
发表于 Post on 2022-6-29 15:28:59 | 只看该作者 Only view this author
milanox 发表于 2022-6-29 15:02
你好,我在使用您给的测试文件计算时。出现图中的报错, 显示 “Cannot open file 'Multiwfn_in.txt'”。  ...

在什么目录下运行的?看起来是没有权限,换个普通点的目录试试。

3

帖子

0

威望

27

eV
积分
30

Level 2 能力者

36#
发表于 Post on 2022-9-11 14:06:12 | 只看该作者 Only view this author
大佬你好,我用了您发的Multiwfn计算转移积分的程序,发现计算我的4个结构单元的寡聚物体系不收敛出现报错,提示(4700,4700)不收敛,请问下是超过4700个轨道的分子就无法计算嘛?谢谢~

3

帖子

0

威望

27

eV
积分
30

Level 2 能力者

37#
发表于 Post on 2022-9-11 14:18:26 | 只看该作者 Only view this author
报错如图所示,麻烦大佬了!

202209111417563792..png (15.62 KB, 下载次数 Times of downloads: 59)

202209111417563792..png

6

帖子

0

威望

61

eV
积分
67

Level 2 能力者

38#
发表于 Post on 2022-9-12 13:37:48 | 只看该作者 Only view this author
用Fortran的原因除了性能还有其它吗?

83

帖子

0

威望

647

eV
积分
730

Level 4 (黑子)

39#
发表于 Post on 2022-9-14 15:27:27 | 只看该作者 Only view this author
楼主 想请教一下 这样读取数值对吗?

202209141527102877..png (183.46 KB, 下载次数 Times of downloads: 65)

202209141527102877..png

4

帖子

0

威望

159

eV
积分
163

Level 3 能力者

40#
发表于 Post on 2022-10-17 17:56:41 | 只看该作者 Only view this author
老师您好,我计算c12h26 c14h28等几种有机小分子电荷迁移模型,在cp2k中已经进行了结构优化用了,pbe泛函和tzv2p基组,从优化了的结构中我选择了几个二聚体,用cp2k molden文件计算二聚体转移积分,已经添加了盒子信息和有效电荷,出现overlap matrix 可能在哪里有问题呢

4

帖子

0

威望

159

eV
积分
163

Level 3 能力者

41#
发表于 Post on 2022-10-17 18:15:29 | 只看该作者 Only view this author
这个是报错信息

202210171815166839..png (14.11 KB, 下载次数 Times of downloads: 57)

202210171815166839..png

4

帖子

0

威望

159

eV
积分
163

Level 3 能力者

42#
发表于 Post on 2022-10-17 21:59:27 | 只看该作者 Only view this author
不好意思已经解决了,感谢

4

帖子

0

威望

159

eV
积分
163

Level 3 能力者

43#
发表于 Post on 2022-10-17 22:02:10 | 只看该作者 Only view this author
未名|梧桐雨 发表于 2022-9-11 14:18
报错如图所示,麻烦大佬了!

这个我也遇到了解决了,你控制打印让他打印出的和基函数一致就行了

419

帖子

1

威望

2549

eV
积分
2988

Level 5 (御坂)

44#
发表于 Post on 2022-10-17 22:10:50 | 只看该作者 Only view this author
兰兰爱学习 发表于 2022-9-14 15:27
楼主 想请教一下 这样读取数值对吗?

读J_eff,后两组数据平时一般不用。读第一行和最后一行没问题。

419

帖子

1

威望

2549

eV
积分
2988

Level 5 (御坂)

45#
发表于 Post on 2022-10-17 22:13:22 | 只看该作者 Only view this author
seahurt 发表于 2022-9-12 13:37
用Fortran的原因除了性能还有其它吗?

读Gaussian里用Fortran输出的科学记数法双精度浮点数比较省事,就是那个带D的格式,因为有时候甚至不含那个D,比如指数超过两位,1.234567E+102会被输出成1.234567+102而不是1.234567D+102(这个虽然可以写程序的时候控制),别的语言读这个太麻烦了。另外就是Fortran调LAPACK比C方便一些,数组运算也方便很多,不用自己写一堆辅助函数。大概是这些原因。

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 20:41 , Processed in 0.274648 second(s), 27 queries , Gzip On.

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