“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办，这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会！请点击此链接查看详情和报名方式，欢迎参加！

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会！报名正在进行中，请点击此链接查看详情，欢迎参加！

捐赠本论坛 (Donate this site)
开启辅助访问

计算化学公社

 找回密码 Forget password
 注册 Register
计算化学公社»论坛首页 › 大杂烩 (Miscellaneous) › 编程、软件开发 (Developement) › 做了个新的反应活性位点的预测程序，欢迎大家试用

做了个新的反应活性位点的预测程序，欢迎大家试用

查看数: 831 | 评论数: 9 | 收藏 Add to favorites 16
关灯 | 提示：支持键盘翻页<-左 右->
帖子模式
    组图打开中，请稍候......
yjcmwgk
发布时间: 2026-2-20 16:12

正文摘要:

本帖最后由 yjcmwgk 于 2026-2-25 18:50 编辑 这半年多，做了个新的反应活性预测算法，带程序的。 我自认为它能准确预测亲电/亲核反应活性的具体位点和强度。 程序的exe文件在这儿下载 https://www.rsc.org/ ...

回复 Reply

wxyhgk 发表于 Post on 2026-3-28 17:47:43
yjcmwgk 发表于 2026-3-27 09:26
确实没测试100原子以上的大分子啊……这样，你提几个，我来试试？

公司实际开发分子，你可以试试

