<div dir="ltr"><div>Hi <span style="color:rgb(31,31,31);font-size:0.875rem;font-family:"Google Sans",Roboto,RobotoDraft,Helvetica,Arial,sans-serif;white-space:nowrap">Jürg, </span></div><div><span style="color:rgb(31,31,31);font-size:0.875rem;font-family:"Google Sans",Roboto,RobotoDraft,Helvetica,Arial,sans-serif;white-space:nowrap"><br></span></div><div><font color="#1f1f1f" face="Google Sans, Roboto, RobotoDraft, Helvetica, Arial, sans-serif"><span style="font-size:14px;white-space:nowrap">Okay, I use a larger one, say 80Ry. That is 0.025 for PW_CUTOFF here in CP2K. However, as I tried to run it, it stopped with unknown errors</span></font></div><div><font color="#1f1f1f" face="Google Sans, Roboto, RobotoDraft, Helvetica, Arial, sans-serif"><span style="font-size:14px;white-space:nowrap"><br></span></font></div><div><font color="#1f1f1f" face="Google Sans, Roboto, RobotoDraft, Helvetica, Arial, sans-serif"><span style="font-size:14px;white-space:nowrap">... ... </span></font></div><div><br></div>XC functionals<br>==============<br>1) XC_GGA_X_PBE: Perdew, Burke & Ernzerhof<br>J. P. Perdew, K. Burke, and M. Ernzerhof., Â Phys. Rev. Lett. 77, 3865 (1996) (10.1103/PhysRevLett.77.3865)<br>J. P. Perdew, K. Burke, and M. Ernzerhof., Â Phys. Rev. Lett. 78, 1396 (1997) (10.1103/PhysRevLett.78.1396)<br><br>2) XC_GGA_C_PBE: Perdew, Burke & Ernzerhof<br>J. P. Perdew, K. Burke, and M. Ernzerhof., Â Phys. Rev. Lett. 77, 3865 (1996) (10.1103/PhysRevLett.77.3865)<br>J. P. Perdew, K. Burke, and M. Ernzerhof., Â Phys. Rev. Lett. 78, 1396 (1997) (10.1103/PhysRevLett.78.1396)<br><br><br>memory consumption<br>==================<br>approximate number of G+k vectors Â  Â  Â  Â : 0<br>approximate number of G vectors Â  Â  Â  Â  Â : 0<br>approximate number of coarse G vectors Â  : 0<br>approximate size of wave-functions for each k-point: 0 Mb, 0 Mb/rank<br>approximate memory consumption of Davidson solver: 0 Mb/rank<br>approximate memory consumption of FFT transforms: 0 Mb/rank<br><br>===================================================================================<br>= Â  BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES<br>= Â  RANK 0 PID 334186 RUNNING AT ubuntua<br>= Â  KILLED BY SIGNAL: 9 (Killed)<br>===================================================================================<br><br><div><br></div><div>&GLOBAL<br>  PROJECT FCC_Copper<br>  RUN_TYPE CELL_OPT<br>  PRINT_LEVEL LOW<br>&END GLOBAL<br><br>&FORCE_EVAL<br>  METHOD SIRIUS<br>  STRESS_TENSOR ANALYTICAL<br>  &SUBSYS<br>  Â  &CELL<br>  Â  Â  A 3.61 0.00 0.00<br>  Â  Â  B 0.00 3.61 0.00<br>  Â  Â  C 0.00 0.00 3.61<br>  Â  Â  SYMMETRY CUBIC<br>  Â  Â  PERIODIC XYZ<br>  Â  &END CELL<br>  Â  &COORD<br>  Â  Â  Cu Â  0.000 0.000 0.000<br>  Â  Â  Cu Â  0.000 1.805 1.805<br>  Â  Â  Cu Â  1.805 0.000 1.805<br>  Â  Â  Cu Â  1.805 1.805 0.000<br>  Â  &END COORD<br>  Â  &KIND Cu<br>  Â  Â  Â POTENTIAL GTH-PBE-q11<br>  Â  &END KIND<br>  &END SUBSYS<br>  &PW_DFT<br>  Â  &CONTROL<br>  Â  Â  Â VERBOSITY 2<br>  Â  &END CONTROL<br>  Â  &PARAMETERS<br>  Â  Â  Â ELECTRONIC_STRUCTURE_METHOD Â pseudopotential<br>  Â  Â  Â GK_CUTOFF 0.012<br>  Â  Â  Â PW_CUTOFF 0.025<br>  Â  Â  Â NUM_DFT_ITER 399<br>  Â  Â  Â NGRIDK 10 10 10<br>  Â  &END PARAMETERS<br>  Â  &ITERATIVE_SOLVER<br>  Â  Â  Â ENERGY_TOLERANCE 1e-5<br>  Â  Â  Â NUM_STEPS 20<br>  Â  Â  Â SUBSPACE_SIZE 4<br>  Â  &END ITERATIVE_SOLVER<br>  Â  &MIXER<br>  Â  Â  Â BETA 0.3<br>  Â  Â  Â TYPE BROYDEN2<br>  Â  &END MIXER<br>  &END PW_DFT<br>  &DFT<br>  Â  &XC<br>  Â  Â  &XC_FUNCTIONAL<br>  Â  Â  Â  &GGA_X_PBE<br>  Â  Â  Â  &END GGA_X_PBE<br>  Â  Â  Â  &GGA_C_PBE<br>  Â  Â  Â  &END GGA_C_PBE<br>  Â  Â  &END XC_FUNCTIONAL<br>  Â  &END XC<br>  &END DFT<br>&END FORCE_EVAL<br><br>&MOTION<br>  &CELL_OPT<br>  Â  MAX_ITER 200<br>  Â  OPTIMIZER BFGS<br>  Â  KEEP_SYMMETRY .TRUE.<br>  &END CELL_OPT<br>&END MOTION<br></div><div><br></div><div>Regards</div><div><br></div><div>Jibiao Li</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 29 Jun 2023 at 19:39, Jürg Hutter <<a href="mailto:hutter@chem.uzh.ch">hutter@chem.uzh.ch</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi<br>
<br>
where are these numbers from? To me 41 Ry for a 11 electron Cu PP seems very low.<br>
According to the manual units for PW_CUTOFF and GK_CUTOFF are [a.u.^-1]<br>
<br>
JH<br>
<br>
________________________________________<br>
From: <a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a> <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> on behalf of Jibiao Li <<a href="mailto:science.sweden@gmail.com" target="_blank">science.sweden@gmail.com</a>><br>
Sent: Thursday, June 29, 2023 1:32 PM<br>
To: <a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a><br>
Subject: Re: [CP2K:19033] Terrible lattice constant of FCC copper from cell optimization<br>
<br>
Hi Jürg,<br>
<br>
For cell optimization for lattice constants, 500 eV or 41 Ry is sufficient. The problem is how should  I convert eV or Ry to the unit of PW_CUTOFF.<br>
<br>
May I get  an answer for it?<br>
<br>
Best<br>
<br>
Jibiao Li<br>
<br>
<br>
On Thu, 29 Jun 2023 at 19:12, Jürg Hutter <<a href="mailto:hutter@chem.uzh.ch" target="_blank">hutter@chem.uzh.ch</a><mailto:<a href="mailto:hutter@chem.uzh.ch" target="_blank">hutter@chem.uzh.ch</a>>> wrote:<br>
Hi<br>
<br>
removing the basis doesn't change the calculation (in this case).<br>
The SIRIUS output is in Bohr (6.82) the input was in Angstrom (3.61).<br>
<br>
Have you investigated if these parameters are adequate for your system/pseudopotential?<br>
<br>
  Â  Â  Â GK_CUTOFF 5<br>
  Â  Â  Â PW_CUTOFF 20<br>
  Â  Â  Â NGRIDK 10 10 10<br>
