<div dir="ltr">Hi Matthias<div> </div><img src="cid:70579cd7-420f-4751-8a5f-af302dd89bb9" alt="TIM截图20190320230817.png" width="320" height="97" style="margin-left: 1em; margin-right: 1em;"> <div>I am having some trouble when using the BS section: non-compatible end of section BETA , File: 'Ni_surface.inp', Line: 128, Column: 1, Chunk: <BEAT> . I don't know how to solve it.</div><div> </div><div>In example of UO2 (<a href="https://pc2.uni-paderborn.de/teaching/trainings/cp2k-tutorial/">https://pc2.uni-paderborn.de/teaching/trainings/cp2k-tutorial/</a>), there are two U of spin up and two U of spin down in UO2 cell, but we only use KIND Ua and KIND Ub to distinguish spin-up and spin-down U atoms. My question is how to let CP2K know which two U belong to KIND Ua， which two U belong to KIND Ub?</div><div> </div><div>Thinks!</div><div> </div><div> </div><div> 在 2018年11月8日星期四 UTC+8下午11:37:18，Matthias Krack写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"> <div lang="EN-US" link="blue" vlink="purple"> <div> Dear Katarina you have indeed to approach each spin channel in the BS section from a “closed shell” point of view. The slides <a href="https://pc2.uni-paderborn.de/fileadmin/pc2/cp2k_workshop/Implicit_Solvent_Methods_and_DFT_U.pdf#page=22" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpc2.uni-paderborn.de%2Ffileadmin%2Fpc2%2Fcp2k_workshop%2FImplicit_Solvent_Methods_and_DFT_U.pdf%23page%3D22\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH4gU3JP2wHtUZrL--31alEmbPwnw';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpc2.uni-paderborn.de%2Ffileadmin%2Fpc2%2Fcp2k_workshop%2FImplicit_Solvent_Methods_and_DFT_U.pdf%23page%3D22\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNH4gU3JP2wHtUZrL--31alEmbPwnw';return true;"> 22 and 23</a> from the last <a href="https://pc2.uni-paderborn.de/teaching/trainings/cp2k-tutorial" target="_blank" rel="nofollow" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpc2.uni-paderborn.de%2Fteaching%2Ftrainings%2Fcp2k-tutorial\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFhi-TdnDSnjAM5mfZIIUwbpdmWyA';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fpc2.uni-paderborn.de%2Fteaching%2Ftrainings%2Fcp2k-tutorial\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNFhi-TdnDSnjAM5mfZIIUwbpdmWyA';return true;"> CP2K tutorial in Paderborn</a> show an example for a BS section for UO2. The BS section is not perfectly intuitive. HTH Matthias From: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="I3bAyOMVBQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="I3bAyOMVBQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>> On Behalf Of katarína stanciaková Sent: Donnerstag, 8. November 2018 14:36 To: cp2k <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="I3bAyOMVBQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>> Subject: Re: [CP2K:10917] Re: Initial guess of the density matrix <div> <div> Dear Matthias, </div> <div> thank you for the suggestion. Could you, please, explain why NEL in alpha channel is 1? </div> <div> Is it becasue neutral Cu has 4s1, so we add one electron and then divide by 2? What if the result would be even number? </div> <div> </div> <div> Thanks! </div> Dňa štvrtok, 8. novembra 2018 11:49:26 UTC+1 Matthias Krack napísal(-a): <blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> <div> Dear Katarina You may try for Cu1 alpha NEL 1 0 L 0 2 N 4 3 instead of NEL 0 0 L 0 2 N 4 3 Hope that works. Matthias From: <a>cp...@googlegroups.com</a> <<a>cp...@googlegroups.com</a>> On Behalf Of katarína stanciaková Sent: Donnerstag, 8. November 2018 11:13 To: cp2k <<a>cp...@googlegroups.com</a>> Subject: Re: [CP2K:10913] Re: Initial guess of the density matrix <div> Dear all, I am trying to set up calculation for my Cu+ system, however, even after reading this thread I cannot understand how the electrons are assigned to channels in CP2K. I wanna test two possible spin states (just for testing purposes) of Cu+: Cu+ with electron configuration 3d9 4s1 (Cu1) and 3d10 4s0 (Cu2). I tried to build my input sections: &KIND Cu1 ELEMENT Cu BASIS_SET DZVP-MOLOPT-GTH POTENTIAL GTH-BLYP-q11 &BS &ALPHA &ALPHA NEL 0 0 L 0 2 N 4 3 &END &BETA NEL -1 -2 L 0 2 N 4 3 &END &END and &KIND Cu2 ELEMENT Cu BASIS_SET DZVP-MOLOPT-GTH POTENTIAL GTH-BLYP-q11 &BS &ALPHA &ALPHA NEL -1 0 L 0 2 N 4 3 &END &BETA NEL -1 0 L 0 2 N 4 3 &END &END Is this correct? <div> </div> <div> Also,regarding previous example on Cl1 - what is the meaning of half-electrons in alpha/beta channel? </div> <div> </div> Thank you very much for the help Katarina Dňa streda, 10. februára 2016 11:53:26 UTC+1 S Ling napísal(-a): <blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> You cannot have half an electron on a p orbital? <div> </div> <div> In your case of Cl1, you have one extra electron in the alpha spin channel, and by doing that, you break the symmetry of alpha and beta spin channels, i.e. you MAY end up with different occupations in the alpha and beta spin channels after the calculation. Of course, this also depends on your specific system, that a symmetry-broken solution is energetically more stable. </div> <div> </div> <div> SL </div> <div> </div> </div> <div> <div> On 9 February 2016 at 07:49, tao liu <<a>liu...@gmail.com</a>> wrote: <blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt"> <div> Dear All, I can understand the &BS setting for Cu1, but still confused on Cl and Cl1. &KIND Cl BASIS_SET DZVP-MOLOPT-GTH POTENTIAL GTH-BLYP-q7 &BS &ALPHA NEL 2 L 1 N 3 &END &BETA NEL 2 L 1 N 3 &END &END &END &KIND Cl1 ELEMENT Cl BASIS_SET DZVP-MOLOPT-GTH POTENTIAL GTH-BLYP-q7 &BS &ALPHA NEL 2 L 1 N 3 &END &BETA NEL 0 L 1 N 3 &END &END &END As suggested by Marcella, for 'Cl' in the input file, Alpha is (5+2)/2=3.5, and same for Beta, in the end 'Cl' should have 3s-2 3p-7 for 'Cl1' in the input file, Alpha is (5+2)/2=3.5, and (5+0)/2 =2.5 for Beta, in the end 'Cl1' should have 3s-2 3p-6 Could anyone tell me if I understand correctly or not ? Thanks, Tao </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a>cp2k+...@googlegroups.com</a>. To post to this group, send email to <a>cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/optout</a>. </blockquote> </div> </div> </blockquote> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a>cp2k+...@googlegroups.com</a>. To post to this group, send email to <a>cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;"> https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;"> https://groups.google.com/d/optout</a>. </div> </div> </blockquote> </div> -- You received this message because you are subscribed to the Google Groups "cp2k" group. To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="I3bAyOMVBQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.com</a>. To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="I3bAyOMVBQAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>. Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/group/cp2k</a>. For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/optout</a>. </div> </div> </blockquote></div></div>