Static screening

From The Yambo Project
Revision as of 19:19, 16 April 2017 by Myrta (talk | contribs)
Jump to navigation Jump to search

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

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.

BSE1-Eq2.png


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

Links