As a final demonstration of mode decomposition, we compute the diffraction spectrum of a liquid-crystal polarization grating. These types of beam splitters use birefringence to produce diffraction orders which are circularly polarized. We will investigate two kinds of polarization gratings: (1) a homogeneous uniaxial grating (commonly known as a circular-polarization grating), and (2) a twisted-nematic bilayer grating as described in Optics Letters, Vol. 33, No. 20, pp. 2287-9, 2008 (pdf). The homogeneous uniaxial grating is just a special case of the twisted-nematic grating with a nematic director rotation angle of φ=0°.
A schematic of the grating geometry is shown below. The grating is a 2d slab in the xy-plane with two parameters: birefringence (Δn) and thickness (d). The twisted-nematic grating consists of two layers of thickness d each with equal and opposite rotation angles of φ=70° for the nematic director. Both gratings contain only three diffraction orders: m=0, ±1. The m=0 order is linearly polarized and the m=±1 orders are circularly polarized with opposite chirality. For the uniaxial grating, the diffraction efficiencies for a mode with wavelength λ can be computed analytically: η0=cos2(πΔnd/λ), η±1=0.5sin2(πΔnd/λ). The derivation of these formulas is presented in Optics Letters, Vol. 24, No. 9, pp. 584-6, 1999. We will verify these analytic results and also demonstrate that the twisted-nematic grating produces a broader bandwidth response for the ±1 orders than the homogeneous uniaxial grating. An important property of these polarization gratings for e.g. display applications is that for a circular-polarized input planewave and phase delay (Δnd/λ) of nearly 0.5, there is only a single diffraction order (+1 or -1) with opposite chiraity to that of the input. This is also demonstrated below.
In this example, the input is a linear-polarized planewave pulse at normal incidence with center wavelength of λ=0.54 μm. The linear polarization is in the yz-plane with a rotation angle of 45° counter clockwise around the x axis. Two sets of mode coefficients are computed in the air region adjacent to the grating for each orthogonal polarization: ODD_Z+EVEN_Y
and EVEN_Z+ODD_Y
, which correspond to +ky + -ky (cosine) and +ky - -ky (sine) modes. From these coefficients for linear-polarized modes, the power in the circular-polarized modes can be computed: |ODD_Z+EVEN_Y|2+|EVEN_Z+ODD_Y|2. The power is identical for the two circular-polarized modes with opposite chiralities since the input is linearly polarized and at normal incidence. The transmittance for the diffraction orders are computed from the mode coefficients. As usual, this requires a separate normalization run to compute the power of the input planewave.
The main part of the script is the function pol_grating
which computes the mode coefficients for a grating with thickness d
, twisted-nematic rotation angle ph
, and periodicity gp
. The anisotropic permittivity of the grating is specified using the material function lc_mat
which involves a position-dependent rotation of the diagonal ε tensor about the x-axis. For φ=0°, the nematic director is oriented along the z-axis: Ez has a larger permittivity than Ey where the birefringence (Δn) is 0.159. The grating has a periodicity of Λ=6.5 μm in the y direction.
import meep as mp
import math
import numpy as np
import matplotlib.pyplot as plt
resolution = 30 # pixels/μm
dpml = 1.0 # PML thickness
dsub = 1.0 # substrate thickness
dpad = 1.0 # padding thickness
k_point = mp.Vector3(0, 0, 0)
pml_layers = [mp.PML(thickness=dpml, direction=mp.X)]
n_0 = 1.55
delta_n = 0.159
epsilon_diag = mp.Matrix(
mp.Vector3(n_0**2, 0, 0),
mp.Vector3(0, n_0**2, 0),
mp.Vector3(0, 0, (n_0 + delta_n) ** 2),
)
wvl = 0.54 # center wavelength
fcen = 1 / wvl # center frequency
def pol_grating(d, ph, gp, nmode):
sx = dpml + dsub + d + d + dpad + dpml
sy = gp
cell_size = mp.Vector3(sx, sy, 0)
# twist angle of nematic director; from equation 1b
def phi(p):
xx = p.x - (-0.5 * sx + dpml + dsub)
if (xx >= 0) and (xx <= d):
return math.pi * p.y / gp + ph * xx / d
else:
return math.pi * p.y / gp - ph * xx / d + 2 * ph
# return the anisotropic permittivity tensor for a uniaxial, twisted nematic liquid crystal
def lc_mat(p):
# rotation matrix for rotation around x axis
Rx = mp.Matrix(
mp.Vector3(1, 0, 0),
mp.Vector3(0, math.cos(phi(p)), math.sin(phi(p))),
mp.Vector3(0, -math.sin(phi(p)), math.cos(phi(p))),
)
lc_epsilon = Rx * epsilon_diag * Rx.transpose()
lc_epsilon_diag = mp.Vector3(lc_epsilon[0].x, lc_epsilon[1].y, lc_epsilon[2].z)
lc_epsilon_offdiag = mp.Vector3(
lc_epsilon[1].x, lc_epsilon[2].x, lc_epsilon[2].y
)
return mp.Medium(
epsilon_diag=lc_epsilon_diag, epsilon_offdiag=lc_epsilon_offdiag
)
geometry = [
mp.Block(
center=mp.Vector3(-0.5 * sx + 0.5 * (dpml + dsub)),
size=mp.Vector3(dpml + dsub, mp.inf, mp.inf),
material=mp.Medium(index=n_0),
),
mp.Block(
center=mp.Vector3(-0.5 * sx + dpml + dsub + d),
size=mp.Vector3(2 * d, mp.inf, mp.inf),
material=lc_mat,
),
]
# linear-polarized planewave pulse source
src_pt = mp.Vector3(-0.5 * sx + dpml + 0.3 * dsub, 0, 0)
sources = [
mp.Source(
mp.GaussianSource(fcen, fwidth=0.05 * fcen),
component=mp.Ez,
center=src_pt,
size=mp.Vector3(0, sy, 0),
),
mp.Source(
mp.GaussianSource(fcen, fwidth=0.05 * fcen),
component=mp.Ey,
center=src_pt,
size=mp.Vector3(0, sy, 0),
),
]
sim = mp.Simulation(
resolution=resolution,
cell_size=cell_size,
boundary_layers=pml_layers,
k_point=k_point,
sources=sources,
default_material=mp.Medium(index=n_0),
)
tran_pt = mp.Vector3(0.5 * sx - dpml - 0.5 * dpad, 0, 0)
tran_flux = sim.add_flux(
fcen, 0, 1, mp.FluxRegion(center=tran_pt, size=mp.Vector3(0, sy, 0))
)
sim.run(until_after_sources=100)
input_flux = mp.get_fluxes(tran_flux)
sim.reset_meep()
sim = mp.Simulation(
resolution=resolution,
cell_size=cell_size,
boundary_layers=pml_layers,
k_point=k_point,
sources=sources,
geometry=geometry,
)
tran_flux = sim.add_flux(
fcen, 0, 1, mp.FluxRegion(center=tran_pt, size=mp.Vector3(0, sy, 0))
)
sim.run(until_after_sources=300)
res1 = sim.get_eigenmode_coefficients(
tran_flux, range(1, nmode + 1), eig_parity=mp.ODD_Z + mp.EVEN_Y
)
res2 = sim.get_eigenmode_coefficients(
tran_flux, range(1, nmode + 1), eig_parity=mp.EVEN_Z + mp.ODD_Y
)
angles = [math.degrees(math.acos(kdom.x / fcen)) for kdom in res1.kdom]
return input_flux[0], angles, res1.alpha[:, 0, 0], res2.alpha[:, 0, 0];
The properties of the two gratings are computed over the thickness range of 0.1 to 3.4 μm corresponding to phase delays (Δnd/λ) of approximately 0 to 1.
ph_uniaxial = 0 # chiral layer twist angle for uniaxial grating
ph_twisted = 70 # chiral layer twist angle for bilayer grating
gp = 6.5 # grating period
nmode = 5 # number of mode coefficients to compute
dd = np.arange(0.2, 3.5, 0.2) # chiral layer thickness
m0_uniaxial = np.zeros(dd.size)
m1_uniaxial = np.zeros(dd.size)
ang_uniaxial = np.zeros(dd.size)
m0_twisted = np.zeros(dd.size)
m1_twisted = np.zeros(dd.size)
ang_twisted = np.zeros(dd.size)
for k in range(len(dd)):
input_flux, angles, coeffs1, coeffs2 = pol_grating(
0.5 * dd[k], math.radians(ph_uniaxial), gp, nmode
)
tran = (abs(coeffs1) ** 2 + abs(coeffs2) ** 2) / input_flux
for m in range(nmode):
print("tran (uniaxial):, {}, {:.2f}, {:.5f}".format(m, angles[m], tran[m]))
m0_uniaxial[k] = tran[0]
m1_uniaxial[k] = tran[1]
ang_uniaxial[k] = angles[1]
input_flux, angles, coeffs1, coeffs2 = pol_grating(
dd[k], math.radians(ph_twisted), gp, nmode
)
tran = (abs(coeffs1) ** 2 + abs(coeffs2) ** 2) / input_flux
for m in range(nmode):
print("tran (twisted):, {}, {:.2f}, {:.5f}".format(m, angles[m], tran[m]))
m0_twisted[k] = tran[0]
m1_twisted[k] = tran[1]
ang_twisted[k] = angles[1]
----------- Initializing structure... time for choose_chunkdivision = 5.79357e-05 s Working in 2D dimensions. Computational cell is 4.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0339479 s ----------- run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 2.71797e-05 s Working in 2D dimensions. Computational cell is 4.2 x 6.5 x 0 with resolution 30 block, center = (-1.1,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-8.32667e-17,0,0) size (0.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 0.758735 s ----------- Meep progress: 224.88333333333333/408.0 = 55.1% done in 4.0s, 3.3s to go on time step 13493 (time=224.883), 0.000296471 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.89247 tran (uniaxial):, 1, 4.77, 0.01875 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 4.41074e-05 s Working in 2D dimensions. Computational cell is 4.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0652409 s ----------- Meep progress: 193.48333333333332/208.0 = 93.0% done in 4.0s, 0.3s to go on time step 11609 (time=193.483), 0.000344592 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.00815e-05 s Working in 2D dimensions. Computational cell is 4.4 x 6.5 x 0 with resolution 30 block, center = (-1.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-1.66533e-16,0,0) size (0.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 3.35973 s ----------- Meep progress: 100.36666666666666/408.0 = 24.6% done in 4.0s, 12.3s to go on time step 6022 (time=100.367), 0.000664325 s/step Meep progress: 200.35/408.0 = 49.1% done in 8.0s, 8.3s to go on time step 12022 (time=200.367), 0.000666721 s/step Meep progress: 301.4166666666667/408.0 = 73.9% done in 12.0s, 4.2s to go on time step 18088 (time=301.467), 0.000659487 s/step Meep progress: 400.9/408.0 = 98.3% done in 16.0s, 0.3s to go on time step 24058 (time=400.967), 0.000670094 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 29 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 55 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.84732 tran (twisted):, 1, 4.77, 0.03707 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 8.70228e-05 s Working in 2D dimensions. Computational cell is 4.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0973799 s ----------- Meep progress: 104.83333333333333/208.0 = 50.4% done in 4.0s, 3.9s to go on time step 6290 (time=104.833), 0.00063597 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 5.91278e-05 s Working in 2D dimensions. Computational cell is 4.4 x 6.5 x 0 with resolution 30 block, center = (-1.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-1.66533e-16,0,0) size (0.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 3.38153 s ----------- Meep progress: 100.73333333333333/408.0 = 24.7% done in 4.0s, 12.2s to go on time step 6044 (time=100.733), 0.00066186 s/step Meep progress: 201.35/408.0 = 49.4% done in 8.0s, 8.2s to go on time step 12083 (time=201.383), 0.000662443 s/step Meep progress: 301.6333333333333/408.0 = 73.9% done in 12.0s, 4.2s to go on time step 18102 (time=301.7), 0.000664572 s/step Meep progress: 403.3333333333333/408.0 = 98.9% done in 16.0s, 0.2s to go on time step 24205 (time=403.417), 0.000655435 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 66 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.79050 tran (uniaxial):, 1, 4.77, 0.06529 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 6.79493e-05 s Working in 2D dimensions. Computational cell is 4.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0874732 s ----------- Meep progress: 102.21666666666667/208.0 = 49.1% done in 4.0s, 4.1s to go on time step 6133 (time=102.217), 0.000652292 s/step Meep progress: 205.86666666666667/208.0 = 99.0% done in 8.0s, 0.1s to go on time step 12353 (time=205.883), 0.000643173 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.41346e-05 s Working in 2D dimensions. Computational cell is 4.8 x 6.5 x 0 with resolution 30 block, center = (-1.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (1.11022e-16,0,0) size (0.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 6.61127 s ----------- Meep progress: 97.56666666666666/408.0 = 23.9% done in 4.0s, 12.7s to go on time step 5854 (time=97.5667), 0.000683333 s/step Meep progress: 195.35/408.0 = 47.9% done in 8.0s, 8.7s to go on time step 11722 (time=195.367), 0.000681668 s/step Meep progress: 293.75/408.0 = 72.0% done in 12.0s, 4.7s to go on time step 17627 (time=293.783), 0.000677395 s/step Meep progress: 391.4/408.0 = 95.9% done in 16.0s, 0.7s to go on time step 23487 (time=391.45), 0.000682659 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.64693 tran (twisted):, 1, 4.77, 0.14670 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 7.70092e-05 s Working in 2D dimensions. Computational cell is 4.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0854559 s ----------- Meep progress: 102.86666666666666/208.0 = 49.5% done in 4.0s, 4.1s to go on time step 6172 (time=102.867), 0.000648176 s/step Meep progress: 207.0/208.0 = 99.5% done in 8.0s, 0.0s to go on time step 12421 (time=207.017), 0.000640153 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 4.6 x 6.5 x 0 with resolution 30 block, center = (-1.3,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (2.22045e-16,0,0) size (0.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 4.98919 s ----------- Meep progress: 98.01666666666667/408.0 = 24.0% done in 4.0s, 12.7s to go on time step 5881 (time=98.0167), 0.000680194 s/step Meep progress: 197.98333333333332/408.0 = 48.5% done in 8.0s, 8.5s to go on time step 11880 (time=198), 0.00066678 s/step Meep progress: 297.0833333333333/408.0 = 72.8% done in 12.0s, 4.5s to go on time step 17827 (time=297.117), 0.000672654 s/step Meep progress: 395.51666666666665/408.0 = 96.9% done in 16.0s, 0.5s to go on time step 23735 (time=395.583), 0.000677137 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.66119 tran (uniaxial):, 1, 4.77, 0.14065 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 2.69413e-05 s Working in 2D dimensions. Computational cell is 5.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.090482 s ----------- Meep progress: 97.23333333333333/208.0 = 46.7% done in 4.0s, 4.6s to go on time step 5834 (time=97.2333), 0.000685701 s/step Meep progress: 196.66666666666666/208.0 = 94.6% done in 8.0s, 0.5s to go on time step 11801 (time=196.683), 0.000670365 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.22272e-05 s Working in 2D dimensions. Computational cell is 5.2 x 6.5 x 0 with resolution 30 block, center = (-1.6,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (1.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 10.1716 s ----------- Meep progress: 89.91666666666667/408.0 = 22.0% done in 4.0s, 14.2s to go on time step 5395 (time=89.9167), 0.000741585 s/step Meep progress: 183.15/408.0 = 44.9% done in 8.0s, 9.8s to go on time step 10990 (time=183.167), 0.000715001 s/step Meep progress: 274.8333333333333/408.0 = 67.4% done in 12.0s, 5.8s to go on time step 16493 (time=274.883), 0.000726979 s/step Meep progress: 367.5333333333333/408.0 = 90.1% done in 16.0s, 1.8s to go on time step 22056 (time=367.6), 0.000719077 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 36 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 50 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.39057 tran (twisted):, 1, 4.77, 0.27182 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000102043 s Working in 2D dimensions. Computational cell is 4.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0766759 s ----------- Meep progress: 102.18333333333334/208.0 = 49.1% done in 4.0s, 4.1s to go on time step 6131 (time=102.183), 0.000652499 s/step Meep progress: 205.46666666666667/208.0 = 98.8% done in 8.0s, 0.1s to go on time step 12330 (time=205.5), 0.000645285 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.38962e-05 s Working in 2D dimensions. Computational cell is 4.8 x 6.5 x 0 with resolution 30 block, center = (-1.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (1.11022e-16,0,0) size (0.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 6.60763 s ----------- Meep progress: 93.03333333333333/408.0 = 22.8% done in 4.0s, 13.5s to go on time step 5582 (time=93.0333), 0.00071674 s/step Meep progress: 188.96666666666667/408.0 = 46.3% done in 8.0s, 9.3s to go on time step 11340 (time=189), 0.000694738 s/step Meep progress: 285.1666666666667/408.0 = 69.9% done in 12.0s, 5.2s to go on time step 17113 (time=285.217), 0.0006929 s/step Meep progress: 381.0/408.0 = 93.4% done in 16.0s, 1.1s to go on time step 22865 (time=381.083), 0.000695525 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 58 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 54 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.47140 tran (uniaxial):, 1, 4.77, 0.23506 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 8.98838e-05 s Working in 2D dimensions. Computational cell is 5.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.100205 s ----------- Meep progress: 92.38333333333333/208.0 = 44.4% done in 4.0s, 5.0s to go on time step 5543 (time=92.3833), 0.000721698 s/step Meep progress: 186.2/208.0 = 89.5% done in 8.0s, 0.9s to go on time step 11174 (time=186.233), 0.000710356 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.50883e-05 s Working in 2D dimensions. Computational cell is 5.6 x 6.5 x 0 with resolution 30 block, center = (-1.8,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (2.22045e-16,0,0) size (1.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 71.5192% done, 1.70329 s remaining subpixel-averaging is 71.5192% done, 1.68457 s remaining subpixel-averaging is 71.5192% done, 1.68515 s remaining time for set_epsilon = 13.1647 s ----------- Meep progress: 87.73333333333333/408.0 = 21.5% done in 4.0s, 14.6s to go on time step 5264 (time=87.7333), 0.00076 s/step Meep progress: 176.35/408.0 = 43.2% done in 8.0s, 10.5s to go on time step 10583 (time=176.383), 0.000752144 s/step Meep progress: 265.5833333333333/408.0 = 65.1% done in 12.0s, 6.4s to go on time step 15938 (time=265.633), 0.000747004 s/step Meep progress: 353.1166666666667/408.0 = 86.5% done in 16.0s, 2.5s to go on time step 21191 (time=353.183), 0.000761606 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 64 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.17500 tran (twisted):, 1, 4.77, 0.38032 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 7.10487e-05 s Working in 2D dimensions. Computational cell is 5 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0600371 s ----------- Meep progress: 100.0/208.0 = 48.1% done in 4.0s, 4.3s to go on time step 6000 (time=100), 0.000666774 s/step Meep progress: 201.83333333333334/208.0 = 97.0% done in 8.0s, 0.2s to go on time step 12113 (time=201.883), 0.000654451 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.10352e-05 s Working in 2D dimensions. Computational cell is 5 x 6.5 x 0 with resolution 30 block, center = (-1.5,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (1,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 8.54403 s ----------- Meep progress: 93.58333333333333/408.0 = 22.9% done in 4.0s, 13.4s to go on time step 5615 (time=93.5833), 0.000712491 s/step Meep progress: 189.16666666666666/408.0 = 46.4% done in 8.0s, 9.3s to go on time step 11352 (time=189.2), 0.000697323 s/step Meep progress: 283.1666666666667/408.0 = 69.4% done in 12.0s, 5.3s to go on time step 16993 (time=283.217), 0.000709178 s/step Meep progress: 378.1166666666667/408.0 = 92.7% done in 16.0s, 1.3s to go on time step 22691 (time=378.183), 0.000702061 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 35 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 61 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.27041 tran (uniaxial):, 1, 4.77, 0.32724 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 3.60012e-05 s Working in 2D dimensions. Computational cell is 6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.131748 s ----------- Meep progress: 87.66666666666667/208.0 = 42.1% done in 4.0s, 5.5s to go on time step 5260 (time=87.6667), 0.000760564 s/step Meep progress: 179.5/208.0 = 86.3% done in 8.0s, 1.3s to go on time step 10772 (time=179.533), 0.000725698 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.29425e-05 s Working in 2D dimensions. Computational cell is 6 x 6.5 x 0 with resolution 30 block, center = (-2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 68.5946% done, 1.84034 s remaining subpixel-averaging is 64.3072% done, 2.3664 s remaining subpixel-averaging is 64.3072% done, 2.43912 s remaining time for set_epsilon = 15.5353 s ----------- Meep progress: 83.8/408.0 = 20.5% done in 4.0s, 15.5s to go on time step 5028 (time=83.8), 0.000795612 s/step Meep progress: 168.46666666666667/408.0 = 41.3% done in 8.0s, 11.4s to go on time step 10110 (time=168.5), 0.000787236 s/step Meep progress: 253.21666666666667/408.0 = 62.1% done in 12.0s, 7.3s to go on time step 15197 (time=253.283), 0.000786469 s/step Meep progress: 337.46666666666664/408.0 = 82.7% done in 16.0s, 3.3s to go on time step 20253 (time=337.55), 0.00079122 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.04972 tran (twisted):, 1, 4.77, 0.44671 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000128031 s Working in 2D dimensions. Computational cell is 5.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.082103 s ----------- Meep progress: 98.0/208.0 = 47.1% done in 4.0s, 4.5s to go on time step 5880 (time=98), 0.000680371 s/step Meep progress: 197.83333333333334/208.0 = 95.1% done in 8.0s, 0.4s to go on time step 11871 (time=197.85), 0.000667708 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.10352e-05 s Working in 2D dimensions. Computational cell is 5.2 x 6.5 x 0 with resolution 30 block, center = (-1.6,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-1.11022e-16,0,0) size (1.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 9.882 s ----------- Meep progress: 92.15/408.0 = 22.6% done in 4.0s, 13.7s to go on time step 5529 (time=92.15), 0.000723491 s/step Meep progress: 185.61666666666667/408.0 = 45.5% done in 8.0s, 9.6s to go on time step 11139 (time=185.65), 0.00071314 s/step Meep progress: 278.5833333333333/408.0 = 68.3% done in 12.0s, 5.6s to go on time step 16719 (time=278.65), 0.000716963 s/step Meep progress: 371.01666666666665/408.0 = 90.9% done in 16.0s, 1.6s to go on time step 22267 (time=371.117), 0.000721017 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.14816 tran (uniaxial):, 1, 4.77, 0.39120 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000138044 s Working in 2D dimensions. Computational cell is 6.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.109352 s ----------- Meep progress: 85.88333333333333/208.0 = 41.3% done in 4.0s, 5.7s to go on time step 5153 (time=85.8833), 0.00077639 s/step Meep progress: 172.66666666666666/208.0 = 83.0% done in 8.0s, 1.6s to go on time step 10361 (time=172.683), 0.000768107 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 6.4 x 6.5 x 0 with resolution 30 block, center = (-2.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-2.22045e-16,0,0) size (2.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 58.4164% done, 3.10351 s remaining subpixel-averaging is 58.4164% done, 3.12547 s remaining subpixel-averaging is 58.4164% done, 3.0991 s remaining time for set_epsilon = 19.9794 s ----------- Meep progress: 80.01666666666667/408.0 = 19.6% done in 4.0s, 16.4s to go on time step 4801 (time=80.0167), 0.000833327 s/step Meep progress: 160.8/408.0 = 39.4% done in 8.0s, 12.3s to go on time step 9649 (time=160.817), 0.000825088 s/step Meep progress: 241.7/408.0 = 59.2% done in 12.0s, 8.3s to go on time step 14504 (time=241.733), 0.000823907 s/step Meep progress: 320.8333333333333/408.0 = 78.6% done in 16.0s, 4.3s to go on time step 19253 (time=320.883), 0.000842394 s/step Meep progress: 401.3/408.0 = 98.4% done in 20.0s, 0.3s to go on time step 24084 (time=401.4), 0.000828145 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 44 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 29 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.00536 tran (twisted):, 1, 4.77, 0.46488 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 6.10352e-05 s Working in 2D dimensions. Computational cell is 5.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0821059 s ----------- Meep progress: 93.78333333333333/208.0 = 45.1% done in 4.0s, 4.9s to go on time step 5627 (time=93.7833), 0.000710922 s/step Meep progress: 189.56666666666666/208.0 = 91.1% done in 8.0s, 0.8s to go on time step 11376 (time=189.6), 0.000695889 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.29425e-05 s Working in 2D dimensions. Computational cell is 5.4 x 6.5 x 0 with resolution 30 block, center = (-1.7,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-1.11022e-16,0,0) size (1.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) time for set_epsilon = 11.8445 s ----------- Meep progress: 89.43333333333334/408.0 = 21.9% done in 4.0s, 14.2s to go on time step 5366 (time=89.4333), 0.000745554 s/step Meep progress: 178.88333333333333/408.0 = 43.8% done in 8.0s, 10.2s to go on time step 10735 (time=178.917), 0.000745129 s/step Meep progress: 269.46666666666664/408.0 = 66.0% done in 12.0s, 6.2s to go on time step 16172 (time=269.533), 0.000735702 s/step Meep progress: 359.71666666666664/408.0 = 88.2% done in 16.0s, 2.1s to go on time step 21588 (time=359.8), 0.000738555 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 63 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.02627 tran (uniaxial):, 1, 4.77, 0.45435 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 7.79629e-05 s Working in 2D dimensions. Computational cell is 6.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.105287 s ----------- Meep progress: 82.33333333333333/208.0 = 39.6% done in 4.0s, 6.1s to go on time step 4940 (time=82.3333), 0.000809818 s/step Meep progress: 165.56666666666666/208.0 = 79.6% done in 8.0s, 2.1s to go on time step 9936 (time=165.6), 0.000800751 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.00951e-05 s Working in 2D dimensions. Computational cell is 6.8 x 6.5 x 0 with resolution 30 block, center = (-2.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-2.22045e-16,0,0) size (2.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 53.5143% done, 3.8263 s remaining subpixel-averaging is 53.5143% done, 3.81257 s remaining subpixel-averaging is 53.5143% done, 3.81274 s remaining time for set_epsilon = 23.647 s ----------- Meep progress: 74.46666666666667/408.0 = 18.3% done in 4.0s, 17.9s to go on time step 4468 (time=74.4667), 0.000895362 s/step Meep progress: 150.38333333333333/408.0 = 36.9% done in 8.0s, 13.7s to go on time step 9024 (time=150.4), 0.000878078 s/step Meep progress: 226.41666666666666/408.0 = 55.5% done in 12.0s, 9.6s to go on time step 13587 (time=226.45), 0.000876645 s/step Meep progress: 302.23333333333335/408.0 = 74.1% done in 16.0s, 5.6s to go on time step 18137 (time=302.283), 0.000879289 s/step Meep progress: 377.98333333333335/408.0 = 92.6% done in 20.0s, 1.6s to go on time step 22683 (time=378.05), 0.000879927 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.00085 tran (twisted):, 1, 4.77, 0.46955 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000133991 s Working in 2D dimensions. Computational cell is 5.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.096647 s ----------- Meep progress: 93.56666666666666/208.0 = 45.0% done in 4.0s, 4.9s to go on time step 5614 (time=93.5667), 0.000712629 s/step Meep progress: 188.45/208.0 = 90.6% done in 8.0s, 0.8s to go on time step 11308 (time=188.467), 0.000702542 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.29425e-05 s Working in 2D dimensions. Computational cell is 5.6 x 6.5 x 0 with resolution 30 block, center = (-1.8,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (2.22045e-16,0,0) size (1.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 71.5192% done, 1.72158 s remaining subpixel-averaging is 71.5192% done, 1.70206 s remaining subpixel-averaging is 71.5192% done, 1.70269 s remaining time for set_epsilon = 13.3156 s ----------- Meep progress: 87.53333333333333/408.0 = 21.5% done in 4.0s, 14.6s to go on time step 5252 (time=87.5333), 0.00076165 s/step Meep progress: 176.3/408.0 = 43.2% done in 8.0s, 10.5s to go on time step 10580 (time=176.333), 0.000750865 s/step Meep progress: 263.75/408.0 = 64.6% done in 12.0s, 6.6s to go on time step 15829 (time=263.817), 0.000762183 s/step Meep progress: 352.4/408.0 = 86.4% done in 16.0s, 2.5s to go on time step 21149 (time=352.483), 0.000751941 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 42 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 37 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.00005 tran (uniaxial):, 1, 4.77, 0.46429 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 7.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0993609 s ----------- Meep progress: 78.7/208.0 = 37.8% done in 4.0s, 6.6s to go on time step 4722 (time=78.7), 0.000847187 s/step Meep progress: 160.36666666666667/208.0 = 77.1% done in 8.0s, 2.4s to go on time step 9624 (time=160.4), 0.000816089 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.81877e-05 s Working in 2D dimensions. Computational cell is 7.2 x 6.5 x 0 with resolution 30 block, center = (-2.6,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (3.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 49.3713% done, 4.3752 s remaining subpixel-averaging is 78.9961% done, 1.1067 s remaining subpixel-averaging is 49.3713% done, 4.42276 s remaining subpixel-averaging is 78.9961% done, 1.10731 s remaining subpixel-averaging is 49.3713% done, 4.43155 s remaining subpixel-averaging is 78.9961% done, 1.11122 s remaining time for set_epsilon = 26.3941 s ----------- Meep progress: 72.93333333333334/408.0 = 17.9% done in 4.0s, 18.4s to go on time step 4376 (time=72.9333), 0.000914306 s/step Meep progress: 148.9/408.0 = 36.5% done in 8.0s, 13.9s to go on time step 8935 (time=148.917), 0.000877482 s/step Meep progress: 224.29999999999998/408.0 = 55.0% done in 12.0s, 9.8s to go on time step 13460 (time=224.333), 0.000884124 s/step Meep progress: 297.3333333333333/408.0 = 72.9% done in 16.0s, 6.0s to go on time step 17843 (time=297.383), 0.00091272 s/step Meep progress: 372.43333333333334/408.0 = 91.3% done in 20.0s, 1.9s to go on time step 22351 (time=372.517), 0.000887407 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 58 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 52 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 52 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.00061 tran (twisted):, 1, 4.77, 0.46431 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000139952 s Working in 2D dimensions. Computational cell is 5.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.103642 s ----------- Meep progress: 92.38333333333333/208.0 = 44.4% done in 4.0s, 5.0s to go on time step 5543 (time=92.3833), 0.000721725 s/step Meep progress: 186.01666666666665/208.0 = 89.4% done in 8.0s, 0.9s to go on time step 11163 (time=186.05), 0.000711761 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 5.8 x 6.5 x 0 with resolution 30 block, center = (-1.9,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (1.11022e-16,0,0) size (1.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 67.7217% done, 2.0545 s remaining subpixel-averaging is 67.7217% done, 2.04229 s remaining subpixel-averaging is 67.7217% done, 2.04423 s remaining time for set_epsilon = 14.95 s ----------- Meep progress: 86.31666666666666/408.0 = 21.2% done in 4.0s, 14.9s to go on time step 5179 (time=86.3167), 0.00077244 s/step Meep progress: 173.38333333333333/408.0 = 42.5% done in 8.0s, 10.8s to go on time step 10404 (time=173.4), 0.000765589 s/step Meep progress: 261.1/408.0 = 64.0% done in 12.0s, 6.8s to go on time step 15668 (time=261.133), 0.000759921 s/step Meep progress: 347.5333333333333/408.0 = 85.2% done in 16.0s, 2.8s to go on time step 20857 (time=347.617), 0.000770998 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 52 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 40 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.02864 tran (uniaxial):, 1, 4.77, 0.44626 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 5.00679e-05 s Working in 2D dimensions. Computational cell is 7.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.137346 s ----------- Meep progress: 78.08333333333333/208.0 = 37.5% done in 4.0s, 6.7s to go on time step 4685 (time=78.0833), 0.000853827 s/step Meep progress: 156.88333333333333/208.0 = 75.4% done in 8.0s, 2.6s to go on time step 9414 (time=156.9), 0.000845896 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 7.6 x 6.5 x 0 with resolution 30 block, center = (-2.8,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (2.22045e-16,0,0) size (3.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 45.8237% done, 5.1262 s remaining subpixel-averaging is 73.3197% done, 1.54776 s remaining subpixel-averaging is 45.8237% done, 5.0838 s remaining subpixel-averaging is 73.3197% done, 1.54977 s remaining subpixel-averaging is 45.8237% done, 5.10381 s remaining subpixel-averaging is 73.3197% done, 1.55772 s remaining time for set_epsilon = 30.0956 s ----------- Meep progress: 70.91666666666667/408.0 = 17.4% done in 4.0s, 19.0s to go on time step 4255 (time=70.9167), 0.000940257 s/step Meep progress: 142.71666666666667/408.0 = 35.0% done in 8.0s, 14.9s to go on time step 8564 (time=142.733), 0.000928472 s/step Meep progress: 213.79999999999998/408.0 = 52.4% done in 12.0s, 10.9s to go on time step 12830 (time=213.833), 0.000937849 s/step Meep progress: 284.56666666666666/408.0 = 69.7% done in 16.0s, 6.9s to go on time step 17078 (time=284.633), 0.000941825 s/step Meep progress: 355.2833333333333/408.0 = 87.1% done in 20.0s, 3.0s to go on time step 21322 (time=355.367), 0.000942543 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.00230 tran (twisted):, 1, 4.77, 0.45960 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 5.38826e-05 s Working in 2D dimensions. Computational cell is 6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.0905509 s ----------- Meep progress: 90.25/208.0 = 43.4% done in 4.0s, 5.2s to go on time step 5415 (time=90.25), 0.000738743 s/step Meep progress: 181.95/208.0 = 87.5% done in 8.0s, 1.1s to go on time step 10918 (time=181.967), 0.000726975 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.41346e-05 s Working in 2D dimensions. Computational cell is 6 x 6.5 x 0 with resolution 30 block, center = (-2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 64.3072% done, 2.38699 s remaining subpixel-averaging is 64.3072% done, 2.40457 s remaining subpixel-averaging is 64.3072% done, 2.40796 s remaining time for set_epsilon = 16.6695 s ----------- Meep progress: 82.75/408.0 = 20.3% done in 4.0s, 15.7s to go on time step 4965 (time=82.75), 0.000805785 s/step Meep progress: 175.51666666666665/408.0 = 43.0% done in 8.0s, 10.6s to go on time step 10533 (time=175.55), 0.000718519 s/step Meep progress: 259.15/408.0 = 63.5% done in 12.0s, 6.9s to go on time step 15553 (time=259.217), 0.000796936 s/step Meep progress: 342.56666666666666/408.0 = 84.0% done in 16.0s, 3.1s to go on time step 20559 (time=342.65), 0.00079916 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 68 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 55 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.15199 tran (uniaxial):, 1, 4.77, 0.39654 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 8.79765e-05 s Working in 2D dimensions. Computational cell is 8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.125251 s ----------- Meep progress: 75.01666666666667/208.0 = 36.1% done in 4.0s, 7.1s to go on time step 4501 (time=75.0167), 0.000888818 s/step Meep progress: 150.11666666666667/208.0 = 72.2% done in 8.0s, 3.1s to go on time step 9009 (time=150.15), 0.000887502 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.39098e-05 s Working in 2D dimensions. Computational cell is 8 x 6.5 x 0 with resolution 30 block, center = (-3,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 42.7517% done, 5.69995 s remaining subpixel-averaging is 68.4044% done, 1.92679 s remaining subpixel-averaging is 42.7517% done, 5.7258 s remaining subpixel-averaging is 68.4044% done, 1.91778 s remaining subpixel-averaging is 42.7517% done, 5.72558 s remaining subpixel-averaging is 68.4044% done, 1.9235 s remaining time for set_epsilon = 32.8454 s ----------- Meep progress: 67.6/408.0 = 16.6% done in 4.0s, 20.1s to go on time step 4056 (time=67.6), 0.000986321 s/step Meep progress: 136.23333333333332/408.0 = 33.4% done in 8.0s, 16.0s to go on time step 8175 (time=136.25), 0.00097132 s/step Meep progress: 204.93333333333334/408.0 = 50.2% done in 12.0s, 11.9s to go on time step 12298 (time=204.967), 0.000970339 s/step Meep progress: 273.5833333333333/408.0 = 67.1% done in 16.0s, 7.9s to go on time step 16418 (time=273.633), 0.000971036 s/step Meep progress: 341.8666666666667/408.0 = 83.8% done in 20.0s, 3.9s to go on time step 20515 (time=341.917), 0.000976364 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.05159 tran (twisted):, 1, 4.77, 0.43938 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 5.29289e-05 s Working in 2D dimensions. Computational cell is 6.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.104505 s ----------- Meep progress: 86.15/208.0 = 41.4% done in 4.0s, 5.7s to go on time step 5169 (time=86.15), 0.000773898 s/step Meep progress: 174.43333333333334/208.0 = 83.9% done in 8.0s, 1.5s to go on time step 10468 (time=174.467), 0.000754955 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.50883e-05 s Working in 2D dimensions. Computational cell is 6.2 x 6.5 x 0 with resolution 30 block, center = (-2.1,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (2.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 61.2204% done, 2.79341 s remaining subpixel-averaging is 61.2204% done, 2.64891 s remaining subpixel-averaging is 61.2204% done, 2.75618 s remaining time for set_epsilon = 18.3206 s ----------- Meep progress: 81.01666666666667/408.0 = 19.9% done in 4.0s, 16.1s to go on time step 4861 (time=81.0167), 0.000823039 s/step Meep progress: 163.05/408.0 = 40.0% done in 8.0s, 12.0s to go on time step 9784 (time=163.067), 0.00081255 s/step Meep progress: 246.23333333333332/408.0 = 60.4% done in 12.0s, 7.9s to go on time step 14776 (time=246.267), 0.000801294 s/step Meep progress: 327.5333333333333/408.0 = 80.3% done in 16.0s, 3.9s to go on time step 19655 (time=327.583), 0.000819844 s/step Meep progress: 407.8666666666667/408.0 = 100.0% done in 20.0s, 0.0s to go on time step 24477 (time=407.95), 0.000829674 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 51 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 49 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.30148 tran (uniaxial):, 1, 4.77, 0.31081 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 9.58443e-05 s Working in 2D dimensions. Computational cell is 8.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.128375 s ----------- Meep progress: 72.65/208.0 = 34.9% done in 4.0s, 7.5s to go on time step 4359 (time=72.65), 0.000917816 s/step Meep progress: 144.25/208.0 = 69.4% done in 8.0s, 3.5s to go on time step 8656 (time=144.267), 0.000930968 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 8.4 x 6.5 x 0 with resolution 30 block, center = (-3.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (4.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 40.0657% done, 6.53195 s remaining subpixel-averaging is 64.1067% done, 2.36242 s remaining subpixel-averaging is 40.0657% done, 6.48991 s remaining subpixel-averaging is 64.1067% done, 2.36408 s remaining subpixel-averaging is 40.0657% done, 6.52898 s remaining subpixel-averaging is 64.1067% done, 2.35834 s remaining time for set_epsilon = 36.7265 s ----------- Meep progress: 65.65/408.0 = 16.1% done in 4.0s, 20.9s to go on time step 3939 (time=65.65), 0.00101563 s/step Meep progress: 132.1/408.0 = 32.4% done in 8.0s, 16.7s to go on time step 7927 (time=132.117), 0.00100305 s/step Meep progress: 198.91666666666666/408.0 = 48.8% done in 12.0s, 12.6s to go on time step 11937 (time=198.95), 0.000997598 s/step Meep progress: 265.51666666666665/408.0 = 65.1% done in 16.0s, 8.6s to go on time step 15934 (time=265.567), 0.00100085 s/step Meep progress: 331.93333333333334/408.0 = 81.4% done in 20.0s, 4.6s to go on time step 19920 (time=332), 0.00100371 s/step Meep progress: 397.23333333333335/408.0 = 97.4% done in 24.0s, 0.7s to go on time step 23839 (time=397.317), 0.00102074 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.17455 tran (twisted):, 1, 4.77, 0.38170 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 4.50611e-05 s Working in 2D dimensions. Computational cell is 6.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.11933 s ----------- Meep progress: 86.23333333333333/208.0 = 41.5% done in 4.0s, 5.6s to go on time step 5174 (time=86.2333), 0.000773265 s/step Meep progress: 173.26666666666665/208.0 = 83.3% done in 8.0s, 1.6s to go on time step 10398 (time=173.3), 0.000765817 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.50883e-05 s Working in 2D dimensions. Computational cell is 6.4 x 6.5 x 0 with resolution 30 block, center = (-2.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (2.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 58.4164% done, 3.09306 s remaining subpixel-averaging is 58.4164% done, 2.89371 s remaining subpixel-averaging is 58.4164% done, 3.09871 s remaining time for set_epsilon = 19.6879 s ----------- Meep progress: 80.91666666666667/408.0 = 19.8% done in 4.0s, 16.2s to go on time step 4855 (time=80.9167), 0.000824082 s/step Meep progress: 162.7/408.0 = 39.9% done in 8.0s, 12.1s to go on time step 9763 (time=162.717), 0.000815005 s/step Meep progress: 244.65/408.0 = 60.0% done in 12.0s, 8.0s to go on time step 14681 (time=244.683), 0.000813452 s/step Meep progress: 326.3833333333333/408.0 = 80.0% done in 16.0s, 4.0s to go on time step 19586 (time=326.433), 0.00081554 s/step Meep progress: 407.9166666666667/408.0 = 100.0% done in 20.0s, 0.0s to go on time step 24479 (time=407.983), 0.000817639 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 69 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 49 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 37 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 36 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.46820 tran (uniaxial):, 1, 4.77, 0.22878 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000147104 s Working in 2D dimensions. Computational cell is 8.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.139945 s ----------- Meep progress: 69.91666666666667/208.0 = 33.6% done in 4.0s, 7.9s to go on time step 4195 (time=69.9167), 0.000953672 s/step Meep progress: 140.83333333333334/208.0 = 67.7% done in 8.0s, 3.8s to go on time step 8450 (time=140.833), 0.000940074 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.10487e-05 s Working in 2D dimensions. Computational cell is 8.8 x 6.5 x 0 with resolution 30 block, center = (-3.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (4.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 37.6973% done, 7.13631 s remaining subpixel-averaging is 60.3172% done, 2.75956 s remaining subpixel-averaging is 82.9371% done, 0.865693 s remaining subpixel-averaging is 37.6973% done, 7.11049 s remaining subpixel-averaging is 60.3172% done, 2.75036 s remaining subpixel-averaging is 82.9371% done, 0.858843 s remaining subpixel-averaging is 37.6973% done, 7.12251 s remaining subpixel-averaging is 60.3172% done, 2.75345 s remaining subpixel-averaging is 82.9371% done, 0.866829 s remaining time for set_epsilon = 39.6605 s ----------- Meep progress: 61.016666666666666/408.0 = 15.0% done in 4.0s, 22.7s to go on time step 3661 (time=61.0167), 0.00109264 s/step Meep progress: 122.33333333333333/408.0 = 30.0% done in 8.0s, 18.7s to go on time step 7342 (time=122.367), 0.00108692 s/step Meep progress: 185.16666666666666/408.0 = 45.4% done in 12.0s, 14.4s to go on time step 11113 (time=185.217), 0.00106086 s/step Meep progress: 248.03333333333333/408.0 = 60.8% done in 16.0s, 10.3s to go on time step 14886 (time=248.1), 0.00106035 s/step Meep progress: 310.8/408.0 = 76.2% done in 20.0s, 6.3s to go on time step 18653 (time=310.883), 0.00106216 s/step Meep progress: 373.48333333333335/408.0 = 91.5% done in 24.0s, 2.2s to go on time step 22416 (time=373.6), 0.0010633 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 27 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.40540 tran (twisted):, 1, 4.77, 0.26104 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 7.20024e-05 s Working in 2D dimensions. Computational cell is 6.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.098484 s ----------- Meep progress: 84.7/208.0 = 40.7% done in 4.0s, 5.8s to go on time step 5082 (time=84.7), 0.000787265 s/step Meep progress: 168.76666666666665/208.0 = 81.1% done in 8.0s, 1.9s to go on time step 10127 (time=168.783), 0.000792943 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 6.6 x 6.5 x 0 with resolution 30 block, center = (-2.3,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (2.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 55.858% done, 3.41444 s remaining subpixel-averaging is 55.858% done, 3.40405 s remaining subpixel-averaging is 55.858% done, 3.39686 s remaining time for set_epsilon = 21.5757 s ----------- Meep progress: 78.68333333333334/408.0 = 19.3% done in 4.0s, 16.7s to go on time step 4721 (time=78.6833), 0.000847361 s/step Meep progress: 157.51666666666665/408.0 = 38.6% done in 8.0s, 12.7s to go on time step 9452 (time=157.533), 0.000845538 s/step Meep progress: 237.29999999999998/408.0 = 58.2% done in 12.0s, 8.6s to go on time step 14240 (time=237.333), 0.000835482 s/step Meep progress: 315.93333333333334/408.0 = 77.4% done in 16.0s, 4.7s to go on time step 18959 (time=315.983), 0.000847783 s/step Meep progress: 394.68333333333334/408.0 = 96.7% done in 20.0s, 0.7s to go on time step 23686 (time=394.767), 0.000846355 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 53 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 38 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 70 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.66179 tran (uniaxial):, 1, 4.77, 0.14253 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000156164 s Working in 2D dimensions. Computational cell is 9.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.172037 s ----------- Meep progress: 68.06666666666666/208.0 = 32.7% done in 4.0s, 8.2s to go on time step 4084 (time=68.0667), 0.000979654 s/step Meep progress: 137.1/208.0 = 65.9% done in 8.0s, 4.1s to go on time step 8226 (time=137.1), 0.000965725 s/step Meep progress: 205.7/208.0 = 98.9% done in 12.0s, 0.1s to go on time step 12343 (time=205.717), 0.000971618 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.98566e-05 s Working in 2D dimensions. Computational cell is 9.2 x 6.5 x 0 with resolution 30 block, center = (-3.6,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (5.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 35.5933% done, 7.99426 s remaining subpixel-averaging is 56.9506% done, 3.19121 s remaining subpixel-averaging is 78.308% done, 1.1792 s remaining subpixel-averaging is 35.5933% done, 7.82132 s remaining subpixel-averaging is 56.9506% done, 3.19462 s remaining subpixel-averaging is 78.308% done, 1.17005 s remaining subpixel-averaging is 35.5933% done, 7.79931 s remaining subpixel-averaging is 56.9506% done, 3.1716 s remaining subpixel-averaging is 78.308% done, 1.18297 s remaining time for set_epsilon = 43.3607 s ----------- Meep progress: 61.21666666666667/408.0 = 15.0% done in 4.0s, 22.7s to go on time step 3673 (time=61.2167), 0.00108931 s/step Meep progress: 122.93333333333334/408.0 = 30.1% done in 8.0s, 18.6s to go on time step 7377 (time=122.95), 0.00108006 s/step Meep progress: 183.58333333333334/408.0 = 45.0% done in 12.0s, 14.7s to go on time step 11017 (time=183.617), 0.00109908 s/step Meep progress: 244.83333333333334/408.0 = 60.0% done in 16.0s, 10.7s to go on time step 14693 (time=244.883), 0.00108818 s/step Meep progress: 306.65/408.0 = 75.2% done in 20.0s, 6.6s to go on time step 18403 (time=306.717), 0.00107845 s/step Meep progress: 367.8/408.0 = 90.1% done in 24.0s, 2.6s to go on time step 22073 (time=367.883), 0.00109002 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.67195 tran (twisted):, 1, 4.77, 0.13686 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 9.08375e-05 s Working in 2D dimensions. Computational cell is 6.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.099133 s ----------- Meep progress: 83.68333333333334/208.0 = 40.2% done in 4.0s, 5.9s to go on time step 5021 (time=83.6833), 0.000796734 s/step Meep progress: 168.46666666666667/208.0 = 81.0% done in 8.0s, 1.9s to go on time step 10110 (time=168.5), 0.000786127 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 6.8 x 6.5 x 0 with resolution 30 block, center = (-2.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-2.22045e-16,0,0) size (2.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 53.5143% done, 3.79773 s remaining subpixel-averaging is 53.5143% done, 3.81162 s remaining subpixel-averaging is 53.5143% done, 3.8089 s remaining time for set_epsilon = 23.541 s ----------- Meep progress: 74.93333333333334/408.0 = 18.4% done in 4.0s, 17.8s to go on time step 4496 (time=74.9333), 0.000889767 s/step Meep progress: 151.51666666666665/408.0 = 37.1% done in 8.0s, 13.5s to go on time step 9092 (time=151.533), 0.000870396 s/step Meep progress: 229.31666666666666/408.0 = 56.2% done in 12.0s, 9.4s to go on time step 13761 (time=229.35), 0.000856824 s/step Meep progress: 307.1666666666667/408.0 = 75.3% done in 16.0s, 5.3s to go on time step 18434 (time=307.233), 0.00085611 s/step Meep progress: 383.71666666666664/408.0 = 94.0% done in 20.0s, 1.3s to go on time step 23028 (time=383.8), 0.000870738 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 69 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 34 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.80825 tran (uniaxial):, 1, 4.77, 0.06307 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000152111 s Working in 2D dimensions. Computational cell is 9.6 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.157447 s ----------- Meep progress: 64.9/208.0 = 31.2% done in 4.0s, 8.8s to go on time step 3894 (time=64.9), 0.00102741 s/step Meep progress: 130.88333333333333/208.0 = 62.9% done in 8.0s, 4.7s to go on time step 7854 (time=130.9), 0.00101027 s/step Meep progress: 196.55/208.0 = 94.5% done in 12.0s, 0.7s to go on time step 11795 (time=196.583), 0.00101506 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.20024e-05 s Working in 2D dimensions. Computational cell is 9.6 x 6.5 x 0 with resolution 30 block, center = (-3.8,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (-4.44089e-16,0,0) size (5.6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 33.7117% done, 8.53515 s remaining subpixel-averaging is 53.94% done, 3.58281 s remaining subpixel-averaging is 74.1684% done, 1.47227 s remaining subpixel-averaging is 33.7117% done, 8.64079 s remaining subpixel-averaging is 53.94% done, 3.58974 s remaining subpixel-averaging is 74.1684% done, 1.46986 s remaining subpixel-averaging is 33.7117% done, 8.56235 s remaining subpixel-averaging is 53.94% done, 3.61465 s remaining subpixel-averaging is 74.1684% done, 1.47245 s remaining time for set_epsilon = 46.6436 s ----------- Meep progress: 58.0/408.0 = 14.2% done in 4.0s, 24.1s to go on time step 3480 (time=58), 0.0011497 s/step Meep progress: 116.86666666666666/408.0 = 28.6% done in 8.0s, 19.9s to go on time step 7013 (time=116.883), 0.00113224 s/step Meep progress: 175.05/408.0 = 42.9% done in 12.0s, 16.0s to go on time step 10505 (time=175.083), 0.00114558 s/step Meep progress: 233.93333333333334/408.0 = 57.3% done in 16.0s, 11.9s to go on time step 14039 (time=233.983), 0.00113194 s/step Meep progress: 293.5/408.0 = 71.9% done in 20.0s, 7.8s to go on time step 17613 (time=293.55), 0.00111922 s/step Meep progress: 352.01666666666665/408.0 = 86.3% done in 24.0s, 3.8s to go on time step 21125 (time=352.083), 0.00113908 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 47 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 35 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 54 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 38 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 31 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.86580 tran (twisted):, 1, 4.77, 0.02804 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000117064 s Working in 2D dimensions. Computational cell is 7 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.103678 s ----------- Meep progress: 80.11666666666666/208.0 = 38.5% done in 4.0s, 6.4s to go on time step 4807 (time=80.1167), 0.000832301 s/step Meep progress: 161.7/208.0 = 77.7% done in 8.0s, 2.3s to go on time step 9704 (time=161.733), 0.000816961 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.39098e-05 s Working in 2D dimensions. Computational cell is 7 x 6.5 x 0 with resolution 30 block, center = (-2.5,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (2.22045e-16,0,0) size (3,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 51.3594% done, 4.1118 s remaining subpixel-averaging is 51.3594% done, 4.09325 s remaining subpixel-averaging is 51.3594% done, 4.09493 s remaining time for set_epsilon = 24.965 s ----------- Meep progress: 75.21666666666667/408.0 = 18.4% done in 4.0s, 17.7s to go on time step 4513 (time=75.2167), 0.000886419 s/step Meep progress: 150.08333333333334/408.0 = 36.8% done in 8.0s, 13.7s to go on time step 9006 (time=150.1), 0.00089035 s/step Meep progress: 222.83333333333334/408.0 = 54.6% done in 12.0s, 10.0s to go on time step 13373 (time=222.883), 0.000916091 s/step Meep progress: 297.3666666666667/408.0 = 72.9% done in 16.0s, 6.0s to go on time step 17846 (time=297.433), 0.000894295 s/step Meep progress: 371.5/408.0 = 91.1% done in 20.0s, 2.0s to go on time step 22295 (time=371.583), 0.000899168 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 64 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.88675 tran (uniaxial):, 1, 4.77, 0.01655 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 7.20024e-05 s Working in 2D dimensions. Computational cell is 10 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.155606 s ----------- Meep progress: 63.06666666666666/208.0 = 30.3% done in 4.0s, 9.2s to go on time step 3784 (time=63.0667), 0.0010573 s/step Meep progress: 127.85/208.0 = 61.5% done in 8.0s, 5.0s to go on time step 7672 (time=127.867), 0.001029 s/step Meep progress: 191.65/208.0 = 92.1% done in 12.0s, 1.0s to go on time step 11501 (time=191.683), 0.00104491 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.00951e-05 s Working in 2D dimensions. Computational cell is 10 x 6.5 x 0 with resolution 30 block, center = (-4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (4.44089e-16,0,0) size (6,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 32.019% done, 9.22186 s remaining subpixel-averaging is 51.2317% done, 3.98466 s remaining subpixel-averaging is 70.4445% done, 1.77788 s remaining subpixel-averaging is 32.019% done, 9.30279 s remaining subpixel-averaging is 51.2317% done, 4.03405 s remaining subpixel-averaging is 70.4445% done, 1.79632 s remaining subpixel-averaging is 32.019% done, 9.17258 s remaining subpixel-averaging is 51.2317% done, 4.02283 s remaining subpixel-averaging is 70.4445% done, 1.78612 s remaining time for set_epsilon = 50.0774 s ----------- Meep progress: 55.21666666666667/408.0 = 13.5% done in 4.0s, 25.6s to go on time step 3313 (time=55.2167), 0.00120767 s/step Meep progress: 112.33333333333333/408.0 = 27.5% done in 8.0s, 21.1s to go on time step 6741 (time=112.35), 0.00116717 s/step Meep progress: 168.81666666666666/408.0 = 41.4% done in 12.0s, 17.0s to go on time step 10131 (time=168.85), 0.00118012 s/step Meep progress: 225.75/408.0 = 55.3% done in 16.0s, 12.9s to go on time step 13549 (time=225.817), 0.00117055 s/step Meep progress: 281.48333333333335/408.0 = 69.0% done in 20.0s, 9.0s to go on time step 16894 (time=281.567), 0.00119605 s/step Meep progress: 338.5333333333333/408.0 = 83.0% done in 24.0s, 4.9s to go on time step 20317 (time=338.617), 0.00116858 s/step Meep progress: 395.55/408.0 = 96.9% done in 28.0s, 0.9s to go on time step 23739 (time=395.65), 0.00116892 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 48 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 32 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.94211 tran (twisted):, 1, 4.77, 0.00255 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.00011301 s Working in 2D dimensions. Computational cell is 7.2 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.111866 s ----------- Meep progress: 80.61666666666666/208.0 = 38.8% done in 4.0s, 6.3s to go on time step 4837 (time=80.6167), 0.000827066 s/step Meep progress: 162.21666666666667/208.0 = 78.0% done in 8.0s, 2.3s to go on time step 9734 (time=162.233), 0.000816884 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.91414e-05 s Working in 2D dimensions. Computational cell is 7.2 x 6.5 x 0 with resolution 30 block, center = (-2.6,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (3.2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 49.3713% done, 4.4399 s remaining subpixel-averaging is 78.9961% done, 1.12404 s remaining subpixel-averaging is 49.3713% done, 4.44133 s remaining subpixel-averaging is 78.9961% done, 1.12126 s remaining subpixel-averaging is 49.3713% done, 4.42219 s remaining subpixel-averaging is 78.9961% done, 1.12413 s remaining time for set_epsilon = 26.6699 s ----------- Meep progress: 73.56666666666666/408.0 = 18.0% done in 4.0s, 18.2s to go on time step 4414 (time=73.5667), 0.00090629 s/step Meep progress: 147.21666666666667/408.0 = 36.1% done in 8.0s, 14.2s to go on time step 8834 (time=147.233), 0.000905054 s/step Meep progress: 221.16666666666666/408.0 = 54.2% done in 12.0s, 10.1s to go on time step 13272 (time=221.2), 0.000901458 s/step Meep progress: 294.68333333333334/408.0 = 72.2% done in 16.0s, 6.2s to go on time step 17684 (time=294.733), 0.000906649 s/step Meep progress: 368.93333333333334/408.0 = 90.4% done in 20.0s, 2.1s to go on time step 22140 (time=369), 0.000897794 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 60 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 44 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 45 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 28 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 65 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 43 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 37 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.93860 tran (uniaxial):, 1, 4.77, 0.00001 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 0.000154018 s Working in 2D dimensions. Computational cell is 10.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.161312 s ----------- Meep progress: 60.9/208.0 = 29.3% done in 4.0s, 9.7s to go on time step 3654 (time=60.9), 0.00109495 s/step Meep progress: 122.28333333333333/208.0 = 58.8% done in 8.0s, 5.6s to go on time step 7338 (time=122.3), 0.00108598 s/step Meep progress: 183.23333333333332/208.0 = 88.1% done in 12.0s, 1.6s to go on time step 10996 (time=183.267), 0.00109353 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.10487e-05 s Working in 2D dimensions. Computational cell is 10.4 x 6.5 x 0 with resolution 30 block, center = (-4.2,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (6.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 30.4883% done, 9.98227 s remaining subpixel-averaging is 48.7824% done, 4.42182 s remaining subpixel-averaging is 67.0766% done, 2.07037 s remaining subpixel-averaging is 85.3708% done, 0.733933 s remaining subpixel-averaging is 30.4883% done, 9.99655 s remaining subpixel-averaging is 48.7824% done, 4.44853 s remaining subpixel-averaging is 67.0766% done, 2.11095 s remaining subpixel-averaging is 85.3708% done, 0.733249 s remaining subpixel-averaging is 30.4883% done, 9.9553 s remaining subpixel-averaging is 48.7824% done, 4.42917 s remaining subpixel-averaging is 67.0766% done, 2.08722 s remaining subpixel-averaging is 85.3708% done, 0.729907 s remaining time for set_epsilon = 53.5962 s ----------- Meep progress: 54.233333333333334/408.0 = 13.3% done in 4.0s, 26.1s to go on time step 3254 (time=54.2333), 0.00122939 s/step Meep progress: 108.7/408.0 = 26.6% done in 8.0s, 22.0s to go on time step 6523 (time=108.717), 0.0012238 s/step Meep progress: 163.1/408.0 = 40.0% done in 12.0s, 18.0s to go on time step 9787 (time=163.117), 0.00122558 s/step Meep progress: 217.76666666666665/408.0 = 53.4% done in 16.0s, 14.0s to go on time step 13068 (time=217.8), 0.00121939 s/step Meep progress: 272.51666666666665/408.0 = 66.8% done in 20.0s, 9.9s to go on time step 16354 (time=272.567), 0.0012174 s/step Meep progress: 327.21666666666664/408.0 = 80.2% done in 24.0s, 5.9s to go on time step 19637 (time=327.283), 0.00121851 s/step Meep progress: 380.9166666666667/408.0 = 93.4% done in 28.0s, 2.0s to go on time step 22860 (time=381), 0.00124152 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 46 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 39 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 29 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 62 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 42 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 30 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.77362 tran (twisted):, 1, 4.77, 0.07215 tran (twisted):, 2, 9.56, 0.00000 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 3.60012e-05 s Working in 2D dimensions. Computational cell is 7.4 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.107628 s ----------- Meep progress: 78.06666666666666/208.0 = 37.5% done in 4.0s, 6.7s to go on time step 4684 (time=78.0667), 0.000854101 s/step Meep progress: 157.66666666666666/208.0 = 75.8% done in 8.0s, 2.6s to go on time step 9461 (time=157.683), 0.000837429 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 6.69956e-05 s Working in 2D dimensions. Computational cell is 7.4 x 6.5 x 0 with resolution 30 block, center = (-2.7,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (3.4,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 47.5314% done, 4.80929 s remaining subpixel-averaging is 76.0521% done, 1.33452 s remaining subpixel-averaging is 47.5314% done, 4.83031 s remaining subpixel-averaging is 76.0521% done, 1.34121 s remaining subpixel-averaging is 47.5314% done, 4.85467 s remaining subpixel-averaging is 76.0521% done, 1.34843 s remaining time for set_epsilon = 28.6162 s ----------- Meep progress: 69.88333333333333/408.0 = 17.1% done in 4.0s, 19.4s to go on time step 4193 (time=69.8833), 0.000954222 s/step Meep progress: 141.2/408.0 = 34.6% done in 8.0s, 15.1s to go on time step 8473 (time=141.217), 0.000934654 s/step Meep progress: 212.45/408.0 = 52.1% done in 12.0s, 11.0s to go on time step 12749 (time=212.483), 0.000935672 s/step Meep progress: 283.9166666666667/408.0 = 69.6% done in 16.0s, 7.0s to go on time step 17038 (time=283.967), 0.000932807 s/step Meep progress: 355.21666666666664/408.0 = 87.1% done in 20.0s, 3.0s to go on time step 21317 (time=355.283), 0.000934988 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 50 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 34 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 31 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 53 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 40 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 33 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 33 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (uniaxial):, 0, 0.00, 0.90581 tran (uniaxial):, 1, 4.77, 0.01950 tran (uniaxial):, 2, 9.56, 0.00000 tran (uniaxial):, 3, 14.43, 0.00000 tran (uniaxial):, 4, 19.41, 0.00000 ----------- Initializing structure... time for choose_chunkdivision = 6.91414e-05 s Working in 2D dimensions. Computational cell is 10.8 x 6.5 x 0 with resolution 30 time for set_epsilon = 0.160437 s ----------- Meep progress: 59.15/208.0 = 28.4% done in 4.0s, 10.1s to go on time step 3549 (time=59.15), 0.00112713 s/step Meep progress: 119.0/208.0 = 57.2% done in 8.0s, 6.0s to go on time step 7141 (time=119.017), 0.00111385 s/step Meep progress: 179.25/208.0 = 86.2% done in 12.0s, 1.9s to go on time step 10757 (time=179.283), 0.00110643 s/step run 0 finished at t = 208.0 (12480 timesteps) ----------- Initializing structure... time for choose_chunkdivision = 7.39098e-05 s Working in 2D dimensions. Computational cell is 10.8 x 6.5 x 0 with resolution 30 block, center = (-4.4,0,0) size (2,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (2.4025,2.4025,2.4025) block, center = (0,0,0) size (6.8,1e+20,1e+20) axes (1,0,0), (0,1,0), (0,0,1) subpixel-averaging is 29.0972% done, 10.5851 s remaining subpixel-averaging is 46.5566% done, 4.83 s remaining subpixel-averaging is 64.0161% done, 2.40549 s remaining subpixel-averaging is 81.4755% done, 0.964889 s remaining subpixel-averaging is 29.0972% done, 10.5591 s remaining subpixel-averaging is 46.5566% done, 4.83987 s remaining subpixel-averaging is 64.0161% done, 2.39451 s remaining subpixel-averaging is 81.4755% done, 0.972159 s remaining subpixel-averaging is 29.0972% done, 10.6613 s remaining subpixel-averaging is 46.5566% done, 4.82302 s remaining subpixel-averaging is 64.0161% done, 2.39088 s remaining subpixel-averaging is 81.4755% done, 0.967715 s remaining time for set_epsilon = 56.7605 s ----------- Meep progress: 51.55/408.0 = 12.6% done in 4.0s, 27.7s to go on time step 3093 (time=51.55), 0.00129338 s/step Meep progress: 103.8/408.0 = 25.4% done in 8.0s, 23.4s to go on time step 6229 (time=103.817), 0.00127578 s/step Meep progress: 157.21666666666667/408.0 = 38.5% done in 12.0s, 19.1s to go on time step 9435 (time=157.25), 0.0012479 s/step Meep progress: 210.46666666666667/408.0 = 51.6% done in 16.0s, 15.0s to go on time step 12631 (time=210.517), 0.00125168 s/step Meep progress: 263.68333333333334/408.0 = 64.6% done in 20.0s, 10.9s to go on time step 15825 (time=263.75), 0.0012525 s/step Meep progress: 316.7833333333333/408.0 = 77.6% done in 24.0s, 6.9s to go on time step 19011 (time=316.85), 0.0012555 s/step Meep progress: 370.1166666666667/408.0 = 90.7% done in 28.0s, 2.9s to go on time step 22212 (time=370.2), 0.00124999 s/step run 0 finished at t = 408.0 (24480 timesteps) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 61 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 45 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 43 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 35 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) MPB solved for omega_1(1.85185,0,0) = 1.85185 after 56 iters Dominant planewave for band 1: (1.851852,-0.000000,0.000000) MPB solved for omega_2(1.85185,0,0) = 1.85823 after 51 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters MPB solved for omega_2(1.84545,0,0) = 1.85185 after 1 iters Dominant planewave for band 2: (1.845450,-0.153846,0.000000) MPB solved for omega_3(1.85185,0,0) = 1.87724 after 41 iters MPB solved for omega_3(1.82612,0,0) = 1.85186 after 1 iters MPB solved for omega_3(1.82611,0,0) = 1.85185 after 1 iters Dominant planewave for band 3: (1.826111,-0.307692,0.000000) MPB solved for omega_4(1.85185,0,0) = 1.9085 after 30 iters MPB solved for omega_4(1.79347,0,0) = 1.85191 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters MPB solved for omega_4(1.79342,0,0) = 1.85185 after 1 iters Dominant planewave for band 4: (1.793415,-0.461538,0.000000) MPB solved for omega_5(1.85185,0,0) = 1.95142 after 32 iters MPB solved for omega_5(1.74693,0,0) = 1.85215 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters MPB solved for omega_5(1.74661,0,0) = 1.85185 after 1 iters Dominant planewave for band 5: (1.746613,-0.615385,0.000000) tran (twisted):, 0, 0.00, 0.49695 tran (twisted):, 1, 4.77, 0.22184 tran (twisted):, 2, 9.56, 0.00001 tran (twisted):, 3, 14.43, 0.00000 tran (twisted):, 4, 19.41, 0.00000
The diffraction spectra is plotted using the script below and shown in the accompanying figure.
cos_angles = [math.cos(math.radians(t)) for t in ang_uniaxial]
tran = m0_uniaxial + 2 * m1_uniaxial
eff_m0 = m0_uniaxial / tran
eff_m1 = (2 * m1_uniaxial / tran) / cos_angles
phase = delta_n * dd / wvl
eff_m0_analytic = [math.cos(math.pi * p) ** 2 for p in phase]
eff_m1_analytic = [math.sin(math.pi * p) ** 2 for p in phase]
plt.figure(dpi=150)
plt.subplot(1, 2, 1)
plt.plot(phase, eff_m0, "bo-", clip_on=False, label="0th order (meep)")
plt.plot(phase, eff_m0_analytic, "b--", clip_on=False, label="0th order (analytic)")
plt.plot(phase, eff_m1, "ro-", clip_on=False, label="±1 orders (meep)")
plt.plot(phase, eff_m1_analytic, "r--", clip_on=False, label="±1 orders (analytic)")
plt.axis([0, 1.0, 0, 1])
plt.xticks([t for t in np.arange(0, 1.2, 0.2)])
plt.xlabel("phase delay Δnd/λ")
plt.ylabel("diffraction efficiency @ λ = 0.54 μm")
plt.legend(loc="center")
plt.title("homogeneous uniaxial grating")
cos_angles = [math.cos(math.radians(t)) for t in ang_twisted]
tran = m0_twisted + 2 * m1_twisted
eff_m0 = m0_twisted / tran
eff_m1 = (2 * m1_twisted / tran) / cos_angles
plt.subplot(1, 2, 2)
plt.plot(phase, eff_m0, "bo-", clip_on=False, label="0th order (meep)")
plt.plot(phase, eff_m1, "ro-", clip_on=False, label="±1 orders (meep)")
plt.axis([0, 1.0, 0, 1])
plt.xticks([t for t in np.arange(0, 1.2, 0.2)])
plt.xlabel("phase delay Δnd/λ")
plt.ylabel("diffraction efficiency @ λ = 0.54 μm")
plt.legend(loc="center")
plt.title("bilayer twisted-nematic grating")
plt.tight_layout()
plt.show()
The left figure shows good agreement between the simulation results and analytic theory for the homogeneous uniaxial grating. Approximately 6% of the power in the input planewave is lost due to reflection from the grating. This value is an average over all phase delays. The total transmittance is therefore around 94%. The twisted-nematic grating, with results shown in the right figure, produces ±1 diffraction orders with nearly-constant peak transmittance over a broader bandwidth around Δnd/λ=0.5 than the homogeneous uniaxial polarization grating. This is consistent with results from the reference. The average reflectance and transmittance for the twisted-nematic grating are similar to those for the homogeneous uniaxial grating.
Finally, we demonstrate that when Δnd/λ=0.5 a circular-polarized planewave input produces just a single ±1 diffraction order. To specify a Ez+iEy circular-polarized planewave requires setting the amplitude
of the Ey source to an imaginary number (from its default of 1):
sources = [mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ez, center=src_pt, size=mp.Vector3(0,sy,0)),
mp.Source(mp.GaussianSource(fcen,fwidth=0.05*fcen), component=mp.Ey, center=src_pt, size=mp.Vector3(0,sy,0), amplitude=1j)]
Note that even though the Jy current amplitude is complex in this example, only its real part is used and the resulting fields are therefore still real (the default).
The figure below shows a snapshot of Ez within the cell for four different cases: phase delays (Δnd/λ) of 0.5 and 1.0, and planewave circular polarization of Ez+iEy and Ez-iEy. The empty regions on the cell sides are PMLs. The thin solid black line denotes the boundary between the grating (on the left) and air. As expected, for Δnd/λ=0.5 there is just a single ±1 diffraction order which depends on the chirality of the input planewave (this is not the case for a linear-polarized planewave). The angle of this diffracted order (±4.8°) agrees with the analytic result. Snapshots of Ey are similar.