Thanks, It worked.<br><br>在 2012年10月21日星期日UTC+8下午2时15分21秒，sandeep写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div>Hi,</div><div> Few comments that might help you in getting right structure. </div><div><br></div><div>1. Your cell parameters are very small. You need to take bigger cell and repeat the calculations. </div><div><br>
</div><div>take at least 2 x 2 x 3 cell. (or better choice is 3 x 3 x 4)</div><div><br></div><div>2. Your cutoff value is less. Try increase to 280 Ry or above. or If you are not sure about cut off value, pl. do set of single point calculations by varying density cutoff. Check for energy convergence. </div>
<div><br></div><div><br></div><div>Wishes,</div><div>Sandeep</div><div><br></div><br><div class="gmail_quote">On Sun, Oct 21, 2012 at 9:52 AM, <a href="javascript:" target="_blank" gdf-obfuscated-mailto="7pVo4vIZGQwJ">jjf...@yahoo.com.cn</a> <span dir="ltr"><<a href="javascript:" target="_blank" gdf-obfuscated-mailto="7pVo4vIZGQwJ">jjf...@yahoo.com.cn</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi, energyone <br><br>I want use cp2k to optimzie the cell of TiO2, The following is the input file:<br><br>&FORCE_EVAL<br>
STRESS_TENSOR ANALYTICAL<br> METHOD Quickstep<br> &DFT<br> BASIS_SET_FILE_NAME /scratch/cp2k/BASIS_MOLOPT<br> POTENTIAL_FILE_NAME /scratch/cp2k/GTH_POTENTIALS<br> &MGRID<br> CUTOFF 200<br> &END MGRID<br>
&QS<br> METHOD GPW<br> &END QS<br> &SCF<br> EPS_SCF 5.E-7<br> MAX_SCF 500<br> &OT<br> MINIMIZER CG<br> &END<br> &END SCF<br> &XC<br> &XC_FUNCTIONAL PBE<br>
&END XC_FUNCTIONAL<br> &END XC<br> &END DFT<br> &SUBSYS<br> &CELL<br> ABC 4.594 4.594 2.959<br> PERIODIC XYZ<br> &END CELL<br> &TOPOLOGY<br> &END TOPOLOGY<br> &COORD<br>
Ti -0.0000 -0.0000 0.0000<br>Ti 2.2970 2.2970 1.4795<br>O 1.4003 1.4003 0.0000<br>O 3.1937 3.1937 0.0000<br>O 0.8967 3.6973 1.4795<br>O 3.6973 0.8967 1.4795<br> &END COORD<br>
&KIND TI<br> BASIS_SET DZVP-MOLOPT-SR-GTH<br> POTENTIAL GTH-PBE-q12<br> &END KIND<br> &KIND O<br> BASIS_SET DZVP-MOLOPT-SR-GTH<br> POTENTIAL GTH-PBE-q6<br> &END KIND<br> &END SUBSYS<br>
&END FORCE_EVAL<br> &MOTION<br> &CELL_OPT<br> KEEP_ANGLES .TRUE.<br> OPTIMIZER CG<br> &CG<br> &LINE_SEARCH<br> TYPE 2PNT<br> &END LINE_SEARCH<br> &END CG<br> &END CELL_OPT<br>
&GEO_OPT<br> OPTIMIZER BFGS<br> &END GEO_OPT<br> &END MOTION<br>&GLOBAL<br> PROJECT TiO2<br> RUN_TYPE CELL_OPT<br> PRINT_LEVEL LOW<br>&END GLOBAL<br><br><br>When completing, It give the cell parameters:<br>
5.54 5.58 2.96<br>and the Coordinate:<br>Ti 0.4095385769 0.2898267710 -0.0352843895<br>Ti 3.1801930298 2.1516147238 1.4539613109<br> O 1.9175831463 1.4494611116 -0.0449395166<br>
O 3.4936489849 3.0632060642 -0.0346941550<br> O 0.7216799169 4.9591730884 1.4533505831<br> O 4.6832547080 0.9868972431 1.4635018888<br><br>The a and b is remarkable larger than exp. value. More importantly, Some Ti-O bond are broken, see the fig. below:<br>
<img alt=""><br>Anyone can give suggestion about how to do it?<span><font color="#888888"><br>
<p></p>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msg/cp2k/-/svvIEFrsLBMJ" target="_blank">https://groups.google.com/d/<wbr>msg/cp2k/-/svvIEFrsLBMJ</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="7pVo4vIZGQwJ">cp...@googlegroups.com</a>.<br>
To unsubscribe from this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="7pVo4vIZGQwJ">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
For more options, visit this group at <a href="http://groups.google.com/group/cp2k?hl=en" target="_blank">http://groups.google.com/<wbr>group/cp2k?hl=en</a>.<br>
</font></span></blockquote></div><br>
</blockquote>