Bethe-Salpeter

UNDER CONSTRUCTION (MP)

In this tutorial you will learn how to:

• generate input files to run a BSE calculation for a two-dimensional material
• run the code and control inputs and outputs
• understand the content of the outputs

File:HBN-2D.png

Atomic structure of two-dimensional hBN

Material: two-dimensional hexagonal BN:

• exagonal lattice
• 2 atoms per cell, B and N (8 electrons, )
• Lattice constants: a = 4.716 [a.u.]
• Plane wave cutoff 40 Ry

You have 2 options now :

Option A: Start from the beginning (DFT runs also): In this case you will need:

• PWSCF input files and pseudopotentials for hBN-2D
• pw.x executable, version 5.0 or later
• p2y executable

Unpack the TARFILE:

$ tar -xcvf hBN-2D.tar
$ cd hBN-2D/PWSCF
$ ls
hbn-2D_nscf.in hbn-2D_scf.in References tmp

Run a DFT self-consistent (scf) and non-self-consistent (nscf) run:

$ pw.x < ./Inputs/hbn_2d_scf.in >  hbn_2d_scf.out
$ pw.x < ./Inputs/hbn_2d_scf.in >  hbn_2d_nscf.out

Generate the SAVE database

$ cd ./tmp/hbn_2d.save
$ p2y 

Move the SAVE database in the YAMBO directory

$ mv  SAVE ../../../YAMBO

Option B): Read directly a prepared SAVE database: $ tar -xcvf hBN-2D.tar $ cd hBN-2D/YAMBO $ mv ./References/SAVE SAVE

Unpack the TARFILE:

$ tar -xcvf hBN-2D.tar
$ cd hBN-2D/PWSCF
$ ls
hbn-2D_bands.in hbn-2D_nscf.in hbn-2D_scf.in	hbn-2D_scf_b.in REFERENCES

• Complete the Generating the Yambo databases tutorial for 2D-hBN
• SAVE folder for hBN-2D.
• yambo executable

Initialization

Every Yambo run must start with this step. Go to the folder containing the hBN-2D SAVE directory:

$ cd TUTORIALS/hBN-2D/YAMBO
$ ls
SAVE

and simply launch the code

yambo 

This will run the initialization (setup) runlevel.
TIP: do not run yambo from inside the SAVE folder! It will complain that "databases not found".

BSE Tutorials at present:

• SiH_4

• Fantastic dimension

• LiF

• you Hydrogen

• GaSb (SOC)

• hexagonal BN (on slides only)

The material

The h-BN is a layered material similar to graphite. It is an indirect wide band-gap material. Its optical spectrum is characterized by large excitonic effects. 5.95 eV is the minimum electronic indirect gap. 6.47 eV is the minimum direct gap. See the literature [1]

Bulk hexagonal BN [2] hexagonal lattice Two atoms per cell (16 electrons) Lattice parameters a = 4.72 [a.u.] c/a= 2.58 Plane waves cutoff 40 Rydberg

[[File:|Band Structure]]

Calculate the screening

A key ingredient to construct the BS kernel in the Screeened EXchange approximation (sex) is the screened electron-electron dielectric matrix W which is normally evaluated in the static approximation.

You have two main options to generate the input :

 1) yambo -b -F  yambo_Ws.in 

This will create the input to calculate W (database ndb.em1s)

 2) yambo -p p -F  yambo_Wp.in 

This will create the input to calculate (or read) W in the PPA approximation (database ndb.pp) This second option is useful if you want to read the static part of W from a previous database ndb.pp generated in a GW-PPA run

Calculate the BS kernel in the SEX (Screened Exchange approximation) in the transitions space

Here we learn how to create the BS kernel in the screened exchange (SEX) approximation which includes both exchange (V) and correlation (-W) terms Runlevel to be used yambo -o b -k sex

Solve the BSE by haydock solver and and calculate spectrum

Generate the input file for solving the BSE using the Haydock solver:

yambo  -o-y h

Solve the BSE by diagonalizing the excitonic matrix and calculate spectrum

Runlevel to be used yambo -y d

Steps:

-Calculate screening

-Calculate the BS Kernel

-Diagonalize the BS Matrix and calculate spectrum

-Visualize/Analyze excitons (with ypp)

-How to Converge

-How to work with SOC