已经优化好的 xyz 坐标
  1. 130
  2. name: Optimized Structure; charge: 0; multiplicity: 1; energy: -205.81913700
  3. C    -1.21340016   0.19899277  -0.20164603
  4. C    -1.22971496  -1.20460695  -0.19100389
  5. C    -0.02928417  -1.88766343  -0.12910457
  6. C     1.19606474  -1.24814589  -0.11651527
  7. C     1.22874960   0.15519299  -0.12721538
  8. C     0.01526204   0.84744012  -0.01509700
  9. N    -2.42566899   0.86738864  -0.27294917
  10. N     2.46581214   0.78017222  -0.12065623
  11. C    -2.78913688   2.17998016  -0.52876360
  12. C    -4.14517664   2.36502965  -0.15595567
  13. C    -4.64575089   1.07898391   0.26122421
  14. C    -3.56763942   0.19358015   0.11129751
  15. C     3.55266235   0.07071485   0.34887615
  16. C     4.64552030   0.92036203   0.57665452
  17. C     4.22188064   2.21940762   0.11641331
  18. C     2.89032427   2.07852938  -0.35149310
  19. C    -5.87178116   0.57404119   0.66709886
  20. C    -5.99581739  -0.79449652   0.89198084
  21. C    -4.91269292  -1.65729518   0.67086131
  22. C    -3.68270884  -1.17364804   0.26458380
  23. C    -2.10538235   3.19673633  -1.17901359
  24. C    -2.72938923   4.41946305  -1.32199920
  25. C    -4.03226508   4.64970902  -0.86486820
  26. C    -4.75209017   3.60066795  -0.30544618
  27. C    -4.63839319   5.98027289  -1.01954182
  28. C    -7.28105672  -1.35058967   1.34300795
  29. C     3.60826798  -1.29882602   0.51432521
  30. C     4.78438901  -1.82083800   1.02002998
  31. C     5.87458387  -0.99277011   1.32509208
  32. C     5.81683671   0.37718454   1.08212167
  33. C     4.88017528   3.43223031   0.00335549
  34. C     4.23656504   4.50180437  -0.60770551
  35. C     2.96006812   4.31491976  -1.15091459
  36. C     2.28659653   3.11475047  -1.04826929
  37. Cl    0.02320400   2.44173346   0.69431401
  38. C    -2.50675191  -2.02456190  -0.13474006
  39. C     2.44051590  -2.10858907   0.01618283
  40. C    -2.32594664  -3.22719204   0.76627994
  41. C    -2.41612461  -4.41850177   0.03310247
  42. C    -2.68493803  -4.08422481  -1.35748238
  43. C    -2.75324517  -2.68939455  -1.47742210
  44. C     2.76778678  -2.75384840  -1.31920167
  45. C     2.65229721  -4.14811893  -1.23447015
  46. C     2.27370421  -4.50123816   0.12553741
  47. C     2.16225667  -3.32235115   0.87570817
  48. C    -2.86507411  -4.88885916  -2.46982575
  49. C    -3.11294466  -4.28524660  -3.69242813
  50. C    -3.18093572  -2.90517172  -3.80437710
  51. C    -2.99985926  -2.09219072  -2.69317071
  52. C    -2.07577252  -3.24553909   2.12059569
  53. C    -1.90781371  -4.47334347   2.74622887
  54. C    -1.98578426  -5.65388239   2.02384511
  55. C    -2.24294705  -5.63993853   0.66227870
  56. C     3.12328686  -2.14110144  -2.49980823
  57. C     3.36617221  -2.93790049  -3.61099117
  58. C     3.25027523  -4.31721091  -3.53424719
  59. C     2.89338980  -4.93646747  -2.34694930
  60. C     2.02502679  -5.72950523   0.71478886
  61. C     1.67185862  -5.76291089   2.05432053
  62. C     1.57322671  -4.59493849   2.79445791
  63. C     1.81626922  -3.36017814   2.20848672
  64. C    -5.96253258   6.12103730  -1.42829121
  65. C    -6.52357976   7.37514916  -1.58092998
  66. C    -5.77402241   8.51050596  -1.32256286
  67. C    -4.45849358   8.38294950  -0.90944078
  68. C    -3.89449776   7.12976889  -0.76147132
  69. C    -7.31641331  -2.33526440   2.32770066
  70. C    -8.52214415  -2.84807219   2.76785097
  71. C    -9.71285942  -2.39051346   2.22886298
  72. C    -9.68945185  -1.41702628   1.24438300
  73. C    -8.48471181  -0.89993919   0.80600870
  74. C     4.89675352   5.80994962  -0.72456551
  75. C     7.10034575  -1.58781953   1.87933099
  76. C     4.17510814   6.98467177  -0.52324334
  77. C     4.78970223   8.21749917  -0.63518594
  78. C     6.13451237   8.29877014  -0.95526648
  79. C     6.86206704   7.13787565  -1.15686397
  80. C     6.25020639   5.90391202  -1.03990466
  81. C     8.35647041  -1.16868715   1.44694913
  82. C     9.50579837  -1.72533964   1.97576144
  83. C     9.42064228  -2.70785591   2.94780461
  84. C     8.17728253  -3.13197385   3.38582880
  85. C     7.02681178  -2.57903499   2.85546362
  86. H    -0.05038904  -2.96507719  -0.08609126
  87. H    -6.71573173   1.22966987   0.81676729
  88. H    -5.05309916  -2.72233049   0.78264101
  89. H    -1.12393241   3.03965644  -1.59588086
  90. H    -2.20847249   5.21586203  -1.83263780
  91. H    -5.77245736   3.74539673   0.01580089
  92. H     4.87960472  -2.88919416   1.14630929
  93. H     6.66797534   1.00519896   1.29555540
  94. H     5.88096515   3.54328787   0.39232264
  95. H     2.50190581   5.12681571  -1.69611985
  96. H     1.33032615   2.98846665  -1.52930847
  97. H    -2.81478667  -5.96435398  -2.38741522
  98. H    -3.25592478  -4.89753707  -4.57086804
  99. H    -3.37667359  -2.45655897  -4.76711866
  100. H    -3.05392473  -1.01680138  -2.78242854
  101. H    -2.00759068  -2.32458288   2.68106830
  102. H    -1.71235276  -4.50894176   3.80765095
  103. H    -1.84540575  -6.59751833   2.53033237
  104. H    -2.30441436  -6.56233380   0.10428638
  105. H     3.21484196  -1.06630162  -2.56207215
  106. H     3.64883825  -2.47710524  -4.54597362
  107. H     3.44258558  -4.91666552  -4.41215303
  108. H     2.80718304  -6.01141720  -2.29127153
  109. H     2.10321495  -6.64216227   0.14305098
  110. H     1.47225761  -6.71208189   2.52966562
  111. H     1.30279389  -4.64608789   3.83859720
  112. H     1.73193463  -2.44877367   2.78231367
  113. H    -6.54792158   5.24005122  -1.64846784
  114. H    -7.54976833   7.46841225  -1.90631081
  115. H    -6.21272598   9.49014320  -1.44265877
  116. H    -3.87042442   9.26489235  -0.69970874
  117. H    -2.87254095   7.03663382  -0.42332562
  118. H    -6.39101868  -2.68444264   2.76288214
  119. H    -8.53391773  -3.60623752   3.53753409
  120. H   -10.65460504  -2.79106464   2.57416519
  121. H   -10.61502402  -1.06069941   0.81545459
  122. H    -8.47287798  -0.15345351   0.02503266
  123. H     3.12930335   6.92726314  -0.25772057
  124. H     4.21800629   9.11953613  -0.47060317
  125. H     6.61312661   9.26255973  -1.04731983
  126. H     7.91088217   7.19515188  -1.41039917
  127. H     6.81981512   5.00278159  -1.21568758
  128. H     8.43000084  -0.41559374   0.67574412
  129. H    10.47300025  -1.39353628   1.62651587
  130. H    10.31911980  -3.14117206   3.36190002
  131. H     8.10459432  -3.89597974   4.14644376
  132. H     6.05956823  -2.90253982   3.21240997
