<div dir="ltr">Hi,<div><br></div><div>I am a bit confused by the convention in the &BS section in CP2K. As in the example of magnetic system (slide 35) in https://www.cp2k.org/_media/events:2015_cecam_tutorial:ling_hybrids.pdf , in order to go from Fe to Fe3+ (high spin with alpha spin more electrons), the following is used for the input:</div><div><br></div><div>&KIND Fe1 </div><div>ELEMENT Fe </div><div>BASIS_SET DZVP-MOLOPT-SR-GTH </div><div>POTENTIAL GTH-PBE-q16 </div><div>&DFT_PLUS_U </div><div>L 2 </div><div>U_MINUS_J [eV] 5.0 </div><div>&END DFT_PLUS_U </div><div><b>&BS </b></div><div><b>&ALPHA </b></div><div><b>NEL +4 -2 </b></div><div><b>L 2 0 </b></div><div><b>N 3 4 </b></div><div><b>&END ALPHA </b></div><div><b>&BETA </b></div><div><b>NEL -6 -2 </b></div><div><b>L 2 0 </b></div><div><b>N 3 4 </b></div><div><b>&END BETA </b></div><div><b>&END BS </b></div><div>&END KIND<br></div><div><br></div><div>In the neutral atom Fe, shouldn't there only be only 2 4s electrons, one alpha and one beta, and 6 3d electrons, 3 alpha and 3 beta? It seems the lines NEL just have numbers double as large as what they should be, e.g. NEL +4 -2 should be NEL +2 -1. If the numbers for NEL is added for both spin channels, one gets -6 instead of -3 (since Fe loses 3 electrons to become Fe3+). Could you please tell me where I get it wrong? Thanks.</div><div><br></div><div>BR.</div><div>Chang Liu</div></div>