24小时热门版块排行榜    

查看: 879  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

胭脂。花开

铁虫 (初入文坛)

[交流] 【求助】高斯计算单点能想要电荷用了PBC,报错。MDV<1

各位高手。我现在在算晶体材料的电荷,不能很好的选取CLUSTER,因为Unit Cell原子很少,考虑用整个胞算能量得到电荷。
下面是我的输入文件和报错信息,还请各位高手帮忙。。
%chk=cpl-2.chk
%mem=170MW
%rwf=1,2000mb,2,2000mb,3,2000mb,4,-1
%nproc=2
# b3lyp/6-31+g(d) pbc=nkpoints=1 pop=(readradii,chelpg)

cpl-2

0 1
Cu                 1.06366137    5.16263199    6.19161811
O                  1.77562608    4.78164751    8.06027023
O                  2.89993688    5.33712233    9.91906309
O                 -0.23081526    8.20549859   10.05419931
O                  3.90084132    8.42061985    0.36747193
N                  1.51086689    7.05224825    6.98989704
N                  2.40126325    9.31644824    8.17440194
N                  0.61339682    3.32811795    5.81117910
C                  2.28237737    5.54417307    8.75513405
C                  2.19911802    7.04028161    8.14739612
C                  1.26860301    8.27034109    6.38226628
C                  1.73919692    9.35513110    7.04508749
C                  2.65343205    8.19603659    8.76102817
C                  3.30896627    8.25447833   10.08742075
C                 -0.27458808    2.78349659    6.67815138
C                  0.90148271    2.63683565    4.66964773
C                  0.39127115    1.31382597    4.47428427
C                  3.93242064    1.47551479    6.48118042
C                 -0.41062865    0.73135663    5.40405577
H                 -0.60749080    3.32283036    7.41695322
H                  1.47952341    3.04731930    3.95549405
H                  0.66843848    0.80983554    3.68329262
H                  3.27228782    1.12709070    7.08795385
H                  0.83250285    8.29316120    5.57262776
H                  1.56646779   10.29687800    6.66572013
Cu                 1.95928176   19.07733199    9.88326967
Cu                 2.52426739   22.66676801    4.52497374
Cu                 1.62864699    8.75206801    0.83332218
O                  1.24731706   18.69634751    8.01461755
O                  1.81230268   23.04775249    2.65632162
O                  2.34061170    9.13305249    2.70197430
O                  0.12300625   19.25182233    6.15582469
O                  0.68799188   22.49227767    0.79752877
O                  3.46492250    8.57757767    4.56076716
O                  3.25375839   22.12019859    6.02068847
O                  3.81874402   19.62390141    0.66239254
O                  0.33417036    5.70920141    4.69590338
O                  0.25207306   22.33531985    4.99082399
O                 -0.31291257   19.40878015   10.34911992
O                  3.33585570    5.49408015    5.72576786
N                  1.51207624   20.96694825    9.08499075
N                  2.07706187   20.77715175    3.72669482
N                  2.07585251    6.86245175    1.63160111
N                  0.62167989   23.23114824    7.90048585
N                  1.18666551   18.51295176    2.54218992
N                  2.96624887    4.59825176    2.81610601
N                  2.40954631   17.24281795   10.26370868
N                  2.97453193   24.50128205    4.90541276
N                  1.17838245   10.58658205    0.45288317
C                  0.74056577   19.45887307    7.31975373
C                  1.30555139   22.28522693    1.96145781
C                  2.84736299    8.37052693    3.39683812
C                  0.82382511   20.95498161    7.92749166
C                  1.38881074   20.78911839    2.56919573
C                  2.76410364    6.87441839    2.78910020
C                  1.75434012   22.18504109    9.69262150
C                  2.31932574   19.55905891    4.33432558
C                  1.83358863    5.64435891    1.02397035
C                  1.28374622   23.26983110    9.02980030
C                  1.84873184   18.47426890    3.67150437
C                  2.30418254    4.55956890    1.68679156
C                  0.36951109   22.11073659    7.31385961
C                  0.93449671   19.63336341    1.95556368
C                  3.21841767    5.71866341    3.40273225
C                 -0.28602313   22.16917833    5.98746704
C                  0.27896249   19.57492167    0.62917111
C                  3.87395189    5.66022167    4.72912482
C                  3.29753122   16.69819659    9.39673640
C                  3.86251684   25.04590341    4.03844047
C                  0.29039754   11.13120341    1.31985545
C                  3.25143167   16.55153565    0.68864819
C                  2.68644605   25.19256435    6.04694412
C                  0.33649708   11.27786435   10.02794366
C                  3.76164323   15.22852597    0.88401166
C                  3.19665760   26.51557403    6.24230759
C                 -0.17371447   12.60087403    9.83258019
C                 -0.90947750   15.39021479    9.59370736
C                 -0.34449188   26.35388521    4.23541143
C                  4.49740626   12.43918521    1.12288449
C                  3.43357179   14.64605663   10.67083201
C                  3.99855741   27.09804337    5.31253608
C                  0.15435697   13.18334337    0.04575985
H                  3.63043393   17.23753036    8.65793456
H                  4.19541955   24.50656964    3.29963863
H                 -0.04250518   10.59186964    2.05865730
H                  2.67339096   16.96201930    1.40280187
H                  2.10840534   24.78208070    6.76109780
H                  0.91453779   10.86738070    9.31378998
H                  3.48447589   14.72453554    1.67500331
H                  2.91949027   27.01956446    7.03329923
H                  0.10345286   13.10486446    9.04158855
H                 -0.24934468   15.04179070    8.98693393
H                  0.31564094   26.70230930    3.62863800
H                  3.83727344   12.78760930    1.72965793
H                  2.19044028   22.20786120   10.50226002
H                  2.75542590   19.53623880    5.14396409
H                  1.39748848    5.62153880    0.21433184
H                  1.45647535   24.21157800    9.40916765
H                  2.02146097   17.53252200    4.05087172
H                  2.13145341    3.61782200    1.30742421
Tv                 4.71790000    0.00000000    0.00000000
Tv                 0.00000000   27.82940000    0.00000000
Tv                -1.12997124    0.00000000   10.71659186

