<div dir="ltr">Dear QE users,<div>I am trying to reproduce the vibrational density of states of ice. As a first attempt I simulated a single water molecule in a large cubic cell with pw.x . Then I caluculated the vibrational modes with ph.x and then, after q2r.x, I used matdyn.w. The vibrational density seems consistent as there are delta functions in correspondence of the bending and antisimmetric stretching frequencies (as it is a single molecule and not crystal). The main problem is that in the vibrational density of states there is not a trace of the symmetric stretching that is present even in the normal mode calculation.</div><div>The structure has been also relaxed and different energy precision has been used for dos calculation.</div><div>here the input file used and a plot where the symmetric stretching at ca. 3800 cm-1 is missing and the bending and asymmetric stretching present.</div><div>-----</div><div>pw.x</div><div>------</div><div>prefix = 'h2o'<br>pseudo_dir = '../pseudo'<br>outdir= './out'<br>/<br>&SYSTEM<br>ecutwfc = 120<br>ibrav = 1<br>celldm(1) = 20<br>nat = 3<br>ntyp = 2<br>/<br>&ELECTRONS<br>conv_thr = 1.D-9<br>/<br>&IONS<br>/<br>ATOMIC_SPECIES<br>H 1 H_ONCV_PBE-1.0.upf<br>O 16 O_ONCV_PBE-1.0.upf<br>ATOMIC_POSITIONS angstrom<br>O 4.7571 4.0000 3.902317<br>H 5.5142 4.0000 4.488417<br>H 4.0000 4.0000 4.488417<br><br>K_POINTS (automatic)<br>1 1 1 0 0 0<br></div><div>-------</div><div>ph.x</div><div>-------</div><div>phonon calc<br>&inputph<br>outdir = './out'<br>prefix = 'h2o'<br>ldisp = .true.<br>amass(1) = 1<br>amass(2) = 16<br>nq1= 1, nq1 =1, nq3=1<br>fildyn = 'h2o.dyn'<br>/<br></div><div>---</div><div>q2r.x</div><div>------</div><div>&input<br>fildyn = 'h2o.dyn'<br>zasr = 'simple'<br>flfrc = 'h20.fc'<br>/<br></div><div>-------</div><div>matdyn.x</div><div>-----------</div><div>&input<br>asr = 'simple'<br>dos = .true.<br>amass(1) = 1<br>amass(2) = 16<br>flfrc = 'h2o.fc'<br>fldos = 'h2o.phdos'<br>nk1= 2, nk2 = 2, nk3 = 2<br>/<br></div><div><div><img src="cid:ii_k5qqyfqe0" alt="image.png" width="452" height="216" class="gmail-CToWUd gmail-a6T" tabindex="0" style="cursor: pointer; outline: 0px;"></div></div></div>