Periodic RESP charges for paracetamol form I crystal
mbellucci
mbell... at gmail.com
Sat Mar 1 20:02:12 UTC 2014
Dear CP2K Users,
I am trying to calculate the RESP charges for a unit cell of the form I
crystal of paracetamol using periodic boundary conditions using version 2.5
of CP2K. I have tried many different things but in every calculation, the
RESP charges seem to be unphysical, i.e. the magnitudes of the charges seem
to very large and most have the wrong sign. The things I have tried in
order to fix this include:
increasing the cutoff
changing the functional
increasing the size of the basis set
changing the values of the LENGTH and RANGE keywords in the PERIODIC_SYS
section of RESP
changing the SURF_DIRECTION in the PERIODIC_SYS section of RESP
I've also experimented with the keywords WIDTH, RESTRAIN_HEAVIES_STRENGTH,
and RESTRAIN_HEAVIES_TO_ZERO in RESP
Nothing seems to work for the RESP charges. However, when I turn off the
periodic boundary conditions the calculations produce RESP charges that
make physical sense. Also, for calculations with or without periodic
boundary conditions, the DDAP charges make physical sense, but I'm hesitant
to use them in MD simulations since they are less frequently used in MD
simulations in the literature. If anyone could help identify an error I am
making when calculating these RESP charges, I would greatly appreciate it!
Also, if anyone could comment on the reliability of DDAP charges in
general, I would appreciate it as well.
For the calculation, the coordinates of the form I crystal as well as the
dimensions of the unit cell were obtained from the cambridge structural
database. The input file I'm using, the outputted RESP charges, and the
outputted DDAP charges are the following:
&GLOBAL
PROJECT form_1
RUN_TYPE ENERGY_FORCE
PRINT_LEVEL MEDIUM
WALLTIME 86400
&END GLOBAL
&FORCE_EVAL
METHOD Quickstep
&DFT
BASIS_SET_FILE_NAME EMSL_BASIS_SETS
POTENTIAL_FILE_NAME POTENTIAL
CHARGE 0
MULTIPLICITY 1
&MGRID
CUTOFF 300
REL_CUTOFF 50
&END MGRID
&QS
METHOD GAPW
EPS_DEFAULT 1.0E-12
MAP_CONSISTENT TRUE
EXTRAPOLATION ASPC
EXTRAPOLATION_ORDER 3
&END QS
&SCF
MAX_SCF 500
SCF_GUESS ATOMIC
EPS_SCF 1.0E-7
&OT ON
MINIMIZER DIIS
PRECONDITIONER FULL_ALL
ENERGY_GAP 0.001
&END OT
&OUTER_SCF
EPS_SCF 1.0E-5
MAX_SCF 100
&END OUTER_SCF
&END SCF
&XC
&XC_FUNCTIONAL PBE
&END XC_FUNCTIONAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 12.93 9.40 7.10
ALPHA_BETA_GAMMA 90.00 115.90 90.00
PERIODIC XYZ
&END CELL
&COORD
C 9.510000 3.334000 5.417000
C 10.770000 3.348000 4.827000
C 11.771000 2.497000 5.288000
C 11.518000 1.636000 6.343000
C 10.266000 1.639000 6.953000
C 9.279000 2.483000 6.489000
C 8.231000 4.763000 3.830000
C 6.944000 5.542000 3.720000
H 10.930000 3.957000 4.145000
H 12.675000 2.557000 4.918000
H 10.101000 1.062000 7.690000
H 8.407000 2.482000 6.891000
H 12.348000 0.517000 7.562000
H 7.724000 4.239000 5.582000
H 6.614000 5.480000 2.836000
H 6.512000 5.574000 4.484000
H 7.181000 6.589000 3.692000
N 8.426000 4.170000 5.017000
O 12.536000 0.800000 6.745000
O 9.035000 4.678000 2.901000
C 6.784000 8.034000 0.970000
C 5.524000 8.048000 1.560000
C 4.523000 7.197000 1.099000
C 4.776000 6.336000 0.044000
C 6.028000 6.339000 -0.567000
C 7.015000 7.183000 -0.102000
C 8.063000 9.463000 2.557000
C 9.350000 10.242000 2.667000
H 5.364000 8.657000 2.242000
H 3.618000 7.257000 1.469000
H 6.193000 5.762000 -1.303000
H 7.887000 7.182000 -0.505000
H 3.945000 5.217000 -1.175000
H 8.570000 8.939000 0.805000
H 9.680000 10.180000 3.551000
H 9.782000 10.274000 1.903000
H 9.113000 11.289000 2.695000
N 7.867000 8.870000 1.370000
O 3.758000 5.500000 -0.358000
O 7.259000 9.378000 3.486000
C 0.319000 6.066000 0.970000
C -0.941000 6.052000 1.560000
C -1.942000 6.903000 1.099000
C -1.689000 7.764000 0.044000
C -0.437000 7.761000 -0.567000
C 0.550000 6.917000 -0.102000
C 1.598000 4.637000 2.557000
C 2.885000 3.858000 2.667000
H -1.101000 5.443000 2.242000
H -2.847000 6.843000 1.469000
H -0.272000 8.338000 -1.303000
H 1.422000 6.918000 -0.505000
H -2.520000 8.883000 -1.175000
H 2.105000 5.161000 0.805000
H 3.215000 3.920000 3.551000
H 3.317000 3.826000 1.903000
H 2.648000 2.811000 2.695000
N 1.402000 5.230000 1.370000
O -2.707000 8.600000 -0.358000
O 0.794000 4.722000 3.486000
C 3.045000 1.366000 5.417000
C 4.305000 1.352000 4.827000
C 5.306000 2.203000 5.288000
C 5.053000 3.064000 6.343000
C 3.801000 3.061000 6.953000
C 2.814000 2.217000 6.489000
C 1.766000 -0.063000 3.830000
C 0.479000 -0.842000 3.720000
H 4.465000 0.743000 4.145000
H 6.210000 2.143000 4.918000
H 3.636000 3.638000 7.690000
H 1.942000 2.218000 6.891000
H 5.883000 4.183000 7.562000
H 1.259000 0.461000 5.582000
H 0.149000 -0.780000 2.836000
H 0.047000 -0.874000 4.484000
H 0.716000 -1.889000 3.692000
N 1.961000 0.530000 5.017000
O 6.071000 3.900000 6.745000
O 2.570000 0.022000 2.901000
&END COORD
&KIND H
BASIS_SET 6-31++Gxx
POTENTIAL GTH-PBE-q1
&END KIND
&KIND C
BASIS_SET 6-31++Gxx
POTENTIAL GTH-PBE-q4
&END KIND
&KIND N
BASIS_SET 6-31++Gxx
POTENTIAL GTH-PBE-q5
&END KIND
&KIND O
BASIS_SET 6-31++Gxx
POTENTIAL GTH-PBE-q6
&END KIND
&END SUBSYS
&PROPERTIES
&RESP
&PERIODIC_SYS
ATOM_LIST 1..80
LENGTH 3.0
RANGE 3.0 5.0
SURF_DIRECTION Z
&END PERIODIC_SYS
&PRINT
&V_RESP_CUBE
&END V_RESP_CUBE
&RESP_CHARGES_TO_FILE
&END RESP_CHARGES_TO_FILE
&PROGRAM_RUN_INFO
&END
&COORD_FIT_POINTS
&END
&END PRINT
&END RESP
&FIT_CHARGE
&END FIT_CHARGE
&END PROPERTIES
&END FORCE_EVAL
RESP charges:
Type | Atom | Charge
RESP 1 C 0.070078
RESP 2 C 0.903292
RESP 3 C 1.113722
RESP 4 C -0.120125
RESP 5 C 1.058323
RESP 6 C 0.655164
RESP 7 C 0.200902
RESP 8 C 2.666416
RESP 9 H -1.222912
RESP 10 H -1.170830
RESP 11 H -0.983979
RESP 12 H -0.818745
RESP 13 H -1.214944
RESP 14 H -1.366249
RESP 15 H -1.224385
RESP 16 H -1.119378
RESP 17 H -1.137280
RESP 18 N 1.259042
RESP 19 O 1.713353
RESP 20 O 0.738109
RESP 21 C 0.041413
RESP 22 C 0.956097
RESP 23 C 1.040182
RESP 24 C -0.082712
RESP 25 C 1.051482
RESP 26 C 0.595410
RESP 27 C 0.202160
RESP 28 C 2.662283
RESP 29 H -1.166123
RESP 30 H -1.126516
RESP 31 H -0.984189
RESP 32 H -0.800838
RESP 33 H -1.188989
RESP 34 H -1.401436
RESP 35 H -1.215709
RESP 36 H -1.113423
RESP 37 H -1.128041
RESP 38 N 1.271056
RESP 39 O 1.679915
RESP 40 O 0.705805
RESP 41 C 0.116651
RESP 42 C 0.813652
RESP 43 C 1.108160
RESP 44 C -0.087035
RESP 45 C 1.066794
RESP 46 C 0.598212
RESP 47 C 0.201601
RESP 48 C 2.665972
RESP 49 H -1.193708
RESP 50 H -1.119339
RESP 51 H -0.987625
RESP 52 H -0.794927
RESP 53 H -1.193920
RESP 54 H -1.403616
RESP 55 H -1.227170
RESP 56 H -1.113105
RESP 57 H -1.132964
RESP 58 N 1.272201
RESP 59 O 1.669491
RESP 60 O 0.737908
RESP 61 C 0.017005
RESP 62 C 0.961044
RESP 63 C 1.064429
RESP 64 C -0.116028
RESP 65 C 1.044203
RESP 66 C 0.665757
RESP 67 C 0.203199
RESP 68 C 2.660164
RESP 69 H -1.159143
RESP 70 H -1.166104
RESP 71 H -0.983609
RESP 72 H -0.830205
RESP 73 H -1.208912
RESP 74 H -1.380887
RESP 75 H -1.210223
RESP 76 H -1.114933
RESP 77 H -1.129467
RESP 78 N 1.264612
RESP 79 O 1.718065
RESP 80 O 0.706396
Total 0.000000
DDAP FULL DENSITY charges:
Atom | Charge
1 C 0.164768
2 C 0.038789
3 C -0.067304
4 C 0.359095
5 C -0.273509
6 C -0.132861
7 C 0.453357
8 C -0.073176
9 H -0.008341
10 H -0.007149
11 H 0.160840
12 H 0.127025
13 H 0.128967
14 H 0.094082
15 H -0.056650
16 H -0.027392
17 H -0.004180
18 N -0.249484
19 O -0.299466
20 O -0.327283
21 C 0.164803
22 C 0.039171
23 C -0.067954
24 C 0.359126
25 C -0.273851
26 C -0.132660
27 C 0.453316
28 C -0.072910
29 H -0.008512
30 H -0.006834
31 H 0.161121
32 H 0.126815
33 H 0.128910
34 H 0.094007
35 H -0.056723
36 H -0.027522
37 H -0.004288
38 N -0.249623
39 O -0.299372
40 O -0.327146
41 C 0.164786
42 C 0.039184
43 C -0.067967
44 C 0.359125
45 C -0.273855
46 C -0.132651
47 C 0.453282
48 C -0.072882
49 H -0.008520
50 H -0.006829
51 H 0.161118
52 H 0.126809
53 H 0.128911
54 H 0.094000
55 H -0.056734
56 H -0.027529
57 H -0.004294
58 N -0.249597
59 O -0.299364
60 O -0.327128
61 C 0.164787
62 C 0.038773
63 C -0.067297
64 C 0.359102
65 C -0.273505
66 C -0.132867
67 C 0.453392
68 C -0.073203
69 H -0.008330
70 H -0.007154
71 H 0.160838
72 H 0.127032
73 H 0.128964
74 H 0.094090
75 H -0.056640
76 H -0.027383
77 H -0.004176
78 N -0.249513
79 O -0.299471
80 O -0.327306
Total 0.000000
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20140301/1b67e43f/attachment.htm>
More information about the CP2K-user
mailing list