Page 1 of 1

Imaginary part of SHG calculation

Posted: Sat Nov 28, 2020 6:17 am
by pangrt
Dear developers,
I am running SHG with Yambo and got the SHG chi, and how could I obtain the imaginary part of the dielectric constant at both
ω/2 and ω? (As in attached, the red solid and dashed lines)



Thanks!

Rongtian

Re: Imaginary part of SHG calculation

Posted: Sat Nov 28, 2020 8:25 am
by lyzhao
Do you mean the imaginary part of the dielectric constant or SHG?
The dielectric constant is linear optical paramater and not SHG.

Re: Imaginary part of SHG calculation

Posted: Sat Nov 28, 2020 11:32 am
by pangrt
Dear Zhao,
I think it is the imaginary part of the dielectric. How to obtain the dielectric constant at ω/2?

Thanks!

Rongtian

Re: Imaginary part of SHG calculation

Posted: Sat Nov 28, 2020 11:41 am
by lyzhao
please see the first tutorial at http://www.attaccalite.com/lumen/tutorials.html.

Re: Imaginary part of SHG calculation

Posted: Sat Nov 28, 2020 12:14 pm
by pangrt
I notice that just an abs. of ω instead of ω/2?

Thanks!
Rongtian

Re: Imaginary part of SHG calculation

Posted: Thu Feb 17, 2022 9:27 am
by AmberLee123
I think the abs. of ω/2 is just to extend the abs. of ω in the picture, ie. f(x/2) compared with f(x).
So there is no specially abs. of ω/2 in the output file from yambo. We just read abs. ω and then recalculate it.

Best.

Re: Imaginary part of SHG calculation

Posted: Thu Feb 17, 2022 2:49 pm
by claudio
Dear Rongtian

once you have the imaginary part of the dielectric constant, just plot it with the argument divided by two
for example in python:


plt.plot(EPS[:,0]/2., EPS[:,1],color="red", linewidth=3.0, linestyle="-", label="$\epsilon_2(\omega/2)$")
plt.plot(EPS[:,0], EPS[:,1],color="red", linewidth=3.0, linestyle="--",label="$\epsilon_2(\omega)$")

best
Claudio