<div dir="ltr">Hi,<div><br></div><div>what do you mean when you say that the results don't agree? Massively wrong, or a small difference?</div><div><br></div><div>Does your SCF converge? I don't see any very obvious problems with your input at a quick glance.</div><div><br></div><div>Matt<br><br>On Tuesday, November 14, 2017 at 11:43:16 AM UTC+1, Amin wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all, <div>I'm having problem with my CELL_OPT calculations. I tried to achieve lattice parameters for Anatase (101) unit cell but my results are very different than the other results coming from literature. Specially the results mentioned in the reference paper of this exercise: <a href="https://www.cp2k.org/exercises:2015_pitt:gga" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fexercises%3A2015_pitt%3Agga\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFn9XnyeJDzijO7J0BvqoI7OH69rg';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fwww.cp2k.org%2Fexercises%3A2015_pitt%3Agga\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFn9XnyeJDzijO7J0BvqoI7OH69rg';return true;">https://www.cp2k.org/<wbr>exercises:2015_pitt:gga</a></div><div>I tried with/without the repetition of the unit cell but the final lattice parameters was not consistent with those mentioned in literature.</div><div>I would be grateful if someone takes a look to my input file and reminds me what's wrong with my calculations. thanks.</div><div><br></div><div><div>&GLOBAL</div><div>  PROJECT cell3</div><div>  RUN_TYPE CELL_OPT</div><div>  PRINT_LEVEL LOW</div><div>  EXTENDED_FFT_LENGTHS</div><div>&END GLOBAL</div><div><br></div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  STRESS_TENSOR ANALYTICAL</div><div>  &DFT</div><div>    BASIS_SET_FILE_NAME ./BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME ./GTH_POTENTIALS</div><div>    UKS</div><div>    WFN_RESTART_FILE_NAME s-RESTART.wfn</div><div>    MULTIPLICITY 1</div><div>    &QS</div><div>      METHOD GPW</div><div>    &END QS</div><div>    &MGRID</div><div>      CUTOFF 800</div><div>      REL_CUTOFF 60</div><div>    &END MGRID</div><div>    &SCF</div><div>      EPS_SCF 1.0E-6</div><div>      MAX_SCF 150</div><div>      &OUTER_SCF</div><div>        EPS_SCF 1.0E-6</div><div>        MAX_SCF 40</div><div>      &END</div><div>      &OT</div><div>        PRECONDITIONER FULL_SINGLE</div><div>        MINIMIZER CG</div><div>      &END OT</div><div>      SCF_GUESS ATOMIC</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL PBE</div><div>      &END XC_FUNCTIONAL</div><div>      &vdW_POTENTIAL</div><div>        DISPERSION_FUNCTIONAL PAIR_POTENTIAL</div><div>        &PAIR_POTENTIAL</div><div>          TYPE DFTD3(BJ)</div><div>          REFERENCE_FUNCTIONAL PBE</div><div>          CALCULATE_C9_TERM .FALSE.</div><div>          PARAMETER_FILE_NAME ./dftd3.dat</div><div>          R_CUTOFF 15.0</div><div>        &END PAIR_POTENTIAL</div><div>      &END vdW_POTENTIAL</div><div>    &END XC</div><div>    &PRINT</div><div>      &HIRSHFELD</div><div>        SELF_CONSISTENT .TRUE.</div><div>      &END HIRSHFELD</div><div>    &END PRINT</div><div>    &POISSON</div><div>      PERIODIC XYZ</div><div>      POISSON_SOLVER PERIODIC</div><div>    &END POISSON</div><div>  &END DFT</div><div>    </div><div>  &SUBSYS</div><div>    &CELL</div><div>      ABC 10.2394 3.7845 3.513</div><div>      PERIODIC XYZ</div><div>      MULTIPLE_UNIT_CELL 1 3 3 </div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>            MULTIPLE_UNIT_CELL 1 3 3</div><div>    &END TOPOLOGY</div><div>    &COORD</div><div>      @INCLUDE 'cell3.coord'</div><div>    &END COORD</div><div>    &KIND O</div><div>      BASIS_SET DZVP-MOLOPT-GTH</div><div>      POTENTIAL GTH-PBE-q6</div><div>    &END KIND</div><div>    &KIND Ti</div><div>      BASIS_SET DZVP-MOLOPT-SR-GTH</div><div>      POTENTIAL GTH-PBE-q12</div><div>    &END KIND</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&MOTION</div><div>  &GEO_OPT</div><div>    TYPE MINIMIZATION</div><div>    OPTIMIZER BFGS</div><div>  &END GEO_OPT</div><div>  &CELL_OPT</div><div>          KEEP_SYMMETRY TRUE</div><div>          KEEP_ANGLES   TRUE</div><div>  &END        </div><div>&END MOTION</div></div><div><br></div><div>cell3.coord:</div><div><br></div><div><div>O     7.692165e-01 -3.336155e-17  1.170579e+00    0.07512 -0.00000  0.33321</div><div>O     5.888893e+00  1.892250e+00  1.170579e+00    0.57512  0.50000  0.33321</div><div>Ti    4.049520e+00  1.892250e+00  1.902225e+00    0.39549  0.50000  0.54148</div><div>Ti    9.169196e+00 -3.976756e-16  1.902225e+00    0.89549 -0.00000  0.54148</div><div>O     3.678746e+00 -1.595502e-16  2.049708e+00    0.35928 -0.00000  0.58346</div><div>O     8.798422e+00  1.892250e+00  2.049708e+00    0.85928  0.50000  0.58346</div><div>O     2.210147e+00  1.892250e+00  2.633872e+00    0.21585  0.50000  0.74975</div><div>O     7.329824e+00 -3.179005e-16  2.633872e+00    0.71585 -0.00000  0.74975</div><div>Ti    1.839373e+00 -7.977512e-17  2.781355e+00    0.17964 -0.00000  0.79173</div><div>Ti    6.959049e+00  1.892250e+00  2.781355e+00    0.67964  0.50000  0.79173</div><div>O     0.000000e+00  0.000000e+00  3.513001e+00    0.00000  0.00000  1.00000</div><div>O     5.119676e+00  1.892250e+00  3.513001e+00    0.50000  0.50000  1.00000</div></div><div><br></div></div></blockquote></div></div>