<div dir="ltr"><div><div>Hello everyone,<br><br></div>I'm trying to compute phonon dispersion properties of graphene nanoribbon using Quantum Espresso. My unit cell is a simple armchair structure with 4 atoms per unit cell. It looks like below:<br><br></div><div>       0----0<br></div><div>      /       \ <br>     0        0<br><br></div><div>The above structure is replicated in the y-direction (the periodic direction) while the x-direction is confined. Using the above structure, here is the outline of the problems I'm running into:<br><br></div><div>1) Negative phonon frequencies (not only at gamma point, but at all other k-points)<br></div><div><br></div><div>2) I tried using 'asr' but that doesn't seem to fix the problem.<br><br></div><div>3) I understand that the problem might be due to instability of the structure (I did notice that the forces on atoms in scf calculation is non-zero , forces of the magnitude 0.2 to 0.4).<br><br></div><div>4) When I try to relax the structure using 'ion_dynamics', the forces did go down but the structure is relaxed so much that it is no longer graphene nanoribbon.<br><br></div><div>Can you please give me some pointers on how to approach this problem, i.e. how do I get rid of negative frequencies while relaxing the structure that constrains it to look like a nanoribbon ?<br><br></div><div>Thanks for your time!<br></div><div><br></div><div>Here are my input files:<br></div><div><br>pw.x (this doesn't include the ion_dynamics relaxation).<br><br> &control<br>    calculation='scf'<br>    restart_mode='from_scratch',<br>    tstress = .true.<br>    tprnfor = .true.<br>    prefix='gnr_arm',<br>    pseudo_dir = '/work/02420/pm24229/QE/QE/pseudo/',<br>    outdir = '/home1/02420/pm24229/tmp/',<br> /<br> &system<br>    ibrav=  0, nat=  4, ntyp= 1, nbnd= 20,<br>    ecutwfc =25.0, occupations='smearing', smearing='marzari-vanderbilt',<br>    degauss=0.02<br> /<br> &electrons<br>    diagonalization='cg'<br>    conv_thr =  1.0d-12<br>    mixing_beta = 0.7<br> /<br>ATOMIC_SPECIES<br> C  12.0107  C.pbe-rrkjus.UPF<br><br>CELL_PARAMETERS {angstrom}<br> 2.5546 0.00000000 0.00000000<br> 0.00000000 10.0 0.00000000<br> 0.00000000 0.00000000 2.00000000<br><br>ATOMIC_POSITIONS {angstrom}<br> C 0.0  0.0  0.0<br> C 1.2773 0.7375 0.000<br> C 1.2773 2.2125 0.00<br> C 0.0  2.95 0.0<br><br>K_POINTS automatic<br>20 1 1 1 1 1<br><br></div><div>ph.x<br><br>phonons on kspace<br> &inputph<br>  tr2_ph=1.0d-14,<br>  prefix='gnr_arm',<br>  ldisp=.true.,<br>  nq1=20, nq2=1, nq3=1<br>  amass(1)=12.0107,<br>  outdir='/home1/02420/pm24229/tmp/',<br>  fildyn='gnr_dyn',<br> /<br><br></div><div>Thanks,<br>Prabhakar Marepalli<br></div><div><br><br></div></div>