<div dir="ltr"><div>You need two define two different (but equal) atomic species if you expect some antiferromagnetic ordering</div><div>(thus starting from opposite starting_magnetization on the two).</div><div><br></div><div>Usually, if you seek a magnetic solution, you can just set starting_magnetization to 0.5 or 0.7 and the </div><div>self consistent charge density automatically will produce a solution in which the spin density will be localized only</div><div>on the atoms carrying a magnetic moment. So, while you need to distinguish between atoms with opposite magnetizations,</div><div>you shouldn't need, usually, to distinguish between atoms with no net magnetization and atoms with net magnetization.</div><div><br></div><div>I would, at a first glance, try sto set starting_magnetization = 0.7 on all Si atoms and see what happens.</div><div>Giovanni</div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">Il giorno sab 19 apr 2025 alle ore 19:46 Lorenzo Sponza <<a href="mailto:lorenzo.sponza@onera.fr">lorenzo.sponza@onera.fr</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thanks Giovanni for your quick answer.<br>
<br>
Still, it's not clear to me how to set the starting magnetization. <br>
Putting just starting_magnetization=X will be useless because this will <br>
set a starting magnetization in all Si atoms, that will break no <br>
symmetry.<br>
On the other hand, if I understood well I can define two atomic species, <br>
say Si1 and Si2, and define the starting magnetization for only one of <br>
the two. I think that this is the thing I should do, but I can not <br>
figure out which atoms should belong to one group or the other and <br>
wheather I need more than two groups for this specific problem.<br>
I was hoping to get help from someone who had already done similar <br>
simulations.<br>
Cheers<br>
<br>
Le 2025-04-19 17:25, Giovanni Cantele a écrit :<br>
> Dear Lorenzo,<br>
> <br>
> the variable starting_magnetization is usually used to drive the<br>
> system<br>
> out of the non magnetic solution, in that starting_magnetization(i) =<br>
> 0 for all atomic species<br>
> should always bring to a spin unpolarized system. The starting value<br>
> in some cases is not so important<br>
> (different values will produce the same magnetic ground state), but as<br>
> far as I remember it could be<br>
> in some (not so frequent) cases where different magnetic energy minima<br>
> exist.<br>
> <br>
> For the other questions:<br>
> - no special attention paid to the pseudopotential choice, if you mean<br>
> to choose a pseudopotential "specific"<br>
> for a spin-polarized calculation<br>
> - I would not use the input_DFT variable unless you need it for<br>
> specific purposes (e.g. non local functionals, like vdw-DF).<br>
> <br>
> Ciao,<br>
> Giovanni<br>
> <br>
> --<br>
> <br>
> Giovanni Cantele, PhD<br>
> CNR-SPIN<br>
> c/o Dipartimento di Fisica<br>
> Universita' di Napoli "Federico II"<br>
> Complesso Universitario M. S. Angelo - Ed. 6<br>
> Via Cintia, I-80126, Napoli, Italy<br>
> e-mail: <a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a><br>
> Phone: +39 081 676910<br>
> Skype contact: giocan74<br>
> <br>
> ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" rel="noreferrer" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>
> Web page: <a href="https://sites.google.com/view/giovanni-cantele/home" rel="noreferrer" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br>
> <br>
> Il giorno sab 19 apr 2025 alle ore 16:37 Lorenzo Sponza<br>
> <<a href="mailto:lorenzo.sponza@onera.fr" target="_blank">lorenzo.sponza@onera.fr</a>> ha scritto:<br>
> <br>
>> Dear Quantum ESPRESSO users and developpers,<br>
>> <br>
>> I need some help to set up a calculation of a Si(111) 1x1 surface.<br>
>> This<br>
>> surface is characterized by dangling bonds that spoil the electronic<br>
>> <br>
>> states around Fermi. In literature, two approaches are usually<br>
>> adopted.<br>
>> One consists in saturating the surface with H atoms, the other<br>
>> consists<br>
>> in performing a spin polarized calculation.<br>
>> I want to give both methods a try, but I stuck when it comes to set<br>
>> up<br>
>> the spin polarized input file. I set nspin=2, but I'm perplexed<br>
>> about<br>
>> the value to give to starting_magnetization since I have only Si<br>
>> atoms<br>
>> here, and in this specific case it is not clear to me how this<br>
>> variable<br>
>> is interpreted and used by the code.<br>
>> Finally, is there any attention to be paid in choosing the<br>
>> pseudopotential or in specifying the input_DFT variable? I've used<br>
>> PBE<br>
>> scalar relativistic PAW pseudos so far.<br>
>> <br>
>> Many thanks for your help.<br>
>> <br>
>> --<br>
>> Dr. Lorenzo Sponza<br>
>> <br>
>> Chargé de Recherche au CNRS<br>
>> Université Paris-Saclay, ONERA, CNRS Laboratoire d'étude des<br>
>> microstructures (LEM)<br>
>> 29 Avenue de la division Leclerc, 92322 Châtillon<br>
>> Tel: +33146734464<br>
>> <br>
>> ETSF Research Team Leader<br>
>> European Theoretical Spectroscopy Facility<br>
>> 17 Sart-Tilman B-4000 Liège, Belgium<br>
>> <br>
> _______________________________________________________________________________<br>
>> The Quantum ESPRESSO Foundation stands in solidarity with all<br>
>> civilians worldwide who are victims of terrorism, military<br>
>> aggression, and indiscriminate warfare.<br>
>> <br>
> --------------------------------------------------------------------------------<br>
>> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a> [1])<br>
>> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
> <br>
> <br>
> Links:<br>
> ------<br>
> [1] <a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">http://www.max-centre.eu</a><br>
> _______________________________________________________________________________<br>
> The Quantum ESPRESSO Foundation stands in solidarity with all<br>
> civilians worldwide who are victims of terrorism, military aggression,<br>
> and indiscriminate warfare.<br>
> --------------------------------------------------------------------------------<br>
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br>
<br>
-- <br>
Dr. Lorenzo Sponza<br>
<br>
Chargé de Recherche au CNRS<br>
Université Paris-Saclay, ONERA, CNRS Laboratoire d'étude des <br>
microstructures (LEM)<br>
29 Avenue de la division Leclerc, 92322 Châtillon<br>
Tel: +33146734464<br>
<br>
ETSF Research Team Leader<br>
European Theoretical Spectroscopy Facility<br>
17 Sart-Tilman B-4000 Liège, Belgium<br>
</blockquote></div>