LOCALIZATION method??

lml... at gmail.com lml... at gmail.com
Wed May 16 00:06:42 UTC 2007


As for LOCALIZATION part, there are four kinds of methods: JACOBI|
CRAZY|DIRECT|L1SD. I tested for several system, L1SD converges very
fast for any system, but gives the wrong results. Do I make any
mistake for this method, or is there some bug for L1SD method?

Thanks in advance!

Limin
-------------------------------------------------------------------
Results for H2O Dimer:
Crazy: MOLECULAR_DIPOLES
                           # molecule nr, charge, dipole vector,
dipole (Debye)
     1    0.000000   -1.278306    0.048275   -1.840541    2.241425
     2    0.000000   -1.125128   -0.183263    1.829523    2.155610

L1SD: MOLECULAR_DIPOLES
                           # molecule nr, charge, dipole vector,
dipole (Debye)
     1    0.000000   -6.676929   -0.192657   -0.914185    6.741975
     2    0.000000    4.422305    0.081061    0.928242    4.519401

My input file:

&FORCE_EVAL
  METHOD Quickstep
  &DFT
    &PRINT
      &MOLECULAR_STATES
      &END MOLECULAR_STATES
      &LOCALIZATION
        LOC_TYPE L1SD  #CRAZY
       EPS_LOCALIZATION 1.0E-8
      &END LOCALIZATION
    &END PRINT
    &MGRID
      CUTOFF 280
    &END MGRID
    &QS
      EXTRAPOLATION USE_GUESS
#      EXTRAPOLATION_ORDER 2
    &END QS
    &SCF
      SCF_GUESS ATOMIC
    &END SCF
    &XC
      &XC_FUNCTIONAL BLYP
      &END XC_FUNCTIONAL
      &XC_GRID
        XC_DERIV SPLINE2
        XC_SMOOTH_RHO NONE
      &END XC_GRID
    &END XC
  &END DFT
  &SUBSYS
    &CELL
      ABC 15 15 15
      UNIT ANGSTROM
    &END CELL
    &COORD
  O        -1.3459376534       -0.0354156443        0.1465219667
  H        -1.6005239833        0.7471984224       -0.3684408375
  H        -1.6119945737       -0.7853664624       -0.4092127268
  O         1.5271489754        0.0518925950       -0.1437561542
 H         0.5752834905        0.0079790198        0.0847062186
  H         1.9800438583        0.0028822544        0.7108287421

&END COORD
    &KIND O
      BASIS_SET TZV2P-GTH-q6
      POTENTIAL GTH-BLYP-q6
    &END KIND
    &KIND H
      BASIS_SET TZV2P-GTH-q1
      POTENTIAL GTH-BLYP-q1
    &END KIND
    &TOPOLOGY
       CONN_FILE 2h2o.psf
       CONNECTIVITY UPSF
    &END TOPOLOGY
  &END SUBSYS
&END FORCE_EVAL
&GLOBAL
  FFTLIB FFTSG
  PROJECT water_dipole
  PRINT_LEVEL LOW




More information about the CP2K-user mailing list