[Pw_forum] PWscf can not find out symmetry operations!!!

Yingli Niu niuyingli at gmail.com
Thu Aug 31 19:48:40 CEST 2006


Respected Prof. Baroni
    Sorry for that I didn't browsed in the mailing list's archive carefully!
    I think you didn't understand me due to my english language problems.
1. I did a relax calculation for a slab with 10 layers, 20 C atoms. Than I
got 20 relaxed positions.
2. I did a scf calculation using the 20 positions. And pw.x found "No
symmetry".
3. I inputed the structure into Material Studio, and it find 8 symmetry
operations and 7 inequivalent atoms. The tolerance is 0.01A.
4. I used 8 symmetry operations to act on the 7 positions, then got 56
positions(some of them are superpositions).
5  I selected 20 positions in the primary unit cell(there were only 20
positions in one primary unit cell) and input them into the PWscf's run_scf
script. Then pw.x find 8 symmetry operations.
    This method is equivalent to modify the atom positions manually in order
to make the pw.x to find symmetry operations!
    The manual say that :

>the number of significant figures in the atomic positions is not large
enough. In file
>PW/eqvect.f90, the variable accep is used to decide whether a rotation is a
symmetry operation.
>Its current value (10-5 ) is quite strict: a rotated atom must coincide
with another atom to
>5 significant digits. You may change the value of accep and recompile.

    In the last email I said that I had changed the value of accep and
recompiled. But it made no use to find any symmetry operations. So I have no
method other than to modify the atom positions by symmetry manually.
  Now I think the problem has been solved!
  Many thanks again and forgive me to bring so many troubles!
Best wishes!
Yingli Niu


On Aug 31, 2006, at 04:01:40 +0200, Prof. Baroni wrote:

>> Dear Pro. Baroni

>>* Thanks a lot for your help!
*>>*   I have said that It was no use to set accep=3D1.0D-2
*>>* and recompiled pw.x. The program still could not find
*>>* any symmetry operation.
*>>*   Then I selected 7 inequivalent atoms, then used 8
*>>* symmetry operations:
*>>*   x,y,z
*>>*   -x,-y,z
*>>*   -x+1/2,y+1/2,-z
*>>*   x+1/2,-y+1/2,-z
*>>*   -x+1/2,-y+1/2,-z
*>>*   x+1/2,y+1/2,-z
*>>*   x,-y,z
*>>*   -x,y,z
*>>* to transform the 7 inequivalent atoms to 20 atoms in
*>>* the primary unit cell.
*>
>I can hardly understand what you did. If by applying some symmetry =20
>operations to a bunch of positions ("atoms") you obtain a larger =20
>number of atoms, doesn't this mean that those operations are not =20
>"symmetrry operations"???
>>*   Then I ran scf. I use Norm-Conserving PP and chose
*>>* ecutwfc=3D70,ecutrho=3D600.
*>
>there is ABSOLUTELY no point in using ecutrho>ecutwfc with norm =20
>conserving PP's
>
>>*   There appeared message "warning:symmetry operation #
*>>*  3/4/5/6 not allowed  fractional translation:" in the
*>>* out file as discussed before. Only 4 symmetry
*>>* operations without inversion could be found!
*>>* I set ecutrho =3D 500 and calculated it again. Then
*>>* all 8 symmetry operations was found!
*>>*   The stress tensor did not change too much!
*>>*   stress tensor(kbar) from original atom positions:
*>>*     49.20      0.01      0.07
*>>*      0.01    -16.72      0.16
*>>*      0.07      0.16      2.61
*>>*   stress tensor(kbar) from changed atom positions:
*>>*     49.98      0.00      0.00
*>>*      0.00    -16.13      0.00
*>>*      0.00      0.00      7.30
*>>*   All the nondiagonal elements changed to 0.00.
*>>*   I think changed accpe from 1.0D-5 to 1.0D-2 will be
*>>* harmful for the calculations. So I reset it to 1.0D-5
*>>* and recompiled.
*>>*   I have no other method to make the program find out
*>>* all 8 symmetry operations than this method.
*
>Have you browsed in the mailing list's archive, as I suggested? =20
>Please, let me NOT rediscuss this point that I and others discussed =20
>the last time a few days ago
>
>>*
*>>*   Another question:
*>>*   I use Yahoo mail box. How can I post the reply
*>>* message next to the message which is replied in the
*>>* mailing list. I alway started a new discussion line
*>>* with reply message! Sorry for that!
*>>
>Please, ask Terry Semel. He certainly knows better than me.
>
>SB
>
>---
>Stefano Baroni - SISSA  &  DEMOCRITOS National Simulation Center - =20
>Trieste
>[+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: /pipermail/attachments/20060901/6b1d6746/attachment.htm 


More information about the Pw_forum mailing list