<br>
regards<br>
JH<br>
<br>
________________________________________<br>
From: <a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a><mailto:<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>> <<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a><mailto:<a href="mailto:cp2k@googlegroups.com" target="_blank">cp2k@googlegroups.com</a>>> on behalf of Jibiao Li <<a href="mailto:science.sweden@gmail.com" target="_blank">science.sweden@gmail.com</a><mailto:<a href="mailto:science.sweden@gmail.com" target="_blank">science.sweden@gmail.com</a>>><br>
Sent: Thursday, June 29, 2023 12:36 PM<br>
To: cp2k<br>
Subject: Re: [CP2K:19031] Terrible lattice constant of FCC copper from cell optimization<br>
<br>
Hi  Jürg,<br>
<br>
I have removed basis set and found that the problem remains. A difficult thing hard to understand : Why a different value (6.8219113397) was taken in the very beginning of the calculations? There must be somthing hard to understand.<br>
<br>
The lattice prior to the optimization:<br>
lattice vectors<br>
  a1 :  Â  Â  Â 6.8219113397  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000<br>
  a2 :  Â  Â  Â 0.0000000000  Â  Â  Â 6.8219113397  Â  Â  Â 0.0000000000<br>
  a3 :  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000  Â  Â  Â 6.8219113397<br>