[ Last edited by xuefei06 on 2009-2-26 at 23:18 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胭脂。花开

铁虫 (初入文坛)

各位大侠。还请多帮忙。。
被这个弄晕了。我试了运行。确定应该是材料坐标的问题,那个应该怎么写啊。。有没有高手知道啊。。这边先谢哈
9楼2009-02-27 16:59:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

胭脂。花开

铁虫 (初入文坛)

这是输出文件后面部分。
  30          8             0        0.018338    9.133052   -2.701974
   31          8             0        0.546647   -4.781648   -2.656322
   32          8             0       -1.670958    5.337122    0.797529
   33          8             0       -1.105972    8.577578   -4.560767
   34          8             0        1.670958   -5.337122   -0.797529
   35          8             0        1.459794    8.205499    0.662393
   36          8             0        2.024780    5.709201   -4.695903
   37          8             0       -1.459794   -8.205499   -0.662393
   38          8             0       -1.541891    8.420620   -0.367472
   39          8             0       -2.106877    5.494080    4.990824
   40          8             0        1.541891   -8.420620    0.367472
   41          7             0       -0.281888    7.052248    3.726695
   42          7             0        0.283097    6.862452   -1.631601
   43          7             0        0.281888   -7.052248   -3.726695
   44          7             0       -1.172284    9.316448    2.542190
   45          7             0       -0.607299    4.598252   -2.816106
   46          7             0        1.172284   -9.316448   -2.542190
   47          7             0        0.615582    3.328118    4.905413
   48          7             0        1.180568   10.586582   -0.452883
   49          7             0       -0.615582   -3.328118   -4.905413
   50          6             0       -1.053399    5.544173    1.961458
   51          6             0       -0.488413    8.370527   -3.396838
   52          6             0        1.053399   -5.544173   -1.961458
   53          6             0       -0.970139    7.040282    2.569196
   54          6             0       -0.405154    6.874418   -2.789100
   55          6             0        0.970139   -7.040282   -2.569196
   56          6             0       -0.039624    8.270341    4.334326
   57          6             0        0.525361    5.644359   -1.023970
   58          6             0        0.039624   -8.270341   -4.334326
   59          6             0       -0.510218    9.355131    3.671504
   60          6             0        0.054767    4.559569   -1.686792
   61          6             0        0.510218   -9.355131   -3.671504
   62          6             0       -1.424453    8.196037    1.955564
   63          6             0       -0.859468    5.718663   -3.402732
   64          6             0        1.424453   -8.196037   -1.955564
   65          6             0       -2.079988    8.254478    0.629171
   66          6             0       -1.515002    5.660222   -4.729125
   67          6             0        2.079988   -8.254478   -0.629171
   68          6             0        1.503567    2.783497    4.038440
   69          6             0        2.068552   11.131203   -1.319855
   70          6             0       -1.503567   -2.783497   -4.038440
   71          6             0        1.457467    2.636836   -4.669648
   72          6             0        0.892482   11.277864    0.688648
   73          6             0       -1.457467   -2.636836    4.669648
   74          6             0        1.967679    1.313826   -4.474284
   75          6             0        1.402693   12.600874    0.884012
   76          6             0       -1.967679   -1.313826    4.474284
   77          6             0       -2.703442    1.475515    4.235411
   78          6             0       -2.138456   12.439185   -1.122884
   79          6             0        2.703442   -1.475515   -4.235411
   80          6             0        1.639607    0.731357    5.312536
   81          6             0        2.204593   13.183343   -0.045760
   82          6             0       -1.639607   -0.731357   -5.312536
   83          1             0        1.836470    3.322830    3.299639
   84          1             0        2.401455   10.591870   -2.058657
   85          1             0       -1.836470   -3.322830   -3.299639
   86          1             0        0.879427    3.047319   -3.955494
   87          1             0        0.314441   10.867381    1.402802
   88          1             0       -0.879427   -3.047319    3.955494
   89          1             0        1.690512    0.809836   -3.683293
   90          1             0        1.125526   13.104864    1.675003
   91          1             0       -1.690512   -0.809836    3.683293
   92          1             0       -2.043309    1.127091    3.628638
   93          1             0       -1.478323   12.787609   -1.729658
   94          1             0        2.043309   -1.127091   -3.628638
   95          1             0        0.396476    8.293161    5.143964
   96          1             0        0.961462    5.621539   -0.214332
   97          1             0       -0.396476   -8.293161   -5.143964
   98          1             0       -0.337489   10.296878    4.050872
   99          1             0        0.227497    3.617822   -1.307424
  100          1             0        0.337489  -10.296878   -4.050872
  101         -2             0        4.717900    0.000000    0.000000
  102         -2             0        0.000000   27.829400    0.000000
  103         -2             0       -1.129971    0.000000   10.716592
---------------------------------------------------------------------
Lengths of translation vectors:      4.717900   27.829400   10.776000
  Angles of translation vectors:     90.000000   96.019100   90.000000
---------------------------------------------------------------------
Standard basis: 6-31+G(d) (6D, 7F)
Integral buffers will be    262144 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned off.
  1612 basis functions,  2868 primitive gaussians,  1624 cartesian basis functions
   308 alpha electrons      308 beta electrons
       nuclear repulsion energy     14284.6201146382 Hartrees.
NAtoms=  100 NActive=  100 NUniq=  100 SFac= 7.50D-01 NAtFMM=   80 NAOKFM=T Big=T
FOutLm=  100.00.
Periodicity:                           1               1               1
Max integer dimensions:               12               2               5
   PBC vector 1 X=     8.9155 Y=    0.0000 Z=    0.0000
   PBC vector 2 X=     0.0000 Y=   52.5899 Z=    0.0000
   PBC vector 3 X=    -2.1353 Y=    0.0000 Z=   20.2514
  Recp vector 1 X=     0.1122 Y=    0.0000 Z=    0.0118
  Recp vector 2 X=     0.0000 Y=    0.0190 Z=    0.0000
  Recp vector 3 X=     0.0000 Y=    0.0000 Z=    0.0494
Generated k point mesh (from -Pi to Pi):
K space mesh: X=      32 Y=       6 Z=      14
A half-cell shift: 1
Using k point mesh (from -Pi to Pi):
K space mesh: X=      32 Y=       6 Z=      14
A half-cell shift: 1
CountK=T Total number of k points:       0
CountK=T Total number of k points:    1344
MDV < 1 in IniBc0.
Error termination via Lnk1e in /home/gauss/g03/l302.exe at Thu Jul  3 11:20:24 2008.
Job cpu time:  0 days  0 hours  6 minutes 55.2 seconds.
File lengths (MBytes):  RWF=    778 Int=      0 D2E=      0 Chk=      1 Scr=      1
2楼2009-02-26 15:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

★ ★ ★ ★ ★ ★
胭脂。花开(金币+2,VIP+0):高手。很感谢。还希望再帮忙多看看。谢谢。 2-26 15:24
胭脂。花开(金币+2,VIP+0):高手。很感谢。还希望再帮忙多看看。谢谢。 2-26 15:25
胭脂。花开(金币+2,VIP+0):高手。很感谢。还希望再帮忙多看看。谢谢。 2-26 15:25
b3lyp/6-31+g(d) pbc=nkpoints=1 pop=(readradii,chelpg)


呵呵,好像头一次见这样的命令
3楼2009-02-26 15:03:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胭脂。花开

铁虫 (初入文坛)

引用回帖:
Originally posted by erylingjet at 2009-2-26 15:03:
b3lyp/6-31+g(d) pbc=nkpoints=1 pop=(readradii,chelpg)


呵呵,好像头一次见这样的命令

版主大人见笑了。我现在Gauss 处于小白阶段。
我就是想放一个元胞算电荷,用chelpg方法,b3lyp这个泛函,6-31+g(d)基组。还请帮忙看看怎么写输入文件。
超级感谢了。。
4楼2009-02-26 15:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见