Hi,

I’m trying to calculate the q resolved Eliashberg function on a grid of q-points in the BZ, which is a thing I’ve seen done in various papers with Wannier interpolation in QE+W90, but I was wondering if it is possible to do in Abinit with the new eph code.

I tried to follow the tutorial, from what I understand the q dependent Eliashberg function is stored in the A2F.nc file, but I’m having trouble extracting the data from this file. Usually I read NetCDF files with Matlab, but in this case I’m struggling with understanding the meaning of the quantities contained in the A2F.nc file. I would have expected to just find a vector containing the reduced coordinates of the q-points and another vector of the same size containing the values of alpha^2F in those points, but it doesn’t look like that to me. Can anybody give me some insight?

Thanks,

Luca

Hi Luca,

the file contains several different versions of a2F, but not q resolved. They are calculated on different q grids, with and without interpolation (qintp and qcoarse resp.)

Should be easy to add as an output option, but it’s not done afaik.

In the log/ouput files there are mode and q presentations of the lambda and line width values, for the original input q points. Might be a start for your analysis.

best

Matthieu

Hi Matthieu,

thanks for your reply. So, for example, if I take one version of a2F, it looks like the netcdf file contains these 3 arrays of data

```
a2ftr_mesh_qintp
Size: 473x1
Dimensions: a2ftr_nomega_qintp
Datatype: double
a2ftr_values_qintp
Size: 473x3x3x4x1
Dimensions: a2ftr_nomega_qintp,three,three,natom3p1,number_of_spins
Datatype: double
a2ftr_lambdaw_qintp
Size: 473x3x3x4x1
Dimensions: a2ftr_nomega_qintp,three,three,natom3p1,number_of_spins
Datatype: double
```

could you help me understand what they are? I wasn’t able to find any documentation on the Abinit website about this.

Luca