# 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)
>
>
> [..]
>
>
>
>
>
> > 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-
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
# 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
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
POSITION     7.5000000000000000E+00
DIRECTION  WALL_MINUS
K     4.0000000000000001E-02
&END WALL
&END METAVAR
&PRINT
&TEMPERATURE_COLVAR  SILENT
COMMON_ITERATION_LEVELS               1
&EACH
&END EACH
&END TEMPERATURE_COLVAR
&COLVAR  SILENT
COMMON_ITERATION_LEVELS               1
&EACH
&END EACH
&END COLVAR
&HILLS  SILENT
COMMON_ITERATION_LEVELS               1
&EACH
&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 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
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