<span style="background-color:rgb(255,255,255);color:rgb(0,0,0)">Hi PW list,</span><br style="background-color:rgb(255,255,255);color:rgb(0,0,0)"><br style="background-color:rgb(255,255,255);color:rgb(0,0,0)"><span style="color:rgb(0,0,0);background-color:rgb(255,255,255)">I am working with Quantum Espresso (v4.3.2.1) in the hopes of attaining Raman phonon frequencies of a benzoic acid crystal. I attained the optimized output geometry from a CASTEP calculation in the supplemental material of a paper. This geometry was used by the authors to calculate accurate IR phonon frequencies. I used a 'relax' calculation in Quantum Espresso on the optimized geometry but this led to large shifts in the atomic positions in the lattice and unreasonable structures for the molecules in the lattice after finding the optimized structure. My question is does anyone have a some experience calculating organic crystals in Quantum Espresso that could give me some tips? <br>
<br>Here is my input file minus the atomic posititions:</span><br><br>&CONTROL<br> title = 'benzacidscftest' ,<br> calculation = 'relax' ,<br> restart_mode = 'from_scratch' ,<br>
outdir = '/home/mreish/tmp/' ,<br> pseudo_dir = '/usr/share/espresso/pseudo/' ,<br> prefix = 'batest5' ,<br> disk_io = 'default' ,<br>
nstep = 600 ,<br> /<br> &SYSTEM<br> ibrav = 12,<br> A = 5.4996 ,<br> B = 5.1283 ,<br> C = 21.950 ,<br>
cosAB = 0 ,<br> cosAC = -.128276341 ,<br> cosBC = 0 ,<br> nat = 60,<br> ntyp = 3,<br> ecutwfc = 5 ,<br>
ecutrho = 40 ,<br> /<br> &ELECTRONS<br> electron_maxstep = 200,<br> conv_thr = 1.0d-11 ,<br> /<br> &IONS<br> /<br>ATOMIC_SPECIES<br> C 12.01000 C.pbe-rrkjus.UPF <br>
O 16.00000 O.pbe-rrkjus.UPF <br> H 1.00000 H.pbe-rrkjus.UPF <br>ATOMIC_POSITIONS angstrom <br> <br>K_POINTS automatic <br> 1 1 1 1 1 1<br><br><br>Cheers,<br><br>Matthew Reish<br>University of Otago<br>
Dunedin<br>New Zealand<br>