I would like to compute the dipole transition matrix elements between valence and conduction states. I believe that with LOPTICS, these matrix elements will be calculated but I don't know how to access them. They are not written out in the OUTCAR with LOPTICS on. Can anyone tell me how to get these matrix elements?
Thank you in advance.
Dipole transition matrix elements from LOPTICS
Moderators: Global Moderator, Moderator
-
ssharifz
- Newbie
- Posts: 7
- Joined: Tue Mar 13, 2007 3:18 pm
- License Nr.: 284
Dipole transition matrix elements from LOPTICS
Last edited by ssharifz on Wed Jun 04, 2014 9:23 pm, edited 1 time in total.
-
ccao
- Newbie
- Posts: 2
- Joined: Wed Jun 04, 2014 5:35 am
- License Nr.: 5-1070
Dipole transition matrix elements from LOPTICS
If I understand correctly, you will have to set NPAR=1 and KPAR=1 with LOPTICS on. VASP will then output a binary file called "OPTICS". However, I don't know the detailed format & the unit used in that file.
Last edited by ccao on Thu Jun 05, 2014 12:01 am, edited 1 time in total.
-
5-1669
- Newbie
- Posts: 1
- Joined: Tue Jul 23, 2013 11:30 am
Re: Dipole transition matrix elements from LOPTICS
I had been looking for this and I think I found the answer in the code (v 5.3.5), so posting it here to verify and for others reference if correct.
The elements are written in binary format to the WAVEDER file (and can be written in text format by uncommenting a line in linear_optics.F that calls a function to write the formatted WAVEDERF file)
The elements of the derivatives with respect to k_i, i=x,y,z, <phi_nk|d/dk_i|phi_nk> are equal to - i <phi_nk|r_i|phi_nk>.
The elements are written in binary format to the WAVEDER file (and can be written in text format by uncommenting a line in linear_optics.F that calls a function to write the formatted WAVEDERF file)
The elements of the derivatives with respect to k_i, i=x,y,z, <phi_nk|d/dk_i|phi_nk> are equal to - i <phi_nk|r_i|phi_nk>.