BSSE with ECP


Click here for full thread
Clicked A Few Times
Hi,
I'm trying to perform BSSE on Gd-DOTA complex with lanl2dz and MWB-53 pseudopotential.
I'm using DFT level

After fist SCF cycle, I obtain this error
'*** The number of electrons 53 in the ECP does not
accord with a known set of completely filled shells !!!
------------------------------------------------------------------------
unrecognised pseudopotential 911'


how I could fix it?

regards

Marcel

echo

start Gd-DOTA-1WAT
memory 800 MW
scratch_dir /tmp/tmp_nwchem

title "Gd-DOTA 'one water' lanl2dz - Gd MWB 53 "
charge -1

geometry units angstroms print xyz autosym
N                    -2.06016503     0.47199520     1.13075948
C -2.53653619 -0.82671569 1.72745356
C -1.37948836 -1.70652258 2.20700787
N -0.44288518 -2.10763149 1.10052874
C 0.86184266 -2.58421054 1.68411392
C 1.74187419 -1.42825232 2.16842760
N 2.13985976 -0.48062609 1.06985023
C 2.63158794 0.81380556 1.66232659
C 1.48826638 1.68979734 2.18173698
N 0.51827185 2.09998646 1.10788034
C -0.76829030 2.57113898 1.73401013
C -1.63552060 1.41010971 2.22799849
C -3.14515471 1.07970558 0.29372962
C -2.64243404 1.97820234 -0.85104242
O -1.40577938 1.78940165 -1.24347266
O -3.46244072 2.79358381 -1.38484009
C -1.05115393 -3.18141624 0.24955295
C -1.94520563 -2.66897314 -0.89433342
O -2.77063004 -3.47588183 -1.42834346
O -1.73765759 -1.43379111 -1.29140960
C 3.19981001 -1.08544504 0.19965754
C 2.65795645 -1.96893961 -0.93917273
O 3.44596380 -2.80514215 -1.48808792
O 1.42055228 -1.74408369 -1.31175049
C 1.09781594 3.17885605 0.24359620
C 1.94248252 2.66706040 -0.93745770
O 2.74084831 3.47817045 -1.50808014
O 1.72757362 1.43036446 -1.31903511
H -3.69826881 0.27068157 -0.19581476
H -3.85892168 1.64596346 0.90819275
H -1.62069570 -3.90008276 0.85513042
H -0.24298121 -3.73151249 -0.24458853
H 3.92549632 -1.66190288 0.79005562
H 3.74465507 -0.27523987 -0.29729373
H 0.27362386 3.73669083 -0.21329078
H 1.69551080 3.88919226 0.83143532
H -3.11595558 -1.35462570 0.96836729
H -3.20978755 -0.63376410 2.57706123
H -1.79341763 -2.60365740 2.69307928
H -0.79986599 -1.17464994 2.96539948
H 0.67746885 -3.26404261 2.53034952
H 1.38644703 -3.15716209 0.91813562
H 2.64102769 -1.84514211 2.64844776
H 1.21089172 -0.85610629 2.93320016
H 3.32885118 0.61585727 2.49121097
H 3.18828262 1.34793436 0.89060463
H 0.93150086 1.15198871 2.95318334
H 1.91711641 2.58302270 2.66216288
H -1.31315376 3.15345203 0.98982623
H -0.55879470 3.24004928 2.58315981
H -1.08364633 0.83046978 2.97218757
H -2.52258611 1.82051808 2.73524248
Gd 0.01945991 0.00548459 -0.48501874
H -0.69790624 0.90074484 -3.23989474
O -0.35710650 -0.00479794 -3.07613145
H -1.12118268 -0.62870948 -3.03273133
end
BASIS "ao basis" spherical PRINT
Gd S
7.9250430 -0.1522640
6.2717590 0.4939230
2.5909130 -1.2945530
Gd S
0.6425090 1.0
Gd S
0.3128260 1.0
Gd S
0.0584800 1.0
Gd S
0.0269240 1.0
Gd P
5.2050550 0.1141430
3.1110210 -0.4426100
0.7437650 1.1665320
Gd P
0.3127020 1.0
Gd P
0.0943750 1.0
Gd P
0.0336140 1.0
Gd D
2.5342090 -0.0394870
0.8067320 0.3908310
0.3203550 0.7016280
Gd D
0.1242070 1.0
Gd D
0.0480910 1.0
  1. BASIS SET: (4s) -> [2s]
H S
    19.2384000              0.0328280        
2.8987000 0.2312040
0.6535000 0.8172260
H S
     0.1776000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
C S
  4233.0000000              0.0012200        
634.9000000 0.0093420
146.1000000 0.0454520
42.5000000 0.1546570
14.1900000 0.3588660
5.1480000 0.4386320
1.9670000 0.1459180
C S
     5.1480000             -0.1683670        
0.4962000 1.0600910
C S
     0.1533000              1.0000000        
C P
    18.1600000              0.0185390        
