[Pw_forum] number of kpoints for lambda.x

yelena yelena at ipb.ac.rs
Mon Feb 3 15:58:10 CET 2014


Dear Lorenzo,
Thank you for answer. I'm aware of that fact, I read almost every post 
written of pw_forum about this subject( what are k and k+q points and 
what exactly kpoints.x does, symemetry...)but still I have problem 
making input for lambda.x.
So I generate kpoints.x mesh, there is 19 points and I use it in input 
file. After that goes elph. files. there are 43 of them. And when I 
execute, I get error

At line 157 of file lambda.f90 (unit = 5, file = 'stdin')
Fortran runtime error: Bad real number in item 1 of list input

I checked my elph files, they are good, no NAN, no missing characters.
My input goes like this:
20 0.12  0
    19
      0.0000000  0.0000000  0.0000000   1.00
      0.0833333  0.0481125  0.0000000   6.00
      0.1666667  0.0962250  0.0000000   6.00
      0.2500000  0.1443376  0.0000000   6.00
      0.3333333  0.1924501  0.0000000   6.00
      0.4166667  0.2405626  0.0000000   6.00
      0.5000000  0.2886751  0.0000000   3.00
      0.0833333  0.1443376  0.0000000   6.00
      0.1666667  0.1924501  0.0000000  12.00
      0.2500000  0.2405626  0.0000000  12.00
      0.3333333  0.2886751  0.0000000  12.00
      0.4166667  0.3367877  0.0000000  12.00
      0.1666667  0.2886751  0.0000000   6.00
      0.2500000  0.3367877  0.0000000  12.00
      0.3333333  0.3849002  0.0000000  12.00
      0.4166667  0.4330127  0.0000000   6.00
      0.2500000  0.4330127  0.0000000   6.00
      0.3333333  0.4811252  0.0000000  12.00
      0.3333333  0.5773503  0.0000000   2.00
elph. 0.000000. 0.000000. 0.000000
elph. 0.000000. 0.096225. 0.000000
elph. 0.000000. 0.192450. 0.000000
elph. 0.000000. 0.288675. 0.000000
elph. 0.000000. 0.384900. 0.000000
elph. 0.000000. 0.481125. 0.000000
elph. 0.000000.-0.577350. 0.000000
elph. 0.083333. 0.144338. 0.000000
elph. 0.083333. 0.240563. 0.000000
elph. 0.083333. 0.336788. 0.000000
elph. 0.083333. 0.433013. 0.000000
elph. 0.083333. 0.529238. 0.000000
elph. 0.083333.-0.048113. 0.000000
elph. 0.083333.-0.144338. 0.000000
elph. 0.083333.-0.240563. 0.000000
elph. 0.083333.-0.336788. 0.000000
elph. 0.083333.-0.433013. 0.000000
elph. 0.083333.-0.529238. 0.000000
elph. 0.166667. 0.288675. 0.000000
elph. 0.166667. 0.384900. 0.000000
elph. 0.166667. 0.481125. 0.000000
elph. 0.166667. 0.577350. 0.000000
elph. 0.166667.-0.096225. 0.000000
elph. 0.166667.-0.192450. 0.000000
elph. 0.166667.-0.288675. 0.000000
elph. 0.166667.-0.384900. 0.000000
elph. 0.166667.-0.481125. 0.000000
elph. 0.250000. 0.433013. 0.000000
elph. 0.250000. 0.529238. 0.000000
elph. 0.250000. 0.625463. 0.000000
elph. 0.250000.-0.144338. 0.000000
elph. 0.250000.-0.240563. 0.000000
elph. 0.250000.-0.336788. 0.000000
elph. 0.250000.-0.433013. 0.000000
elph. 0.333333. 0.577350. 0.000000
elph. 0.333333. 0.673575. 0.000000
elph. 0.333333.-0.192450. 0.000000
elph. 0.333333.-0.288675. 0.000000
elph. 0.333333.-0.384900. 0.000000
elph. 0.416667. 0.721688. 0.000000
elph. 0.416667.-0.240563. 0.000000
elph. 0.416667.-0.336788. 0.000000
elph.-0.500000.-0.866025. 0.000000
  0.10

What am I missing?
Best,
Jelena


On 31 Jan 2014 13:01, Lorenzo Paulatto wrote:
> On 01/31/2014 01:02 PM, yelena wrote:
>> Hi!
>> I'm trying to make lambda.x input, but I can't get proper number of
>> kpoints using kpoints.x. Is there any other way to calculate kpoints 
>> for
>> lambda.x?
>> In electron phonon calculation I'm using 12 12 1 grid and get 43 q
>> points and when I do kpoints.x for 12 12 1 I get 19 points. I'm 
>> having
>> hexagonal structure with 7 atoms. What can I do?
>>
> Dear Yelena,
> the electron phonon calculation may be using both the k and k+q
> points, eventually with different symmerty. They are not necessary 
> the
> same produced by kpoints.x
>
> bests




More information about the users mailing list