<div dir="ltr"><div>Dear all CP2K users,</div><div><br></div><div><br></div><div>I'm doing some exercises by following "exercises:2016_summer_school:geometry_and_cell_optimization(<a href="https://www.cp2k.org/exercises:2016_summer_school:geometry_and_cell_optimization">https://www.cp2k.org/exercises:2016_summer_school:geometry_and_cell_optimization</a>)".</div><div><br></div><div>What I want to know from this optimazation calculation is the functional that fits with the experimental data.</div><div><br></div><div>Now here are some questions.</div><div><br></div><div>1.</div><div>On the example page, there are 3steps.</div><div>1) Geometry optimize with Fist(MM Force Field) with each "NREP" for "MULTIPLE_UNIT_CELL ${NREP} ${NREP} ${NREP}"</div><div>2) Energy calculation with QS(DFT) with each "SCALING_FACTOR" for "ABC 5.620*${SCALE_FACTOR} 5.620*${SCALE_FACTOR} 5.620*${SCALE_FACTOR}"</div><div>3) Cell optimization with QS(DFT)</div><div><br></div><div>From the results, 1) shows no big defferences with MULTIPLE_UNIT_CELL. Should I do  this?</div><div>Would doing it with QS(DFT) not with Fist(MM Force Field) works?</div><div><br></div><div>I just used the input on the page and fixed it to do the 3) CELL_OPT calculation. It contains the KEEP_SYMMETRY option.</div><div>But the result shows that the length of the each nearest atoms are not equal.</div><div>What happened? Did I do something wrong? (input and output files are attached)</div><div><br></div><div>From 2) I can get the proper scaleing factor by caculating energy of the unitcell with various lattice parameter.</div><div>But I'm not sure that this scaled structure is same as the optimized structure because of the bad results of 3).</div><div>How can I check it?</div><div><br></div><div><br></div><div><br></div><div>2.</div><div>I tried with 2) Energy calculation and 3) Cell optimization with both optimizer BFGS and LBFGS.</div><div>From the output, BFGS and LBFGS run times were very similar.</div><div>Does optimizer have no effects on calculating energy and cell optimization?</div><div>Followings are the runtime of the Energy calculation.</div><div>----------------------------------------------------------------------------------------------------------------------------</div><div>OPTIMIZER : BFGS, NREP : 1, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 15:28:03.704<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 15:28:13.963</div><div>   0:10<br>OPTIMIZER : BFGS, NREP : 2, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 15:28:21.649<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 15:29:04.595</div><div>   0:37<br>OPTIMIZER : BFGS, NREP : 3, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 15:29:13.594<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 15:32:48.711</div><div>   3:35<br>OPTIMIZER : BFGS, NREP : 4, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 15:33:49.742<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 15:47:38.693</div><div>   13:89<br></div><div>OPTIMIZER : LBFGS, NREP : 1, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 16:15:59.167<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 16:16:11.632</div><div>   0:12<br>OPTIMIZER : LBFGS, NREP : 2, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 16:16:14.969<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 16:17:20.282</div><div>   1:06<br>OPTIMIZER : LBFGS, NREP : 3, SCALE_FACTOR : 1, CUTOFF : 400 <br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 16:20:41.212<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 16:24:30.740</div><div>   3:49<br>OPTIMIZER : LBFGS, NREP : 3, SCALE_FACTOR : 1, CUTOFF : 400<br>  **** **** ******  **  PROGRAM STARTED AT               2017-09-25 16:25:31.836<br>  **** **** ******  **  PROGRAM ENDED AT                 2017-09-25 16:39:20.152</div><div>   13:49</div><div>----------------------------------------------------------------------------------------------------------------------------</div><div><br></div><div>3. Can I calculate the bandgap of the bulk material  from the unitcell information with CP2K?</div></div>