<div><div><div>Dear Quantum Espresso experts, </div><div> </div><div>I am having trouble calculating the phonon bands of graphene.</div><div> </div><div>First, I took the cif-file of the already relaxed graphite cell from the Materials Project site https://next-gen.materialsproject.org/ and made graphene cell from it by hand.</div><div>Then I relaxed the cell obtained with vc-relax and 12x12x1 k-point set.</div><div>I used the same parameters to perform a scf calculation:</div><div> </div><div>&CONTROL</div><div>calculation='scf'</div><div>title='graphene'</div><div>prefix='graphene'</div><div>verbosity='high'</div><div>restart_mode='from_scratch'</div><div>nstep=1000</div><div>iprint=1</div><div>tstress=.true.</div><div>tprnfor=.true.</div><div>disk_io='default'</div><div>pseudo_dir = '/home/mxm2/Desktop/Graphene',</div><div>outdir='/home/mxm2/Desktop/Graphene/temp'</div><div>/</div><div>&SYSTEM</div><div>ibrav = 0,</div><div>celldm(1)=1.889726,</div><div>nat = 2,</div><div>ntyp = 1,</div><div>ecutwfc = 60.0 ,</div><div>ecutrho = 480.0 ,</div><div>/</div><div>&ELECTRONS</div><div>electron_maxstep = 1000,</div><div>conv_thr = 1.0d-12 ,</div><div>mixing_mode = 'plain' ,</div><div>mixing_beta = 0.7d0 ,</div><div>/</div><div>ATOMIC_SPECIES</div><div>C 12.011 C.pz-rrkjus.UPF</div><div>ATOMIC_POSITIONS (crystal)</div><div>C        0.000000000   0.000000000   0.000861580</div><div>C        0.333333000   0.666667000   0.000838420</div><div>CELL_PARAMETERS (alat)</div><div>   2.438515837  -0.000000000   0.000000000</div><div>  -1.219257919   2.111816662   0.000000000</div><div>  -0.000000000  -0.000000000  19.995916314</div><div> </div><div>K_POINTS automatic</div><div>12 12 1 0 0 0</div><div> </div><div>Then, I run ph.in calculation:</div><div> </div><div>&inputph</div><div>recover=.true.</div><div>tr2_ph=1.0d-14,  </div><div>prefix='graphene',</div><div>ldisp=.true. </div><div>nq1=12, </div><div>nq2=12, </div><div>nq3=1,  </div><div>amass(1)=12.011,   </div><div>outdir='/home/mxm2/Desktop/Graphene/temp', </div><div>fildyn='graphene.dyn',    </div><div>/</div><div> </div><div>My q2r:</div><div> </div><div>&input</div><div>fildyn='graphene.dyn',</div><div>zasr='simple',</div><div>flfrc='graphene.12x12x1.fc'</div><div>/</div><div> </div><div>Matdyn input:</div><div> </div><div>&input</div><div>asr='simple',</div><div>amass(1)=12.011,</div><div>flfrc='graphene.12x12x1.fc', </div><div>flfrq='graphene.disp.freq' </div><div>q_in_band_form=.true.,</div><div>/</div><div>4 ! number of q-points</div><div>0.000 0.000 0.000 10 ! gamma point</div><div>0.500 0.000 0.000 10 ! M</div><div>0.333 0.333 0.000 10 ! K</div><div>0.000 0.000 0.000  1 ! gamma point</div><div> </div><div>Plotband:</div><div> </div><div>graphene.disp.freq</div><div>0 1600</div><div>graphene.disp.freq.dat</div><div>graphene.disp.freq.ps</div><div>0</div><div>100 0</div><div> </div><div>Why the phonon dispersion does not match the pictures from the papers and for calculations with e.g. 18x18x1 k-point and nq-point sets the output file of graphene.dyn1 look like:</div><div> </div><div>Dielectric Tensor:</div><div> </div><div>************************        166.332153320312          0.000000000000</div><div>        166.332275390625************************          0.000000000000</div><div>          0.000000000000          0.000000000000          1.116060253916</div><div> </div><div>     Effective Charges E-U: Z_{alpha}{s,beta}</div><div> </div><div>     atom #    1</div><div>  -11117471.575740195811         -0.002632206306          0.000000000000</div><div>         -0.002632205375-11117471.579664016142            0.000000000000</div><div>          0.000000000000          0.000000000000         -0.003734609171</div><div>          </div><div><div>Also, what is the best way to set the lattice parameters, via Ibrav or CELL_PARAMETERS?</div></div><div>Any help is appreciated.</div><div> </div><div>Best regards,</div><div>Evgeniy Sysoev.</div></div></div>