Hi qing zhou, <div><br /></div><div>I checked your input and the output, and the spectrum. </div><div>In spectrum, there is only one state at 6167.202695 eV, which I would say it bad but not so bad. </div><div>but the only one state should be noticed, which it is due to you set ENERGY_RANGE 30. The expected Fe 1s around 7000, is out of this range tooo far. </div><div>and I would suggest: </div><div>1. increase the Grid size.  <b>GRID Fex 500 1000 </b></div><div>2. decrease the number for searching, <b>N_SEARCH 5  </b></div><div>Because you only have one Fe here, and calculate its 1s excitations. The MOs order should be Fe_1s, at the lowest. </div><div>3. increase <b>RI_REGION</b> (you can test, ending up to contain all atoms in the structure) </div><div>4. maybe increase the <b>ENERGY_RANGE</b>, I would suggest 1000 or 2000. since now your first excitation state is 6167, the range should include the expected states (for Fe 1s 7000 eV approximately). After you improve the first excitation energy, it can be reduced. </div><div><br /></div><div>another option, decrease HF portion. PBE0(25%) or B3LYP(25%) might give a better spectrum, which I preferred for Fe 2p spectrum. </div><div><br /></div><div>Best, </div><div>Nanchen </div><div><br /></div><div><br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">在2025年11月13日星期四 UTC+1 06:19:15<qing zhou> 写道：<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear CP2K users and developers,<div><p>I am using CP2K to calculate the XAS of an open-shell system with a spin multiplicity of 2. The OT method was enabled for the SCF calculation to aid convergence, and I have set EXCITATIONS OS_SPIN_CONS in the &XAS_TDP section. For the functional, I am using a PBE-based functional with 45% Hartree-Fock exchange. The def2-TZVP basis set was used for the Fe atom, while the DZVP-MOLOPT basis set was used for all other atoms.</p><p>As shown in the output file, the calculation finished without any errors, and the SCF converged successfully. However, the excitation energies in the resulting spectrum file are obviously incorrect. I recall that the energy for the Fe K-edge should be around 7000 eV. Furthermore, the spectrum file contains only a single data point.</p><p>I would like to ask what potential issues might exist in my calculation setup. How can I adjust my settings to correctly obtain the XAS spectrum for this open-shell system?</p><p>I have attached my input (.inp), output (.out), and spectrum files for your reference. Thank you for your help!</p></div></blockquote></div>

<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 visit <a href="https://groups.google.com/d/msgid/cp2k/3cf8c482-18fa-4ac7-8f98-3a108b7fa9f6n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/3cf8c482-18fa-4ac7-8f98-3a108b7fa9f6n%40googlegroups.com</a>.<br />