[Pw_forum] Problem with high accuracy and parallel ph.x (Mohammad Saghayezhian)

Andrea Dal Corso dalcorso at sissa.it
Fri Jul 1 13:03:34 CEST 2011


Which is the error message, if there is one. Note that the phonon code
at a low symmetry q point requires much more k points than the pw.x code.

Andrea


Quoting Mohammad Saghayezhian <m.saghayezhian at gmail.com>:

> Dear Andrea,
>
> Thanks for reply.
>
> The nkp in my calculation is far less than 40000. As I know, the symmetry
> reduces nkp,
> therefore for a mesh of 30*30*28 the max nkp is 25200, So if the problem
> rooted in max nkp parameter,
> above mesh should not make any problem, there are some
> tests that I did for checking the valid nkp and cut-off energy:
> Cut-off(Ry) , Kpoints: 33, 22*22*18 (OK)
>                                 35, 25*25*22 (OK)
>                                 35, 30*30*26 ()
>                                 38, 22*22*18 (OK)
>                                 38, 30*30*26 (Failed)
>                                 40, 30*30*26 (Failed)
> All of above number kpoints are lower that nkp max (40000). I also run the
> calculation
> without npool option to provide higher memory for large cut-off. but still I
> wonder why
> the last two calculation did not worked!
>
> Could it be due to optimization deficiency in compilation of ph.x files?
> Please take a look at this links:
> http://www.democritos.it/pipermail/pw_forum/2005-February/002070.html
>  http://www.democritos.it/pipermail/pw_forum/2005-February/002074.html
> Thanks a lot
>
>
>
> On Thu, Fri, 01 Jul 2011 10:01:48 +0200,  Andrea Dal Corso wrote:
>> Please note that the maximum number of k points is still a parameter in
>> QE (presently set to 40000), so to make phonon calculations with a very
>> large mesh of k points you have to increase the value of npk in
>> Modules/parameters.f90.
>
>
>> HTH,
>
>> Andrea
>
>
>
> On Thu, 2011-06-30 at 11:58 +0330, Mohammad Saghayezhian wrote:
>> Dear All,
>>
>> When I perform the scf calculation with high kp (45 45 45) and cut-off
>> (42 Ry), the parallel execution ph.x fails and stops (q-points in
>> ph.in: 4 4 3). but lowering the kp (15 15 15) and cut-off (22 Ry)
>> solves the problem and the parallel ph.x goes on without any problem.
>> Does this mean that I can not calculate phonon spectra with this
>> accuracy?
>> Also, the calculations are for metals.
> --
> -------------------------------------------------------------------
> Mohammad Saghayezhian
> Computational Condensed Matter Research Lab
> Physics Department, Isfahan University of Technology, Isfahan, Iran
>
> Tel lab: +98 311 391 3731     Fax Office: +98311 391 3746
> -------------------------------------------------------------------
>



----------------------------------------------------------------
   SISSA Webmail https://webmail.sissa.it/
   Powered by Horde http://www.horde.org/





More information about the users mailing list