<p dir="ltr">OK, I guess I did do MSR1a+charge.</p>
<p dir="ltr">I don't think there is an issue, just switch to MSR1. For insulators MSR1a often converges the charge & forces very well. For metals a -cc of 0.001 is more normal, and you probably have a metal (degenerate semiconductor) due to the charge.</p>
<p dir="ltr">---<br>
Professor Laurence Marks<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</a><br>
Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent">www.cfw.org/100-percent</a><br>
Co-Editor, Acta Cryst A<br>
    </p>
<div class="gmail_quote">On Feb 22, 2016 22:21, "Hu, Wenhao" <<a href="mailto:wenhao-hu@uiowa.edu">wenhao-hu@uiowa.edu</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi, Professor Marks:<br>
<br>
The following is the results from the label of :NEC:<br>
<br>
:NEC03: NUCLEAR AND ELECTRONIC CHARGE    814.00000   813.99598 Valence<br>
:NEC01: NUCLEAR AND ELECTRONIC CHARGE    814.00000   813.96208<br>
:NEC02: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00000<br>
:NEC03: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00092 Valence<br>
:NEC01: NUCLEAR AND ELECTRONIC CHARGE    814.00000   813.96207<br>
:NEC02: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00000<br>
:NEC03: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00091 Valence<br>
:NEC01: NUCLEAR AND ELECTRONIC CHARGE    814.00000   813.96207<br>
:NEC02: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00000<br>
:NEC03: NUCLEAR AND ELECTRONIC CHARGE    814.00000   814.00063 Valence<br>
<br>
About the charged state of defects in diamond, transitions among different charged states are possible to show up depending on the location of the Fermi energy. You can think about the nitrogen vacancy center, which has two stable states, i.e. neutral and -1 charged. The latter is proved to exhibit a spin one ground state.  Even though the supercell calculations have some flaws, my calculations so far qualitatively show a similar result in some charged isolated vacancy depending on various factors. The charged electron should be delocalized, but I can’t confirm the extent without a converged results. The thing I can confirm is that the magnetic moment are all from the defect level. Or maybe I didn’t get your question very well. If any other information about my calculation is needed, please let me know.<br>
<br>
Thank you very much for your help. I really appreciate it.<br>
<br>
Wenhao<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>