回复: [CP2K:4204] SCF Convergence

shoutian sun sunshou... at yahoo.com.cn
Thu Nov 22 03:36:08 UTC 2012


Hi SeokYun,

I have some suggestions for your input file.
 
1. The conformer you used for optimization is very poor, even using atomic unit,
    the C-C distance of the ring is about 3.33 anstrom (kinds be wrong) or 3.33 bohr 
    (about 1.76 angstrom, also a bit larger, isn't it?).  So, my suggestion is rebuilding
    your conformer.
 
2. use the POISSON_SOLVER (I think your systerm is a cluster, right?)
    &POISSON
         PERIODIC none
         POISSON_SOLVER wavelet
    &END POISSON

3. use  MINIMIZER CG for OT sec. (If SCF doesn't converge, sometimes, it works)
 
4. XC_FUNCTIONAL Pade is not a good choice. (BLYP, I always used.)
 
5. sometime, SCF doesn't converge, you have to use smooth method in the &XC section
    , as follows:
      &XC_GRID 
         XC_SMOOTH_RHO NN10
         XC_DERIV SPLINE2_SMOOTH
      &END XC_GRID 

6. cell parameter:      ABC 40.000 40.000 40.000
    I think the default unit is angstrom, right?
    40.0 angstrom for 35 atoms conformer is quite big.
    And very time consuming, Don't you think so?
    (PERIODIC none for clusters.)
 
7. The coordinates of your conformer. I think your comformer
    is build in atomic unit. And you do not modify it and just put
    it in the &COORD.....&END COORD section. Am I right?
    However, the problem is the default unit of atom coordinates
    for cp2k is angstrom, not bohr. So, please modify it, or add
    key word: UNIT Bohr in the following sec.
 
    &COORD
        C    15.02103029        21.43572257     20.18000314
        C    12.13797029        19.86029257     20.10217314
            ..................
       &END COORD
 
8. using POISSON_SOLVER wavelet, you have to use the following key words
    to put the cluster in the center of the box.
 
    &TOPOLOGY
        &CENTER_COORDINATES 
        &END
    &END 

 
I rebuild your comformer and do the optimization with my input file, it converge well.
For detail, please refer to cp2k input reference.
 
Hope, it works for you.
 
Regards,
sun
Chemistry, CUHK, Hong Kong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20121122/27e6f900/attachment.htm>


More information about the CP2K-user mailing list