<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Hi<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Your 2x2x2 model system contains 8 Fe atoms with a high-spin 3d4 configuration which results in 8*4=32 unpaired alpha (spin-up) 3d electrons and
thus the proper multiplicity should be rather 8*4+1=33 than 65 isn’t it?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">HTH<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Matthias<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> cp...@googlegroups.com <...@googlegroups.com>
<b>On Behalf Of </b>Jiapeng Liu<br>
<b>Sent:</b> Donnerstag, 28. November 2019 09:10<br>
<b>To:</b> cp2k <...@googlegroups.com><br>
<b>Subject:</b> [CP2K:12567] BS setting for spin polarized calculation<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Dear CP2K users,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I am performing a spin-polarized calculation of perovskite SrFeO3, which has been optimized using VASP in material project. According to the results shown here <a href="https://materialsproject.org/materials/mp-510624/">https://materialsproject.org/materials/mp-510624/</a>,
the material shows a magnetic moment of 4.17 \muB. I would like to reproduce the results using CP2K and the following are my input files for the calculation. For the BS section, I would like to simulate the Sr2+, Fe4+(3d4,4s0) high spin state, O2-. And the
multiplicity is set as 2*total spin+1=65 as the simulation supercell is 2x2x2.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">However, I noticed that in the output file, the total electrons in spin 1 is 208 and 144 in spin 2 after Re-scaling. <o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> Spin 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Re-scaling the density matrix to get the right number of electrons for spin 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> # Electrons Trace(P) Scaling factor<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 208 192.000 1.083<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"> Spin 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Re-scaling the density matrix to get the right number of electrons for spin 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> # Electrons Trace(P) Scaling factor<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 144 160.000 0.900<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I am wondering why it is not 192 in spin 1 and 160 in spin 2. Am I correct in setting the multiplicity state? I am quite confused about this part. And my other question is how to accelerate the scf convergence, the results shown in the
output seem to oscillate near the final energy for quite many steps. Can you help me with this, thank you for your help.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">&GLOBAL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> PROJECT SrFeO3_bulk<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> RUN_TYPE ENERGY<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> PRINT_LEVEL MEDIUM<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&END GLOBAL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">&FORCE_EVAL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> METHOD Quickstep<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &DFT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> UKS .TRUE.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> CHARGE 0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MULTIPLICITY 65<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">! RELAX_MULTIPLICITY 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> BASIS_SET_FILE_NAME BASIS_MOLOPT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> POTENTIAL_FILE_NAME GTH_POTENTIALS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &POISSON<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> PERIODIC XYZ<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END POISSON<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &QS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> EXTRAPOLATION USE_GUESS ! required for K-Point sampling<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> METHOD GPW<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END QS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &MGRID<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NGRIDS 4<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> CUTOFF 500<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> REL_CUTOFF 60<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END MGRID<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &SCF<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> SCF_GUESS ATOMIC<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> EPS_SCF 1.0E-6<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MAX_SCF 200<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ADDED_MOS 40<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &DIAGONALIZATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ALGORITHM STANDARD<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> EPS_ADAPT 0.01<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END DIAGONALIZATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &OUTER_SCF<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> EPS_SCF 1.0E-5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MAX_SCF 3<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END OUTER_SCF<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &SMEAR ON<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> METHOD FERMI_DIRAC<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ELECTRONIC_TEMPERATURE [K] 500<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END SMEAR<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &MIXING<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> METHOD BROYDEN_MIXING<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ALPHA 0.4<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> BETA 1.5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NBROYDEN 8<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END MIXING<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END SCF<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &XC<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &XC_FUNCTIONAL PBE<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END XC_FUNCTIONAL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END XC<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &KPOINTS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> SCHEME MONKHORST-PACK 4 4 4<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">! WAVEFUNCTIONS REAL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> SYMMETRY .TRUE.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> FULL_GRID .FALSE.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> PARALLEL_GROUP_SIZE -1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END KPOINTS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END DFT<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &SUBSYS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &CELL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ABC [angstrom] 3.91181515 3.91181515 3.91181515<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> PERIODIC XYZ<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MULTIPLE_UNIT_CELL 2 2 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END CELL<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &TOPOLOGY<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">! COORD_FILE_NAME sfo_mp-510624.cif<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">! COORD_FILE_FORMAT CIF<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> MULTIPLE_UNIT_CELL 2 2 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END TOPOLOGY<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &COORD<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> SCALED<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Sr 0.5 0.5 0.5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> Fe 0.0 0.0 0.0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> O 0.5 0.0 0.0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> O 0.0 0.5 0.0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> O 0.0 0.0 0.5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &KIND Sr<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ELEMENT Sr<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> BASIS_SET DZVP-MOLOPT-SR-GTH<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> POTENTIAL GTH-PBE-q10<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL -2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL -2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 0<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 5<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END KIND<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &KIND Fe<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ELEMENT Fe<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> BASIS_SET DZVP-MOLOPT-SR-GTH<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> POTENTIAL GTH-PBE-q16<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">! MAGNETIZATION 4<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL -2 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 0 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 4 3<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL -2 -6<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 0 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 4 3<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END KIND<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &KIND O<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ELEMENT O<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> BASIS_SET DZVP-MOLOPT-SR-GTH<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> POTENTIAL GTH-PBE-q6<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END ALPHA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> NEL 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> L 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> N 2<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BETA<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END BS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END KIND<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> &END SUBSYS<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">&END FORCE_EVAL<o:p></o:p></p>
</div>
</div>
</div>
<p class="MsoNormal">-- <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+unsu...@googlegroups.com">cp2k+unsu...@googlegroups.com</a>.<br>
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/3ae28d6c-88ba-4b12-a3e1-b8ba905e5cea%40googlegroups.com?utm_medium=email&utm_source=footer">
https://groups.google.com/d/msgid/cp2k/3ae28d6c-88ba-4b12-a3e1-b8ba905e5cea%40googlegroups.com</a>.<o:p></o:p></p>
</div>
</body>
</html>