复制代码


yjcmwgk 发表于 Post on 2026-3-27 09:26:37
wxyhgk 发表于 2026-3-26 23:06
在预测性的算法在我的测试里面小分子都不错，但是一旦分子大了，基本都不准，我看你也没测试大分子（100 原 ...

确实没测试100原子以上的大分子啊……这样，你提几个，我来试试？
wxyhgk 发表于 Post on 2026-3-26 23:06:14
在预测性的算法在我的测试里面小分子都不错，但是一旦分子大了，基本都不准，我看你也没测试大分子（100 原子以上的），不知道效果如果
yjcmwgk 发表于 Post on 2026-2-24 23:17:58
牧生 发表于 2026-2-24 17:02
谢谢贴心指导，确实是因为这个原因。

现在已经完全走通了流程

啊 谢谢你的反馈 正好借楼重申一下 如果cube的头定义（格子定义）不一样 那两个cube互相运算将会毫无意义 因为每个空间采样点的坐标都不一样 完全没意义了 所以这个程序一开始就会检测cube的头定义
牧生 发表于 Post on 2026-2-24 17:02:46
yjcmwgk 发表于 2026-2-24 15:38
就是这一句啊
ATTENTION!!! Cubes definitions are different, please check!

谢谢贴心指导，确实是因为这个原因。

现在已经完全走通了流程
牧生 发表于 Post on 2026-2-24 14:22:10
本帖最后由 牧生 于 2026-2-24 14:52 编辑

由于我们没有买gaussian，所以我用了orca来计算了单点能，以期望重复楼主的操作。以我的理解，主要就是要得到某个分子0电荷单点能和静电势，+1和-1电荷激发态的单点能。

但有一步卡住了，请楼主帮忙解答一下。

为了在我的笔记本电脑上操作，使用的计算级别比较低。我用EthyleneKetene进行计算

1. 在b97-3c级别下，优化电中性的EthyleneKetene，并得到了当前级别下的单点能文件EthyleneKetene.gbw，直接用Multiwfn自带的ESPiso.bat （计算ESP时改成了高质量格点的选项）计算得到了ESP1.cub和density1.cub文件，分别改名为ESP0.cub和density0.cub

2. 用优化后的EthyleneKetene，再次生成单点能计算文件，电荷和自旋多重度该为1 2，同样方法得到density+1.cub

3. 改2中的电荷和自旋多重度该为-1 2，同样方法得到density-1.cub

4.把上述得到的4个cub文件和RAI.exe放在同一个文件夹里面，打开RAI.exe, 依次把density-1.cub，density0.cub，density+1.cub，ESP0.cub拖进去    （拖的顺序是带负电，0电荷，正电荷）

---- RAI(+) and RAI(-) Calculating Program (v 1.0) ----
Input the cube file name of negative charge file ( -1 e ) unit of length must be Bohr, the same below:
C:\Users\gxw85\Desktop\g6\1\density-1.cub
Input the cube file name of electron neutral file (  0 e ):
C:\Users\gxw85\Desktop\g6\1\density0.cub
Input the cube file name of positive charge file  ( +1 e ):
C:\Users\gxw85\Desktop\g6\1\density+1.cub
Input the cube file name of electrostatic file    (  esp ):
C:\Users\gxw85\Desktop\g6\1\ESP0.cub
ATTENTION!!! Cubes definitions are different, please check!
PAUSE
To resume execution, type go.  Other input will terminate the job.



在这里就不对了。请问可能错在哪一步呢？






yjcmwgk 发表于 Post on 2026-2-20 21:20:59
wal 发表于 2026-2-20 18:51
程序直接附在论文里好啊 赞
见过发论文挂GitHub仓库假装开源的那种 发完论文仓库直接撤了 令人感到幽默

哈哈哈哈哈哈哈哈 这是我的习惯
如果论文里我用的什么自编代码 我总是随着论文一起发表
明面上的理由：为了让大家复现计算结果更方便！
暗戳戳的心理：你用了我的程序就要老老实实引我文章！

评分 Rate

参与人数
Participants 1		eV +5 收起 理由
Reason
student0618 + 5 赞！

查看全部评分 View all ratings

wal 发表于 Post on 2026-2-20 18:51:34
程序直接附在论文里好啊 赞
见过发论文挂GitHub仓库假装开源的那种 发完论文仓库直接撤了 令人感到幽默

评分 Rate

参与人数
Participants 2		eV +10 收起 理由
Reason
yjcmwgk + 5 233333
ChrisZheng + 5 233333

查看全部评分 View all ratings

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

GMT+8, 2026-4-13 23:29 , Processed in 0.179167 second(s), 26 queries , Gzip On.

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