<div dir="ltr"><div>Funny. Look at the electronic structure: with smearing, you have two degenerate states with up spin at the Fermi energy, atoms n.71, 125, 179 (those close to the vacancy) with the same spin polarization, 0.15 or so (whatever it means: atomic charges and polarizations reprinted by the code should not be taken too seriously). By fixing occupancies, you break the degeneracy and get a sort of "broken-symmetry" solution, with a gap, two of the three atoms with a larger polarization (0.29) and one with opposite polarization (-0.18). I guess that in order to reproduce the former solution with fixed occupancies, you have to set occupancies (from input) as 1 1 1 ... 0.5 0.5 0 0 ... for spin up, 1 1 1 ... 0 0 ... for spin down.</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 13, 2019 at 4:03 AM Brendan Smith <<a href="mailto:bsmith24@buffalo.edu">bsmith24@buffalo.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi QE Experts,</div><div><br></div><div>Please know that before sending this question, I have done my best to extensively read the existing body of past forum messages and other QE related learning materials. <br></div><div><br></div><div>I am currently studying NV centers in bulk diamond, and I have a question regarding a single NV0 (charge neutral) center I have made in a 3x3x3 (216) atom supercell of bulk diamond. In this system, we have an odd number of electrons, leaving a single lone electron 
(spin-up in this case). Therefore, I expect that my total magnetization should be 1.0 
Bohr mag/cell, and I see this consistently in my QE outputs. I am comfortable in thinking that my system is ferromagnetic, based on my deductions from reading the literature. However, I have some discrepancy when it comes to the absolute magnetization, and it is here where my question lies.</div><div><br></div><div>When I fix the 
total magnetization to be = 1.0 
Bohr mag/cell, and use a fixed occupations scheme, the absolute magnetization for some reason converges to a value of 2.11 Bohr mag/cell. I find this to be strange, since it should be around the value of the total magnetization, which is 1 .0
Bohr mag/cell. What this tells me is that for some reason, perhaps I am converging to some anti-ferromagnetic solution. Then again, I don't know if this is even possible given that I am using fixed occupations.</div><div><br></div><div>Next, when I use occupations = "smearing", ("cold") and set only a guess starting magnetization (do not fix total magnetization), I find that I actually get a value for BOTH the total and absolute magnetization to be around 1.0 
Bohr mag/cell, which is more or less what I expect to see, given that I believe that my system is rightfully a simple ferromagnetic system.<br></div><div><br></div><div>So my question is, what the heck is going on in the case of fixed occupations? why am I converging to  solution in which the total magnetization is 1.0 
Bohr mag/cell, and the absolute magnetization is 2.11 
Bohr mag/cell? Is there any apparent misunderstanding I may be having in my analysis? I seek your expertise in this regard. <br></div><div><br></div><div>I have attached my PWscf input and output files. The output file that was used with smearing is labeled "x0.scf_smear.out", and the output file in which fixed occupations were used along with setting tot_magnetization is labeled "x0.scf_fixed.out". The input file is labeled "<a href="http://x0.scf.in" target="_blank">x0.scf.in</a>". To run the calculations with the fixed occupations, I just comment out the smearing part and uncomment the tot_magnetization line. <br></div><div><br></div><div>Thank you for your time in reading my email. I look forward to any suggestions.</div><div><br></div><div>Best,<br></div><div>Brendan A. Smith</div><div>Ph.D candidate at the State University of New York at Buffalo<br></div><div><br></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</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></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>