[CP2K:3128] Re: Metadynamics Bug report in colvar: DISTANCE

Teodoro Laino teodor... at gmail.com
Wed Mar 2 22:01:09 UTC 2011


Dear Ata,

finally the situation is clear : the problem is not the CV itself, rather the extended lagrangian scheme, which does not take into account the PBC.

Thanks for pointing this out, we will fix this problem ASAP.
The correct fix is anyway not to comment the line containing NINT, rather taking care of PBC in the extended lagrangian scheme.

Just for completeness: if you run a non-lagrangian MTD (which is what 90% of the people do nowadays) you will not see any jump and everything will run OK.

Thanks again and regards,
Teo

On Mar 2, 2011, at 10:53 PM, Ata Roudgar wrote:

> 
> 
> 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
> 
> 
> -- 
> You received this message because you are subscribed to the Google Groups "cp2k" group.
> To post to this group, send email to cp... at googlegroups.com.
> To unsubscribe from this group, send email to cp2k+uns... at googlegroups.com.
> For more options, visit this group at http://groups.google.com/group/cp2k?hl=en.
> 




More information about the CP2K-user mailing list