<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I am trying to relax a vacancy in a cubic 512 atom unit cell of diamond. It has<br></div>been seen that the T_d symmetry around the vacancy is lowered to D_{2d}. <br></div>But my system does not show this because all the atoms (first 4 neighbors) around the vacancy are at identical distances from the vacancy and each other.<br><br></div><div>Initially, using a FCC supercell also gave the same relaxed structure. This prompted me to start with a cubic structure.</div><br><div></div><br><div>My question is, is there a trick to do such vacancy calculations?<br clear="all"></div><div><div><div><div><div><div><br></div><div>Like starting by substituting a carbon atom with another atom to distort the lattice<br></div><div>before creating a vacancy at the site?<br><br></div><div>Or just shifting the atoms slightly before starting the relaxation?<br></div><div><br></div><div>I have done few tests to ensure that the Ecut is high (80Ryd), its a Gamma point calculation, LDA potentials and an optimized lattice constant. <br><br></div><div>Thanks and regards,<br><br></div><div>Nandan.<br></div><div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>***********************************************************************************************<br>Nandan Tandon<br>Michigan State University<br></div><div>East Lansing<br>Michigan 48824-1226</div><div>*************************************************************************************************</div></div></div>
</div></div></div></div></div></div></div>