Question about method 'azimuthal_integral_times_sin_beta'

[Jaebum]

Dear Amos

I had a chance to access your paper (Amos Egel et al 2017 J. Opt. 19 025605) and got interested. For few days I attempted to follow your demo of specular/diffuse reflection of nano-strucutred OLED thin film stacks. 

However, there seems to be some misunderstanding of mine about SMUTHI. As far as I'm concerned, the intensity of the reflected ray is at its maximum at specular reflection direction (=incident angle). But the value I_beta_j (azimuthal_integral_times_sin_beta) doesn't agree with the above and was at maximum elsewhere.

The incident beam's polar angle was 3/5 * pi, and the reflected beam's maximum value was detected at 34.5 degree (0.191666 pi). From what I've thought, the maximum should have been at 2/5 *pi but it wasn't. I think I have confused/missed something badly and I need some help.

 For detailed information, I include some snapshots.

(If it is okay for you to provide your code, that would be really thankful.)

Thanks

Jaebum Noh.

[Amos Egel]

Dear Jaebum,

thanks for your interest in Smuthi and in our work on OLED simulation.

For an incident angle of 3/5 Pi, the reflected beam should also be at 3/5 Pi. Please plot the 2D-distribution of totalfarfield to check that this is the case.

The method "azimuthal_integral_times_sin_beta()" returns the signal, integrated over the azimuthal direction and multiplied with the sine of the polar angle. The sin(beta) factor may shift the maximum towards larger angles. I don't understand, however, why it should shift the maximum to smaller angles. If you want me to look into that issue, I would need a fully running minimal script (minimal in the sense that it is a short script to reproduce your observation).

By the way, two years after the 2017 paper, I have also published a dissertation on the topic. It gives more details on OLED simulation with Smuthi.

Best regards, Amos

--
Smuthi project repository: https://gitlab.com/AmosEgel/smuthi
Online documentation https://smuthi.readthedocs.io/en/latest/
---
You received this message because you are subscribed to the Google Groups "Smuthi" group.
To unsubscribe from this group and stop receiving emails from it, send an email to smuthi+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/smuthi/f5923093-3091-4ece-a661-5bec3298db9bn%40googlegroups.com.

[Jaebum]

Dear Amos,

thank you for the fast and sincere reply.

I have set the value 'beam_waist' of the gaussian beam to 16.5nm, which should have been set as 16.5um. After changing the value and removing the sin(beta) term from I_jbeta, the maximum reflection was found at the right place! Everything else seems to work perfectly. All thanks to you.

And also thanks for the recommendation, I'll check the paper later.

Thanks

Jaebum Noh.

 

2022년 8월 23일 화요일 오전 12시 39분 46초 UTC+9에 amos...@gmail.com님이 작성: