homogeneous vs. randomized k/q-mesh in BZ integration
Posted: Sat Aug 20, 2016 7:35 pm
Dear EPW-team,
first of all congratulations to the refurbished, extended version of EPW. I was using EPW2.3
for quite some time as a benchmark for my own clutter of el-ph code, but the new version is
really a pleasure.
My question:
In the example for doped Diamond http://epw.org.uk/Documentation/B-dopedDiamond
in Fig. 3f) the phonon linewidths for the highest phonon mode is shown - randomized 50^3
Sobol k-points and a smearing of 100meV was used. I note \gamma:max~7meV around
the \Gamma point.
In the next Fig. 4, the same calculation is done with an homogeneous/unshifted 50^3 grid.
Here I note \gamma:max~11meV around the \Gamma point.
While I know that convergence with random meshes is much faster and a 50^3 homogenous
grid is far from convergence, I wonder why the discrepancy between both "integration method"
in the absolute value of the selfenergy is so large.
Is there a certain criterion (or empirical knowledge) on the convergence dependency
of Brillouin-zone integrations with randomized/homogenous meshes?
thanks for your effort!
Nicki
CAMd, DTU Lyngby
first of all congratulations to the refurbished, extended version of EPW. I was using EPW2.3
for quite some time as a benchmark for my own clutter of el-ph code, but the new version is
really a pleasure.
My question:
In the example for doped Diamond http://epw.org.uk/Documentation/B-dopedDiamond
in Fig. 3f) the phonon linewidths for the highest phonon mode is shown - randomized 50^3
Sobol k-points and a smearing of 100meV was used. I note \gamma:max~7meV around
the \Gamma point.
In the next Fig. 4, the same calculation is done with an homogeneous/unshifted 50^3 grid.
Here I note \gamma:max~11meV around the \Gamma point.
While I know that convergence with random meshes is much faster and a 50^3 homogenous
grid is far from convergence, I wonder why the discrepancy between both "integration method"
in the absolute value of the selfenergy is so large.
Is there a certain criterion (or empirical knowledge) on the convergence dependency
of Brillouin-zone integrations with randomized/homogenous meshes?
thanks for your effort!
Nicki
CAMd, DTU Lyngby