<div dir="ltr">On Wed, Sep 26, 2018 at 2:08 PM, Carina Backtorp <span dir="ltr"><<a href="mailto:carinabacktorp@hotmail.com" target="_blank">carinabacktorp@hotmail.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"> <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">As an example, if one places ONE  water molecule in a large cubic box and runs a vc-relax calculation in QE with ibrav = 1, i.e. cubic symmetry, then QE does a symmetry analysis concluding that the crystal symmetry is not ibrav 1 and therefore does the calculation with the correct symmetry.  
</div></blockquote><div><br></div><div>exactly: at the end of the calculation, your cell may no longer be cubic </div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
 If so, is there a key word in QE that I can include into the input file so that QE gives the obtained corrected symmetry in the output file?</div></div></blockquote><div><br></div><div>no <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
Also, is it possible to give QE a constraint in the input file so that it keeps the symmetry that was given? For example: even if the crystal structure does not have the cubic symmetry according to QE, can one still keep all three box lengths equal during a
 vc-relax calculation? </div></div></blockquote><div><br></div><div>option cell_dofree='volume' should keep the three axis as they are, varying only the volume<br></div><div><br></div><div>Paolo<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
Thank you for all help!</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
Kind regards,</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
Carina</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif;font-size:12pt">
</div>
<hr style="display:inline-block;width:98%">
<div id="m_926662176841979861divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>><br>
<b>Sent:</b> Wednesday, September 19, 2018 3:44 PM<br>
<b>To:</b> Quantum Espresso users Forum<br>
<b>Cc:</b> Carina Backtorp<br>
<b>Subject:</b> Re: [QE-users] Keep the symmetry during a calculation</font>
<div> </div>
</div>

<div>
<div dir="ltr">On Sat, Sep 15, 2018 at 10:09 AM, Carina Backtorp <span dir="ltr">
<<a href="mailto:carinabacktorp@hotmail.com" target="_blank">carinabacktorp@hotmail.com</a>></span> wrote:<br>
<div class="m_926662176841979861x_gmail_extra">
<div class="m_926662176841979861x_gmail_quote">
<blockquote class="m_926662176841979861x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_926662176841979861x_m_1727868156708052799divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<br>
<div style="margin-top:0;margin-bottom:0">1) When doing an vc-relax calculation in espresso, I expected that once the symmetry was given (ibrav=1) espresso should keep the symmetry during all the calculation</div>
</div>
</div>
</blockquote>
<div><br>
</div>
<div style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
it does. Note however that the symmetry of the lattice is just part of the story.  QE finds the starting symmetry of the crystal, and this is what is kept during the run (although occasionally, numerical noise and poor convergence may lead to the loss of the
 original symmetry). If the starting symmetry is not cubic, the final cell may lose its original cubic aspect.  Also note that QE uses its own criteria for symmetry, that may differ from criteria used by other codes, so something that QE deems cubic may be
 classified as orthorhombic by a pickier code</div>
<div style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<br>
</div>
<div style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
Paolo<br>
</div>
<div> </div>
<blockquote class="m_926662176841979861x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr"></div>
<br>
______________________________<wbr>_________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.o<wbr>rg</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso<wbr>.org/mailman/listinfo/users</a><span class="HOEnZb"><font color="#888888"><br>
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
</font></span></div><span class="HOEnZb"><font color="#888888">
<br>
<br clear="all">
<br>
-- <br>
<div class="m_926662176841979861x_gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. <a href="https://maps.google.com/?q=Udine,+via+delle+Scienze+208,+33100+Udine,+Italy&entry=gmail&source=g">Udine, via delle Scienze 208, 33100 Udine, Italy</a><br>
Phone +39-0432-558216, fax +39-0432-558222<br>
<br>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
</div>
</div>
</div>

</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>