paw+gw in parallel + wannier90

GW, Bethe-Salpeter …

Moderators: maryam.azizi, bruneval

Locked
raul_l
Posts: 74
Joined: Sun Jan 08, 2012 7:45 pm

paw+gw in parallel + wannier90

Post by raul_l » Tue Dec 18, 2012 5:37 pm

(I removed this post from another thread and thought it would be better to create a new one.)

I modified the tests/wannier/Input/t03.in example by using paw and splitting the input file in half, one for calculating GW corrections in parallel and one for calculating MLWFs in serial. I used Si-LDA.paw from Psps_for_tests.
If I run

Code: Select all

# t03.in
# Silicon structure
     acell     10.263  10.263  10.263
     rprim      0.00   0.50   0.50
                0.50   0.00   0.50
                0.50   0.50   0.00
     natom      2
      xred      0.00   0.00   0.00
                0.25   0.25   0.25
    ntypat      1
     typat      1      1
     znucl     14.00
 symmorphi      0

# Parameters common to all runs
      ecut      6.00
 pawecutdg     12
    istwfk      8*1
     ngkpt      4      4      4
     nstep    100
   nshiftk      1
    shiftk      0.00   0.00   0.00
    enunit      2

    ndtset      6

# Self-consistent run to get the density
    toldfe1     1.00d-6

# Non-self-consistent run to get all cg wavefunctions
    getden2     1
    getwfk2     1
      iscf2     -2
    tolwfr2     1.0d-10
     nband2     30
   kssform2     3
  nbandkss2     27
    npwkss2     89

# Calculation of the dielectric matrix - iteration 1
 optdriver3     3
 gwcalctyp3    28
    getwfk3     2
     nband3    10
    npweps3    27
    ecutwfn3     1.49923969E+00 Hartree
      awtr3       0

# Calculation of the model GW corrections - iteration 1
 optdriver4     4
 gwcalctyp4    28
    getwfk4     2
    getscr4     3
     nband4    10
  ecutsigx4     0.0
   npwsigx4    27
   ecutwfn4     1.49923969E+00 Hartree

# Calculation of the dielectric matrix - iteration 2
 optdriver5     3
 gwcalctyp5    28
    getkss5     2
    getqps5     4

     nband5    10
    npweps5    27
   ecutwfn5     1.49923969E+00 Hartree
      awtr5       0

# Calculation of the model GW corrections - iteration 2
 optdriver6     4
 gwcalctyp6    28
    getwfk6     2
    getqps6     4
    getscr6     5
     nband6    10
  ecutsigx6     0.0
   npwsigx6    27
   ecutwfn6     1.49923969E+00 Hartree

#Common to all model GW calculations
   rhoqpmix     0.5
     nkptgw     8
      kptgw     
              0.00000000E+00  0.00000000E+00  0.00000000E+00
              2.50000000E-01  0.00000000E+00  0.00000000E+00
              5.00000000E-01  0.00000000E+00  0.00000000E+00
              2.50000000E-01  2.50000000E-01  0.00000000E+00
              5.00000000E-01  2.50000000E-01  0.00000000E+00
             -2.50000000E-01  2.50000000E-01  0.00000000E+00
              5.00000000E-01  5.00000000E-01  0.00000000E+00
             -2.50000000E-01  5.00000000E-01  2.50000000E-01
       bdgw   
              1 8
              1 8
              1 8
              1 8
              1 8
              1 8
              1 8
              1 8

I get

Code: Select all

