<p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">Dear CP2K users</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">I would like to calculate core electron binding energy using GW100 method in CP2K.</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"><a href="https://www.cp2k.org/howto:gw" style="color: rgb(149, 79, 114); text-decoration: underline;">https://www.cp2k.org/howto:gw</a></span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">I started with a simple test, using a sulfate ion, SO42-, in vacuum to calculate the S2p lines.</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">I first optimized the structure of SO42- in vacuum at PBE level, and then perform a GW100 calculation on the optimized structure to obtain the core levels. See attached the input files I used</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">All the calculations run well, with no problem.</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">What bothers me is the S2p lines. My understanding reading the literature was that GW100 should include spin-orbit coupling. Thus S2p lines should be split into two levels: S2p(1/2) and S2p(3/2) (fine structure). Looking at the output (GW100.out) I see no splitting</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">  Molecular orbital   E_SCF (eV)       Sigc (eV)   Sigx-vxc (eV)       E_GW (eV)</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     1 ( occ )         -2389.519           4.316        -105.931       -2491.133</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     2 ( occ )          -498.313           4.023         -48.840        -543.129</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     3 ( occ )          -498.312           4.022         -48.840        -543.130</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     4 ( occ )          -498.306           4.024         -48.840        -543.122</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     5 ( occ )          -498.299           4.020         -48.840        -543.119</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     6 ( occ )          -201.514           7.495         -34.306        -228.325</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     7 ( occ )          -147.655           7.735         -25.358        <b>-165.279</b></span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     8 ( occ )          -147.655           7.735         -25.359        <b>-165.279</b></span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">     9 ( occ )          -147.654           7.735         -25.359        <b>-165.279</b></span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US" style="font-size: 10pt; font-family: Arial, sans-serif;">…..</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">Is there something wrong in my understating of GW100 and/or in my output?</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">Thank you very much for any possible help or suggestion…the method is new to me.</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US"> </span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">Best regards</span></p><p style="margin: 0cm; font-size: medium; font-family: Calibri, sans-serif; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);"><span lang="EN-US">Ivan</span></p>

<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/49340b35-e503-4e27-bcda-4b07cc145df7n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/49340b35-e503-4e27-bcda-4b07cc145df7n%40googlegroups.com</a>.<br />