<div dir="ltr">Dear Matt,<div><br></div><div>Thank you for your correction! Now the script finally works.</div><div><br></div><div>For this test system, the cut off distance is set very small, due to the short lattice parameter in c direction (2.5A).</div><div>For a real run, I will use a cutoff distance at least 6A (< 0.5*L_cell).</div><div><br></div><div>If dispersion correction is needed in the XC section, I would try to use</div><div><br></div><div>--------------------------------------------------------------------</div><div><div>&VDW_POTENTIAL</div><div>     POTENTIAL_TYPE                             PAIR_POTENTIAL</div><div>     &PAIR_POTENTIAL</div><div>         TYPE                                               DFTD3 </div><div>         PARAMETER_FILE_NAME            dftd3.dat </div><div>         REFERENCE_FUNCTIONAL         B3LYP</div><div>         #CALCULATE_C9_TERM              TRUE</div><div>         #LONG_RANGE_CORRECTION   FALSE</div><div>         R_CUTOFF [angstrom]                  16.0</div><div>     &END PAIR_POTENTIAL </div><div>&END VDW_POTENTIAL</div></div><div>-------------------------------------------------------------</div><div><br></div><div><br></div><div>I wish this piece of script is appropriate for the pe.inp script.</div><div><br></div><div><br></div><div>Besides, I find that for this periodic system, I cannot use</div><div>  POTENTIAL_TYPE   COULOMB</div><div>in the HF section, which results in very long OT optimization process and never converges due to my test.</div><div>Is it the case that periodic Hartree Fock calculation must use truncated potential?<br></div><div><br></div><div>Thank you for your time.</div><div><br></div><div>Regard,</div><div>Yiyang</div><div><br><br>On Tuesday, October 25, 2016 at 12:49:09 PM UTC-7, Matt W wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Hi,</div><div><br></div>you need to set the FIRST_SNAPSHOT variable to 1 (or higher). <div><br></div><div>BTW, your settings for exact exchange in the B3LYP functional are not good - you truncate the coulomb interaction very hard, but don't include a long range correction.</div><div><br></div><div>Matt<br></div></div></blockquote></div></div>