<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=us-ascii">

<meta name="Generator" content="Microsoft Word 14 (filtered medium)">

<style><!--

/* Font Definitions */

@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:0in;

        margin-bottom:.0001pt;

        font-size:11.0pt;

        font-family:"Calibri","sans-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;}

span.EmailStyle17

        {mso-style-type:personal-compose;

        font-family:"Calibri","sans-serif";

        color:windowtext;}

.MsoChpDefault

        {mso-style-type:export-only;

        font-family:"Calibri","sans-serif";}

@page WordSection1

        {size:8.5in 11.0in;

        margin:1.0in 1.0in 1.0in 1.0in;}

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">For 64 atoms cell you can test nearest neighbor(NN) , 2<sup>nd</sup> NN, and so on.  This can give you some local information about local arrangements of the defects.<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">If you reduce the cell to 32 and insert one dopant ion, then you are simulating an ordered array of defects because of the periodic boundary conditions.

<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">In fact the interaction across the cell boundary applies also to the case of 64 atoms above, but one hopes (and can possibly show) that the results represent the local interactions between the two defects within the cell rather than the

 image-image interaction across the cell boundary.<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">In my opinion, a clever way to study concentrated solutions of defects is to resort to model Hamiltonians along the lines of the cluster expansion approach.  This can work nicely if you care about neutral defects in semiconductors (eg:

 the oxygen vacancy in CeO2 as in <a href="http://arxiv.org/abs/1206.5429">http://arxiv.org/abs/1206.5429</a> ). For charged defects where long range interactions become important, there is no clean way (as far as I know) to get a model Hamiltonian for the

 interacting defects system. In the literature (especially that of CeO2) you may find few attempts for concentrated charged defects. See for example:

<a href="http://journals.aps.org/prb/abstract/10.1103/PhysRevB.90.115120">http://journals.aps.org/prb/abstract/10.1103/PhysRevB.90.115120</a><o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">Mostafa Youssef<o:p></o:p></p>

<p class="MsoNormal">MIT<o:p></o:p></p>

</div>

</body>

</html>