Transport calculations

Formalism

The conductivity and the Seebeck coefficient in direction \(\alpha\beta\) are defined as [1]:

\[\sigma_{\alpha\beta} = \beta e^{2} A_{0,\alpha\beta} \ \ \ \text{and} \ \ \ S_{\alpha\beta} = -\frac{k_B}{|e|}\frac{A_{1,\alpha\beta}}{A_{0,\alpha\beta}},\]

in which the kinetic coefficients \(A_{n,\alpha\beta}\) are given by

\[A_{n,\alpha\beta} = N_{sp} \pi \hbar \int{d\omega \left(\beta\omega\right)^n f\left(\omega\right)f\left(-\omega\right)\Gamma_{\alpha\beta}\left(\omega,\omega\right)}.\]

Here \(N_{sp}\) is the spin factor and \(f(\omega)\) is the Fermi function. The transport distribution \(\Gamma_{\alpha\beta}\left(\omega_1,\omega_2\right)\) is defined as

\[\Gamma_{\alpha\beta}\left(\omega_1,\omega_2\right) = \frac{1}{V} \sum_k Tr\left(v_{k,\alpha}A_{k}(\omega_1)v_{k,\beta}A_{k}\left(\omega_2\right)\right),\]

where \(V\) is the unit cell volume. In multi-band systems the velocities \(v_{k}\) and the spectral function \(A(k,\omega)\) are matrices in the band indices \(i\) and \(j\). The frequency depended optical conductivity is given by

\[\sigma(\Omega) = N_{sp} \pi e^2 \hbar \int{d\omega \Gamma_{\alpha\beta}(\omega+\Omega/2,\omega-\Omega/2)\frac{f(\omega-\Omega/2)-f(\omega+\Omega/2)}{\Omega}}.\]

Prerequisites

First perform a standard DFT+DMFT calculation for your desired material and obtain the real-frequency self energy by doing an analytic continuation.

Warning

This package does NOT provide an explicit method to do an analytic continuation of self energies and Green functions from Matsubara frequencies to the real frequency axis! There are methods included e.g. in the ALPS package, which can be used for these purposes. Keep in mind that all these methods have to be used very carefully. Especially for optics calculations it is crucial to perform the analytic continuation in such a way that the obtained real frequency self energy is accurate around the Fermi energy as low energy features strongly influence the final results!

Besides the self energy the Wien2k files read by the transport converter (convert_transport_input) are:
  • .struct: The lattice constants specified in the struct file are used to calculate the unit cell volume.
  • .outputs: In this file the k-point symmetries are given.
  • .oubwin: Contains the indices of the bands within the projected subspace (written by dmftproj) for each k-point.
  • .pmat: This file is the output of the Wien2k optics package and contains the velocity (momentum) matrix elements between all bands in the desired energy window for each k-point. How to use the optics package is described below.
  • .h5: The hdf5 archive has to be present and should contain the dft_input subgroup. Otherwise convert_dft_input needs to be called before convert_transport_input.

Wien2k optics package

The basics steps to calculate the matrix elements of the momentum operator with the Wien2k optics package are:
  1. Perform a standard Wien2k calculation for your material.
  2. Run x kgen to generate a dense k-mesh.
  3. Run x lapw1.
  4. For metals change TETRA to 101.0 in case.in2.
  5. Run x lapw2 -fermi.
  6. Run x optic.

Additionally the input file case.inop is required. A detail description on how to setup this file can be found in the Wien2k user guide [2] on page 166. The optics energy window should be chosen according to the window used for dmftproj. Note that the current version of the transport code uses only the smaller of those two windows. However, keep in mind that the optics energy window has to be specified in absolute values and NOT relative to the Fermi energy! You can read off the Fermi energy from the case.scf2 file. Please do not set the optional parameter NBvalMAX in case.inop. Furthermore it is necessary to set line 6 to “ON” and put a “1” in the following line to enable the printing of the matrix elements to case.pmat.

Using the transport code

First we have to read the Wien2k files and store the relevant information in the hdf5 archive:

from pytriqs.applications.dft.converters.wien2k_converter import *
from pytriqs.applications.dft.sumk_dft_tools import *

Converter = Wien2kConverter(filename='case', repacking=True)
Converter.convert_transport_input()

SK = SumkDFTTools(hdf_file='case.h5', use_dft_blocks=True)

The converter convert_transport_input reads the required data of the Wien2k output and stores it in the dft_transp_input subgroup of your hdf file. Additionally we need to read and set the self energy, the chemical potential and the double counting:

ar = HDFArchive('case.h5', 'a')
SK.set_Sigma([ar['dmft_output']['Sigma_w']])
chemical_potential,dc_imp,dc_energ = SK.load(['chemical_potential','dc_imp','dc_energ'])
SK.set_mu(chemical_potential)
SK.set_dc(dc_imp,dc_energ)
del ar

As next step we can calculate the transport distribution \(\Gamma_{\alpha\beta}(\omega)\):

SK.transport_distribution(directions=['xx'], Om_mesh=[0.0, 0.1], energy_window=[-0.3,0.3],
                                             with_Sigma=True, broadening=0.0, beta=40)

Here the transport distribution is calculated in \(xx\) direction for the frequencies \(\Omega=0.0\) and \(0.1\). To use the previously obtained self energy we set with_Sigma to True and the broadening to \(0.0\). As we also want to calculate the Seebeck coefficient we have to include \(\Omega=0.0\) in the mesh. Note that the current version of the code repines the \(\Omega\) values to the closest values on the self energy mesh. For complete description of the input parameters see the transport_distribution reference.

The resulting transport distribution is not automatically saved, but this can be easily achieved with:

SK.save(['Gamma_w','Om_meshr','omega','directions'])

You can retrieve it from the archive by:

SK.Gamma_w, SK.Om_meshr, SK.omega, SK.directions = SK.load(['Gamma_w','Om_meshr','omega','directions'])

Finally the optical conductivity \(\sigma(\Omega)\) and the Seebeck coefficient \(S\) can be obtained with:

SK.conductivity_and_seebeck(beta=40)
SK.save(['seebeck','optic_cond'])

It is strongly advised to check convergence in the number of k-points!