<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 />