Dear all, <div><br></div><div>I am quite new to CP2K and working with transition metals. </div><div><br></div><div>I am trying to reproduce some results of an O2 molecule dissociating from an heme group (Fe2+). I am struggling to converge the SCF wavefunction even in a singlepoint calculation and I am not sure what I am doing wrong here. </div><div><br></div><div>As an initial guess, I am using PBE with the following combination of basis sets and pseudopotentials: </div><div>- Fe: BASIS_SET DZVP-MOLOPT-SR-GTH & GTH-PBE-q16. </div><div>- All the other atoms: DZVP-MOLOPT-GTH with GTH-PBE-q* . </div><div><br></div><div>I have attached the input file and the initial geometry.</div><div><br></div><div>Can you help me out? </div><div><br></div><div>Thank you very much in advance, </div><div><br></div><div>Salome</div>