def test_transform(self):
e_sus = [mp.LorentzianSusceptibility(sigma_diag=mp.Vector3(1, 2, 3),
sigma_offdiag=mp.Vector3(12, 13, 14)),
mp.DrudeSusceptibility(sigma_diag=mp.Vector3(1, 2, 3),
sigma_offdiag=mp.Vector3(12, 13, 14))]
mat = mp.Medium(epsilon_diag=mp.Vector3(1, 2, 3), epsilon_offdiag=mp.Vector3(12, 13, 14),
E_susceptibilities=e_sus)
rot_angle = math.radians(23.9)
rot_matrix = mp.Matrix(mp.Vector3(math.cos(rot_angle), math.sin(rot_angle), 0),
mp.Vector3(-math.sin(rot_angle), math.cos(rot_angle), 0),
mp.Vector3(0, 0, 1))
mat.transform(rot_matrix)
expected_diag = mp.Vector3(-7.72552, 10.72552, 3)
expected_offdiag = mp.Vector3(7.69024, 6.21332, 18.06640)
self.assertTrue(mat.epsilon_diag.close(expected_diag, tol=4))
self.assertTrue(mat.epsilon_offdiag.close(expected_offdiag, tol=4))
self.assertEqual(mat.mu_diag, mp.Vector3(1, 1, 1))
self.assertEqual(mat.mu_offdiag, mp.Vector3())
self.assertEqual(len(mat.E_susceptibilities), 2)
self.assertTrue(mat.E_susceptibilities[0].sigma_diag.close(expected_diag, tol=4))
self.assertTrue(mat.E_susceptibilities[0].sigma_offdiag.close(expected_offdiag, tol=4))
self.assertTrue(mat.E_susceptibilities[1].sigma_diag.close(expected_diag, tol=4))
self.assertTrue(mat.E_susceptibilities[1].sigma_offdiag.close(expected_offdiag, tol=4))
def calc(varibles,orginal,depth):
um_scale = 1
say = 0
eV_um_scale = um_scale/1.23984193
Au_plasma_frq = 9.03*eV_um_scale
metal_range = mp.FreqRange(min=um_scale/6.1992, max=um_scale/.24797)
out = []
var = [np.linspace(varibles[j][0],varibles[j][1],5) for j in range(6)]
correlation = -1
Au_f0 = 0.760
for i in range(5):
Au_frq0 = var[0][i]
Au_gam0 = 0.053*eV_um_scale
Au_sig0 = Au_f0*Au_plasma_frq**2/Au_frq0**2
Au_f1 = 0.024
for j in range(5):
Au_frq1 = var[1][j]*eV_um_scale # 2.988 um
Au_gam1 = 0.241*eV_um_scale
Au_sig1 = Au_f1*Au_plasma_frq**2/Au_frq1**2
Au_f2 = 0.010
for k in range(5):
Au_frq2 = var[2][k]*eV_um_scale # 1.494 um
Au_gam2 = 0.345*eV_um_scale
Au_sig2 = Au_f2*Au_plasma_frq**2/Au_frq2**2
Au_f3 = 0.071
for l in range(5):
Au_frq3 = var[3][l]*eV_um_scale # 0.418 um
Au_gam3 = 0.870*eV_um_scale
Au_sig3 = Au_f3*Au_plasma_frq**2/Au_frq3**2
Au_f4 = 0.601
for m in range(5):
Au_frq4 = var[4][m]*eV_um_scale # 0.288 um
Au_gam4 = 2.494*eV_um_scale
Au_sig4 = Au_f4*Au_plasma_frq**2/Au_frq4**2
Au_f5 = 4.384
for n in range(5):
Au_frq5 = var[5][n]*eV_um_scale # 0.093 um
Au_gam5 = 2.214*eV_um_scale
Au_sig5 = Au_f5*Au_plasma_frq**2/Au_frq5**2
Au_susc = [mp.DrudeSusceptibility(frequency=Au_frq0, gamma=Au_gam0, sigma=Au_sig0),
mp.LorentzianSusceptibility(frequency=Au_frq1, gamma=Au_gam1, sigma=Au_sig1),
mp.LorentzianSusceptibility(frequency=Au_frq2, gamma=Au_gam2, sigma=Au_sig2),
mp.LorentzianSusceptibility(frequency=Au_frq3, gamma=Au_gam3, sigma=Au_sig3),
mp.LorentzianSusceptibility(frequency=Au_frq4, gamma=Au_gam4, sigma=Au_sig4),
mp.LorentzianSusceptibility(frequency=Au_frq5, gamma=Au_gam5, sigma=Au_sig5)]
Au = mp.Medium(epsilon=1.0, E_susceptibilities=Au_susc, valid_freq_range=metal_range)
freqs = np.linspace(um_scale/0.7,um_scale/0.35,101)
au = np.empty(101,dtype=np.cdouble)
say=say+1
for ii in range (101):
au[ii]= Au.epsilon(freqs[ii])[1][1]
temp = np.corrcoef(au,orginal)[0][1]
if(correlation<temp):
correlation = temp
print(correlation)
out = [i,j,k,l,m,n]
print(out)
def test_1d_with_overlap(self):
# A z=30 cell split into three fragments of size 10 each, with a block
# covering the middle fragment, and half of the two outer fragments.
mat = mp.Medium(H_susceptibilities=[mp.DrudeSusceptibility()])
fs = self.get_fragment_stats(mp.Vector3(z=20),
mp.Vector3(z=30),
1,
def_mat=mat)
self.assertEqual(len(fs), 3)
# Middle fragment is completely covered by the block
self.check_stats(fs[1],
a_eps=100,
a_mu=100,
nonlin=300,
susc=300,
cond=300)
# Outer two fragments are half covered by the block, and half covered by default_material 'mat'
for i in [0, 2]:
self.check_stats(fs[i],
a_eps=50,
a_mu=50,
nonlin=150,
susc=200,
cond=150)
def test_no_geometry(self):
mat = mp.Medium(epsilon=12,
epsilon_offdiag=mp.Vector3(x=1),
mu_offdiag=mp.Vector3(x=20),
E_chi2_diag=mp.Vector3(1, 1),
H_chi3_diag=mp.Vector3(x=1),
E_susceptibilities=[
mp.LorentzianSusceptibility(),
mp.NoisyLorentzianSusceptibility()
],
H_susceptibilities=[mp.DrudeSusceptibility()],
D_conductivity_diag=mp.Vector3(y=1),
B_conductivity_diag=mp.Vector3(x=1, z=1))
fs = self.get_fragment_stats(mp.Vector3(),
mp.Vector3(10, 10),
2,
def_mat=mat,
geom=[])
self.assertEqual(len(fs), 1)
self.check_stats(fs[0],
a_eps=10000,
a_mu=10000,
nonlin=30000,
susc=30000,
cond=30000)
def get_fragment_stats(self, block_size, cell_size, dims, box_center=mp.Vector3(), dft_vecs=None,
def_mat=mp.air, sym=[]):
mat = mp.Medium(
epsilon=12,
epsilon_offdiag=mp.Vector3(z=1),
mu_offdiag=mp.Vector3(x=20),
E_chi2_diag=mp.Vector3(1, 1),
H_chi3_diag=mp.Vector3(z=1),
E_susceptibilities=[mp.LorentzianSusceptibility(), mp.NoisyLorentzianSusceptibility()],
H_susceptibilities=[mp.DrudeSusceptibility()],
D_conductivity_diag=mp.Vector3(y=1),
B_conductivity_diag=mp.Vector3(x=1, z=1)
)
geom = [mp.Block(size=block_size, center=box_center, material=mat)]
sim = mp.Simulation(cell_size=cell_size, resolution=10, geometry=geom, dimensions=dims,
default_material=def_mat, symmetries=sym)
if dft_vecs:
if dft_vecs['flux_regions']:
sim.add_flux(1, 0.5, 5, *dft_vecs['flux_regions'])
if dft_vecs['n2f_regions']:
sim.add_near2far(1, 0.5, 7, *dft_vecs['n2f_regions'])
if dft_vecs['force_regions']:
sim.add_force(1, 0.5, 9, *dft_vecs['force_regions'])
if dft_vecs['fields_components']:
sim.add_dft_fields(dft_vecs['fields_components'], 0, 1, 5, where=dft_vecs['fields_where'],
center=dft_vecs['fields_center'], size=dft_vecs['fields_size'])
gv = sim._create_grid_volume(False)
stats = sim._compute_fragment_stats(gv)
return stats
def test_3d_with_overlap(self):
# A 30 x 30 x 30 cell with a 20 x 20 x 20 block placed at the center, split
# into 27 10 x 10 x 10 fragments
mat = mp.Medium(E_susceptibilities=[mp.DrudeSusceptibility()])
fs = self.get_fragment_stats(mp.Vector3(20, 20, 20), mp.Vector3(30, 30, 30), 3, def_mat=mat)
self.assertEqual(len(fs), 27)
# Middle fragment contains entire block
idx = 13
self.check_stats(fs[idx], a_eps=10000, a_mu=10000, nonlin=30000, susc=30000, cond=30000)
# Six fragments adjacent to the middle fragment faces will be half covered by the block,
# and half covered by default_material 'mat'
for i in [4, 10, 12, 14, 16, 22]:
self.check_stats(fs[i], a_eps=5000, a_mu=5000, nonlin=15000, susc=20000, cond=15000)
# The corners will be 1/8 covered by the block and 7/8 covered by default_material 'mat'
for i in [0, 2, 6, 8, 18, 20, 24, 26]:
self.check_stats(fs[i], a_eps=1250, a_mu=1250, nonlin=3750, susc=12500, cond=3750)
# The rest will be 1/4 covered by the block and 3/4 covered by default_material 'mat'
for i in [1, 3, 5, 7, 9, 11, 15, 17, 19, 21, 23, 25]:
self.check_stats(fs[i], a_eps=2500, a_mu=2500, nonlin=7500, susc=15000, cond=7500)
def main(args):
resolution = 20
cell_size = mp.Vector3(z=10)
dimensions = 1
# conversion factor fro eV to 1/um
eV_um_scale = 1 / 1.23984193
# Al, from Rakic et al., Applied Optics, vol. 32, p. 5274 (1998)
Al_eps_inf = 1
Al_plasma_frq = 14.98 * eV_um_scale
Al_f0 = 0.523
Al_frq0 = 1e-10
Al_gam0 = 0.047 * eV_um_scale
Al_sig0 = (Al_f0 * math.sqrt(Al_plasma_frq)) / (math.sqrt(Al_frq0))
Al_f1 = 0.050
Al_frq1 = 1.544 * eV_um_scale # 803 nm
Al_gam1 = 0.312 * eV_um_scale
Al_sig1 = (Al_f1 * math.sqrt(Al_plasma_frq)) / (math.sqrt(Al_frq1))
E_susceptibilities = [
mp.DrudeSusceptibility(frequency=Al_frq0, gamma=Al_gam0,
sigma=Al_sig0),
mp.LorentzianSusceptibility(frequency=Al_frq1,
gamma=Al_gam1,
sigma=Al_sig1)
]
Al = mp.Medium(epsilon=Al_eps_inf, E_susceptibilities=E_susceptibilities)
pml_layers = [
mp.PML(1, direction=mp.Z) if args.pml else mp.Absorber(1,
direction=mp.Z)
]
sources = [
mp.Source(src=mp.GaussianSource(1 / 0.803, fwidth=0.1),
center=mp.Vector3(),
component=mp.Ex)
]
def print_stuff(sim):
print("ex:, {}, {}".format(sim.meep_time(),
sim.get_field_point(mp.Ex, mp.Vector3())))
sim = mp.Simulation(cell_size=cell_size,
resolution=resolution,
dimensions=dimensions,
default_material=Al,
boundary_layers=pml_layers,
sources=sources)
sim.run(mp.at_every(10, print_stuff),
until_after_sources=mp.stop_when_fields_decayed(
50, mp.Ex, mp.Vector3(), 1e-6))
def lossy_Pd(r):
"""Give a input r between 0 and 1, the function will return a modified
Pd with lower reflection. The absorption coefficient is lower
too. We have to use a thicker layer the prevent any transmission.
"""
um_scale = 1.0
# conversion factor for eV to 1/um [=1/hc]
eV_um_scale = um_scale / 1.23984193
metal_range = mp.FreqRange(min=um_scale / 12.4, max=um_scale / 0.2)
Pd_plasma_frq = 9.72 * eV_um_scale
Pd_f0 = 0.330
Pd_frq0 = 1e-10
Pd_gam0 = 0.008 * eV_um_scale
Pd_sig0 = Pd_f0 * Pd_plasma_frq**2 / Pd_frq0**2
Pd_f1 = 0.649
Pd_frq1 = 0.336 * eV_um_scale # 3.690 um
Pd_gam1 = 2.950 * eV_um_scale
Pd_sig1 = Pd_f1 * Pd_plasma_frq**2 / Pd_frq1**2
Pd_f2 = 0.121
Pd_frq2 = 0.501 * eV_um_scale # 2.475 um
Pd_gam2 = 0.555 * eV_um_scale
Pd_sig2 = Pd_f2 * Pd_plasma_frq**2 / Pd_frq2**2
Pd_f3 = 0.638
Pd_frq3 = 1.659 * eV_um_scale # 0.747 um
Pd_gam3 = 4.621 * eV_um_scale
Pd_sig3 = Pd_f3 * Pd_plasma_frq**2 / Pd_frq3**2
Pd_f4 = 0.453
Pd_frq4 = 5.715 * eV_um_scale # 0.217 um
Pd_gam4 = 3.236 * eV_um_scale
Pd_sig4 = Pd_f4 * Pd_plasma_frq**2 / Pd_frq4**2
Pd_susc = [
mp.DrudeSusceptibility(frequency=Pd_frq0,
gamma=Pd_gam0,
sigma=Pd_sig0 * r),
mp.LorentzianSusceptibility(frequency=Pd_frq1,
gamma=Pd_gam1,
sigma=Pd_sig1 * r),
mp.LorentzianSusceptibility(frequency=Pd_frq2,
gamma=Pd_gam2,
sigma=Pd_sig2 * r),
mp.LorentzianSusceptibility(frequency=Pd_frq3,
gamma=Pd_gam3,
sigma=Pd_sig3 * r),
mp.LorentzianSusceptibility(frequency=Pd_frq4,
gamma=Pd_gam4,
sigma=Pd_sig4 * r)
]
Pd = mp.Medium(epsilon=1.0,
E_susceptibilities=Pd_susc,
valid_freq_range=metal_range)
return Pd
def test_1d_with_overlap(self):
# A z=30 cell, with a block covering the middle 20 units.
mat = mp.Medium(H_susceptibilities=[mp.DrudeSusceptibility()])
fs = self.get_fragment_stats(mp.Vector3(z=20),
mp.Vector3(z=30),
1,
def_mat=mat)
self.check_stats(fs,
a_eps=200,
a_mu=200,
nonlin=600,
susc=700,
cond=600)
def setUp(self):
resolution = 20
cell_size = mp.Vector3(z=10)
dimensions = 1
# conversion factor from eV to 1/um
eV_um_scale = 1 / 1.23984193
# Al, from Rakic et al., Applied Optics, vol. 32, p. 5274 (1998)
Al_eps_inf = 1
Al_plasma_frq = 14.98 * eV_um_scale
Al_f0 = 0.523
Al_frq0 = 1e-10
Al_gam0 = 0.047 * eV_um_scale
Al_sig0 = (Al_f0 * math.sqrt(Al_plasma_frq)) / (math.sqrt(Al_frq0))
Al_f1 = 0.050
Al_frq1 = 1.544 * eV_um_scale # 803 nm
Al_gam1 = 0.312 * eV_um_scale
Al_sig1 = (Al_f1 * math.sqrt(Al_plasma_frq)) / (math.sqrt(Al_frq1))
E_susceptibilities = [
mp.DrudeSusceptibility(frequency=Al_frq0,
gamma=Al_gam0,
sigma=Al_sig0),
mp.LorentzianSusceptibility(frequency=Al_frq1,
gamma=Al_gam1,
sigma=Al_sig1)
]
Al = mp.Medium(epsilon=Al_eps_inf,
E_susceptibilities=E_susceptibilities)
absorber_layers = [mp.Absorber(1, direction=mp.Z)]
sources = [
mp.Source(src=mp.GaussianSource(1 / 0.803, fwidth=0.1),
center=mp.Vector3(),
component=mp.Ex)
]
self.sim = mp.Simulation(cell_size=cell_size,
resolution=resolution,
dimensions=dimensions,
default_material=Al,
boundary_layers=absorber_layers,
sources=sources)
def test_2d_with_overlap(self):
# A 30 x 30 cell, with a 20 x 20 block in the middle, split into 9 10 x 10 fragments.
mat = mp.Medium(H_susceptibilities=[mp.DrudeSusceptibility()])
fs = self.get_fragment_stats(mp.Vector3(20, 20), mp.Vector3(30, 30), 2, def_mat=mat)
self.assertEqual(len(fs), 9)
# Middle fragment contains entire block
idx = 4
self.check_stats(fs[idx], a_eps=1000, a_mu=1000, nonlin=3000, susc=3000, cond=3000)
# Top-middle, bottom-middle, left-middle, and right-middle fragments are half
# covered by the block, and half covered by default_material 'mat'.
for i in [1, 3, 5, 7]:
self.check_stats(fs[i], a_eps=500, a_mu=500, nonlin=1500, susc=2000, cond=1500)
# The four corner fragments are quarter-filled by the block, and 3/4 filled by
# default_material 'mat'
for i in [0, 2, 6, 8]:
self.check_stats(fs[i], a_eps=250, a_mu=250, nonlin=750, susc=1500, cond=750)
frq_max = 1/wvl_min
frq_cen = 0.5*(frq_min+frq_max)
dfrq = frq_max-frq_min
um_scale = alpha# default unit length is 1 um
eV_um_scale = um_scale/1.23984193# conversion factor for eV to 1/um [=1/hc]
h_plankeV=4.135667662*10**-15 # eVs
## Define material susceptibility
# Drude
drude_range = mp.FreqRange(min=um_scale/(wvl_max*2), max=um_scale/(wvl_min/2)) # Valid frequency range um_scale/wavelength
drude_frq = Omega_p*eV_um_scale # converts energy to meep frequency
drude_gam = Gama*eV_um_scale # converts energy to meep frequency
drude_sig = Sigma
drude_sig_diag=Sigma_diag
epsilon = epsilon_inf # epsilon infinity
drude_susc = [mp.DrudeSusceptibility(frequency=drude_frq, gamma=drude_gam, sigma=drude_sig, sigma_diag=drude_sig_diag)]
drude = mp.Medium(epsilon_diag=mp.Vector3(epsilon, epsilon, epsilon), E_susceptibilities=drude_susc, valid_freq_range=drude_range)
c1=mp.Vector3(tr_x, 0.0, 0.0)
c2=mp.Vector3(0.0, tr_y, 0.0)
c3=mp.Vector3(0.0, 0.0, tr_z)
m_tr=mp.Matrix(c1,c2,c3)
drude.transform(m_tr)
material_set=drude
## Create geometry
lx=lx*tr_x
ly=ly*tr_y
lz=lz*tr_z
material_ref=mp.Medium(epsilon_diag=mp.Vector3(1, 1, 1),mu_diag=mp.Vector3(1,1,1))
material_ref.transform(m_tr)
def import_medium(name, from_um_factor=1, paper="R"):
"""Returns Medium instance from string with specified length scale
It's widely based on Meep Materials Library, merely adapted to change
according to the length unit scale used.
Parameters
----------
name: str
Name of the desired material.
from_um_factor=1 : int, float, optional
Factor to transform from SI μm to the chosen length unit. For example,
to work with 100 nm as 1 Meep unit, from_um_factor=.1 must be specified.
paper="R": str
Name of desired source for experimental input data of medium.
Returns
-------
medium : mp.Medium
The mp.Medium instance of the desired material.
Raises
------
"No media called ... is available" : SyntaxError
If no material is found whose name matches the string given.
"""
if "Ag" not in name and "Au" not in name:
raise SyntaxError("No media called {} is available".format(name))
return
# Default unit length is 1 um
eV_from_um_factor = from_um_factor / 1.23984193 # Conversion factor: eV to 1/um [=1/hc]
# from_um_factor used to be um_scale (but that's a shady name)
############ GOLD #################################################
if name == "Au" and paper == "R":
#------------------------------------------------------------------
# Elemental metals from A.D. Rakic et al., Applied Optics, Vol. 37, No. 22, pp. 5271-83, 1998
# Wavelength range: 0.2 - 12.4 um
# Gold (Au)
metal_range = mp.FreqRange(min=from_um_factor * 1e3 / 6199.2,
max=from_um_factor * 1e3 / 247.97)
Au_plasma_frq = 9.03 * eV_from_um_factor
Au_f0 = 0.760
Au_frq0 = 1e-10
Au_gam0 = 0.053 * eV_from_um_factor
Au_sig0 = Au_f0 * Au_plasma_frq**2 / Au_frq0**2
Au_f1 = 0.024
Au_frq1 = 0.415 * eV_from_um_factor # 2.988 um
Au_gam1 = 0.241 * eV_from_um_factor
Au_sig1 = Au_f1 * Au_plasma_frq**2 / Au_frq1**2
Au_f2 = 0.010
Au_frq2 = 0.830 * eV_from_um_factor # 1.494 um
Au_gam2 = 0.345 * eV_from_um_factor
Au_sig2 = Au_f2 * Au_plasma_frq**2 / Au_frq2**2
Au_f3 = 0.071
Au_frq3 = 2.969 * eV_from_um_factor # 0.418 um
Au_gam3 = 0.870 * eV_from_um_factor
Au_sig3 = Au_f3 * Au_plasma_frq**2 / Au_frq3**2
Au_f4 = 0.601
Au_frq4 = 4.304 * eV_from_um_factor # 0.288 um
Au_gam4 = 2.494 * eV_from_um_factor
Au_sig4 = Au_f4 * Au_plasma_frq**2 / Au_frq4**2
Au_f5 = 4.384
Au_frq5 = 13.32 * eV_from_um_factor # 0.093 um
Au_gam5 = 2.214 * eV_from_um_factor
Au_sig5 = Au_f5 * Au_plasma_frq**2 / Au_frq5**2
Au_susc = [
mp.DrudeSusceptibility(frequency=Au_frq0,
gamma=Au_gam0,
sigma=Au_sig0),
mp.LorentzianSusceptibility(frequency=Au_frq1,
gamma=Au_gam1,
sigma=Au_sig1),
mp.LorentzianSusceptibility(frequency=Au_frq2,
gamma=Au_gam2,
sigma=Au_sig2),
mp.LorentzianSusceptibility(frequency=Au_frq3,
gamma=Au_gam3,
sigma=Au_sig3),
mp.LorentzianSusceptibility(frequency=Au_frq4,
gamma=Au_gam4,
sigma=Au_sig4),
mp.LorentzianSusceptibility(frequency=Au_frq5,
gamma=Au_gam5,
sigma=Au_sig5)
]
Au = mp.Medium(epsilon=1.0,
E_susceptibilities=Au_susc,
valid_freq_range=metal_range)
Au.from_um_factor = from_um_factor
return Au
elif name == "Au" and paper == "JC":
#------------------------------------------------------------------
# Metals from D. Barchiesi and T. Grosges, J. Nanophotonics, Vol. 8, 08996, 2015
# Wavelength range: 0.4 - 0.8 um
# Gold (Au)
# Fit to P.B. Johnson and R.W. Christy, Physical Review B, Vol. 6, pp. 4370-9, 1972
# metal_visible_range = mp.FreqRange(min=from_um_factor*1e3/800,
# max=from_um_factor*1e3/400)
# Au_JC_visible_frq0 = 1*from_um_factor/0.139779231751333
# Au_JC_visible_gam0 = 1*from_um_factor/26.1269913352870
# Au_JC_visible_sig0 = 1
# Au_JC_visible_frq1 = 1*from_um_factor/0.404064525036786
# Au_JC_visible_gam1 = 1*from_um_factor/1.12834046202759
# Au_JC_visible_sig1 = 2.07118534879440
# Au_JC_visible_susc = [mp.DrudeSusceptibility(frequency=Au_JC_visible_frq0, gamma=Au_JC_visible_gam0, sigma=Au_JC_visible_sig0),
# mp.LorentzianSusceptibility(frequency=Au_JC_visible_frq1, gamma=Au_JC_visible_gam1, sigma=Au_JC_visible_sig1)]
# Au_JC_visible = mp.Medium(epsilon=6.1599,
# E_susceptibilities=Au_JC_visible_susc,
# valid_freq_range=metal_visible_range)
# Au_JC_visible.from_um_factor = from_um_factor
# return Au_JC_visible
#------------------------------------------------------------------
# Metal from my own fit
# Wavelength range: 0.1879 - 1.937 um
# Gold (Au)
# Fit to P.B. Johnson and R.W. Christy, Physical Review B, Vol. 6, pp. 4370-9, 1972
Au_JC_range = mp.FreqRange(min=from_um_factor / 1.937,
max=from_um_factor / .1879)
freq_0 = 1.000000082740371e-10 * from_um_factor
gamma_0 = 4.4142682842363e-09 * from_um_factor
sigma_0 = 3.3002903009040977e+21
freq_1 = 0.3260039379724786 * from_um_factor
gamma_1 = 0.03601307014052124 * from_um_factor
sigma_1 = 103.74591029640469
freq_2 = 0.47387215165339414 * from_um_factor
gamma_2 = 0.34294093699162054 * from_um_factor
sigma_2 = 14.168079504545002
freq_3 = 2.3910144662345445 * from_um_factor
gamma_3 = 0.5900378254265015 * from_um_factor
sigma_3 = 0.8155991478264435
freq_4 = 3.3577076082530537 * from_um_factor
gamma_4 = 1.6689250252226686 * from_um_factor
sigma_4 = 2.038193481751111
freq_5 = 8.915719663759013 * from_um_factor
gamma_5 = 7.539763092679264 * from_um_factor
sigma_5 = 3.74409654935571
Au_JC_susc = [
mp.DrudeSusceptibility(frequency=freq_0,
gamma=gamma_0,
sigma=sigma_0),
mp.LorentzianSusceptibility(frequency=freq_1,
gamma=gamma_1,
sigma=sigma_1),
mp.LorentzianSusceptibility(frequency=freq_2,
gamma=gamma_2,
sigma=sigma_2),
mp.LorentzianSusceptibility(frequency=freq_3,
gamma=gamma_3,
sigma=sigma_3),
mp.LorentzianSusceptibility(frequency=freq_4,
gamma=gamma_4,
sigma=sigma_4),
mp.LorentzianSusceptibility(frequency=freq_5,
gamma=gamma_5,
sigma=sigma_5)
]
Au_JC = mp.Medium(
epsilon=1.0, #6.1599,
E_susceptibilities=Au_JC_susc,
valid_freq_range=Au_JC_range)
Au_JC.from_um_factor = from_um_factor
return Au_JC
elif name == "Au" and paper == "P":
#------------------------------------------------------------------
# Metals from D. Barchiesi and T. Grosges, J. Nanophotonics, Vol. 8, 08996, 2015
# Wavelength range: 0.4 - 0.8 um
# Gold (Au)
# Fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985
metal_visible_range = mp.FreqRange(min=from_um_factor * 1e3 / 800,
max=from_um_factor * 1e3 / 400)
Au_visible_frq0 = 1 * from_um_factor / 0.0473629248511456
Au_visible_gam0 = 1 * from_um_factor / 0.255476199605166
Au_visible_sig0 = 1
Au_visible_frq1 = 1 * from_um_factor / 0.800619321082804
Au_visible_gam1 = 1 * from_um_factor / 0.381870287531951
Au_visible_sig1 = -169.060953137985
Au_visible_susc = [
mp.DrudeSusceptibility(frequency=Au_visible_frq0,
gamma=Au_visible_gam0,
sigma=Au_visible_sig0),
mp.LorentzianSusceptibility(frequency=Au_visible_frq1,
gamma=Au_visible_gam1,
sigma=Au_visible_sig1)
]
Au_visible = mp.Medium(epsilon=0.6888,
E_susceptibilities=Au_visible_susc,
valid_freq_range=metal_visible_range)
Au_visible.from_um_factor = from_um_factor
return Au_visible
elif name == "Au":
raise ValueError("No source found for Au with that name")
############ SILVER ###############################################
if name == "Ag" and paper == "R":
#------------------------------------------------------------------
# Elemental metals from A.D. Rakic et al., Applied Optics, Vol. 37, No. 22, pp. 5271-83, 1998
# Wavelength range: 0.2 - 12.4 um
# Silver (Ag)
metal_range = mp.FreqRange(min=from_um_factor * 1e3 / 12398,
max=from_um_factor * 1e3 / 247.97)
Ag_plasma_frq = 9.01 * eV_from_um_factor
Ag_f0 = 0.845
Ag_frq0 = 1e-10
Ag_gam0 = 0.048 * eV_from_um_factor
Ag_sig0 = Ag_f0 * Ag_plasma_frq**2 / Ag_frq0**2
Ag_f1 = 0.065
Ag_frq1 = 0.816 * eV_from_um_factor # 1.519 um
Ag_gam1 = 3.886 * eV_from_um_factor
Ag_sig1 = Ag_f1 * Ag_plasma_frq**2 / Ag_frq1**2
Ag_f2 = 0.124
Ag_frq2 = 4.481 * eV_from_um_factor # 0.273 um
Ag_gam2 = 0.452 * eV_from_um_factor
Ag_sig2 = Ag_f2 * Ag_plasma_frq**2 / Ag_frq2**2
Ag_f3 = 0.011
Ag_frq3 = 8.185 * eV_from_um_factor # 0.152 um
Ag_gam3 = 0.065 * eV_from_um_factor
Ag_sig3 = Ag_f3 * Ag_plasma_frq**2 / Ag_frq3**2
Ag_f4 = 0.840
Ag_frq4 = 9.083 * eV_from_um_factor # 0.137 um
Ag_gam4 = 0.916 * eV_from_um_factor
Ag_sig4 = Ag_f4 * Ag_plasma_frq**2 / Ag_frq4**2
Ag_f5 = 5.646
Ag_frq5 = 20.29 * eV_from_um_factor # 0.061 um
Ag_gam5 = 2.419 * eV_from_um_factor
Ag_sig5 = Ag_f5 * Ag_plasma_frq**2 / Ag_frq5**2
Ag_susc = [
mp.DrudeSusceptibility(frequency=Ag_frq0,
gamma=Ag_gam0,
sigma=Ag_sig0),
mp.LorentzianSusceptibility(frequency=Ag_frq1,
gamma=Ag_gam1,
sigma=Ag_sig1),
mp.LorentzianSusceptibility(frequency=Ag_frq2,
gamma=Ag_gam2,
sigma=Ag_sig2),
mp.LorentzianSusceptibility(frequency=Ag_frq3,
gamma=Ag_gam3,
sigma=Ag_sig3),
mp.LorentzianSusceptibility(frequency=Ag_frq4,
gamma=Ag_gam4,
sigma=Ag_sig4),
mp.LorentzianSusceptibility(frequency=Ag_frq5,
gamma=Ag_gam5,
sigma=Ag_sig5)
]
Ag = mp.Medium(epsilon=1.0,
E_susceptibilities=Ag_susc,
valid_freq_range=metal_range)
Ag.from_um_factor = from_um_factor
return Ag
elif name == "Ag" and paper == "P":
#------------------------------------------------------------------
# Metals from D. Barchiesi and T. Grosges, J. Nanophotonics, Vol. 8, 08996, 2015
# Wavelength range: 0.4 - 0.8 um
## WARNING: unstable; field divergence may occur
# Silver (Au)
# Fit to E.D. Palik, Handbook of Optical Constants, Academic Press, 1985
metal_visible_range = mp.FreqRange(min=from_um_factor * 1e3 / 800,
max=from_um_factor * 1e3 / 400)
Ag_visible_frq0 = 1 * from_um_factor / 0.142050162130618
Ag_visible_gam0 = 1 * from_um_factor / 18.0357292925015
Ag_visible_sig0 = 1
Ag_visible_frq1 = 1 * from_um_factor / 0.115692151792108
Ag_visible_gam1 = 1 * from_um_factor / 0.257794324096575
Ag_visible_sig1 = 3.74465275944019
Ag_visible_susc = [
mp.DrudeSusceptibility(frequency=Ag_visible_frq0,
gamma=Ag_visible_gam0,
sigma=Ag_visible_sig0),
mp.LorentzianSusceptibility(frequency=Ag_visible_frq1,
gamma=Ag_visible_gam1,
sigma=Ag_visible_sig1)
]
Ag_visible = mp.Medium(epsilon=0.0067526,
E_susceptibilities=Ag_visible_susc,
valid_freq_range=metal_visible_range)
Ag_visible.from_um_factor = from_um_factor
return Ag_visible
elif name == "Ag":
raise ValueError("No source found for Ag with that name")
else:
raise ValueError("No source found for that material")
# find the best performing set of parameters
idx_opt = np.where(np.min(mins) == mins)
print(f"optimal:, {ps[idx_opt]}, {mins[idx_opt]}")
# define a Meep `Medium` class object using the optimal fitting parameters
E_susceptibilities = []
for n in range(num_lorentzians):
mymaterial_freq = ps[idx_opt][0][3 * n + 1]
mymaterial_gamma = ps[idx_opt][0][3 * n + 2]
if mymaterial_freq == 0:
mymaterial_sigma = ps[idx_opt][0][3 * n + 0]
E_susceptibilities.append(
mp.DrudeSusceptibility(frequency=1.0,
gamma=mymaterial_gamma,
sigma=mymaterial_sigma))
else:
mymaterial_sigma = ps[idx_opt][0][3 * n + 0] / mymaterial_freq**2
E_susceptibilities.append(
mp.LorentzianSusceptibility(frequency=mymaterial_freq,
gamma=mymaterial_gamma,
sigma=mymaterial_sigma))
mymaterial = mp.Medium(epsilon=eps_inf, E_susceptibilities=E_susceptibilities)
# plot the fit and the actual data for comparison
mymaterial_eps = [mymaterial.epsilon(f)[0][0] for f in freqs_reduced]
plt.subplot(1, 2, 1)
Ag_gam1 = 3.886*eV_um_scale
Ag_sig1 = Ag_f1*Ag_plasma_frq**2/Ag_frq1**2
Ag_f2 = 0.124
Ag_frq2 = 4.481*eV_um_scale # 0.273 um
Ag_gam2 = 0.452*eV_um_scale
Ag_sig2 = Ag_f2*Ag_plasma_frq**2/Ag_frq2**2
Ag_f3 = 0.011
Ag_frq3 = 8.185*eV_um_scale # 0.152 um
Ag_gam3 = 0.065*eV_um_scale
Ag_sig3 = Ag_f3*Ag_plasma_frq**2/Ag_frq3**2
Ag_f4 = 0.840
Ag_frq4 = 9.083*eV_um_scale # 0.137 um
Ag_gam4 = 0.916*eV_um_scale
Ag_sig4 = Ag_f4*Ag_plasma_frq**2/Ag_frq4**2
Ag_susc = [ mp.DrudeSusceptibility(frequency=Ag_frq0, gamma=Ag_gam0, sigma=Ag_sig0),
mp.LorentzianSusceptibility(frequency=Ag_frq1, gamma=Ag_gam1, sigma=Ag_sig1),
mp.LorentzianSusceptibility(frequency=Ag_frq2, gamma=Ag_gam2, sigma=Ag_sig2),
mp.LorentzianSusceptibility(frequency=Ag_frq3, gamma=Ag_gam3, sigma=Ag_sig3),
mp.LorentzianSusceptibility(frequency=Ag_frq4, gamma=Ag_gam4, sigma=Ag_sig4) ]
Ag = mp.Medium(epsilon=1.0, E_susceptibilities=Ag_susc, valid_freq_range=metal_range)
#------------------------------------------------------------------
Au_plasma_frq = 9.03*eV_um_scale
Au_f0 = 0.760
Au_frq0 = 1e-10
Au_gam0 = 0.053*eV_um_scale
Au_sig0 = Au_f0*Au_plasma_frq**2/Au_frq0**2
Au_f1 = 0.024
Ag_sig1 = Ag_f1 * Ag_plasma_frq**2 / Ag_frq1**2
Ag_f2 = 0.124
Ag_frq2 = 4.481 * eV_um_scale # 0.273 um
Ag_gam2 = 0.452 * eV_um_scale
Ag_sig2 = Ag_f2 * Ag_plasma_frq**2 / Ag_frq2**2
Ag_f3 = 0.011
Ag_frq3 = 8.185 * eV_um_scale # 0.152 um
Ag_gam3 = 0.065 * eV_um_scale
Ag_sig3 = Ag_f3 * Ag_plasma_frq**2 / Ag_frq3**2
Ag_f4 = 0.840
Ag_frq4 = 9.083 * eV_um_scale # 0.137 um
Ag_gam4 = 0.916 * eV_um_scale
Ag_sig4 = Ag_f4 * Ag_plasma_frq**2 / Ag_frq4**2
Ag_susc = [
mp.DrudeSusceptibility(frequency=Ag_frq0, gamma=Ag_gam0, sigma=Ag_sig0),
mp.LorentzianSusceptibility(frequency=Ag_frq1,
gamma=Ag_gam1,
sigma=Ag_sig1),
mp.LorentzianSusceptibility(frequency=Ag_frq2,
gamma=Ag_gam2,
sigma=Ag_sig2),
mp.LorentzianSusceptibility(frequency=Ag_frq3,
gamma=Ag_gam3,
sigma=Ag_sig3),
mp.LorentzianSusceptibility(frequency=Ag_frq4,
gamma=Ag_gam4,
sigma=Ag_sig4)
]
Ag = mp.Medium(epsilon=1.0, E_susceptibilities=Ag_susc)
um_scale = 1.0
eV_um_scale = um_scale / 1.23984193
metal_range = mp.FreqRange(min=um_scale / 12.4, max=um_scale / 0.2)
############################
# Indium 295K
############################
metal_range = mp.FreqRange(min=um_scale / 10, max=um_scale / 0.1)
In295_plasma_frq = 7.462 * um_scale
In295_frq0 = 1e-10
In295_gam0 = 0.147 * um_scale
In295_sig0 = In295_plasma_frq**2 / In295_frq0**2
In295_susc = [
mp.DrudeSusceptibility(frequency=In295_frq0,
gamma=In295_gam0,
sigma=In295_sig0),
]
#:Personal material, Indium at 295K
In295 = mp.Medium(epsilon=3.363,
E_susceptibilities=In295_susc,
valid_freq_range=metal_range)
############################
# Indium 4.2K
############################
metal_range = mp.FreqRange(min=um_scale / 10, max=um_scale / 0.1)
In4p2_plasma_frq = 6.8353 * um_scale
In4p2_frq0 = 1e-10
In4p2_gam0 = 0.0471 * um_scale
# ---------------------------------------------------------------------------- #
# Import Libraries
# ---------------------------------------------------------------------------- #
import meep as mp
import numpy as np
import matplotlib.pyplot as plt
# ---------------------------------------------------------------------------- #
# Important constants and materials
# ---------------------------------------------------------------------------- #
eps0 = 7.40523989e+02
susceptibilities = [
mp.DrudeSusceptibility(frequency=5.28366673e-01,
gamma=1.62316031,
sigma=1.17034609e-03),
mp.LorentzianSusceptibility(frequency=4.19962269e+02,
gamma=1.84154007e+03,
sigma=7.44917949e+02)
]
FR4 = mp.Medium(epsilon=eps0, E_susceptibilities=susceptibilities)
# ---------------------------------------------------------------------------- #
# Create Geometry
# ---------------------------------------------------------------------------- #
ww = 0.5 # waveguide width (microns)
wl = 10 # waveguie height (microns)
waveguide = mp.Block(material=FR4,
center=mp.Vector3(0, 0, 0),
Ag_frq3 = 8.185 * eV_um_scale # 0.152 um
Ag_gam3 = 0.065 * eV_um_scale
Ag_sig3 = Ag_f3 * Ag_plasma_frq**2 / Ag_frq3**2
Ag_f4 = 0.840
Ag_frq4 = 9.083 * eV_um_scale # 0.137 um
Ag_gam4 = 0.916 * eV_um_scale
Ag_sig4 = Ag_f4 * Ag_plasma_frq**2 / Ag_frq4**2
Ag_f5 = 5.646
Ag_frq5 = 20.29 * eV_um_scale # 0.137 um
Ag_gam5 = 2.419 * eV_um_scale
Ag_sig5 = Ag_f5 * Ag_plasma_frq**2 / Ag_frq5**2
Ag_susc = [
mp.DrudeSusceptibility(frequency=Ag_frq0, gamma=Ag_gam0, sigma=Ag_sig0),
mp.LorentzianSusceptibility(frequency=Ag_frq1,
gamma=Ag_gam1,
sigma=Ag_sig1),
mp.LorentzianSusceptibility(frequency=Ag_frq2,
gamma=Ag_gam2,
sigma=Ag_sig2),
mp.LorentzianSusceptibility(frequency=Ag_frq3,
gamma=Ag_gam3,
sigma=Ag_sig3),
mp.LorentzianSusceptibility(frequency=Ag_frq4,
gamma=Ag_gam4,
sigma=Ag_sig4),
mp.LorentzianSusceptibility(frequency=Ag_frq5,
gamma=Ag_gam5,
sigma=Ag_sig5)
Au_gam2 = 0.345*eV_um_scale
Au_sig2 = Au_f2*Au_plasma_frq**2/Au_frq2**2
Au_f3 = 0.071
Au_frq3 = var[0]*eV_um_scale # 0.418 um
Au_gam3 = 0.870*eV_um_scale
Au_sig3 = Au_f3*Au_plasma_frq**2/Au_frq3**2
Au_f4 = 0.601
Au_frq4 = var[0]*eV_um_scale # 0.288 um
Au_gam4 = 2.494*eV_um_scale
Au_sig4 = Au_f4*Au_plasma_frq**2/Au_frq4**2
Au_f5 = 4.384
Au_frq5 = var[4]*eV_um_scale # 0.093 um
Au_gam5 = 2.214*eV_um_scale
Au_sig5 = Au_f5*Au_plasma_frq**2/Au_frq5**2
Au_susc = [mp.DrudeSusceptibility(frequency=Au_frq0, gamma=Au_gam0, sigma=Au_sig0),
mp.LorentzianSusceptibility(frequency=Au_frq1, gamma=Au_gam1, sigma=Au_sig1),
mp.LorentzianSusceptibility(frequency=Au_frq2, gamma=Au_gam2, sigma=Au_sig2),
mp.LorentzianSusceptibility(frequency=Au_frq3, gamma=Au_gam3, sigma=Au_sig3),
mp.LorentzianSusceptibility(frequency=Au_frq4, gamma=Au_gam4, sigma=Au_sig4),
mp.LorentzianSusceptibility(frequency=Au_frq5, gamma=Au_gam5, sigma=Au_sig5)]
Auour = mp.Medium(epsilon=1.0, E_susceptibilities=Au_susc, valid_freq_range=metal_range)
freqs = np.linspace(um_scale/0.7,um_scale/0.35,101)
au = np.empty(101,dtype=np.cdouble)
freqs = np.linspace(um_scale/0.7,um_scale/0.35,101)
au = np.empty(101,dtype=np.cdouble)
auour = np.empty(101,dtype=np.cdouble)
for i in range (101):
auour[i]= Auour.epsilon(freqs[i])[1][1]
au[i]= Au.epsilon(freqs[i])[1][1]
def main(args):
resolution = 30
nSiO2 = 1.4
SiO2 = mp.Medium(index=nSiO2)
# conversion factor for eV to 1/um
eV_um_scale = 1 / 1.23984193
# W, from Rakic et al., Applied Optics, vol. 32, p. 5274 (1998)
W_eps_inf = 1
W_plasma_frq = 13.22 * eV_um_scale
W_f0 = 0.206
W_frq0 = 1e-10
W_gam0 = 0.064 * eV_um_scale
W_sig0 = W_f0 * W_plasma_frq**2 / W_frq0**2
W_f1 = 0.054
W_frq1 = 1.004 * eV_um_scale # 1.235 um
W_gam1 = 0.530 * eV_um_scale
W_sig1 = W_f1 * W_plasma_frq**2 / W_frq1**2
W_f2 = 0.166
W_frq2 = 1.917 * eV_um_scale # 0.647
W_gam2 = 1.281 * eV_um_scale
W_sig2 = W_f2 * W_plasma_frq**2 / W_frq2**2
W_susc = [
mp.DrudeSusceptibility(frequency=W_frq0, gamma=W_gam0, sigma=W_sig0),
mp.LorentzianSusceptibility(frequency=W_frq1,
gamma=W_gam1,
sigma=W_sig1),
mp.LorentzianSusceptibility(frequency=W_frq2,
gamma=W_gam2,
sigma=W_sig2)
]
W = mp.Medium(epsilon=W_eps_inf, E_susceptibilities=W_susc)
# crystalline Si, from M.A. Green, Solar Energy Materials and Solar Cells, vol. 92, pp. 1305-1310 (2008)
# fitted Lorentzian parameters, only for 600nm-1100nm
Si_eps_inf = 9.14
Si_eps_imag = -0.0334
Si_eps_imag_frq = 1 / 1.55
Si_frq = 2.2384
Si_gam = 4.3645e-02
Si_sig = 14.797 / Si_frq**2
Si = mp.Medium(epsilon=Si_eps_inf,
D_conductivity=2 * math.pi * Si_eps_imag_frq * Si_eps_imag /
Si_eps_inf,
E_susceptibilities=[
mp.LorentzianSusceptibility(frequency=Si_frq,
gamma=Si_gam,
sigma=Si_sig)
])
a = args.a # lattice periodicity
cone_r = args.cone_r # cone radius
cone_h = args.cone_h # cone height
wire_w = args.wire_w # metal-grid wire width
wire_h = args.wire_h # metal-grid wire height
trench_w = args.trench_w # trench width
trench_h = args.trench_h # trench height
Np = args.Np # number of periods in supercell
dair = 1.0 # air gap thickness
dmcl = 1.7 # micro lens thickness
dsub = 3000.0 # substrate thickness
dpml = 1.0 # PML thickness
sxy = Np * a
sz = dpml + dair + dmcl + dsub + dpml
cell_size = mp.Vector3(sxy, sxy, sz)
boundary_layers = [
mp.PML(dpml, direction=mp.Z, side=mp.High),
mp.Absorber(dpml, direction=mp.Z, side=mp.Low)
]
geometry = []
if args.substrate:
geometry = [
mp.Sphere(material=SiO2,
radius=dmcl,
center=mp.Vector3(0, 0, 0.5 * sz - dpml - dair - dmcl)),
mp.Block(material=Si,
size=mp.Vector3(mp.inf, mp.inf, dsub + dpml),
center=mp.Vector3(0, 0, -0.5 * sz + 0.5 * (dsub + dpml))),
mp.Block(
material=W,
size=mp.Vector3(mp.inf, wire_w, wire_h),
center=mp.Vector3(0, -0.5 * sxy + 0.5 * wire_w,
-0.5 * sz + dpml + dsub + 0.5 * wire_h)),
mp.Block(
material=W,
size=mp.Vector3(mp.inf, wire_w, wire_h),
center=mp.Vector3(0, +0.5 * sxy - 0.5 * wire_w,
-0.5 * sz + dpml + dsub + 0.5 * wire_h)),
mp.Block(
material=W,
size=mp.Vector3(wire_w, mp.inf, wire_h),
center=mp.Vector3(-0.5 * sxy + 0.5 * wire_w, 0,
-0.5 * sz + dpml + dsub + 0.5 * wire_h)),
mp.Block(material=W,
size=mp.Vector3(wire_w, mp.inf, wire_h),
center=mp.Vector3(+0.5 * sxy - 0.5 * wire_w, 0,
-0.5 * sz + dpml + dsub + 0.5 * wire_h))
]
if args.substrate and args.texture:
for nx in range(Np):
for ny in range(Np):
cx = -0.5 * sxy + (nx + 0.5) * a
cy = -0.5 * sxy + (ny + 0.5) * a
geometry.append(
mp.Cone(material=SiO2,
radius=0,
radius2=cone_r,
height=cone_h,
center=mp.Vector3(
cx, cy,
0.5 * sz - dpml - dair - dmcl - 0.5 * cone_h)))
if args.substrate:
geometry.append(
mp.Block(material=SiO2,
size=mp.Vector3(mp.inf, trench_w, trench_h),
center=mp.Vector3(
0, -0.5 * sxy + 0.5 * trench_w,
0.5 * sz - dpml - dair - dmcl - 0.5 * trench_h)))
geometry.append(
mp.Block(material=SiO2,
size=mp.Vector3(mp.inf, trench_w, trench_h),
center=mp.Vector3(
0, +0.5 * sxy - 0.5 * trench_w,
0.5 * sz - dpml - dair - dmcl - 0.5 * trench_h)))
geometry.append(
mp.Block(material=SiO2,
size=mp.Vector3(trench_w, mp.inf, trench_h),
center=mp.Vector3(
-0.5 * sxy + 0.5 * trench_w, 0,
0.5 * sz - dpml - dair - dmcl - 0.5 * trench_h)))
geometry.append(
mp.Block(material=SiO2,
size=mp.Vector3(trench_w, mp.inf, trench_h),
center=mp.Vector3(
+0.5 * sxy - 0.5 * trench_w, 0,
0.5 * sz - dpml - dair - dmcl - 0.5 * trench_h)))
k_point = mp.Vector3(0, 0, 0)
lambda_min = 0.7 # minimum source wavelength
lambda_max = 1.0 # maximum source wavelength
fmin = 1 / lambda_max
fmax = 1 / lambda_min
fcen = 0.5 * (fmin + fmax)
df = fmax - fmin
sources = [
mp.Source(mp.GaussianSource(fcen, fwidth=df),
component=mp.Ex,
center=mp.Vector3(0, 0, 0.5 * sz - dpml - 0.5 * dair),
size=mp.Vector3(sxy, sxy, 0))
]
sim = mp.Simulation(resolution=resolution,
cell_size=cell_size,
boundary_layers=boundary_layers,
geometry=geometry,
dimensions=3,
k_point=k_point,
sources=sources)
nfreq = 50
refl = sim.add_flux(
fcen, df, nfreq,
mp.FluxRegion(center=mp.Vector3(0, 0, 0.5 * sz - dpml),
size=mp.Vector3(sxy, sxy, 0)))
trans_grid = sim.add_flux(
fcen, df, nfreq,
mp.FluxRegion(center=mp.Vector3(0, 0,
-0.5 * sz + dpml + dsub + wire_h),
size=mp.Vector3(sxy, sxy, 0)))
trans_sub_top = sim.add_flux(
fcen, df, nfreq,
mp.FluxRegion(center=mp.Vector3(0, 0, -0.5 * sz + dpml + dsub),
size=mp.Vector3(sxy, sxy, 0)))
trans_sub_bot = sim.add_flux(
fcen, df, nfreq,
mp.FluxRegion(center=mp.Vector3(0, 0, -0.5 * sz + dpml),
size=mp.Vector3(sxy, sxy, 0)))
sim.run(mp.at_beginning(mp.output_epsilon), until=0)
if args.substrate:
sim.load_minus_flux('refl-flux', refl)
sim.run(until_after_sources=mp.stop_when_fields_decayed(
50, mp.Ex, mp.Vector3(0, 0, -0.5 * sz + dpml + 0.5 * dsub), 1e-9))
if not args.substrate:
sim.save_flux('refl-flux', refl)
sim.display_fluxes(refl, trans_grid, trans_sub_top, trans_sub_bot)