*** glibc detected *** abinit: free(): invalid next size (fast): 0x00000000067f3220 ***
======= Backtrace: =========
/lib64/libc.so.6[0x3ce8275916]
/lib64/libc.so.6[0x3ce8278443]
/opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/intel64/libifcore.so.5(for_deallocate+0xb9)[0x7fbb5d3b7749]
/opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/intel64/libifcore.so.5(for_dealloc_allocatable+0x70)[0x7fbb5d3b7650]
abinit[0x92612d]
abinit[0x48fa99]
abinit[0x412b41]
abinit[0x40a43b]
abinit[0x40852c]
/lib64/libc.so.6(__libc_start_main+0xfd)[0x3ce821ecdd]
abinit[0x408429]
======= Memory map: ========
00400000-027fe000 r-xp 00000000 08:01 791275                             /usr/local/bin/abinit
029fd000-02b9c000 rw-p 023fd000 08:01 791275                             /usr/local/bin/abinit
02b9c000-05a9d000 rw-p 00000000 00:00 0
06689000-078ed000 rw-p 00000000 00:00 0                                  [heap]
3ce7e00000-3ce7e20000 r-xp 00000000 08:01 1177370                        /lib64/ld-2.12.so
3ce801f000-3ce8020000 r--p 0001f000 08:01 1177370                        /lib64/ld-2.12.so
3ce8020000-3ce8021000 rw-p 00020000 08:01 1177370                        /lib64/ld-2.12.so
3ce8021000-3ce8022000 rw-p 00000000 00:00 0
3ce8200000-3ce8389000 r-xp 00000000 08:01 1180960                        /lib64/libc-2.12.so
3ce8389000-3ce8589000 ---p 00189000 08:01 1180960                        /lib64/libc-2.12.so
3ce8589000-3ce858d000 r--p 00189000 08:01 1180960                        /lib64/libc-2.12.so
3ce858d000-3ce858e000 rw-p 0018d000 08:01 1180960                        /lib64/libc-2.12.so
3ce858e000-3ce8593000 rw-p 00000000 00:00 0
3ce8600000-3ce8683000 r-xp 00000000 08:01 1186763                        /lib64/libm-2.12.so
3ce8683000-3ce8882000 ---p 00083000 08:01 1186763                        /lib64/libm-2.12.so
3ce8882000-3ce8883000 r--p 00082000 08:01 1186763                        /lib64/libm-2.12.so
3ce8883000-3ce8884000 rw-p 00083000 08:01 1186763                        /lib64/libm-2.12.so
3ce8a00000-3ce8a02000 r-xp 00000000 08:01 1205524                        /lib64/libdl-2.12.so
3ce8a02000-3ce8c02000 ---p 00002000 08:01 1205524                        /lib64/libdl-2.12.so
3ce8c02000-3ce8c03000 r--p 00002000 08:01 1205524                        /lib64/libdl-2.12.so
3ce8c03000-3ce8c04000 rw-p 00003000 08:01 1205524                        /lib64/libdl-2.12.so
3ce8e00000-3ce8e17000 r-xp 00000000 08:01 1205516                        /lib64/libpthread-2.12.so
3ce8e17000-3ce9017000 ---p 00017000 08:01 1205516                        /lib64/libpthread-2.12.so
3ce9017000-3ce9018000 r--p 00017000 08:01 1205516                        /lib64/libpthread-2.12.so
3ce9018000-3ce9019000 rw-p 00018000 08:01 1205516                        /lib64/libpthread-2.12.so
3ce9019000-3ce901d000 rw-p 00000000 00:00 0
3ce9600000-3ce9607000 r-xp 00000000 08:01 1206259                        /lib64/librt-2.12.so
3ce9607000-3ce9806000 ---p 00007000 08:01 1206259                        /lib64/librt-2.12.so
3ce9806000-3ce9807000 r--p 00006000 08:01 1206259                        /lib64/librt-2.12.so
3ce9807000-3ce9808000 rw-p 00007000 08:01 1206259                        /lib64/librt-2.12.so
3cea200000-3cea216000 r-xp 00000000 08:01 1206256                        /lib64/libresolv-2.12.so
3cea216000-3cea416000 ---p 00016000 08:01 1206256                        /lib64/libresolv-2.12.so
3cea416000-3cea417000 r--p 00016000 08:01 1206256                        /lib64/libresolv-2.12.so
3cea417000-3cea418000 rw-p 00017000 08:01 1206256                        /lib64/libresolv-2.12.so
3cea418000-3cea41a000 rw-p 00000000 00:00 0
3cec200000-3cec216000 r-xp 00000000 08:01 1205515                        /lib64/libgcc_s-4.4.6-20110824.so.1
3cec216000-3cec415000 ---p 00016000 08:01 1205515                        /lib64/libgcc_s-4.4.6-20110824.so.1
3cec415000-3cec416000 rw-p 00015000 08:01 1205515                        /lib64/libgcc_s-4.4.6-20110824.so.1
7fbb58afe000-7fbb59010000 rw-p 00000000 00:00 0
7fbb59010000-7fbb59051000 rw-s 00000000 00:10 3410897                    /dev/shm/mpich_shar_tmp7gX3Ie (deleted)
7fbb59051000-7fbb59092000 rw-s 00000000 00:10 3410895                    /dev/shm/mpich_shar_tmpdHQUHe (deleted)
7fbb59092000-7fbb5ac26000 r-xp 00000000 08:01 544834                     /opt/intel/composer_xe_2011_sp1.9.293/mkl/lib/intel64/libmkl_mc3.so
7fbb5ac26000-7fbb5ae26000 ---p 01b94000 08:01 544834                     /opt/intel/composer_xe_2011_sp1.9.293/mkl/lib/intel64/libmkl_mc3.so
7fbb5ae26000-7fbb5ae7b000 rw-p 01b94000 08:01 544834                     /opt/intel/composer_xe_2011_sp1.9.293/mkl/lib/intel64/libmkl_mc3.so
7fbb5ae7b000-7fbb5ae7d000 rw-p 00000000 00:00 0
7fbb5ae7d000-7fbb5ae82000 r-xp 00000000 08:01 1177369                    /lib64/libnss_dns-2.12.so
7fbb5ae82000-7fbb5b081000 ---p 00005000 08:01 1177369                    /lib64/libnss_dns-2.12.so
7fbb5b081000-7fbb5b082000 r--p 00004000 08:01 1177369                    /lib64/libnss_dns-2.12.so
7fbb5b082000-7fbb5b083000 rw-p 00005000 08:01 1177369                    /lib64/libnss_dns-2.12.so
7fbb5b083000-7fbb5b085000 r-xp 00000000 08:01 1205065                    /lib64/libnss_mdns4_minimal.so.2
7fbb5b085000-7fbb5b284000 ---p 00002000 08:01 1205065                    /lib64/libnss_mdns4_minimal.so.2
7fbb5b284000-7fbb5b285000 rw-p 00001000 08:01 1205065                    /lib64/libnss_mdns4_minimal.so.2
7fbb5b285000-7fbb5b291000 r-xp 00000000 08:01 1199272                    /lib64/libnss_files-2.12.so
7fbb5b291000-7fbb5b491000 ---p 0000c000 08:01 1199272                    /lib64/libnss_files-2.12.so
7fbb5b491000-7fbb5b492000 r--p 0000c000 08:01 1199272                    /lib64/libnss_files-2.12.so
7fbb5b492000-7fbb5b493000 rw-p 0000d000 08:01 1199272                    /lib64/libnss_files-2.12.so
7fbb5b493000-7fbb5c554000 rw-s 00000000 00:10 3409789                    /dev/shm/mpich_shar_tmpMyV9sN (deleted)
7fbb5c554000-7fbb5c557000 rw-p 00000000 00:00 0
7fbb5c557000-7fbb5c5a0000 r-xp 00000000 08:01 525280                     /opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/intel64/libintlc.so.5
7fbb5c5a0000-7fbb5c69f000 ---p 00049000 08:01 525280                     /opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/intel64/libintlc.so.5

