The transition dipole moment calculation

I am trying to get the transition dipole moment from the ABINIT calculation. These values are used in the calculation of linear and nonlinear susceptibility, so I believe the transition dipole moment should be calculated in the program. Where can I find this data? If it is not printed out, is there any method to pick it up?

These values are used in the calculation of linear and nonlinear susceptibility, so I believe the transition dipole moment should be calculated in the program.

The matrix elements of the velocity operator

<nk| \dfrac{\partial{H}}{\partial{k_i}} |mk>

are computed in the DFPT code when we perform a DDK calculation.

See e.g. this section of the optics tutorial.
These matrix elements are stored in the 1WFK files produced by the DFPT run (one file for each reduced direction).

To understand how to read these matrix elements from file, you may want to look at
the optic Fortran tool (https://github.com/abinit/abinit/blob/9670f8e7ff0b3310581270c073eb60e28f2fc347/src/98_main/optic.F90#L287)

Note that optic is also able to read the same matrix elements from the EVK.nc netcdf files that are automatically produced by the DDK runs.
The EVK.files are much smaller that the 1WFK files as they do not contain the first-order wavefunctions, only the matrix elements of dH/dk

1 Like