<br>
<br>
<br>
The head of The output file:<br>
<br>
 DBCSR| CPU Multiplication driver  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â XSMM (U)<br>
 DBCSR| Multrec recursion limit  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  512 (U)<br>
 DBCSR| Multiplication stack size  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 1000 (D)<br>
 DBCSR| Maximum elements for images  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  UNLIMITED (U)<br>
 DBCSR| Multiplicative factor virtual images  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 1 (U)<br>
 DBCSR| Use multiplication densification  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â T (D)<br>
 DBCSR| Multiplication size stacks  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 3 (U)<br>
 DBCSR| Use memory pool for CPU allocation  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â F (U)<br>
 DBCSR| Number of 3D layers  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â SINGLE (U)<br>
 DBCSR| Use MPI memory allocation  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  F (U)<br>
 DBCSR| Use RMA algorithm  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  F (U)<br>
 DBCSR| Use Communication thread  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â T (U)<br>
 DBCSR| Communication thread load  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 42 (D)<br>
 DBCSR| MPI: My process id  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 0<br>
 DBCSR| MPI: Number of processes  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 8<br>
 DBCSR| OMP: Current number of threads  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  13<br>
 DBCSR| OMP: Max number of threads  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  13<br>
 DBCSR| Split modifier for TAS multiplication algorithm  Â  Â  Â  Â  Â  Â  Â  Â  1.0E+00 (U)<br>
<br>
SIRIUS 7.3.2, git hash: <a href="https://api.github.com/repos/electronic-structure/SIRIUS/git/ref/tags/v7.3.2" rel="noreferrer" target="_blank">https://api.github.com/repos/electronic-structure/SIRIUS/git/ref/tags/v7.3.2</a><br>
<br>
  **** **** ******  **  PROGRAM STARTED AT  Â  Â  Â  Â  Â  Â  Â 2023-06-29 18:25:42.062<br>
 ***** ** ***  *** **  Â PROGRAM STARTED ON  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ubuntua<br>
 **  Â  ****  Â ******  Â  PROGRAM STARTED BY  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  jibiaoli<br>
 ***** **  Â  ** ** **  Â PROGRAM PROCESS ID  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  333548<br>
  **** **  *******  **  PROGRAM STARTED IN  Â  Â  Â  Â  Â  Â  /home/jibiaoli/calc/test<br>
