Dear devellopers, dear all, I would like to perform a vc-relax that concerves the point group of each atom and the symmetry of the cell. I know it is possible to impose symmetries on the BZ's k-points or to fix some atoms coordinates, but this does not give the expected result. For example, in a trigonal cell, ibrav 5 imposes the symmetries: v1 = a(tx,-ty,tz) v2 = a(0,2ty,tz) v3 = a(-tx,-ty,tz) but these symmetries seem to be broken by the relaxation, so that at the end of the relaxation, v2(1) is not exactelly null and v1(1) is not exactelly -v3(1), or for the atomic positions, a coordinate (x, x, x) alat becomes (x, 0.99 x, 1.001x ).
?? Is it possible to impose that these constraints on the unit-cell and on the atomic positions stay fix during the relaxation process ?? 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. 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. 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... Best regards, Antoine Jay Save the date le 14 octobre 2017 : Venez d��couvrir l'ISAE-SUPAERO de 10h �� 18h pour les portes ouvertes du leader mondial de la formation en ing��nierie a��rospatiale. Rencontre avec Thomas Pesquet, dipl��m�� Supa��ro 2001, au programme. Plus d'infos sur isae-supaero.fr
_______________________________________________ Pw_forum mailing list [email protected] http://pwscf.org/mailman/listinfo/pw_forum
