<div dir="ltr">Maybe this information is useful: <a href="https://www.cp2k.org/howto:running_qe_computation">https://www.cp2k.org/howto:running_qe_computation</a><div><br></div><div>in short terms, the energy (E) of a plane wave with a wave vector G is: E = (h/2pi)²G²/(2m_e) which is reduced to (if a.u. is used): E = G²/2</div><div>so, E in Ha (a.u.) and G in a.u.</div><div>The G value is the input for SIRIUS, so you need: G [a.u.] = sqrt(2*E [Ha]) , and in QE the energy is used by default in Rydberg, so E_QE [Ry] = 2E [Ha]</div><div>Finally, G [a.u.] = sqrt(E_QE [Ry])</div><div><br></div><div>Regards - Lucas Lodeiro<br><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 4 jul 2023 a las 9:24, Jibiao Li (<<a href="mailto:science.sweden@gmail.com">science.sweden@gmail.com</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><font face="tahoma, sans-serif"><span lang="EN-US">Hi all,</span><span lang="EN-US"></span></font></p>

<p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><span lang="EN-US"><font face="tahoma, sans-serif"> </font></span></p>

<p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><font face="tahoma, sans-serif"><span lang="EN-US">According to the manual, units for
PW_CUTOFF and GK_CUTOFF are [a.u.^-1]. I am still puzzled by the question of
how I should convert a cutoff energy of 80 Ry to that can be accepted by
PW_CUTOFF in CP2K. </span><span lang="EN-US"></span></font></p>

<p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><span lang="EN-US"><font face="tahoma, sans-serif"> </font></span></p>

<p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><font face="tahoma, sans-serif"><span lang="EN-US">I discussed it with </span><span lang="EN-US" style="color:rgb(31,31,31)">Dr. </span><span lang="EN-US">Jürg Hutter, </span><span lang="EN-US" style="color:rgb(31,31,31)">but no definite answer was received.
Could anyone provide a clear answer?</span><span lang="EN-US"></span></font></p><p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><span lang="EN-US" style="color:rgb(31,31,31)"><font face="tahoma, sans-serif"><br></font></span></p><p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><span lang="EN-US" style="color:rgb(31,31,31)"><font face="tahoma, sans-serif">Regards</font></span></p><p class="MsoNormal" align="left" style="margin:0cm 0cm 0.0001pt"><font face="tahoma, sans-serif">Jibiao Li</font></p><table cellpadding="0" style="border-collapse:collapse;margin-top:0px;width:auto;font-size:14px;display:block"><tbody style="display:block"></tbody></table><div><font face="tahoma, sans-serif"><br></font></div><div><font face="tahoma, sans-serif"><span style="color:rgb(80,0,80)">&GLOBAL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  PROJECT FCC_Copper</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  RUN_TYPE CELL_OPT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  PRINT_LEVEL LOW</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">&END GLOBAL</span><br style="color:rgb(80,0,80)"><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">&FORCE_EVAL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  METHOD SIRIUS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  STRESS_TENSOR ANALYTICAL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &SUBSYS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &CELL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      A 3.61 0.00 0.00</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      B 0.00 3.61 0.00</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      C 0.00 0.00 3.61</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      SYMMETRY CUBIC</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      PERIODIC XYZ</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END CELL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &COORD</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      Cu   0.000 0.000 0.000</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      Cu   0.000 1.805 1.805</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      Cu   1.805 0.000 1.805</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      Cu   1.805 1.805 0.000</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END COORD</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &KIND Cu</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       POTENTIAL GTH-PBE-q11</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END KIND</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &END SUBSYS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &PW_DFT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &CONTROL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       VERBOSITY 2</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END CONTROL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &PARAMETERS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       ELECTRONIC_STRUCTURE_METHOD  pseudopotential</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       GK_CUTOFF ???</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       PW_CUTOFF ???</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       NUM_DFT_ITER 399</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       NGRIDK 10 10 10</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END PARAMETERS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &ITERATIVE_SOLVER</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       ENERGY_TOLERANCE 1e-5</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       NUM_STEPS 20</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       SUBSPACE_SIZE 4</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END ITERATIVE_SOLVER</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &MIXER</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       BETA 0.3</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">       TYPE BROYDEN2</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END MIXER</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &END PW_DFT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &DFT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &XC</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      &XC_FUNCTIONAL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">        &GGA_X_PBE</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">        &END GGA_X_PBE</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">        &GGA_C_PBE</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">        &END GGA_C_PBE</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">      &END XC_FUNCTIONAL</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    &END XC</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &END DFT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">&END FORCE_EVAL</span><br style="color:rgb(80,0,80)"><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">&MOTION</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &CELL_OPT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    MAX_ITER 200</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    OPTIMIZER BFGS</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">    KEEP_SYMMETRY .TRUE.</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">  &END CELL_OPT</span><br style="color:rgb(80,0,80)"><span style="color:rgb(80,0,80)">&END MOTION</span><br></font><div><br></div></div></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" 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/CAAk9cEgu9zMKLgO32U3dHjUJwiR9CQ2xAkr-BtB7F9x9CXjzaQ%40mail.gmail.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/CAAk9cEgu9zMKLgO32U3dHjUJwiR9CQ2xAkr-BtB7F9x9CXjzaQ%40mail.gmail.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/CAOFT4PLMRmsw%3DysRumFTj8htOfYHTPTyyCA79KSXdn3ruSMyNQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAOFT4PLMRmsw%3DysRumFTj8htOfYHTPTyyCA79KSXdn3ruSMyNQ%40mail.gmail.com</a>.<br />