<div dir="ltr">Hi, you can try the following settings based on your input Â Â &SCF Â Â Â ADDED_MOS 3 Â Â Â &SMEAR Â Â Â Â #Â Â Â Â Â LIST 1-3 1/3 Â Â Â Â METHOD ENERGY_WINDOW Â Â Â Â WINDOW_SIZE 0.1 Â Â Â Â FIXED_MAGNETIC_MOMENT -1 Â Â Â &END SMEAR Â Â Â EPS_SCF 2.0E-5 Â Â &END SCF which should result in an almost degenerate 1/3 occupation for that specific case of an oxygen atom. The degeneracy will improve for a setup with a really isolated oxygen atom, e.g.by using a larger cubic cell (> (10A)^3). In your (cell) setup the oxygen atom is still interacting with its images and its environment is not isotropic. Matthias On Thursday, 4 August 2016 14:46:38 UTC+2, Michael Trumm wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Ahoi, i am trying to perform some atomistic calculations using fractional occupation numbers. (in this example it would be the neutral oxygen atom) after i initialized them via &SMEAR the occupation looks kind of ok, but during the run cp2k chooses one of the tripletts. hence my question. is there any way to fix the fractional occupation during the whole run to get a fractional triplett with the beta p-orbitals occupied p_x 1/3Â p_y 1/3 and p_z 1/3 ? thanks for any help, Michael. OUTPUT: Â ALPHA MO EIGENVALUES AND MO OCCUPATION NUMBERS AFTER SCF STEP 20 # MO indexÂ Â Â Â Â Â Â Â Â MO eigenvalue [a.u.]Â Â Â Â Â Â Â Â Â Â Â MO occupation Â Â Â Â Â Â Â Â 1Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.917069Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 2Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.394950Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 3Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.394539Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 4Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.312729Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 5Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.468611Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.000000 Â Â Â Â Â Â Â Â 6Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.493158Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.000000 # SumÂ Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 4.000000 Â Fermi energy:Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.241121 Â HOMO-LUMO gap:Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.781340 =Â 21.26 eV Â BETA MO EIGENVALUES AND MO OCCUPATION NUMBERS AFTER SCF STEP 20 # MO indexÂ Â Â Â Â Â Â Â Â MO eigenvalue [a.u.]Â Â Â Â Â Â Â Â Â Â Â MO occupation Â Â Â Â Â Â Â Â 1Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.778066Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 2Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.268824Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 1.000000 Â Â Â Â Â Â Â Â 3Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.213094Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.000000 Â Â Â Â Â Â Â Â 4Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â -0.211852Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 0.000000 # SumÂ Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â 2.000000 INPUT: &GLOBAL Â PROJECT test Â RUN_TYPE ENERGY &END GLOBAL &FORCE_EVAL Â METHOD QS Â &DFT Â Â LSD Â Â MULTIPLICITY 3 Â Â &PRINT Â Â Â Â &MO Â Â Â Â Â Â OCCUPATION_NUMBERS Â Â Â Â &END MO Â Â &END PRINT Â Â BASIS_SET_FILE_NAME BASIS_MOLOPT Â Â POTENTIAL_FILE_NAME GTH_POTENTIALS Â Â CHARGEÂ 0 Â Â &SCF Â Â Â ADDED_MOS 2 Â Â Â &SMEAR Â Â Â Â Â LIST 1-3 1/3 Â Â Â &END SMEAR Â Â Â EPS_SCF 2.0E-5 Â Â &END SCF Â Â &XC Â Â Â &XC_FUNCTIONAL PBE Â Â Â &END XC_FUNCTIONAL Â Â &END XC Â &END DFT Â &SUBSYS Â Â &CELL Â Â Â Â ABC 5.189 8.953 20.129 Â Â Â Â ALPHA_BETA_GAMMA 90 101.11 90 Â Â Â Â PERIODIC XYZ Â Â &END CELL Â Â &KIND O Â Â Â Â BASIS_SET DZVP-MOLOPT-GTH Â Â Â Â POTENTIAL GTH-PBE-q6 Â Â &END KIND Â Â &COORD Â OÂ Â Â Â Â Â Â -0.9968591421Â Â Â Â Â Â Â 2.9869391388Â Â Â Â Â Â Â 4.8752526707 Â Â &END COORD Â &END SUBSYS &END FORCE_EVAL </div></blockquote></div>