<div dir="ltr"><div>Thanks Giuseppe.<br><br></div>Nandan.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 7, 2018 at 5:49 AM, Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> Or just shifting the atoms slightly before starting the relaxation?<br>
<br>
</span>Yes! If the code finds a starting symmetry, it does not break it.<br>
<br>
HTH<br>
Giuseppe<br>
<div><div class="h5"><br>
<br>
Quoting Nandan Tandon <<a href="mailto:nandan.tandon@gmail.com">nandan.tandon@gmail.com</a>>:<br>
<br>
> Hi,<br>
><br>
> I am trying to relax a vacancy in a cubic 512 atom unit cell of diamond. It<br>
> has<br>
> been seen that the T_d symmetry around the vacancy is lowered to D_{2d}.<br>
> But my system does not show this because all the atoms (first 4 neighbors)<br>
> around the vacancy are at identical distances from the vacancy and each<br>
> other.<br>
><br>
> Initially, using a FCC supercell also gave the same relaxed structure. This<br>
> prompted me to start with a cubic structure.<br>
><br>
><br>
> My question is, is there a trick to do such vacancy calculations?<br>
><br>
> Like starting by substituting a carbon atom with another atom to distort<br>
> the lattice<br>
> before creating a vacancy at the site?<br>
><br>
> Or just shifting the atoms slightly before starting the relaxation?<br>
><br>
> I have done few tests to ensure that the Ecut is high (80Ryd), its a Gamma<br>
> point calculation, LDA potentials and an optimized lattice constant.<br>
><br>
> Thanks and regards,<br>
><br>
> Nandan.<br>
> --<br>
> ******************************<wbr>******************************<wbr>******************************<wbr>*****<br>
> Nandan Tandon<br>
> Michigan State University<br>
> East Lansing<br>
> Michigan 48824-1226<br>
> ******************************<wbr>******************************<wbr>******************************<wbr>*******<br>
<br>
<br>
<br>
</div></div>GIUSEPPE MATTIOLI<br>
CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
Via Salaria Km 29,300 - C.P. 10<br>
I-00015 - Monterotondo Scalo (RM)<br>
Mob (*preferred*) <a href="tel:%2B39%20373%207305625" value="+393737305625">+39 373 7305625</a><br>
Tel <a href="tel:%2B%2039%2006%2090672342" value="+390690672342">+ 39 06 90672342</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316">+39 06 90672316</a><br>
E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>><br>
<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <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>