oldtests/c_Mie_tmatrices.py

   1 #!/usr/bin/env python
   2 # coding: utf-8
   3
   4 # In[1]:
   5
   6
   7 from qpms import *
   8
   9
  10 # In[6]:
  11
  12
  13 R = 40e-9
  14 ω_p = 9*eV/ℏ #9*eV/ℏ
  15 ε_inf = 4.6
  16 γ_p = 0.1*eV/ℏ
  17 ε_b = 2.13
  18 lMax = 3
  19
  20 ω = 1.5*eV/ℏ
  21
  22
  23 # In[7]:
  24
  25
  26 ε_m = ε_drude(ε_inf, ω_p, γ_p, ω)
  27
  28
  29 # In[9]:
  30
  31
  32 k_i = cmath.sqrt(ε_m)*ω/c
  33 k_e = cmath.sqrt(ε_b)*ω/c
  34 RH, RV, TH, TV = mie_coefficients(a=R, nmax=lMax, k_i=k_i, k_e=k_e, J_ext=1, J_scat=3)
  35
  36
  37 # In[11]:
  38
  39
  40 spec = BaseSpec(lMax=lMax)
  41 cT = CTMatrix.spherical(spec, R, k_i, k_e, 1, 1)
  42
  43
  44 # In[16]:
  45
  46
  47 print(np.diag(cT[...]))
  48
  49
  50 # In[18]:
  51
  52
  53 print(RV)
  54 print(RH)
  55
  56
  57 # In[ ]:
  58
  59
  60
  61