Dear CP2K users,<div><br /></div><div>I'm encountering convergence problems when trying to optimize the bulk iron supercell consisting from 64 Fe atoms, i.e. 4x4x4 supercell. I wanted to simulate adsorption of different organic molecules on Fe slab later, so this is my starting point.</div><div><br /></div><div>SCF can't converge no matter what MIXING and DIAGONALIZATION method I apply. I always get something like this:</div><div><br /></div><div>...</div><div>    33 Broy./Diag. 0.70E+00    2.7     1.02510462     -7902.0344477276  2.42E+01<br />    34 Broy./Diag. 0.70E+00    2.7     1.02694146     -7882.2903511652  1.97E+01<br />    35 Broy./Diag. 0.70E+00    2.7     1.02572163     -7833.4164745080  4.89E+01<br />    36 Broy./Diag. 0.70E+00    2.7     1.02859994     -7878.1467068004 -4.47E+01<br />    37 Broy./Diag. 0.70E+00    2.7     1.02924599     -7812.2513589335  6.59E+01<br />    38 Broy./Diag. 0.70E+00    2.7     1.02822047     -7899.1170081476 -8.69E+01<br />    39 Broy./Diag. 0.70E+00    2.7     1.02715598     -7876.2960251819  2.28E+01<br />    40 Broy./Diag. 0.70E+00    2.7     1.02871900     -7939.8877789157 -6.36E+01<br />    41 Broy./Diag. 0.70E+00    2.7     1.02766599     -7827.2772337753  1.13E+02<br />    42 Broy./Diag. 0.70E+00    2.7     1.02780452     -7823.8921903136  3.39E+00<br />    43 Broy./Diag. 0.70E+00    2.7     1.02878646     -7862.7815187639 -3.89E+01<br />    44 Broy./Diag. 0.70E+00    2.7     1.02806441     -7910.5114238974 -4.77E+01<br /></div><div>...</div><div><br /></div><div>I followed the instructions from the previous discussions here in group and assembled the attached input file. Also, I looked for discussions regarding metallic systems including Fe such as:</div><div>https://groups.google.com/g/cp2k/c/ugkJLz4PSSI/m/C-ProBR6BwAJ</div><div><br /></div><div>I have an experience with the CP2K and molecular systems but not for solid state calculations.</div><div><br /></div><div>I also tried the following:</div><div>1) To use the smaller cell and MONKHORST-PACK Scheme (somewhere I saw that for cp2k it is recommended to perform calculations on the supercells)</div><div>2) To introduce a small random displacement in the system</div><div>3) To change ALPHA, BETA and other parameters within the MIXING section.</div><div>4) To try different MULTIPLICITY and initial MAGNETIZATION for Fe atoms.</div><div>but problem persists.</div><div><br /></div><div>I'm not sure what RELAX_MULTIPLICITY constant to use and if MAGNETIZATION and MULTIPLICITY are correct in this case. Need clarification on this as well.</div><div><br /></div><div>I attached my input and coord file, not sure where I'm making mistake.</div><div><br /></div><div>Thanks in advance!</div><div><br /></div><div>All best,</div><div>Branislav</div><div><br /></div><div> </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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d9d14e45-f243-4c4a-aa1f-988018aa889fn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d9d14e45-f243-4c4a-aa1f-988018aa889fn%40googlegroups.com</a>.<br />