3.9860000 0.1154360
1.1430000 0.3861880
0.3594000 0.6401140
C P
     0.1146000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
N S
  5909.0000000              0.0011900        
887.5000000 0.0090990
204.7000000 0.0441450
59.8400000 0.1504640
20.0000000 0.3567410
7.1930000 0.4465330
2.6860000 0.1456030
N S
     7.1930000             -0.1604050        
0.7000000 1.0582150
N S
     0.2133000              1.0000000        
N P
    26.7900000              0.0182540        
5.9560000 0.1164610
1.7070000 0.3901780
0.5314000 0.6371020
N P
     0.1654000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
O S
  7817.0000000              0.0011760        
1176.0000000 0.0089680
273.2000000 0.0428680
81.1700000 0.1439300
27.1800000 0.3556300
9.5320000 0.4612480
3.4140000 0.1402060
O S
     9.5320000             -0.1541530        
0.9398000 1.0569140
O S
     0.2846000              1.0000000        
O P
    35.1800000              0.0195800        
7.9040000 0.1242000
2.3050000 0.3947140
0.7171000 0.6273760
O P
     0.2137000              1.0000000 
bqGd S
7.9250430 -0.1522640
6.2717590 0.4939230
2.5909130 -1.2945530
bqGd S
0.6425090 1.0
bqGd S
0.3128260 1.0
bqGd S
0.0584800 1.0
bqGd S
0.0269240 1.0
bqGd P
5.2050550 0.1141430
3.1110210 -0.4426100
0.7437650 1.1665320
bqGd P
0.3127020 1.0
bqGd P
0.0943750 1.0
bqGd P
0.0336140 1.0
bqGd D
2.5342090 -0.0394870
0.8067320 0.3908310
0.3203550 0.7016280
bqGd D
0.1242070 1.0
bqGd D
0.0480910 1.0
  1. BASIS SET: (4s) -> [2s]
bqH S
    19.2384000              0.0328280        
2.8987000 0.2312040
0.6535000 0.8172260
bqH S
     0.1776000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
bqC S
  4233.0000000              0.0012200        
634.9000000 0.0093420
146.1000000 0.0454520
42.5000000 0.1546570
14.1900000 0.3588660
5.1480000 0.4386320
1.9670000 0.1459180
bqC S
     5.1480000             -0.1683670        
0.4962000 1.0600910
bqC S
     0.1533000              1.0000000        
bqC P
    18.1600000              0.0185390        
3.9860000 0.1154360
1.1430000 0.3861880
0.3594000 0.6401140
bqC P
     0.1146000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
bqN S
  5909.0000000              0.0011900        
887.5000000 0.0090990
204.7000000 0.0441450
59.8400000 0.1504640
20.0000000 0.3567410
7.1930000 0.4465330
2.6860000 0.1456030
bqN S
     7.1930000             -0.1604050        
0.7000000 1.0582150
bqN S
     0.2133000              1.0000000        
bqN P
    26.7900000              0.0182540        
5.9560000 0.1164610
1.7070000 0.3901780
0.5314000 0.6371020
bqN P
     0.1654000              1.0000000        
  1. BASIS SET: (10s,5p) -> [3s,2p]
bqO S
  7817.0000000              0.0011760        
1176.0000000 0.0089680
273.2000000 0.0428680
81.1700000 0.1439300
27.1800000 0.3556300
9.5320000 0.4612480
3.4140000 0.1402060
bqO S
     9.5320000             -0.1541530        
0.9398000 1.0569140
bqO S
     0.2846000              1.0000000        
bqO P
    35.1800000              0.0195800        
7.9040000 0.1242000
2.3050000 0.3947140
0.7171000 0.6273760
bqO P
     0.2137000              1.0000000  
end
ecp
Gd nelec 53
Gd ul
2 1.000000000 0.000000000
Gd S
2 5.0260000 139.60161900
2 2.5130000 -6.87969400
Gd P
2 4.3022000 79.87382400
2 2.1511000 0.93841900
Gd D
2 3.0011000 32.70972900
2 1.5055000 1.24195200
Gd F
2 6.2049000 -74.24879000
bqGd nelec 53
bqGd ul
2 1.000000000 0.000000000
bqGd S
2 5.0260000 139.60161900
2 2.5130000 -6.87969400
bqGd P
2 4.3022000 79.87382400
2 2.1511000 0.93841900
bqGd D
2 3.0011000 32.70972900
2 1.5055000 1.24195200
bqGd F
2 6.2049000 -74.24879000
end
cosmo
solvent water
parameters cosmo.par
end
dft
vectors input hcore output Gd-DOTA-1WAT.movecs
mult 1
xc b3lyp
disp vdw 3
maxiter 600
convergence energy 1d-7
convergence damp 33 lshift .3
grid fine
end
bsse
mon Gd 53
charge 3 mult 1
mon water 54 55 56
charge 0 mult 1
mon DOTA 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
charge -4 mult 1
end

task dft energy