The exchange term in BSE

Deals with issues related to computation of optical spectra in reciprocal space: RPA, TDDFT, local field effects.

Moderators: Davide Sangalli, andrea.ferretti, myrta gruning, andrea marini, Daniele Varsano, Conor Hogan

Post Reply
sdwang
Posts: 287
Joined: Fri Apr 09, 2010 12:30 pm

The exchange term in BSE

Post by sdwang » Fri Jun 30, 2023 3:18 am

Dear all,
When dealing with the system with SOC, how can we obtain the exchange interaction energy after BSE calculation?
Thanks!

Shudong
S. D. Wang
IMU,HOHHOT,CHINA
E-mail: sdwang@imu.edu.cn

User avatar
Davide Sangalli
Posts: 614
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: The exchange term in BSE

Post by Davide Sangalli » Fri Jun 30, 2023 9:41 am

Dear Shudong,
can you better specify what you would like to compute ?

Best,
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/

sdwang
Posts: 287
Joined: Fri Apr 09, 2010 12:30 pm

Re: The exchange term in BSE

Post by sdwang » Fri Jun 30, 2023 11:19 am

Dear Davide,
Thanks for your reply!
I want to obtain the magnitude of the exchange term in the BSE kernel with full spinorial calculations, like in TMDCs system with strong SOC.

Thanks!

Shudong
S. D. Wang
IMU,HOHHOT,CHINA
E-mail: sdwang@imu.edu.cn

User avatar
Davide Sangalli
Posts: 614
Joined: Tue May 29, 2012 4:49 pm
Location: Via Salaria Km 29.3, CP 10, 00016, Monterotondo Stazione, Italy
Contact:

Re: The exchange term in BSE

Post by Davide Sangalli » Fri Jun 30, 2023 11:28 am

Can you define "the magnitude of the exchange term" with an equation.

In general, when solving the BSE, the exchange kernel enters the excitonic hamiltonian with many matrix elements, K^x_{cvk,c'v'k'} ...
It is not just one number that you can quantify.

One option (but it is just an idea) could be to compare the BSE solution with and without K^x, and see how much the excitonic poles shift in energy.
You can do this with two calculations, one converging the energy cutoff on the exchange and the other setting it to 0.

With exchange (where E_converged can be 10000, 20000, etc ... depending on the material

Code: Select all

BSENGexx = "E_convergend" mHa 
No exchange

Code: Select all

BSENGexx = 0 mHa 
D.
Davide Sangalli, PhD
CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre
https://sites.google.com/view/davidesangalli
http://www.max-centre.eu/

Post Reply