<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7653.38">
<TITLE>Phonons at generic q-point</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<P><FONT SIZE=2>Dear All,<BR>
<BR>
I have been trying to run a parallel phonon calculation at generic q-point for a SrTiO3 multilayer (16 atoms) using espresso 4.0.3. I run an initial scf calculation with wf_collect=true and then a ph.x calculation with lnscf=true. The calculation always crashes with the following error (full input pasted below):<BR>
<BR>
! ... rest of modes<BR>
<BR>
Representation 47 1 modes - To be done<BR>
<BR>
Representation 48 1 modes - To be done<BR>
<BR>
from openfilq : error # 1<BR>
file sto_mag.wfc not found<BR>
stopping ...<BR>
<BR>
The wfc files are present and written to the same directory from which the calculation is running and this is not a 'recovered' calculation.<BR>
<BR>
So far I have tried the following tests:<BR>
<BR>
* The examples for Si at <A HREF="http://people.sissa.it/~degironc/QE-Tutorial/">http://people.sissa.it/~degironc/QE-Tutorial/</A> (all run successfully);<BR>
* Phonon dispersion and q=0 calculations for bulk SrTiO3 (all run successfully);<BR>
* Phonon calculation at M for bulk SrTiO3, with lnscf=true (run successfully);<BR>
* Phonon calculation at q=0 for SrTiO3 multilayer (run successfully).<BR>
<BR>
I also noticed the message 'Found additional translation' printed at the top of the output. I had a look through the source code (sgam_at.f90 and sgam_at_mag.f90), so I know this message is printed if the identity is found to have fractional translations. I found one previous post about this message in the archive (<A HREF="http://tinyurl.com/5k926j">http://tinyurl.com/5k926j</A>) but it seemed to concern how parallel jobs were re-started before the wf_collect keyword existed (I run all my calculations with wf_collect=true and my phonon calculation uses the same number of processors as the initial scf run). Is this message the source of my problems or is there something else going on? Any hints would be appreciated.<BR>
<BR>
Cheers,<BR>
<BR>
Nicole<BR>
<BR>
Nicole Benedek<BR>
Research Associate<BR>
Department of Materials<BR>
Imperial College London<BR>
<BR>
! INPUT FILE FOR PHONON CALC<BR>
<BR>
Phonons of SrTiO3 multilayer at X<BR>
&inputph<BR>
tr2_ph=1.0d-14,<BR>
amass(1)=87.62,<BR>
amass(2)=47.88,<BR>
amass(3)=15.9994,<BR>
prefix='sto_mag',<BR>
outdir='./'<BR>
fildyn='sto_mag.dynX',<BR>
lnscf=.true.<BR>
/<BR>
0.0 0.5 0.0 ! X point in BZ for orthorhombic primitive systems<BR>
</FONT>
</P>
</BODY>
</HTML>