Metadynamics Bug report in colvar: DISTANCE
Ata Roudgar
rou... at yahoo.com
Wed Mar 2 21:53:40 UTC 2011
On Mar 2, 12:10 pm, Teodoro Laino <teodor... at gmail.com> wrote:
> Dear Ata,
>
> > No Teo, my xj is a fixed point and never change its sign and therefore
> > only xi when it goes far enough it becomes xi --> -xi and therefore
> > (xi-xj) --> (-x1-xj)
>
> sorry this is not the case, since the PBC are applied to the difference xi-xj and not
> to xi and xj independently. See:
>
> > ss=MATMUL(cell%h_inv,xpi-xpj)
>
> anyway.. I try to follow your thoughts:
>
> [..]
>
>
>
>
>
> > now, ss is the transformed version of (Ri-Rj) where Ri=(xi,yi,xi) and
> > Rj=(xj,yj,zj). This is needed because now all ss elements varied from
> > 0 to 1
> > Now after line 2267
> > ss = ss- NINT (ss)
>
> > the ss remains unchanged if its value changes between 0 to 0.5 that
> > means if the distance between two atoms in x direction does not go
> > betound of half of the length of the unitcell in x direction
> > everything is still fine. The smae for y and z direction. The problem
> > occours when the distance between two atoms (in any direction) goes
> > beyond the half of the length of unit cell (in that direction) in that
> > case
>
> > fist, second or third element of ss lets say ss = 0.5+\epsilon after
> > line 2267 the value of ss jump from 0.5+\epsilon to \epsilon because
> > NINT(0.5+x)=x
>
> This is wrong - In fact
> NINT(0.5+\epsilon) = 1
> NINT(0.5-\epsilon)=0
> NINT rounds its argument to the nearest whole number.
yes you are right, but it does not change anything you still have jump
at around 0.5 and that is what I mean.
> ======
>
> Let's make a simple 1D example using the equations you have posted. Let's assume you have a 1D cell of 10.
> Your atom (j) is fixed to position 2.
> Your atom (i) is in position 6.
> The difference ss is equal to 4.
>
> After a short time atom (i) continuously moves to position 7.2 -> we know that the distance is now 5.2 if you do not consider periodic boundary conditions.
> Instead, if you consider PBC the real distance is a bit different - In fact, using your notation:
>
exactly, but PBC is not implied on fictitious CV.
The following input file is my metadynamics input file right before
the jump occurs. You can run this input and and look at the file Test-
metadyn-colvar-1.metadynLog. The first column shows the time, second
and third shows the fictitious and actual CV. They should be around
12.4 Bohr. They will always follow each other because the spring
constant in my metadynamics is relatively strong. After few time step
you will see that the third column jump from 12.875 to 6.469 and
simultaneously the total energy will jump from 463.164 to around
-457.857. If there is no bug in this part of the program please
explain to me why this happens. The ss value right before and after
jump is 0.499 and 0.501. You can choose any cutoff you wish. Please
consider the following input is working for the past 18 ps and there
is no such jump in total energy at all.
Cheers,
Ata
# Version information for this restart file
# current date 2011-02-27 17:10:05.108
# current working dir /global/scratch/aroudgar/metadynamics/at-dcc/
metadynamics/NEW/k-0.3-m-80-w-0.001-YZ/tmp4
# Program compiled at Fri Dec 17
10:47:57 PST 2010
# Program compiled
on bugaboo
# Program compiled for Linux-
x86-64-intel
# Last CVS
entry /lib//
&GLOBAL
PRINT_LEVEL LOW
PROJECT_NAME Test
RUN_TYPE MD
&END GLOBAL
&MOTION
&MD
ENSEMBLE NVT
STEPS 20000
TIMESTEP 4.9999999999999989E-01
STEP_START_VAL 6700
TIME_START_VAL 3.3499999999996744E+03
ECONS_START_VAL 0.0000000000000000E+00
TEMPERATURE 3.0000000000000000E+02
&THERMOSTAT
&NOSE
LENGTH 3
YOSHIDA 3
TIMECON 9.9999999999999986E+01
MTS 2
&COORD
-7.7932446603633299E-01 3.5299788030455950E+01
9.5120091330717464E+00
&END COORD
&VELOCITY
-6.4090980006235532E-06 2.6182433926763283E-04
1.6574644406097236E-04
&END VELOCITY
&MASS
1.5587807263267867E+06 1.6237299232570693E+04
1.6237299232570693E+04
&END MASS
&FORCE
-8.1209715875513483E-09 -5.4566666280671773E-08
1.0041970782827435E-08
&END FORCE
&END NOSE
&END THERMOSTAT
&AVERAGES T
&RESTART_AVERAGES
ITIMES_START 1
AVECPU 2.0259987895522375E+01
AVEHUGONIOT 0.0000000000000000E+00
AVETEMP_BARO 0.0000000000000000E+00
AVEPOT -4.6317693590085298E+02
AVEKIN 4.4643391049105370E-02
AVETEMP 2.9369273103593309E+02
AVEKIN_QM 0.0000000000000000E+00
AVETEMP_QM 0.0000000000000000E+00
AVEVOL 1.2658589352056042E+04
AVECELL_A 2.1987101119477416E+01
AVECELL_B 2.1987101119477416E+01
AVECELL_C 3.0235618126171140E+01
AVEALPHA 9.0000000000000000E+01
AVEBETA 9.0000000000000000E+01
AVEGAMMA 1.2000002726760935E+02
AVE_ECONS -4.0621845777571574E+06
AVE_PRESS 0.0000000000000000E+00
AVE_PXX 0.0000000000000000E+00
&END RESTART_AVERAGES
&END AVERAGES
&END MD
&FREE_ENERGY
&METADYN
NT_HILLS 300
TEMPERATURE 1.0000000000000001E+02
MIN_DISP 1.4999999999999999E-01
OLD_HILL_NUMBER 58
OLD_HILL_STEP 6546
TEMP_TOL 1.0000000000000001E+02
WW 1.0000000000000000E-03
DO_HILLS T
LAGRANGE T
STEP_START_VAL 6701
NHILLS_START_VAL 58
COLVAR_AVG_TEMPERATURE_RESTART 2.3466829732275312E+02
&METAVAR
LAMBDA 2.9999999999999999E-01
MASS 7.9999999999999957E+01
SCALE 1.0000000000000001E-01
COLVAR 1
&WALL
TYPE QUADRATIC
POSITION 7.5000000000000000E+00
&QUADRATIC
DIRECTION WALL_MINUS
K 4.0000000000000001E-02
&END QUADRATIC
&END WALL
&END METAVAR
&PRINT
&TEMPERATURE_COLVAR SILENT
COMMON_ITERATION_LEVELS 1
FILENAME metadyn-temp
&EACH
METADYNAMICS 1
&END EACH
&END TEMPERATURE_COLVAR
&COLVAR SILENT
COMMON_ITERATION_LEVELS 1
FILENAME metadyn-colvar
&EACH
METADYNAMICS 1
&END EACH
&END COLVAR
&HILLS SILENT
COMMON_ITERATION_LEVELS 1
FILENAME metadyn-hills
&EACH
METADYNAMICS 1
&END EACH
&END HILLS
&END PRINT
&SPAWNED_HILLS_POS
9.2037463026715205E+00
9.0519432714782315E+00
8.9009797493567149E+00
8.7490974701248465E+00
8.5978967955333303E+00
8.4466245242519555E+00
8.2956553448546320E+00
8.1444067631395605E+00
8.0319675945167024E+00
7.8806786036045287E+00
7.7305325928669166E+00
7.5804497214034283E+00
7.4769074601723338E+00
7.6269637744000445E+00
7.7787827950570323E+00
7.9290458392419287E+00
8.0808476468441839E+00
8.2312918599376310E+00
8.3830667872535809E+00
8.5347061374312485E+00
8.6856587011131499E+00
8.8374884037431336E+00
8.9881435138506465E+00
9.1386386197760636E+00
9.2898200450113091E+00
9.4406410076979412E+00
9.5914553917365826E+00
9.7416144302800323E+00
9.8920564815469874E+00
1.0043075198551801E+01
1.0092677659270272E+01
9.9411314876890184E+00
9.7899715981862929E+00
9.6389013548177207E+00
9.4870387486160901E+00
9.3362273927757382E+00
9.3119149415735247E+00
9.4621593501258658E+00
9.6136460740207141E+00
9.7656668091063921E+00
9.8638139593661585E+00
1.0006146292427402E+01
1.0156965827472412E+01
1.0307126881087143E+01
1.0457676656080393E+01
1.0607970323009074E+01
1.0751773377754027E+01
1.0902170103710018E+01
1.1053539747576522E+01
1.1204008916441760E+01
1.1354638011511200E+01
1.1505433476918373E+01
1.1656217066386494E+01
1.1808018323310415E+01
1.1958048062272351E+01
1.2108946522669605E+01
1.2259425772459918E+01
1.2410347793704762E+01
&END SPAWNED_HILLS_POS
&SPAWNED_HILLS_SCALE
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
1.0000000000000001E-01
&END SPAWNED_HILLS_SCALE
&SPAWNED_HILLS_HEIGHT
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
1.0000000000000000E-03
&END SPAWNED_HILLS_HEIGHT
&EXT_LAGRANGE_SS0
1.2420202916836748E+01
&END EXT_LAGRANGE_SS0
&EXT_LAGRANGE_VVP
-5.4061182864883049E-05
&END EXT_LAGRANGE_VVP
&END METADYN
&END FREE_ENERGY
&CONSTRAINT
&FIXED_ATOMS
COMPONENTS_TO_FIX XYZ
LIST 1 2 3
&END FIXED_ATOMS
&END CONSTRAINT
&PRINT
&TRAJECTORY SILENT
&EACH
GEO_OPT 1
&END EACH
&END TRAJECTORY
&END PRINT
&PRINT
&TRAJECTORY SILENT
COMMON_ITERATION_LEVELS 1
&END TRAJECTORY
&END PRINT
&END MOTION
&FORCE_EVAL
&DFT
BASIS_SET_FILE_NAME BASIS_SETS
POTENTIAL_FILE_NAME GTH_POTENTIALS
&SCF
MAX_SCF 600
EPS_SCF 1.0000000000000002E-06
EPS_DIIS 1.0000000000000000E-02
SCF_GUESS RESTART
&OT T
MINIMIZER DIIS
&END OT
&END SCF
&QS
METHOD GPW
&END QS
&MGRID
CUTOFF 2.0000000000000000E+02
&END MGRID
&XC
DENSITY_CUTOFF 1.0000000000000000E-10
GRADIENT_CUTOFF 1.0000000000000000E-10
TAU_CUTOFF 1.0000000000000000E-10
&XC_GRID
XC_DERIV SPLINE3
&END XC_GRID
&XC_FUNCTIONAL NO_SHORTCUT
&BECKE88 T
&END BECKE88
&LYP T
&END LYP
&END XC_FUNCTIONAL
&END XC
&END DFT
&SUBSYS
&CELL
A 1.0076269999999999E+01 5.8175340000000002E+00
0.0000000000000000E+00
B -1.0076269999999999E+01 5.8175340000000002E+00
0.0000000000000000E+00
C 0.0000000000000000E+00 0.0000000000000000E+00
1.6000000000000000E+01
MULTIPLE_UNIT_CELL 1 1
1
&END CELL
&COORD
C 7.9880000000000004E+00 5.2590000000000003E+00
1.3000000000000005E+00
C 4.6300010000000000E+00 1.1076000000000004E+01
1.3000000000000005E+00
C 1.2709999999999999E+00 5.2580010000000019E+00
1.3000000000000005E+00
F 9.2226966724628188E+00 4.9805556716761368E+00
7.2225416728690561E-01
F 6.9674989329956896E+00 4.5355675995163027E+00
7.2190957971451930E-01
F 7.6855271444365592E+00 6.5880950230004709E+00
1.1711922900662335E+00
F 4.3602269386555568E+00 1.2333486551590330E+01
8.2818695783094887E-01
F 5.8078259767072753E+00 1.0721620783573288E+01
8.2062919627650444E-01
F 3.6202630737543768E+00 1.0267144632602010E+01
8.1734437544587601E-01
F 1.5765258028677382E+00 6.4757671032161213E+00
8.3030409668983396E-01
F 1.6948176536128772E-02 4.9134308114130709E+00
8.4604301494802581E-01
F 2.1811543014320347E+00 4.2861345228635885E+00
8.4000639789940057E-01
S 4.7482479036855407E+00 1.1172571521586979E+01
3.2373959417324238E+00
S 8.2168129341714291E+00 4.7426222141775050E+00
3.1820580020745837E+00
S 1.4075781203157094E+00 4.9465035650337406E+00
3.1462466817624319E+00
O 9.5502788745047287E+00 5.0650055831565357E+00
3.5903977852428222E+00
O 7.2586332966178713E+00 5.5651200807526164E+00
3.9681177145791113E+00
O 7.8793262610787460E+00 3.2995534520953269E+00
3.0512047957238302E+00
O 4.9503145460018629E+00 5.4524519188509037E+00
4.7169484938001052E+00
O 5.1028457331102652E+00 9.8643354220120276E+00
3.7885051702102679E+00
O 5.8349629599663517E+00 1.2143073205124171E+01
3.3886131073450492E+00
O 3.4264484931362627E+00 1.1709124423423557E+01
3.7288397235166864E+00
O 2.8650557053711059E+00 4.9834964235814736E+00
3.2931624147470857E+00
O 5.5366125873012528E-01 5.8200590903279155E+00
3.9697907251739135E+00
O 8.4976053600938040E-01 3.5187563412120424E+00
3.0901485606405310E+00
O 8.2173321282708933E+00 1.0798627268870031E+00
4.1369435317102159E+00
O 1.3872011714026042E+00 1.2880654699664318E+00
4.1772344335217486E+00
H 7.3048919624611370E+00 7.2626256084752339E-01
3.8069079153820975E+00
H 9.0191423125896630E+00 5.1243655322028336E-01
3.9581205980776644E+00
H 8.2722243544827450E+00 2.0396344711614449E+00
3.7161950202092711E+00
H 1.1953653672980988E+00 2.2810853320956426E+00
3.7450653204402449E+00
H 2.2560194695515978E+00 9.7697117060340422E-01
3.7654416329260143E+00
H 6.5054157983547256E-01 7.2416868324037487E-01
3.8522611258593900E+00
H 5.9869397426689233E+00 5.3312507999716070E+00
4.4115839456702446E+00
H 4.7173070408957134E+00 6.4349959311705529E+00
4.6308949004629243E+00
H 4.1908134051272343E+00 5.0560032567003201E+00
4.1646868827053511E+00
&END COORD
&VELOCITY
0.0000000000000000E+00 0.0000000000000000E+00
0.0000000000000000E+00
0.0000000000000000E+00 0.0000000000000000E+00
0.0000000000000000E+00
0.0000000000000000E+00 0.0000000000000000E+00
0.0000000000000000E+00
-1.6146662355549744E-04 1.4779505784251750E-04
1.1586716996194748E-04
3.1694329641808057E-05 -9.6438148311259813E-05
-3.5242244672038439E-04
-1.5999810385909565E-04 -1.6770896067401615E-05
-6.7935984377535884E-05
9.2976580415781232E-05 3.3231062073675962E-04
-9.1537212833910816E-05
3.6912259551732420E-04 -1.1394763916765328E-04
1.6818294566017949E-04
-2.5085868036771888E-05 1.0898212390356668E-04
-1.0874548263747586E-04
3.7291040039120902E-04 -8.9955246135243807E-05
1.0841076666913192E-04
1.1785214959345341E-04 -1.9716100476346654E-04
-6.8687660868577989E-05
-1.0233905834575182E-04 -2.6183149730414840E-04
-1.6499052868116797E-04
7.3614835446323020E-05 -8.6520116317353961E-05
1.2904544193026454E-04
-2.2527268840117085E-05 1.0043801287136054E-04
7.9153198622836617E-07
4.2256160468752235E-05 -7.0601107896640125E-06
-4.6162931657152544E-05
2.0309810255459651E-04 1.5229331759206605E-04
9.6448719432691273E-05
3.1825779427769679E-04 -2.6541889805121972E-04
4.6799071373241176E-05
1.7861569023351365E-04 3.2065500715914869E-04
-6.8816180449052429E-05
1.6938117540878583E-04 -2.2437336346908790E-04
4.3741094899649310E-05
1.4729869709469105E-04 6.6574364942962287E-05
-9.5928882202402514E-05
1.5618817531899164E-04 -5.3948336682535934E-06
-2.3957296350727427E-05
-9.7923731444223285E-06 2.3869688087500951E-04
-1.9347058118150511E-04
1.4415419054878415E-04 3.0311163055023101E-04
-2.3117486924459675E-04
5.6730620998313021E-05 7.4953007711596917E-05
-6.5822705353396033E-05
1.3407973555425313E-05 -3.1831703531943634E-05
-1.6040793417828106E-04
1.0187809188170929E-04 1.5433000776563738E-04
7.7189800748697637E-05
1.1417116474838422E-04 -6.7124074570798338E-05
1.2896444411374729E-04
-9.5247826723750728E-04 5.5430871122229546E-04
7.2510275465565853E-04
-3.4337906125219547E-04 -8.8968325052729330E-04
-6.9946345217617801E-04
-5.5081071135718105E-04 3.0332929002212776E-04
7.2953690334395316E-04
-5.3118328760707744E-04 -7.6129923838194223E-04
5.5564267204983241E-04
-3.2481493868575227E-04 2.6218631802235148E-04
-1.3863672276590006E-04
2.7133512771225222E-04 4.6531407622418586E-04
5.3971655304983583E-04
-4.5926517713953160E-04 -5.8339528364184396E-04
1.1106336743747319E-04
2.5557125698911061E-04 -1.7545768242913388E-03
-2.3471250294436382E-04
1.6202917735624094E-04 -5.8897834071576808E-04
5.5027483922722451E-04
&END VELOCITY
&KIND C
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q4
&BASIS
1
2 0 2 7 2 2 1
6.132624767898 -0.105576563700 0.024850587600 0.035098108400
-0.056712999000 0.037469425800
2.625196064782 -0.174866621100 0.147440856000 0.114197930900
-0.076743951500 0.071047574800
1.045456957247 0.064464981400 -0.421418379200 0.215908137300
-0.528056757000 0.619330561200
0.478316330874 0.830447035300 0.823051336800 0.355929151000
0.394594045700 0.710730857000
0.178617414302 0.371911373800 0.184880712100 0.271541842000
0.721247475800 0.389148235000
0.075144725465 0.006174464000 -0.742086051600 0.056256557700
0.733902531700 0.346374124100
0.030286753006 0.008966477300 0.083382192000 0.004998059700
0.248172175500 -0.005142430900
&END BASIS
&POTENTIAL
2 2
0.3380660900000000E+00 2 -0.9136268710000000E+01 0.1429259560000000E
+01
2
0.3023222300000000E+00 1 0.9665512280000000E+01
0.2863791200000000E+00 0
&END POTENTIAL
&END KIND
&KIND F
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q7
&BASIS
1
2 0 2 7 2 2 1
15.025672765127 -0.076135699000 0.060516471000 0.043030876500
-0.027382614400 0.006095862800
6.156800581515 -0.145492650600 0.078091919900 0.138591822100
-0.119463647500 0.057880784500
2.411456205564 0.195144281500 -0.051537351700 0.283251574300
-0.152089668600 0.143657436500
0.963986369687 0.553654166300 -0.448874166200 0.357678749700
-0.460859710900 0.782158193600
0.396229721431 0.454888028100 0.143111336600 0.300622862100
0.260540362300 0.349204565600
0.157658905678 0.078338198300 0.433345205900 0.138839114100
0.410677921300 0.106634189000
0.048130255107 0.002175323600 -0.111759214900 0.008740274300
0.051656760200 -0.008599668300
&END BASIS
&POTENTIAL
2 5
0.2138401400000000E+00 2 -0.2190241518000000E+02 0.3275975720000000E
+01
2
0.1945902800000000E+00 1 0.2379658295000000E+02
0.1878632600000000E+00 0
&END POTENTIAL
&END KIND
&KIND S
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q6
&BASIS
1
2 0 2 6 2 2 1
3.817590656735 0.075558928000 -0.045713227100 -0.019219238100
0.001714015800 0.007498037400
2.362751931093 -0.351216488200 0.179094296600 -0.035664835200
0.064608303200 0.104757185100
0.861004180298 -0.099160602500 0.206133423700 0.160834397000
-0.170453754100 0.344818199900
0.417524512433 0.862753116400 -0.737640763800 0.507608206800
-0.536548069500 0.472319101700
0.181513540673 0.307053614900 0.419230078700 0.249563613200
0.786543640200 0.204745694800
0.070570776865 0.001592872500 0.359607313200 0.025974500500
0.740434318200 -0.019952573900
&END BASIS
&POTENTIAL
2 4
0.4200000000000000E+00 1 -0.5988800450000000E+01
2
0.3697326600000000E+00 2 0.1255648450000000E+02 -0.3912904590000000E
+01
0.5051538100000000E+01
0.4128757100000000E+00 1 0.3587721140000000E+01
&END POTENTIAL
&END KIND
&KIND O
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q6
&BASIS
1
2 0 2 7 2 2 1
12.015954705512 -0.060190841200 0.065738617900 0.036543638800
-0.034210557400 0.014807054400
5.108150287385 -0.129597923300 0.110885902200 0.120927648700
-0.120619770900 0.068186159300
2.048398039874 0.118175889400 -0.053732406400 0.251093670300
-0.213719464600 0.290576499200
0.832381575582 0.462964485000 -0.572670666200 0.352639910300
-0.473674858400 1.063344189500
0.352316246455 0.450353782600 0.186760006700 0.294708645200
0.484848376400 0.307656114200
0.142977330880 0.092715833600 0.387201458600 0.173039869300
0.717465919700 0.318346834400
0.046760918300 -0.000255945800 0.003825849600 0.009726110600
0.032498979400 -0.005771736600
&END BASIS
&POTENTIAL
2 4
0.2434202600000000E+00 2 -0.1699189235000000E+02 0.2566142060000000E
+01
2
0.2208314000000000E+00 1 0.1838885102000000E+02
0.2172007000000000E+00 0
&END POTENTIAL
&END KIND
&KIND H
BASIS_SET DZVP-MOLOPT-GTH
POTENTIAL GTH-BLYP-q1
&BASIS
1
2 0 1 7 2 1
11.478000339908 0.024916243200 -0.012512421400 0.024510918200
3.700758562763 0.079825490000 -0.056449071100 0.058140794100
1.446884268432 0.128862675300 0.011242684700 0.444709498500
0.716814589696 0.379448894600 -0.418587548300 0.646207973100
0.247918564176 0.324552432600 0.590363216700 0.803385018200
0.066918004004 0.037148121400 0.438703133000 0.892971208700
0.021708243634 -0.001125195500 -0.059693171300 0.120101316500
&END BASIS
&POTENTIAL
1
0.2000000000000000E+00 2 -0.4195961470000000E+01 0.7304982100000000E
+00
0
&END POTENTIAL
&END KIND
&TOPOLOGY
MULTIPLE_UNIT_CELL 1 1
1
&END TOPOLOGY
&COLVAR
&DISTANCE
ATOMS 2 19
AXIS YZ
&END DISTANCE
&END COLVAR
&END SUBSYS
&END FORCE_EVAL
More information about the CP2K-user
mailing list