Dynamic Dielectric Matrix [05] Free DIP_iR stops code
Posted: Wed Jun 17, 2015 12:06 pm
Hello,
I am running an MoS2 system in order work my way through GW and BSE calculations (see viewtopic.php?f=14&t=1019 ). I am increasing the number of k points with which the calculation is performed (9x9, 18x18, 27x27...) in order to get accurate band structure values for the system. Recently, I ran a 54x54 grid and began the convergence process on the same variables under the configuration conditions and the system stops running. Where it stops running is at
[05] Dynamic Dielectric Matrix (PPA)
Alloc DIP_iR
Free DIP_iR
I watch the system memory on the nodes running the calculation and there is no undue stress on the either the RAM, CPU or SWAP at this point in the process and the program does not write any information on stop. Now, I know that anyone here probably needs more information about the system to help with the crash, but my question is more about what is happening in the code when it crashes.
When I look at all other output files (l_em1d_ppa_HF_and_locXC_gw0_##) for all of my other runs, there is no mention of DIP_iR. A novice reading of this seems to suggest that increased MP grid triggers a different routine and that routine is causing me problems. What is DIP_iR and why have I not encountered it in my previous output files?
Additional information:
total number of k-points in run: 1460
Input file:
HF_and_locXC # [R XX] Hartree-Fock Self-energy and Vxc
ppa # [R Xp] Plasmon Pole Approximation
gw0 # [R GW] GoWo Quasiparticle energy levels
em1d # [R Xd] Dynamical Inverse Dielectric Matrix
FFTGvecs=100
EXXRLvcs=100 RL # [XX] Exchange RL components
Chimod= "Hartree" # [X] IP/Hartree/ALDA/LRC/BSfxc
% BndsRnXp
1 | 400 | # [Xp] Polarization function bands
%
NGsBlkXp= 25 RL # [Xp] Response block size
% LongDrXp
1.000000 | 0.000000 | 0.000000 | # [Xp] [cc] Electric Field
%
PPAPntXp= 27.21138 eV # [Xp] PPA imaginary energy
% GbndRnge
1 | 400 | # [GW] G[W] bands range
GDamping= 0.10000 eV # [GW] G[W] damping
dScStep= 0.10000 eV # [GW] Energy step to evalute Z factors
DysSolver= "n" # [GW] Dyson Equation solver (`n`,`s`,`g`)
%QPkrange # [GW] QP generalized Kpoint/Band indices
623|623| 9|10|
965|965| 9|10|
I am running an MoS2 system in order work my way through GW and BSE calculations (see viewtopic.php?f=14&t=1019 ). I am increasing the number of k points with which the calculation is performed (9x9, 18x18, 27x27...) in order to get accurate band structure values for the system. Recently, I ran a 54x54 grid and began the convergence process on the same variables under the configuration conditions and the system stops running. Where it stops running is at
[05] Dynamic Dielectric Matrix (PPA)
Alloc DIP_iR
Free DIP_iR
I watch the system memory on the nodes running the calculation and there is no undue stress on the either the RAM, CPU or SWAP at this point in the process and the program does not write any information on stop. Now, I know that anyone here probably needs more information about the system to help with the crash, but my question is more about what is happening in the code when it crashes.
When I look at all other output files (l_em1d_ppa_HF_and_locXC_gw0_##) for all of my other runs, there is no mention of DIP_iR. A novice reading of this seems to suggest that increased MP grid triggers a different routine and that routine is causing me problems. What is DIP_iR and why have I not encountered it in my previous output files?
Additional information:
total number of k-points in run: 1460
Input file:
HF_and_locXC # [R XX] Hartree-Fock Self-energy and Vxc
ppa # [R Xp] Plasmon Pole Approximation
gw0 # [R GW] GoWo Quasiparticle energy levels
em1d # [R Xd] Dynamical Inverse Dielectric Matrix
FFTGvecs=100
EXXRLvcs=100 RL # [XX] Exchange RL components
Chimod= "Hartree" # [X] IP/Hartree/ALDA/LRC/BSfxc
% BndsRnXp
1 | 400 | # [Xp] Polarization function bands
%
NGsBlkXp= 25 RL # [Xp] Response block size
% LongDrXp
1.000000 | 0.000000 | 0.000000 | # [Xp] [cc] Electric Field
%
PPAPntXp= 27.21138 eV # [Xp] PPA imaginary energy
% GbndRnge
1 | 400 | # [GW] G[W] bands range
GDamping= 0.10000 eV # [GW] G[W] damping
dScStep= 0.10000 eV # [GW] Energy step to evalute Z factors
DysSolver= "n" # [GW] Dyson Equation solver (`n`,`s`,`g`)
%QPkrange # [GW] QP generalized Kpoint/Band indices
623|623| 9|10|
965|965| 9|10|