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
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