Hi everyone, <br>
<br>
I am calculating the el-ph interactions of graphene with Pwscf 3.1.1. There are two questions here:<br>
<br>
1) First thing is about the symmetry obtained with pw.x. For the same
graphene system, when I use ibrav=0 with cell parameters, I got only 8
symmetry opertions( with inversion). There are six modes in the phonon
spectra. However, since the graphene system is hexagonal, I rotate the
coordinate system and use ibrav=4 to define the system, then I got 24
sym operations (with inversion). There are also six modes (with 2
degenerate modes) obtained. Why does the pw.x obtained different
symmetry for the same system? Can I adjust it manually? Here is my
input:<br>
<br>
<a href="http://C.scf.fit.in">C.scf.fit.in</a>:<br>
<br>
 &control<br>
    calculation='scf'<br>
    restart_mode='from_scratch',<br>
    prefix='C',<br>
    pseudo_dir = '/dsgpfs/jayan/celph',<br>
    outdir='/dsgpfs/jayan/scratch',<br>
    wf_collect=.true.<br>
 /<br>
 &system<br>
    ibrav=  0, celldm(1) =4.628,nat=2,ntyp= 1,<br>
    ecutwfc =30.0,nr1=24,nr2=24,<br>
    occupations='smearing', smearing='methfessel-paxton', degauss=0.02,<br>
    la2F=.true.,<br>
 /<br>
 &electrons<br>
    conv_thr =  1.0d-10<br>
    mixing_beta = 0.7<br>
 /<br>
ATOMIC_SPECIES<br>
 C  12.0107  C_nc<br>
ATOMIC_POSITIONS (alat)<br>
C  0.0000000000    0.0000000000     0.0000000000<br>
C  0.5773502693    0.0000000000     0.0000000000<br>
K_POINTS {automatic}<br>
 20 20 1 0 0 0<br>
CELL_PARAMETERS<br>
  0.866025404   -0.500000000    0.000000000<br>
  0.000000000    1.000000000    0.000000000<br>
  0.000000000    0.000000000    6.000000000<br>
<br>
<br>
2) The second question is about the output of the ph.x. I notice that the C.elph.out output following lines (2 2 1 1):<br>
<br>
     There are   4 irreducible representations<br>
<br>
     Representation     1      2 modes - To be done<br>
<br>
     Representation     2      2 modes - To be done<br>
<br>
     Representation     3      1 modes - To be done<br>
<br>
     Representation     4      1 modes - To be done<br>
     PHONON       :    16.94s CPU time<br>
<br>
but at the end of the calculation, I got the omega as follows:<br>
<br>
     omega( 1) =      16.131001 [THz] =     538.075842 [cm-1]<br>
     omega( 2) =      16.131001 [THz] =     538.075842 [cm-1]<br>
     omega( 3) =      30.116408 [THz] =    1004.581911 [cm-1]<br>
     omega( 4) =      36.922309 [THz] =    1231.603836 [cm-1]<br>
     omega( 5) =      36.922309 [THz] =    1231.603836 [cm-1]<br>
     omega( 6) =      39.058597 [THz] =    1302.863199 [cm-1]<br>
<br>
Note that the order of the omega is different from the calculation
order (2 1 2 1). Then I am quite confused at this, why are the output
omega not the same order as the real calculation? Did the progam sort
them after completion? I am interested in some data for each
representation, so I need to get the omega for the corresponding
representation. <br>
<br>
Thanks for your help in advance!<br>
<br>
Jia-An Yan<br>