<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Paul,<div class=""><br class=""></div><div class="">Just to clarify this point:</div><div class=""><br class=""></div><div class="">> While Vasp can carry out Raman calculations using phon3py […]</div><div class=""><br class=""></div><div class="">Phono3py can interface Quantum ESPRESSO as well, since this code simply needs forces of supercells with finite displacements to evaluate third order force constants. (See also its <a href="https://phonopy.github.io/phono3py/" class="">doc</a>). These are needed in case you would like to have the phonon line widths at different temperatures.</div><div class=""><br class=""></div><div class="">While the computation of phonon line widths is very onerous indeed, it is usually not needed for a comparison with experiments. Moreover, you don’t need them for the computation of the Raman intensities. For this, you need to evaluate susceptibility derivatives in respect to atomic displacements, often called Raman tensors. </div><div class=""><br class=""></div><div class="">These tensors are third order derivatives of the total energy in respect to 2 electric fields and 1 atomic displacement. Thus, there are multiple way to compute them: 2n+1 theorem, derivatives of the charge density and finite differences.</div><div class=""><br class=""></div><div class="">Usually, the implementation of the former 2 are much more efficient, but also very involved in the implementation. This is the reason why in plane-waves codes their implementations are limited to LDA(+U). In Quantum ESPRESSO we are currently limited to LDA and NC pseudos.</div><div class=""><br class=""></div><div class="">To overcome this limitation on the functional (and pseudo), finite differences can still be employed, using e.g. the electric enthalpy approach (see Umari and Pasquarello, Diamond & Related Materials 14 (2005) 1255–1261). There are already some routines in Quantum ESPRESSO, in the PHonon/FD folder, to carry out these calculations semi-automatically. </div><div class=""><br class=""></div><div class="">Since the approach requires many calculations (and there are some pitfalls as well), one year ago I and my supervisor (Prof. Nicola Marzari) decided to implement in the <a href="https://www.aiida.net/" class="">AiiDA</a> framework a workflow for doing all these calculations automatically. The package is a python plugin which extends the functionalities of AiiDA and its aiida-quantuemespresso plugin; the latter already manages many of the usual QE executables). </div><div class=""><br class=""></div><div class="">This plugin will be made open source very soon at the link <a href="https://github.com/bastonero/aiida-vibroscopy" class="">https://github.com/bastonero/aiida-vibroscopy</a>, and also on the AiiDA registry as aiida-vibroscopy.</div><div class=""><br class=""></div><div class="">I hope this is and will be useful.</div><div class=""><br class=""></div><div class="">Best regards,</div><div class="">Lorenzo Bastonero</div><div class=""><br class=""><div class="">
<meta charset="UTF-8" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div><div class="">********************************</div><div class="">Lorenzo Bastonero</div><div class=""><br class=""></div><div class="">PhD Student</div><div class="">U Bremen Excellence Chair,<br class="">Bremen Center for Computational Materials Science,<br class="">and MAPEX Center for Materials and Processes</div><div class=""><br class=""></div><div class="">University of Bremen</div><div class="">Faculty of Production Engineering</div><div class="">TAB-Building, Room 3.32</div><div class="">Am Fallturm 1</div><div class="">28359 Bremen, DE</div><br class="Apple-interchange-newline"><a href="http://www.hmi.uni-bremen.de/" class="">http://www.hmi.uni-bremen.de/</a></div></div></div></div></div></div>
</div>
<br class=""></div></body></html>