<div dir="ltr">Thanks a lot, Dr. Mattioli. If this is the case, then I would just go ahead with a fixed cell relax calculation.<div><br></div><div>Thanks</div><div>Mitul</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Thu, Feb 6, 2014 at 2:08 PM, 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">
<br>
Dear Mitul<br>
If I remember well you cannot perform NEB calculations in a variable-cell framework. You should choose a reasonable set of lattice parameters along<br>
the whole reaction path.<br>
HTH<br>
Giuseppe<br>
<div class=""><br>
On Thursday 06 February 2014 07:36:18 Mitul Mundra wrote:<br>
> Hello,<br>
><br>
> I need to do NEB calculations for my system. My system consists of Si 2*2*2<br>
> supercell with some Li impurities. Prior to NEB, I am doing a vc-relax<br>
> calculation to get the minimum energy states for the first and the last<br>
</div>> images of the NEB. For the vc-relax, I set *ibrav = 1* for the symmetric<br>
<div class="">> cubic system. In the results obtained from vc-relax calculation, the cubic<br>
> symmetry of the system is lost and I get the following call parameters for<br>
> the initial and final state in the output file.<br>
><br>
> CELL_PARAMETERS (alat= 20.51108588)<br>
>    1.011059885  -0.002882944  -0.000008465<br>
>   -0.002882944   1.011059885   0.000008465<br>
>   -0.000008492   0.000008492   1.007297719<br>
><br>
> CELL_PARAMETERS (alat= 20.51108588)<br>
>    1.004932042   0.000000442  -0.000000442<br>
>    0.000000442   1.004688370   0.000104290<br>
>   -0.000000442   0.000104290   1.004688370<br>
><br>
> Now, I am confused about my input file for the NEB. I can not just directly<br>
> change the lattice constant by taking the cube root of the new volume<br>
</div>> obtained. I tried to set* ibrav = 0* and give the cell parameters obtained<br>
<div class="">> but I can only enter them once (correct me if I am wrong) and thus, can not<br>
> incorporate both the set of cell parameters. I would be extremely thankful<br>
> if anyone could suggest me a way to tackle this difficulty.<br>
><br>
><br>
> Thanks<br>
> Mitul Mundra<br>
> Final year Dual Degree Student<br>
> Department of Chemical Engineering<br>
> IIT Kanpur<br>
<br>
<br>
</div>********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l'utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l'oppression.<br>
********************************************************<br>
<br>
   Giuseppe Mattioli<br>
   CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
   v. Salaria Km 29,300 - C.P. 10<br>
   I 00015 - Monterotondo Stazione (RM)<br>
   Tel <a href="tel:%2B%2039%2006%2090672836" value="+390690672836">+ 39 06 90672836</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>
   ResearcherID: F-6308-2012<br>
<br>
_______________________________________________<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" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br></div>