<br>
 CP2K| version string:  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â CP2K version 2023.1<br>
 CP2K| source code revision number:  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  git:b888bd8<br>
 CP2K| cp2kflags: omp fftw3 libxc pexsi elpa parallel scalapack cosma xsmm max_c<br>
 CP2K|  Â  Â  Â  Â  Â  ontr=4 plumed2 spglib mkl sirius libvori libbqb libvdwxc<br>
 CP2K| is freely available from  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  <a href="https://www.cp2k.org/" rel="noreferrer" target="_blank">https://www.cp2k.org/</a><br>
 CP2K| Program compiled at  Â  Â  Â  Â  Â  2023å¹´ 06月 28æ—¥ æ˜ŸæœŸä¸‰ 23:59:21 CST<br>
 CP2K| Program compiled on  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ubuntua<br>
 CP2K| Program compiled for  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â Linux-intel-x86_64<br>
 CP2K| Data directory path  Â  Â  Â  Â  Â  Â  Â  Â /home/jibiaoli/codes/cp2k-2023.1/data<br>
 CP2K| Input file name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  copper.inp<br>
<br>
 GLOBAL| Force Environment number  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  1<br>
 GLOBAL| Basis set file name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â BASIS_SET<br>
 GLOBAL| Potential file name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â POTENTIAL<br>
 GLOBAL| MM Potential file name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â MM_POTENTIAL<br>
 GLOBAL| Coordinate file name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  __STD_INPUT__<br>
 GLOBAL| Method name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  CP2K<br>
 GLOBAL| Project name  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â FCC_Copper<br>
 GLOBAL| Run type  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â CELL_OPT<br>
 GLOBAL| FFT library  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â FFTW3<br>
 GLOBAL| Diagonalization library  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ELPA<br>
 GLOBAL| DGEMM library  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  BLAS<br>
 GLOBAL| Minimum number of eigenvectors for ELPA usage  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  16<br>
 GLOBAL| Orthonormality check for eigenvectors  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  DISABLED<br>
 GLOBAL| Matrix multiplication library  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â COSMA<br>
 GLOBAL| All-to-all communication in single precision  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  F<br>
 GLOBAL| FFTs using library dependent lengths  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  F<br>
 GLOBAL| Grid backend  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â AUTO<br>
 GLOBAL| Global print level  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  LOW<br>
 GLOBAL| MPI I/O enabled  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â T<br>
 GLOBAL| Total number of message passing processes  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 8<br>
 GLOBAL| Number of threads for this process  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 13<br>
 GLOBAL| This output is from process  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 0<br>
 GLOBAL| Stack size for threads created by OpenMP (OMP_STACKSIZE)  Â  Â  Â  default<br>
 GLOBAL| CPU model name  Â  Â  Â  Â  Â  Intel(R) Xeon(R) Platinum 8171M CPU @ 2.30GHz<br>
 GLOBAL| CPUID  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  1003<br>
<br>
 MEMORY| system memory details [Kb]<br>
 MEMORY|  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  rank 0  Â  Â  Â  Â  Â min  Â  Â  Â  Â  Â max  Â  Â  Â average<br>
 MEMORY| MemTotal  Â  Â  Â  Â  Â  Â 98585440  Â  Â  98585440  Â  Â  98585440  Â  Â  98585440<br>
 MEMORY| MemFree  Â  Â  Â  Â  Â  Â  57901124  Â  Â  57901124  Â  Â  57901124  Â  Â  57901124<br>
 MEMORY| Buffers  Â  Â  Â  Â  Â  Â  Â  521896  Â  Â  Â  521896  Â  Â  Â  521896  Â  Â  Â  521896<br>
 MEMORY| Cached  Â  Â  Â  Â  Â  Â  Â 29522272  Â  Â  29522272  Â  Â  29522272  Â  Â  29522272<br>
 MEMORY| Slab  Â  Â  Â  Â  Â  Â  Â  Â  2499004  Â  Â  Â 2499004  Â  Â  Â 2499004  Â  Â  Â 2499004<br>
 MEMORY| SReclaimable  Â  Â  Â  Â  1413276  Â  Â  Â 1413276  Â  Â  Â 1413276  Â  Â  Â 1413276<br>
 MEMORY| MemLikelyFree  Â  Â  Â  89358568  Â  Â  89358568  Â  Â  89358568  Â  Â  89358568<br>
<br>
<br>
 GENERATE|  Preliminary Number of Bonds generated:  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â 0<br>
 GENERATE|  Achieved consistency in connectivity generation.<br>
