Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Backscattering vs scattering #22

Open
nicodemoor opened this issue May 4, 2022 · 1 comment
Open

Backscattering vs scattering #22

nicodemoor opened this issue May 4, 2022 · 1 comment

Comments

@nicodemoor
Copy link

Hello,

I've found something weird when I compare the backscattering efficiency and the scattering efficiency of gold particle (50nm=R).

The former is higher than the last which I think could not be possible!

What's wrong?

comp

Here the code :

`
import PyMieScatt as ps
from Mie2 import Miedata
import matplotlib.pyplot as plt
from numpy import array
md = Miedata("Au")

wvl,n,k = md.indice()

m = n + k*1j

Qext_list = []
Qsca_list = []
Qabs_list = []
g_list = []
Qpr_list = []
Qback_list = []
Qratio_list = []

print(len(wvl))

for i in range(len(wvl)):
[Qext,Qsca,Qabs,g,Qpr,Qback,Qratio] = ps.AutoMieQ(m[i],wvl[i],50,1.33,asCrossSection=False)

Qext_list.append(Qext)
Qsca_list.append(Qsca)
Qabs_list.append(Qabs)
g_list.append(g)
Qpr_list.append(Qpr)
Qback_list.append(Qback)
Qratio_list.append(Qratio)

plt.plot(wvl,Qratio_list,label="Qratio")
plt.plot(wvl,array(Qback_list)/array(Qsca_list),"--",label="calculate")
plt.plot(wvl,Qback_list,label="Qback")
plt.plot(wvl,array(Qsca_list)array(Qratio_list),"--",label="QscaQratio")
plt.plot(wvl,array(Qsca_list),".",label="Qsca")
plt.ylabel("Q")
plt.xlabel("Wavelength (nm)")
plt.legend()
plt.show()

`

@isilber
Copy link

isilber commented Dec 26, 2024

@nicodemoor, unless I missed something, the code appears solid to me. This apparent inconsistency results from a missing Qback normalization (1 / 4 * pi) to render it in proper units of m-1 sr-1. I therefore suspect that this issue could be closed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants