<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"></head><body>
    
<div>Hello Lachlan,</div><div><br></div><div>Yes this seems a good approach. I asked about dispersion because I don't know if the phonon tool is available in the dftb+ release yet. I think not. On the other hand my tool is based on supercell replica but if the supercell is large the computation of the Hessian becomes quite long and has to be broken for safety by selecting movable atoms. If interested let us know.</div><div>With 'modes' you only get the gamma point frequencies. </div><div><br></div><div>Kind regards,</div><div><br></div><div>Alessandro</div><div><br></div><div><br></div><div><br></div><div id="composer_signature"><div style="font-size:88%;color:#364f67" dir="auto">Inviato dal mio dispositivo Samsung</div></div><br><br>-------- Messaggio originale --------<br>Da: Lachlan Oberg <u5805496@anu.edu.au> <br>Data: 26/07/17  03:56  (GMT+01:00) <br>A: User list for DFTB+ related questions <dftb-plus-user@mailman.zfn.uni-bremen.de> <br>Oggetto: Re: [DFTB-Plus-User] R:  Phonon Mode Calculations <br><br>
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Hello Alessandro,</p>
<p><br>
</p>
<p>Thank you for your reply. Yes, we are looking to produce phonon dispersion relations in periodic structures. Afterwards, we were looking to produce phonon dispersion relations in nanowire structures (approximately 50x50 Angstrom). Is this possible if we
 were to obtain accurate skf files (i.e., if we first optimised to known phonon dispersion relations, say, in the bulk material)?</p>
<p><br>
</p>
<p>Kind regards</p>
<p>Lachlan Oberg<br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> DFTB-Plus-User <dftb-plus-user-bounces@mailman.zfn.uni-bremen.de> on behalf of alessandro.pecchia <alessandro.pecchia@ismn.cnr.it><br>
<b>Sent:</b> Wednesday, July 26, 2017 1:18:07 AM<br>
<b>To:</b> User list for DFTB+ related questions<br>
<b>Subject:</b> [DFTB-Plus-User] R: Phonon Mode Calculations</font>
<div> </div>
</div>
<div>
<div>Dear Lachlan,</div>
<div><br>
</div>
<div>you need sk-files with repulsive part included otherwise you do not get the right forces. Then in the input you set SecondDerivatives as diver. The Hessian is currently computed with finite differences.</div>
<div>The 'modes' code uses a LAPACK eigensolver for real symmetric matrices, if I remember correctly. </div>
<div>Good relaxed structures and electronic bandstructures do not necessarily implay accurate vibrations as that's related to the curvature around the minimum, harder to get in many cases. It depends if you can live with about 10-20% accurate frequencies. Do
 you need phonon dispersions ?</div>
<div><br>
</div>
<div>Kind Regards,</div>
<div><br>
</div>
<div>Alessandro Pecchia</div>
<div><br>
</div>
<div><br>
</div>
<div id="composer_signature">
<div style="font-size:88%;color:#364f67" dir="auto">Inviato dal mio dispositivo Samsung</div>
</div>
<br>
<br>
-------- Messaggio originale --------<br>
Da: Lachlan Oberg <u5805496@anu.edu.au> <br>
Data: 24/07/17 07:31 (GMT+01:00) <br>
A: dftb-plus-user@mailman.zfn.uni-bremen.de <br>
Oggetto: [DFTB-Plus-User] Phonon Mode Calculations <br>
<br>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Hello,</p>
<p><br>
</p>
<p>I have two questions regarding the MODES algorithm used to calculate phonon modes. Firstly, what method does it use to calculate the eigenmodes given the Hessian? Secondly, what part of the .skf files is required to determine the Hessian? Are the integral
 tables sufficient, or do I explicitly require .skf files whose geometry can be optimized (i.e., they contain the
<span>repulsive spline interaction)? For example, does a faithful reproduction of the conduction/valence bands necessarily mean
</span>a faithful reproduction of the modes?<br>
</p>
<p><br>
</p>
<p>Kind regards</p>
<p>Lachlan Oberg<br>
</p>
</div>
</div>
</body></html>