<br>
 *******************************************************************************<br>
 *******************************************************************************<br>
 **  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â **<br>
 **  Â #### ##  Â ## ####  ##### #####  Â  Â  Â  Â  Â #### ## ####  ## ## ##  ####  Â **<br>
 **  Â ## ####  Â ## ## ## ###  Â #####  Â  Â  Â  Â  ##  Â  ## ##  # ## ## ## ##  Â  Â  **<br>
 **  Â #### ## # ## ## ## #####  ###  Â  Â  Â  Â  Â ##### ## ####  ## ## ## #####  Â **<br>
 **  Â ##  Â ####### ## ## ##  Â  Â ###  Â  Â  Â  Â  Â  Â  ## ## ## #  ## #####  Â  ##  Â **<br>
 **  Â ##  Â ## # ## ####  ##  Â  Â ###  Â  Â  Â  Â  Â ####  ## ##  # ## ##### ####  Â  **<br>
 **  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â **<br>
 **  Â SIRIUS: Domain specific library for electronic structure calculations  Â **<br>
 **  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â **<br>
 *******************************************************************************<br>
<br>
SIRIUS version : 7.3.2<br>
git hash  Â  Â  Â : <a href="https://api.github.com/repos/electronic-structure/SIRIUS/git/ref/tags/v7.3.2" rel="noreferrer" target="_blank">https://api.github.com/repos/electronic-structure/SIRIUS/git/ref/tags/v7.3.2</a><br>
git branch  Â  Â : release v7.3.2<br>
build time  Â  Â : 2023-06-28 23:11:48<br>
start time  Â  Â : Thu, 29 Jun 2023 18:25:42<br>
<br>
number of MPI ranks  Â  Â  Â  Â  Â : 8<br>
MPI grid  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  : 1 1<br>
maximum number of OMP threads : 13<br>
number of MPI ranks per node  : 8<br>
page size (Kb)  Â  Â  Â  Â  Â  Â  Â  : 4<br>
number of pages  Â  Â  Â  Â  Â  Â  Â : 24646360<br>
available memory (GB)  Â  Â  Â  Â : 94<br>
<br>
FFT context for density and potential<br>
=====================================<br>
  comm size  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 8<br>
  plane wave cutoff  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 20.000000<br>
  grid size  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 48 48 48  Â total : 110592<br>
  grid limits  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : -23 24  Â -23 24  Â -23 24<br>
  number of G-vectors within the cutoff : 21466<br>
  local number of G-vectors  Â  Â  Â  Â  Â  Â : 2672<br>
  number of G-shells  Â  Â  Â  Â  Â  Â  Â  Â  Â  : 395<br>
<br>
FFT context for coarse grid<br>
=====================================<br>
  comm size  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 1<br>
  plane wave cutoff  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 10.000000<br>
  grid size  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : 24 24 24  Â total : 13824<br>
  grid limits  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â : -11 12  Â -11 12  Â -11 12<br>
  number of G-vectors within the cutoff : 2689<br>
  local number of G-vectors  Â  Â  Â  Â  Â  Â : 337<br>
  number of G-shells  Â  Â  Â  Â  Â  Â  Â  Â  Â  : 100<br>
<br>
number of local G-vector blocks: 1<br>
<br>
Unit cell<br>
--------------------------------------------------------------------------------<br>
lattice vectors<br>
  a1 :  Â  Â  Â 6.8219113397  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000<br>
  a2 :  Â  Â  Â 0.0000000000  Â  Â  Â 6.8219113397  Â  Â  Â 0.0000000000<br>
  a3 :  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000  Â  Â  Â 6.8219113397<br>
reciprocal lattice vectors<br>
  b1 :  Â  Â  Â 0.9210300448  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000<br>
  b2 :  Â  Â  Â 0.0000000000  Â  Â  Â 0.9210300448  Â  Â  Â 0.0000000000<br>
  b3 :  Â  Â  Â 0.0000000000  Â  Â  Â 0.0000000000  Â  Â  Â 0.9210300448<br>
