<div dir="ltr">Hi Matt,<div><br></div><div>Thanks for the reply. The way it is specified in CP2K is quite confusing to me, to be honest. I am used to unit cell + translation vectors + # repetitions combination. Incidentally, this setup was actually suggested by another CP2K user from here, and it yielded per-atom nickel cohesive energy quite close to the experimental value (-4.422 eV from my calculations). The same for aluminum yielded -2.9 eV, which isn't great, but also close.</div><div>Now I wonder how that's possible, if the setup is wildly incorrect...<br><br>On Monday, September 29, 2014 2:20:06 AM UTC-6, 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">Hi,<div><br></div><div>you aren't running bullk Ni, but a Ni atom in a large box. See the definition of multiple unit cell in the CELL section</div><div><br></div><div><div>MULTIPLE_UNIT_CELL {Integer} {Integer} {Integer}</div><div>Specifies the numbers of repetition in space (X, Y, Z) of the defined cell, assuming it as a unit cell. This keyword affects only the CELL specification. The same keyword in SUBSYS%TOPOLOGY%MULTIPLE_UNIT_<wbr>CELL should be modified in order to affect the coordinates specification.  [Edit]</div><div><br></div><div>Why you aren't getting any output, I'm not sure. Maybe the memory cost is too high and your system is swapping?</div><div><br></div><div>Matt</div><br>On Monday, September 22, 2014 11:26:18 PM UTC+1, Alex wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi everyone. I'm trying to optimize a periodic (6x6x6) lattice of nickel using either CELL_OPT or GEO_OPT (input attached). The simulation's been running for several days on 8 cores, no SCF updates, it just hangs at<div><br></div><div><div><font face="courier new, monospace"> *****************************<wbr>******************************<wbr>********************</font></div><div><font face="courier new, monospace"> ***                     STARTING GEOMETRY OPTIMIZATION                      ***</font></div><div><font face="courier new, monospace"> ***                           CONJUGATE GRADIENTS                           ***</font></div><div><font face="courier new, monospace"> *****************************<wbr>******************************<wbr>********************</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Spin 1</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Number of electrons:                                                          9</font></div><div><font face="courier new, monospace"> Number of occupied orbitals:                                                  9</font></div><div><font face="courier new, monospace"> Number of molecular orbitals:                                                26</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Spin 2</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Number of electrons:                                                          9</font></div><div><font face="courier new, monospace"> Number of occupied orbitals:                                                  9</font></div><div><font face="courier new, monospace"> Number of molecular orbitals:                                                26</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Number of orbital functions:                                                 26</font></div><div><font face="courier new, monospace"> Number of independent orbital functions:                                     26</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> Extrapolation method: initial_guess</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"> SCF WAVEFUNCTION OPTIMIZATION</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">  Step     Update method      Time    Convergence         Total energy    Change</font></div><div><font face="courier new, monospace">  ------------------------------<wbr>------------------------------<wbr>------------------</font></div></div><div><br></div><div><br></div><div> BFGS minimizer doesn't seem like an option with this system size, so I used CG... Am I doing something seriously wrong?</div><div><br></div><div>Thank you.</div></div></blockquote></div></div></blockquote></div></div>