with abinit-7.0.4. The second input file is

Code: Select all

# t04.in
# Silicon structure, same as above
    ...

# Parameters common to all runs, same as above, except no ndtset
   ...

# Calculation of the quasiparticle Wannier functions
    irdden     1
    irdwfk     1
    irdqps     1
    kptopt     3
    istwfk     64*1
      iscf     -2
     nstep     0
    tolwfr     1.0d-10
     nband     10
   prtwant     3
 w90iniprj     2
 w90prtunk     0

but I didn't get a change to run that because the first one failed.
Could somebody share the input of some simple PAW+GW+WANNIER run where the PAW+GW part is done in parallel? Unless there's a bug I don't see what I'm doing wrong.
Raul Laasner
Netherlands Institute for Space Research

raul_l
Posts: 74
Joined: Sun Jan 08, 2012 7:45 pm

Re: paw+gw in parallel + wannier90

Post by raul_l » Wed Dec 19, 2012 11:14 am

Can anyone tell me whether it is possible in principle - a paw+gw run in parallel followed by wannier90? If it is not implemented in abinit then at least I know I shouldn't spend more time trying to find a mistake in my input.
Raul Laasner
Netherlands Institute for Space Research

raul_l
Posts: 74
Joined: Sun Jan 08, 2012 7:45 pm

Re: paw+gw in parallel + wannier90

Post by raul_l » Tue Jan 01, 2013 11:53 am

I have managed to do paw+gw+wannier90 using abinit-6.12.3, an unshifted k-mesh and the single-shot PP model (not for silicon but my own system). The resulting _GW file has to be renamed "g0w0" for Wannier interpolation. With abinit-7.0.4 I get the same error as above. With self-consistency I also get the same error.
Raul Laasner
Netherlands Institute for Space Research

Locked