<br>
unit cell volume :  Â  Â  Â 317.48134574 [a.u.^3]<br>
1/sqrt(omega)  Â  :  Â  Â  Â  Â 0.05612300<br>
MT volume  Â  Â  Â  : 0.000000 ( 0.00%)<br>
IT volume  Â  Â  Â  : 317.481346 (100.00%)<br>
<br>
<br>
&GLOBAL<br>
  PROJECT FCC_Copper<br>
  RUN_TYPE CELL_OPT<br>
  PRINT_LEVEL LOW<br>
&END GLOBAL<br>
<br>
&FORCE_EVAL<br>
  METHOD SIRIUS<br>
  STRESS_TENSOR ANALYTICAL<br>
  &SUBSYS<br>
  Â  &CELL<br>
  Â  Â  A 3.61 0.00 0.00<br>
  Â  Â  B 0.00 3.61 0.00<br>
  Â  Â  C 0.00 0.00 3.61<br>
  Â  Â  SYMMETRY CUBIC<br>
  Â  Â  PERIODIC XYZ<br>
  Â  &END CELL<br>
  Â  &COORD<br>
  Â  Â  Cu  Â 0.000 0.000 0.000<br>
  Â  Â  Cu  Â 0.000 1.805 1.805<br>
  Â  Â  Cu  Â 1.805 0.000 1.805<br>
  Â  Â  Cu  Â 1.805 1.805 0.000<br>
  Â  &END COORD<br>
  Â  &KIND Cu<br>
  Â  Â  Â POTENTIAL GTH-PBE-q11<br>
  Â  &END KIND<br>
  &END SUBSYS<br>
  &PW_DFT<br>
  Â  &CONTROL<br>
  Â  Â  Â VERBOSITY 2<br>
  Â  &END CONTROL<br>
  Â  &PARAMETERS<br>
  Â  Â  Â ELECTRONIC_STRUCTURE_METHOD  pseudopotential<br>
  Â  Â  Â GK_CUTOFF 5<br>
  Â  Â  Â PW_CUTOFF 20<br>
  Â  Â  Â NUM_DFT_ITER 399<br>
  Â  Â  Â NGRIDK 10 10 10<br>
  Â  &END PARAMETERS<br>
  Â  &ITERATIVE_SOLVER<br>
  Â  Â  Â ENERGY_TOLERANCE 1e-5<br>
  Â  Â  Â NUM_STEPS 20<br>
  Â  Â  Â SUBSPACE_SIZE 4<br>
  Â  &END ITERATIVE_SOLVER<br>
  Â  &MIXER<br>
  Â  Â  Â BETA 0.3<br>
  Â  Â  Â TYPE BROYDEN2<br>
  Â  &END MIXER<br>
  &END PW_DFT<br>
  &DFT<br>
  Â  &XC<br>
  Â  Â  &XC_FUNCTIONAL<br>
  Â  Â  Â  &GGA_X_PBE<br>
  Â  Â  Â  &END GGA_X_PBE<br>
  Â  Â  Â  &GGA_C_PBE<br>
  Â  Â  Â  &END GGA_C_PBE<br>
  Â  Â  &END XC_FUNCTIONAL<br>
  Â  &END XC<br>
  &END DFT<br>
&END FORCE_EVAL<br>
<br>
&MOTION<br>
  &CELL_OPT<br>
  Â  MAX_ITER 200<br>
  Â  OPTIMIZER BFGS<br>
  Â  KEEP_SYMMETRY .TRUE.<br>
  &END CELL_OPT<br>
