[Pw_forum] Error in phonon calculation using matdyn.x
Stefano de Gironcoli
degironc at sissa.it
Tue Mar 29 15:14:46 CEST 2011
increasing the threshold to 1.d-7 will probably not help..
I think the total weight you have is 7.5 instead of 8.0...
is your cell very anisotropic ?
try to increase the VERY VERY SAFE RANGE to 3 times nr1,nr2,nr3 ...
stefano
On 03/29/2011 03:00 PM, Eyvaz Isaev wrote:
> Hi,
>
> Find subroutine frc_blk(...) inside matdyn.f90 the next line:
>
> IF (ABS(total_weight-nr1*nr2*nr3).GT.1.0d-8) THEN
>
> Might be due to some numerical noise this condition is not fulfilled. Then try
> 1.d0d-7 (recompile matdyn.f90 using "make ph") and see what happens. If not,
> then one has to analyze your input files.
>
> Bests,
> Eyvaz.
>
> -------------------------------------------------------------------
> Prof. Eyvaz Isaev,
> Department of Physics, Chemistry, and Biology (IFM), Linkoping University,
> Sweden
>
> Theoretical Physics Department, Moscow State Institute of Steel& Alloys,
> Russia,
>
> isaev at ifm.liu.se, eyvaz_isaev at yahoo.com
>
>
>
> ----- Original Message ----
> From: mayank gupta<mayankaditya at gmail.com>
> To: pw_forum at pwscf.org
> Sent: Tue, March 29, 2011 6:00:15 AM
> Subject: [Pw_forum] Error in phonon calculation using matdyn.x
>
> Dear All I am doing pwscf and phonon calculation for Y2O3. In 2x2x2 Q
> mesh it generates 6 irreducible q points (yo.dyn* files). The phonon
> frequency calculation using matdyn.x run I got error message
>
> 7.5000000000000000
> from frc_blk : error # 1
> wrong total_weight
> stopping ...
>
> What could be the possible error.
>
>
> Thanks
>
>
>
>
>
>
>
More information about the users
mailing list