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