lifetime effects in IR reflectivity  [SOLVED]

MULTIBINIT, aTDEP, ANADDB, optics, cut3d, mrg* postprocessors

Moderators: MMNSchmitt, gonze

Locked
mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

lifetime effects in IR reflectivity

Post by mverstra » Thu Feb 16, 2017 10:54 am

Dear Matthieu

I am one of user abinit code. I want to quality factor by abinit code. I know quality factor is equal inverse dielectric loss tangent (tan δ = ε2/ε1). ε2 , ε1 are real and imaginary part dielectric.
I should used equation 8 of attached 33.pdf for obtaining ε2 , ε1. this relation said that dielectric tensor along the Cartesian coordinates is obtained by j th mode oscillator strength along both a and b direction , eigenfrequency, and damping factor of the jth vibrational mode. I calculated j th mode oscillator strength along both a and b direction by abinit code . Now need to calculated damping frequency. Unfortunately, i do not know how i calculate it by abinit code. could you please explain for me about how can i obtain damping frequency for all modes? In other hand, how do i calculate the IR reflectivity with damping from the frequency dependence of the dielectric constant ( I can calculate IR reflectivity without damping from the frequency dependence of the dielectric constant) by abinit code . should i use iovmov=2 and vis, ...?
You do favor if you provid me by guiding. I would appreciate your consideration and I am anxiously looking forward to your reply.

Best Regards,
Narges
Matthieu Verstraete
University of Liege, Belgium

mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

Re: lifetime effects in IR reflectivity  [SOLVED]

Post by mverstra » Thu Feb 16, 2017 11:09 am

Dear Narges,

I am not sure I understand your question entirely, but here goes. Your quality factor has to do with the zero frequency limit of the ratio involving the real and imaginary parts of epsilon_0, correct? Typically the intrinsic IR reflectivity spectrum is actually quite good without the broadening, but obviously not perfect wrt experiment.


Your reference J. Am. Ceram. Soc., 96 [9] 2898–2905 (2013) fits experimental data with a parametrized curve for the Raman and IR intensities. Their parameters are the frequencies of phonons (those you have) and the lifetimes gamma_j (those you do not have). The lifetimes are due to two main sources of scattering for the phonons: electrons and other phonons. These do not come out of a simple phonon DFPT run, but you can fit them.

The former can be calculated using the elphon tutorial and anaddb, though this is quite heavy because you need a lot of k points. We are working on a much faster interpolated version for the epc, but this is under development, and we'll have to update the tutorial as well.

The latter can be calculated as well from MD or other finite displacements, using external software such as phono3py or ShengBTE. You can set the "prtposcar" flag in abinit to get VASP style outputs, and import the results for the anharmonic forces from abinit. You run phono3 or sheng, and get back the phonon lifetimes due to anharmonic phonon scattering. Again, for big systems this will be quite heavy as you have many perturbations to do for the 3rd interatomic force constants.

best
Matthieu Verstraete
University of Liege, Belgium

mverstra
Posts: 655
Joined: Wed Aug 19, 2009 12:01 pm

Re: lifetime effects in IR reflectivity

Post by mverstra » Thu Feb 16, 2017 1:23 pm

PS: Atsushi Togo has just confirmed he will also make a direct abinit interface for phono3py (should be soon).

See https://atztogo.github.io/phono3py/
https://atztogo.github.io/phonopy/
Matthieu Verstraete
University of Liege, Belgium

Locked