<div dir="ltr">On Mon, Oct 2, 2017 at 4:25 PM, JAY Antoine <span dir="ltr"><<a href="mailto:Antoine.JAY@isae-supaero.fr" target="_blank">Antoine.JAY@isae-supaero.fr</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">I would like to perform a vc-relax that conserves the point group of each atom and the symmetry of the cell.<br></span></blockquote><div><br></div><div>vc-relax conserves the starting symmetry of the crystal. There is no way to conserve a symmetry that is not there at the beginning of the structural optimization.<br>  <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">For example,  in a trigonal cell, ibrav 5 imposes the symmetries:<br>v1 = a(tx,-ty,tz)<br>v2 = a(0,2ty,tz)<br>v3 = a(-tx,-ty,tz)<br>but these symmetries seem to be broken by the relaxation,<br></span></blockquote><div><br></div><div>see above: symmetries that were present in the starting structure are not broken. If a "symmetry" is "broken by structural relaxation", it wasn't a symmetry.<br><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">so that at the end of the relaxation, v2(1) is not exactly null and v1(1) is not exactly -v3(1),<br></span>or for the atomic positions, a coordinate (x, x, x) alat becomes (x, 0.99 x, 1.001x ).<div class="HOEnZb"><div class="h5"><br><br>?? Is it possible to impose that these constraints on the unit-cell and on the atomic positions stay fix during the relaxation process ??<br><br>This is a major issue for the use of these relaxed values in a second step for phonons for example, because it will then not detect the symmetries.<br>I think it could be a good idea to add into the vc-relax output file a lign that gives the corresponding celldm(i) and not only the CELL_PARAMETERS.<br></div></div></blockquote><div><br></div><div>I also think that it is a good idea, but it takes a bit of time and effort, that becomes a byte, or even a word (32 or 64 bits) if one wants to write a piece of code that works perfectly in all cases<br><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 class="HOEnZb"><div class="h5">In fact, in the scf before the phonon use, we would like to have ibrav =/0 and celldm in state of ibrav 0 and quasi symmetric CELL_PARAMETERS...<br><br>Best regards,<br><br>Antoine Jay<br><br><br><br><br> 
<br><img title="" alt="" width="970" height="120">

<br>
</div></div><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><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>