应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】我想算甘油分子的静电势能
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 532  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

castep

木虫 (著名写手)

[交流] 【求助】我想算甘油分子的静电势能

我的输入文件如下, 不知道对不对 ,,,,,,,,,,
大家帮我改一下

%chk=glycerolesp.chk
%mem=6MW
%nproc=1
# prop=efg b3lyp/6-311g(d,p) geom=connectivity

Title Card Required

0 1
C
C                  1              B1
C                  2              B2    1              A1
O                  2              B3    1              A2    3              D1
H                  1              B4    2              A3    4              D2
O                  1              B5    2              A4    4              D3
H                  2              B6    1              A5    6              D4
H                  3              B7    2              A6    1              D5
H                  4              B8    2              A7    1              D6
O                  3              B9    2              A8    1              D7
H                  3             B10    2              A9    1              D8
H                  1             B11    6             A10    2              D9
H                  6             B12    1             A11    2             D10
H                 10             B13    3             A12    2             D11

   B1             1.53366451
   B2             1.53097031
   B3             1.42942357
   B4             1.09215034
   B5             1.42716077
   B6             1.09701650
   B7             1.09275710
   B8             0.96772356
   B9             1.43358443
   B10            1.09621313
   B11            1.09801814
   B12            0.96875826
   B13            0.96484924
   A1           112.78171469
   A2           109.45636475
   A3           109.89838628
   A4           110.60139482
   A5           109.10044882
   A6           110.92778414
   A7           105.29367163
   A8           111.25706595
   A9           109.47014596
   A10          111.71151753
   A11          106.05552054
   A12          105.64191555
   D1           121.20852547
   D2           -68.68095782
   D3            47.99266303
   D4           166.14968661
   D5           -55.65143407
   D6           -20.77720217
   D7            62.62795496
   D8          -175.14840762
   D9          -123.37230752
   D10           58.45971151
   D11           65.24787716

1 2 1.0 5 1.0 6 1.0 12 1.0
2 3 1.0 4 1.0 7 1.0
3 8 1.0 10 1.0 11 1.0
4 9 1.0
5
6 13 1.0
7
8
9
10 14 1.0
11
12
13
14
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
多多交流！
1楼 2009-12-13 16:12:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
分子的静电势值是一个分子空间各处、各点有一个值的概念。

因此要用cubegen命令或cube关键词并指定计算potential（静电势）来进行，这样得到分子空间一个格点文件，每一个格点有一个具体的静电势值，该文件可以作图观察。
赞一下(1人) 回复此楼
2楼2009-12-13 16:42:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
恩 这样算只能得到每个原子处的静电势吧。不过输入文件没有问题
**********************************************************************

            Electrostatic Properties Using The SCF Density

**********************************************************************

       Atomic Center    1 is at  -1.242854 -0.084919  0.685933
       Atomic Center    2 is at  -0.013426  0.801205  0.450523
       Atomic Center    3 is at   1.305318  0.040445  0.611852
       Atomic Center    4 is at  -0.068330  1.351233 -0.867698
       Atomic Center    5 is at  -2.142463  0.532517  0.733552
       Atomic Center    6 is at  -1.436225 -0.977176 -0.411005
       Atomic Center    7 is at  -0.026437  1.631110  1.167833
       Atomic Center    8 is at   1.363226 -0.436538  1.593306
       Atomic Center    9 is at  -0.665585  0.771601 -1.361460
       Atomic Center   10 is at   1.421357 -1.009817 -0.356987
       Atomic Center   11 is at   2.144343  0.739705  0.518280
       Atomic Center   12 is at  -1.148282 -0.625154  1.637167
       Atomic Center   13 is at  -0.622498 -1.498438 -0.479038
       Atomic Center   14 is at   1.469044 -0.569096 -1.213973
-----------------------------------------------------------------

              Electrostatic Properties (Atomic Units)

-----------------------------------------------------------------
    Center     Electric         -------- Electric Field --------
               Potential          X             Y             Z
-----------------------------------------------------------------
    1 Atom    -14.720978      0.001965      0.006001      0.005424
    2 Atom    -14.717241      0.000346     -0.004423      0.007294
    3 Atom    -14.720812     -0.002909      0.007754      0.004963
    4 Atom    -22.376723      0.021767      0.039209     -0.014433
    5 Atom     -1.113387     -0.015706      0.009899      0.000600
    6 Atom    -22.373226     -0.039391     -0.003780     -0.026441
    7 Atom     -1.109289     -0.000307      0.014407      0.011147
    8 Atom     -1.105460      0.001359     -0.008438      0.016328
    9 Atom     -1.006973     -0.015003     -0.012322     -0.013604
   10 Atom    -22.360869      0.004139     -0.045929      0.010639
   11 Atom     -1.112890      0.014450      0.010731     -0.002104
   12 Atom     -1.111041      0.001083     -0.009813      0.014728
   13 Atom     -1.006894      0.019259     -0.013558     -0.002951
   14 Atom     -0.994775      0.002428      0.008844     -0.021994
-----------------------------------------------------------------
赞一下(1人) 回复此楼
3楼2009-12-13 16:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
Electrostatic potential    静电势

Electrostatic Properties  静电性能

Electric  Potential         电势
赞一下(1人) 回复此楼
4楼2009-12-13 19:52:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

castep

木虫 (著名写手)
引用回帖:
Originally posted by zhou2009 at 2009-12-13 16:42:
分子的静电势值是一个分子空间各处、各点有一个值的概念。

因此要用cubegen命令或cube关键词并指定计算potential（静电势）来进行，这样得到分子空间一个格点文件，每一个格点有一个具体的静电势值，该文件可以 ...

具体怎么做，我的帖子已经给了输入文件了，你能不能帮我改成一个能算得输入文件 ，用你说的那个关键词cube， 我不知道怎么做
赞一下 回复此楼
多多交流！
5楼2009-12-14 09:52:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 castep 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂：085600材料与化工，考材科基，总分319 +20 678lucky 2026-03-31 25/1250 2026-04-04 14:56 by 聪明的大松鼠
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 320分人工智能调剂 +6 振—TZ 2026-04-03 6/300 2026-04-04 01:05 by Dream007008
[考研] 295求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 21:01 by zhq0425
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 275求调剂 +14 waltzh 2026-04-01 14/700 2026-04-03 11:24 by qqq850810
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 289求调剂 +4 Acesczlo 2026-03-29 5/250 2026-04-03 10:09 by 不168
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 349求调剂 +10 zwjjjjjj 2026-03-31 10/500 2026-04-02 20:13 by dongzh2009
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 09:15 by olim
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
[考研] 291求调剂 +3 迷蒙木木 2026-04-01 4/200 2026-04-01 11:07 by 逆水乘风
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk￥
[考研] 求收留 +8 1943443204 2026-03-28 8/400 2026-03-31 15:00 by -迷了路啊路
信息提示
关闭
请填处理意见
关闭