<div class="gmail_quote">On Tue, Jul 10, 2012 at 4:44 PM, Torstein Fjermestad <span dir="ltr"><<a href="mailto:torstein.fjermestad@kjemi.uio.no" target="_blank">torstein.fjermestad@kjemi.uio.no</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> The problem is that I have so far only managed to complete step 1, to<br>
 compute the phonons. I am currently doing test calculations on a water<br>
 molecule and this is the input file of the phonon calculation:<br>
<br>... </blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> This calculation crashes presumably because it expects a file called<br>
 h2o.dynG0. But from where do I get this file? As far as I can see the<br>
 phonon calculation only produces the file h2o.dynG, not a h2o.dynG0<br>
 file.<br><br></blockquote><div><br></div><div> Dear Torstein,</div><div>the problem in your case is that q2r.x expects to find a phonon *dispersion* calculation, i.e. the phonon computed over a regular grid of q-points n reciprocal space. On the other hand, you only computed the phonons at Gamma. It is true that Gamma is just a 1x1x1 grid, and if you are sure it is enough to describe your system you can manage to convince q2r to use it. You need to do the following:</div>
<div>1. move h2o.dynG to h2o.dynG1</div><div>2. create a file called h2o.dynG0 which this inside:</div><div> 1 1 1</div><div> 1</div><div> 0.0 0.0 0.0</div><div>3. try again with q2r.x</div><div><br></div><div>Alternatively, you can just repeat the phonon calculation with ldisp = .true.; if you want your grid to be just the Gamma point you specify nq1=1, nq2=1, nq3=1.</div>
<div><br></div><div>I stress that a single-point sampling may not be enough to describe accurately the properties of your system (unless the supercell is very large) but this you need to check for yourself.</div><div><br>
</div><div>good work</div></div><br clear="all"><div><br></div>-- <br><div>Lorenzo Paulatto IdR @ IMPMC/CNRS & Université Paris 6</div><div>phone: +33 (0)1 44275 084 / skype: paulatz</div><div>www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" target="_blank">http://www-int.impmc.upmc.fr/~paulatto/</a></div>
<div>mail:  23-24/4é16 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05</div><br>