<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dear<font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"> Krishnendu</font>,</p>
<p><br>
</p>
<p>in order to perform the so-called Fourier-interpolation you need first to compute phonon frequencies with DPTF on a proper grid of q-points using ph.x. From the dynamical matrix produced by the phonon code you can then compute the interatomic force constants
(IFC) in real space using the q2r.x code. From there, it is possible to recompute the phonon frequencies at any arbitrary q-point by means of the matdyn.x code. At this stage, you can specify the q-mesh (or a list of specific q-points) that is required for
calculating converged values of the quantities you are interested in (in our case a 21x21x21 mesh was sufficient). For more specific information, try to have a look at this very useful document by Paolo Giannozzi
<a href="http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_disp.html">http://www.fisica.uniud.it/~giannozz/QE-Tutorial/tutorial_disp.html</a><br>
</p>
<p><br>
</p>
<p>Both q2r.x and matdyn.x should be automatically compiled by typing "make ph" from the main quantum Espresso directory, as stated (for example) here
<a href="http://www.afs.enea.it/software/qu_esp/Doc/ph_6.0.0.pdf">http://www.afs.enea.it/software/qu_esp/Doc/ph_6.0.0.pdf</a></p>
<p><br>
</p>
<p>Examples of q2r and matdyn input files can be found here</p>
<p><a href="http://www.fisica.uniud.it/~giannozz/QE-Tutorial/examples_disp.tgz">http://www.fisica.uniud.it/~giannozz/QE-Tutorial/examples_disp.tgz</a><br>
</p>
<p>Have also a look at this tutorial (Exercise 2 in particular)</p>
<p><a href="http://www.tcl.t.u-tokyo.ac.jp/pdf/2014_lm/3rd_day_pm.pdf">http://www.tcl.t.u-tokyo.ac.jp/pdf/2014_lm/3rd_day_pm.pdf</a></p>
<p><br>
</p>
<p>In particular you can read:<br>
</p>
<p>##q2r input template</p>
<p>&INPUT<br>
fildyn = ***,<br>
flfrc = ***,<br>
zasr = ***, (we use crystal)<br>
/<br>
</p>
<p>##<br>
</p>
<p><br>
</p>
<p>## matdyn input template<br>
</p>
<p>&INPUT<br>
ntyp = ***, (we have 1 atom type)<br>
flvec = ***,<br>
flfrc = ***,<br>
asr = ***, (we use crystal)<br>
amass(1) = ***,<br>
/<br>
*** (integer # of q points for which you want compute phonon frequencies)<br>
x1 x2 x3 0 (here goes a list of q points coordinates followed by a growing integer number N counting the q-points. In our case, this list should correspond to a 21x21x21 MP grid)</p>
x1 x2 x3 1<br>
x1 x2 x3 2
<p>...</p>
<p>...<br>
</p>
<p>##</p>
<p>Replace *** with the correct input card. An exhaustive description of the input cards available for the two codes can be found in the header of the source codes q2r.f90 and matdyn.f90 in /yourQEfolder/PHonon/PH directory.<br>
</p>
<p><br>
</p>
<p>Hope this helps!<br>
</p>
<p>Best regards, </p>
<p> Daniele <br>
</p>
<div id="Signature">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div>
<div>
<div><font class="Apple-style-span" color="#333333" size="2">__________________________________________</font></div>
<div><font class="Apple-style-span" color="#333333" size="2"><br>
</font></div>
<div><font class="Apple-style-span" color="#333333" size="2">Daniele Francesco Dragoni</font></div>
<div><font class="Apple-style-span" color="#333333" size="2">Laboratoire de théorie et simulation des matériaux</font></div>
<div><font class="Apple-style-span" color="#333333" size="2">Swiss Federal Institute of Technology </font></div>
<div><font class="Apple-style-span" color="#333333" size="2">EPFL STI IMX THEOS</font></div>
<div><font class="Apple-style-span" color="#333333" size="2">ME D2 1020<br>
Station 9<br>
</font></div>
<div><font class="Apple-style-span" color="#333333" size="2">CH-1015, Lausanne<br>
+41 21 69 31956<br>
</font><a href="tel:+41216931956" id="NoLP"><span class="local-color-text"></span></a><br>
</div>
<div><font class="Apple-style-span" color="#333333" size="2"><br>
</font></div>
</div>
</div>
</div>
</div>
</div>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>From:</b> pw_forum-bounces@pwscf.org <pw_forum-bounces@pwscf.org> on behalf of Krishnendu Mukherjee <krishnendu.mukherjee789@gmail.com><br>
<b>Sent:</b> Tuesday, January 2, 2018 10:31 AM<br>
<b>To:</b> pw_forum@pwscf.org<br>
<b>Subject:</b> [Pw_forum] Fourier interpolation for DFPT</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div><br clear="all">
</div>
Dear experts,<br>
<br>
</div>
I am interested in calculating the phonon calculation of bcc iron ferrite. I am following the article titled "Thermoelastic properties of alpha-iron from first-principles" authored by Daniele Dragoni, Davide Ceresoli and Nicola Marzari, published in PRB 91,
104105. In that article it is mentioned that, <br>
<br>
</div>
"Phonon calculations were carried out for each deformation within DFPT. The dynamical matrix and its eigenvalues are calculated on an 4x4x4 mesh of special points in the BZ and Fourier-interpolated on an extended 21x21x21 grid for integration of thermodynamic
quantities."<br>
<br>
</div>
Now I gave the nq1=4, nq2=4, nq3=4 at <a href="http://fe.ph.in">fe.ph.in</a> according to the values mentioned in the article.
<br>
<br>
</div>
But I do not know the parameter for which I have to set the values 21x21x21. Can you kindly let me know the parameter where I have to put the values to set the Fourier-interpolation extension grid for integration of thermodynamic quantities.<br>
<br>
</div>
Thanking you,<br>
</div>
Yours sincerely,<br>
</div>
Krishnendu<br>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div><br>
-- <br>
<div class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div>Dr. Krishnendu Mukherjee,</div>
</div>
<div><br>
</div>
<div>Principal Scientist,</div>
<div>CSIR-NML,</div>
<div>Jamshedpur.</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>