&END MOTION<br>
<br>
Best<br>
<br>
Jibiao Li<br>
On Thursday, 29 June 2023 at 16:24:04 UTC+8 Jürg Hutter wrote:<br>
Hi<br>
<br>
you are doing a plane wave calculation, so you don't have to specify a Gaussian basis set.<br>
However, the crucial parameter is now the PW cutoff. Make sure you have a reasonable value,<br>
either by testing or by comparing to literature.<br>
<br>
regards<br>
JH<br>
<br>
________________________________________<br>
From: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><mailto:<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>> <<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><mailto:<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>>> on behalf of Jibiao Li <<a href="mailto:science...@gmail.com" target="_blank">science...@gmail.com</a><mailto:<a href="mailto:science...@gmail.com" target="_blank">science...@gmail.com</a>>><br>
Sent: Thursday, June 29, 2023 3:23 AM<br>
To: <a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a><mailto:<a href="mailto:cp...@googlegroups.com" target="_blank">cp...@googlegroups.com</a>><br>
Subject: [CP2K:19024] Terrible lattice constant of FCC copper from cell optimization<br>
<br>
Hi all,<br>
<br>
I got a terrible lattice constant (2.71) that is in significant disagreement with the experimental value (3.61). Something wrong with my input file ?<br>
<br>
*******************************************************************************<br>
*** GEOMETRY OPTIMIZATION COMPLETED ***<br>
*******************************************************************************<br>
<br>
Reevaluating energy at the minimum<br>
<br>
CELL| Volume [angstrom^3]: 19.962031<br>
CELL| Vector a [angstrom]: 2.713 0.000 0.000 |a| = 2.712699<br>
CELL| Vector b [angstrom]: 0.000 2.713 0.000 |b| = 2.712699<br>
CELL| Vector c [angstrom]: 0.000 0.000 2.713 |c| = 2.712699<br>
CELL| Angle (b,c), alpha [degree]: 90.000000<br>
CELL| Angle (a,c), beta [degree]: 90.000000<br>
CELL| Angle (a,b), gamma [degree]: 90.000000<br>
CELL| Requested initial symmetry: CUBIC<br>
CELL| Numerically orthorhombic: YES<br>
CELL| Periodicity XYZ<br>
<br>
<br>
&GLOBAL<br>
PROJECT FCC_Copper<br>
RUN_TYPE CELL_OPT<br>
PRINT_LEVEL LOW<br>
&END GLOBAL<br>
<br>
&FORCE_EVAL<br>
METHOD SIRIUS<br>
STRESS_TENSOR ANALYTICAL<br>
&SUBSYS<br>
&CELL<br>
A 3.61 0.00 0.00<br>
B 0.00 3.61 0.00<br>
C 0.00 0.00 3.61<br>
SYMMETRY CUBIC<br>
PERIODIC XYZ<br>
&END CELL<br>
&COORD<br>
Cu 0.000 0.000 0.000<br>
Cu 0.000 1.805 1.805<br>
Cu 1.805 0.000 1.805<br>
Cu 1.805 1.805 0.000<br>
&END COORD<br>
&KIND Cu<br>
BASIS_SET DZVP-MOLOPT-SR-GTH<br>
POTENTIAL GTH-PBE-q11<br>
&END KIND<br>
&END SUBSYS<br>
&PW_DFT<br>
&CONTROL<br>
VERBOSITY 2<br>
&END CONTROL<br>
&PARAMETERS<br>
ELECTRONIC_STRUCTURE_METHOD pseudopotential<br>
GK_CUTOFF 5.0<br>
PW_CUTOFF 20.00<br>
NUM_DFT_ITER 399<br>
NGRIDK 10 10 10<br>
&END PARAMETERS<br>
&ITERATIVE_SOLVER<br>
ENERGY_TOLERANCE 1e-5<br>
NUM_STEPS 20<br>
SUBSPACE_SIZE 4<br>
&END ITERATIVE_SOLVER<br>
&MIXER<br>
BETA 0.3<br>
TYPE BROYDEN2<br>
&END MIXER<br>
&END PW_DFT<br>
&DFT<br>
&XC<br>
&XC_FUNCTIONAL<br>
&GGA_X_PBE<br>
&END GGA_X_PBE<br>
&GGA_C_PBE<br>
&END GGA_C_PBE<br>
&END XC_FUNCTIONAL<br>
&END XC<br>
&END DFT<br>
&END FORCE_EVAL<br>
<br>
&MOTION<br>
&CELL_OPT<br>
MAX_ITER 200<br>
OPTIMIZER BFGS<br>
KEEP_SYMMETRY .TRUE.<br>
&END CELL_OPT<br>
&END MOTION<br>
<br>
Best Regards<br>
<br>
Jibiao Li<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Buns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a><mailto:<a href="mailto:cp2k%252Buns...@googlegroups.com" target="_blank">cp2k%2Buns...@googlegroups.com</a>><mailto:<a href="mailto:cp2k%2Buns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a><mailto:<a href="mailto:cp2k%252Buns...@googlegroups.com" target="_blank">cp2k%2Buns...@googlegroups.com</a>>>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAAk9cEhyfbBYjAthWCTS7Gx_sR9%3DF74R36Q-GHtbwTLzHD40Zg%40mail.gmail.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAAk9cEhyfbBYjAthWCTS7Gx_sR9%3DF74R36Q-GHtbwTLzHD40Zg%40mail.gmail.com</a><<a href="https://groups.google.com/d/msgid/cp2k/CAAk9cEhyfbBYjAthWCTS7Gx_sR9%3DF74R36Q-GHtbwTLzHD40Zg%40mail.gmail.com?utm_medium=email&utm_source=footer" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAAk9cEhyfbBYjAthWCTS7Gx_sR9%3DF74R36Q-GHtbwTLzHD40Zg%40mail.gmail.com?utm_medium=email&utm_source=footer</a>>.<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a><mailto:<a href="mailto:cp2k%252Bunsubscribe@googlegroups.com" target="_blank">cp2k%2Bunsubscribe@googlegroups.com</a>><mailto:<a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a><mailto:<a href="mailto:cp2k%252Bunsubscribe@googlegroups.com" target="_blank">cp2k%2Bunsubscribe@googlegroups.com</a>>>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/dea9b6e7-9721-4663-ab47-53744272bb3fn%40googlegroups.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/dea9b6e7-9721-4663-ab47-53744272bb3fn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/dea9b6e7-9721-4663-ab47-53744272bb3fn%40googlegroups.com?utm_medium=email&utm_source=footer" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/dea9b6e7-9721-4663-ab47-53744272bb3fn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a><mailto:<a href="mailto:cp2k%252Bunsubscribe@googlegroups.com" target="_blank">cp2k%2Bunsubscribe@googlegroups.com</a>>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZR0P278MB07597A89D3E8D6E7EF4D01C19F25A%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/ZR0P278MB07597A89D3E8D6E7EF4D01C19F25A%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM</a>.<br>
<br>
--<br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a><mailto:<a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAAk9cEgvx3Z3aAYfuYFr%3D6Q7oy7-H31FeQyFS453yXAkwkQrCQ%40mail.gmail.com" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAAk9cEgvx3Z3aAYfuYFr%3D6Q7oy7-H31FeQyFS453yXAkwkQrCQ%40mail.gmail.com</a><<a href="https://groups.google.com/d/msgid/cp2k/CAAk9cEgvx3Z3aAYfuYFr%3D6Q7oy7-H31FeQyFS453yXAkwkQrCQ%40mail.gmail.com?utm_medium=email&utm_source=footer" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAAk9cEgvx3Z3aAYfuYFr%3D6Q7oy7-H31FeQyFS453yXAkwkQrCQ%40mail.gmail.com?utm_medium=email&utm_source=footer</a>>.<br>
<br>
-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bunsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759B0541B5B6054C76C532C9F25A%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM" rel="noreferrer" target="_blank">https://groups.google.com/d/msgid/cp2k/ZR0P278MB0759B0541B5B6054C76C532C9F25A%40ZR0P278MB0759.CHEP278.PROD.OUTLOOK.COM</a>.<br>
</blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/CAAk9cEjCqR2k-ENQfsURfrHq5rTAFLyt42kiYYrzqapnTVDq3Q%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAAk9cEjCqR2k-ENQfsURfrHq5rTAFLyt42kiYYrzqapnTVDq3Q%40mail.gmail.com</a>.<br />