<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">As part of a larger problem, I need to simulate a Nb-doped SrTiO3 surface, where 1/16 Ti sites is Nb. To first order this is 1/16 extra electrons. There are two ways I can think of to do this:</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">a) Charged cell calculation, with extra electrons + background charge.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">b) Virtual atom, where the Ti's have a Z of 22.0625.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Does anyone have any experience with something similar? I can think of reasons why both have issues, e.g. the extra potential in the vacuum in a).</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">N.B., doing a super-duper surface cell with 1/16 Ti changed is too vast.</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> <br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Györgyi</div></div></div>