Static screening

In this module you will learn how to calculate the static screening. This is the first step in the calculation of optical spectra within the Bethe-Salpeter framework.

Prerequisites

• You must first complete the "How to use Yambo" modules

You will need:

• The SAVE databases for 3D hBN
• The yambo executable
• gnuplot for plotting spectra

Background

To calculate the correlation part of the kernel W we need the static dielectric screening. This is a calculation of the linear response of the system analogous to the calculation of the RPA/IP dielectric function. One important difference is that here we consider the static dielectric function.

Choosing input parameters

Start by generating the input by invoking yambo with the option "-b" from the command line:

$ yambo -b -F 01_3D_BSE_screening.in

The input opens in the standard editor. Similarly to the other linear response calculations the relevant input variables to be changed are:

% BndsRnXs
  1 | 40 |                 
%
 NGsBlkXs= 4 Ry

The first variable gives how many bands are included in the sum to calculate the static response function. The second is a cutoff for the dimension of the static dielectric matrix.

In the next tutorial you will see how to choose these two parameters. Another relevant input parameter to change is

% LongDrXs
 1.000 | 1.000| 1.000
%

so that the perturbing electric field has component in each direction.

Static screening runlevel

Run the calculation by invoking yambo:

$ yambo -F 01_3D_BSE_screening.in -J 3D_BSE 

In the log (l_em1s) of the calculation you can see that after calculating the dipole matrix elements, for each q vector yambo calculates the IP response function and by inversion the RPA response function

<02s> Xo@q[1] |########################################| [100%] --(E) --(X)
<02s> X@q[1] |########################################| [100%] --(E) --(X)

In the report, r_em1s, the details of the calculations are reported under the 5th section

[05] Static Dielectric Matrix
=============================

This calculation does not produce any human readable output, but both the dipole matrix elements and the static screening dielectric function are saved in a database in the 3D_BSE directory:

3D_BSE/ndb.dip_iR_and_P
3D_BSE/ndb.em1s

which are needed in the BS kernel runlevel.

Summary

From this tutorial you've learned:

• How to compute the static screening as first step in a BSE calculation