# as we are using autograd to perform gradient calculations, we need to import numpy from autograd's wrapper
import autograd as ag
import autograd.numpy as np
import matplotlib.pyplot as plt
import tidy3d as td
import tidy3d.web as web
# we will use Tidy3D's built-in Adam helper to perform the optimization
from tidy3d.plugins.autograd import adam, apply_updates
TIDY3D
LEARNING CENTER
Example Library
Inverse design of an apodized grating coupler through shape optimization¶
A grating coupler (GC) is a key photonic device used to couple light between optical fibers and on-chip waveguides by diffracting light at a specific angle. These devices often employ apodization, a technique that modifies the grating tooth parameters along the device length to improve efficiency and minimize reflection losses.
In this notebook, we will use gradient-based optimization to determine the optimal grating tooth distribution. We start from a uniform grating coupler and modify the tooth and gap sizes simultaneously using gradient based optimization. Feature size constraints on the parameters enable the device to remain fabricable over the course of optimization. In just a handful of iterations, we achieve a design that outperforms an apodized grating coupler designed using semi-analytical methods.
While the simulations are performed using a 2D cross section, the same approach can be modified to 3D without much further complication.
If you are interested in other inverse design examples using tidy3d, you can find many of them here.
If you are unfamiliar with inverse design, we also recommend our intro to inverse design tutorials and our primer on automatic differentiation with tidy3d.
Basic Parameters
First, we define the basic parameters used to define a uniform grating coupler for operation primarily in the C-band.
lda0 = 1.55 # central wavelength
n_ldas = 101 # number of wavelength points
ldas = np.linspace(1.5, 1.6, n_ldas) # wavelength range
freq0 = td.C_0 / lda0 # central frequency
freqs = td.C_0 / ldas # frequency range
num_freqs_objective = 5
ldas_objective = np.linspace(1.5, 1.6, num_freqs_objective)
freqs_objective = td.C_0 / ldas_objective
fwidth = 0.5 * (np.max(freqs) - np.min(freqs)) # width of the source frequency rangeThe material platform is the common silicon on insulator (SOI) with a 260 nm silicon thickness and 2 µm BOX thickness. The top cladding layer is 680 nm oxide. The grating will be partially etched with an etching depth of 160 nm.
# define materials from the material library
si = td.material_library["cSi"]["Palik_LowLoss"]
sio2 = td.material_library["SiO2"]["Palik_LowLoss"]
t_si = 0.26 # thickness of the silicon layer
etch_depth = 0.16 # etching depth
t_tox = 0.68 # top oxide layer thickness
t_box = 2 # bottom oxide layer thicknessWe aim to design the GC for the standard single-mode fiber (SMF) with a mode field diameter (MFD) of 10.8 µm. The fiber is tilted at a 14.5 degree angle.
theta = np.deg2rad(14.5) # fiber tilt angle
mfd = 10.8 # mode field diameter
source_x = 4.5 # x position of the fiber
N = 16 # number of grating teeth to create
inf_eff = 1e3 # effective infinity
buffer = 1.1 * lda0 # buffer spacing to pad the simulation domainNext we compute the uniform tooth parameters (periodicity and fill fraction) to maximize efficiency.
neff_unetch = 2.9 # effective index of the slab mode of the unetched waveguide
neff_etch = 2.2 # effective index of the slab mode of the etched waveguide
n_c = 1.44 # refractive index of the cladding (SiO2)
theta_c = np.sin(theta) / n_c # incident angle in the cladding
def get_periodicity(fill_fraction: float) -> float:
"""periodicity (bragg condition) as function of fill fraction, angle, wavelength, and etching parameters."""
return lda0 / (
fill_fraction * neff_unetch + (1 - fill_fraction) * neff_etch - n_c * np.sin(theta_c)
)
f0 = 0.80
periodicity = get_periodicity(f0)
p_list = N * [periodicity]
f_list = N * [f0]After obtaining \(p_{\text{i}}\) and \(f_{\text{i}}\), we will write a function to return the width of each air gap and each tooth. We will use these gap widths as the parameters for our design. The advantage is that we can impose bounds on these values to ensure the minimal feature size is above the fabrication constraint.
Here we aim to have a minimal feature size above 60 nm. The current design has a minimal feature size right above it.
def get_widths(p_list, f_list):
# calculate the widths of air gaps and silicon teeth
return np.array([item for p, f in zip(p_list, f_list) for item in [p * (1 - f), p * f]])
widths = get_widths(p_list, f_list)
l_grating = np.sum(widths) + 3 * lda0 # total length of the GC
# print the current minimal feature size
print(f"The minimal feature size is {1e3 * np.min(widths):.2f} nm.")The minimal feature size is 123.46 nm.To ensure the minimal feature size is maintained during the optimization, we will use a design parameter in the range from -inf to inf and project it to a tanh function that is bounded by the minimal and maximal feature sizes.
min_width = 0.08 # minimal feature size we want to maintain
max_width = 1 # maximal feature size
# function to project a design parameter between -inf to inf to between min_width and max_width
def project(x):
return 0.5 * (max_width - min_width) * np.tanh(x) + 0.5 * (max_width + min_width)
# function to inversely project a design parameter between min_width and max_width to between -inf to inf
def inverse_project(y):
return np.arctanh((2 * (y - 0.5 * (max_width + min_width))) / (max_width - min_width))
# project the widths to parameters between -inf to inf
params0 = inverse_project(widths)Define Static Components of the Simulation
Next we will define the components that don’t change during the optimization. These include the cladding layer, the output waveguide, the unetched layer, the bottom oxide layer, and the silicon substrate.
mask = td.Structure(
geometry=td.Box(size=(td.inf, td.inf, td.inf), center=(0, 0, 0)),
medium=td.Medium(),
)
# create the top oxide layer
tox = td.Structure(
geometry=td.Box.from_bounds(rmin=(-inf_eff, -inf_eff, 0), rmax=(inf_eff, inf_eff, t_tox)),
medium=sio2,
)
# create the output waveguide
slab_waveguide = td.Structure(
geometry=td.Box.from_bounds(rmin=(-inf_eff, -inf_eff, 0), rmax=(0, inf_eff, t_si)),
medium=si,
)
# create the unetched waveguide
unetched_waveguide = td.Structure(
geometry=td.Box.from_bounds(rmin=(0, -inf_eff, 0), rmax=(inf_eff, inf_eff, t_si - etch_depth)),
medium=si,
)
# create the bottom oxide layer
box = td.Structure(
geometry=td.Box.from_bounds(rmin=(-inf_eff, -inf_eff, -t_box), rmax=(inf_eff, inf_eff, 0)),
medium=sio2,
)
# create the silicon substrate layer
substrate = td.Structure(
geometry=td.Box.from_bounds(
rmin=(-inf_eff, -inf_eff, -inf_eff), rmax=(inf_eff, inf_eff, -t_box)
),
medium=si,
)The source and monitor are not changed either. We use a GaussianBeam to represent the incident fiber mode. A ModeMonitor is placed at the output waveguide to measure the coupling efficiency.
gap = 1 # gap size between the source plane and the top surface of the cladding
# define a gaussian beam source
source = td.GaussianBeam(
size=(2 * mfd, td.inf, 0),
center=[source_x, 0, t_tox + gap],
source_time=td.GaussianPulse(freq0=freq0, fwidth=fwidth),
angle_theta=theta,
direction="-",
waist_radius=mfd / 2,
pol_angle=np.pi / 2, # 90 degree polarization angle for TE polarization
)
# define a mode monitor
mode_monitor = td.ModeMonitor(
center=(-buffer / 2, 0, t_si / 2),
size=(0, td.inf, 6 * t_si),
freqs=freqs,
mode_spec=td.ModeSpec(num_modes=1, target_neff=3),
name="mode",
)Simulate the Designed Apodized GC
To simulate the designed GC and facility future optimization, we define a function that takes in an array of design parameters and returns a Simulation object. Note again that each design parameter has a range from -inf to inf. It will then be projected to a range between min_width and max_width.
def make_2d_sim(design_parameters):
# calculate the widths of air gaps and silicon teeth
widths_si = project(design_parameters[1::2])
widths_air = project(design_parameters[::2])
# initialize the center and size of each silicon teeth
center = 0
size = 0
# create the grating geometries from the given widths
gratings = 0
for width_si, width_air in zip(widths_si, widths_air):
center += width_air + width_si / 2
size = width_si
gratings += td.Box(
center=(center, 0, t_si - etch_depth / 2), size=(size, td.inf, etch_depth)
)
center += width_si / 2
# create the grating structure
gratings = td.Structure(geometry=gratings, medium=si)
# create a box to represent the simulation domain box
sim_box = td.Box.from_bounds(
rmin=(-buffer, 0, -t_box - buffer / 2),
rmax=(l_grating + buffer, 0, t_si + buffer),
)
run_time = 1e-12 # simulation run time
# construct simulation
sim = td.Simulation(
center=sim_box.center,
size=sim_box.size,
grid_spec=td.GridSpec.auto(
min_steps_per_wvl=30, wavelength=lda0
), # use a fine grid to ensure the small features are well resolved
structures=[
mask,
tox,
gratings,
unetched_waveguide,
slab_waveguide,
box,
substrate,
],
sources=[source],
monitors=[mode_monitor],
medium=sio2,
run_time=run_time,
boundary_spec=td.BoundarySpec(
x=td.Boundary.pml(),
y=td.Boundary.periodic(), # set the boundary to periodic in y since it's a 2D simulation
z=td.Boundary.pml(),
),
)
return simCreate the simulation for the designed apodized GC and visualize it.
sim0 = make_2d_sim(params0)
sim0.plot_eps(y=0, freq=freq0)
plt.show()
Also visualize the grid to ensure it resolves the small features well.
ax = sim0.plot(y=0)
sim0.plot_grid(y=0, ax=ax)
ax.set_xlim(-0.1, 1)
ax.set_ylim(0, 0.3)
plt.show()
Run the simulation.
sim_data0 = web.run(simulation=sim0, task_name="initial design")08:08:31 UTC Created task 'initial design' with resource_id 'fdve-3e5ce965-8f28-466d-a9b2-b3d07c101238' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-3e5ce965-8f2 8-466d-a9b2-b3d07c101238'.
Task folder: 'default'.
08:08:32 UTC Estimated FlexCredit cost: 0.025. This assumes the FDTD solver runs for the full simulation time; if early shutoff is reached, the billed cost can be lower. Use 'web.real_cost(task_id)' to get the billed FlexCredit cost after a simulation run.
08:08:35 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:08:47 UTC starting up solver
08:08:48 UTC running solver
08:08:53 UTC early shutoff detected at 47%, exiting.
status = success
View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-3e5ce965-8f2 8-466d-a9b2-b3d07c101238'.
08:08:54 UTC Loading results from simulation_data.hdf5
Result Visualization
After the simulation is done, we can extract the coupling efficiency from the mode monitor data. In addition, we also calculate two important metrics: the maximum coupling efficiency and the 1dB bandwidth. The maximal coupling efficiency is about 2.7dB and the 1dB bandwidth is about 36 nm.
# function to calculate the coupling efficiency from simulation data
def compute_ce(sim_data):
amp = sim_data["mode"].amps.sel(mode_index=0, direction="-").values
ce = np.abs(amp) ** 2 # transmission to the top waveguide
return ce
def dB(ce):
return 10 * np.log10(ce)
# calculate the coupling efficiency for the designed GC
ce0 = compute_ce(sim_data0)
# plot the coupling efficiency
plt.plot(ldas, dB(ce0), c="red", linewidth=2)
plt.xlim(min(ldas), max(ldas))
plt.ylim(-10, 0)
plt.xlabel("Wavelength (µm)")
plt.ylabel("Coupling efficiency (dB)")
plt.grid()
plt.show()
# function to calculate the 1dB bandwidth from the coupling efficiency
def bandwidth(ldas, ce):
max_ce = np.max(dB(ce))
threshold = max_ce - 1
within_1db = np.where(dB(ce) >= threshold)[0]
lambda_min = ldas[within_1db[0]]
lambda_max = ldas[within_1db[-1]]
bandwidth = lambda_max - lambda_min
return 1e3 * bandwidth
# print the maximum coupling efficiency and 1dB bandwidth
print(f"The 1dB bandwidth is {bandwidth(ldas, ce0):.1f} nm")
print(f"The maximum coupling efficiency is {np.max(dB(ce0)):.2f} dB.")
The 1dB bandwidth is 36.0 nm
The maximum coupling efficiency is -2.74 dB.Optimize GC with Inverse Design to Maximize the Efficiency and Bandwidth
In this part, we will apply inverse design to maximize the coupling efficiency of our coupling over the wavelength range from 1500 to 1600 nm.
To do this, we need to express our objective as a function of our parameters returning a single float value to maximize. This function will involve constructing, running, and postprocessing our simulation.
def J(design_parameters: np.ndarray) -> float:
sim = make_2d_sim(design_parameters)
sim_data = web.run(sim, task_name="GC_invdes", verbose=True)
ce = np.sum(compute_ce(sim_data)) / n_ldas
return cedJ = ag.value_and_grad(J)Before running the optimization, we can check the value and gradient of the initial design to ensure the gradient tracking is working properly.
%%time
val, grad = dJ(params0)
print(val)
print(grad)08:08:55 UTC Created task 'GC_invdes' with resource_id 'fdve-0aa7ca50-bcef-42c7-92e6-e759db293886' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-0aa7ca50-bce f-42c7-92e6-e759db293886'.
Task folder: 'default'.
08:08:59 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:09:13 UTC starting up solver
running solver
08:09:19 UTC early shutoff detected at 47%, exiting.
status = postprocess
08:09:26 UTC status = success
08:09:28 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-0aa7ca50-bce f-42c7-92e6-e759db293886'.
08:09:29 UTC Loading results from simulation_data.hdf5
08:09:30 UTC Started working on Batch containing 1 tasks.
08:09:33 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:09:54 UTC Batch complete.
0.3390465892159393
[-1.17321182e-02 6.28505441e-01 2.69946648e-02 5.84308310e-01
-1.64972445e-03 1.32803981e-01 -2.03079973e-03 -2.04572945e-02
9.64456146e-04 6.00574356e-02 1.18689844e-02 5.72078254e-02
1.55931666e-02 -2.16742048e-02 4.13397433e-03 -7.84230979e-02
7.88709567e-03 -7.10540198e-02 1.24795132e-03 -8.56533503e-02
4.50736483e-03 -7.83018856e-02 -1.96388528e-03 -7.23079247e-02
2.53005638e-04 -5.06262475e-02 -3.73559248e-04 -3.75730912e-02
-6.11800676e-04 -2.56284733e-02 -5.57500031e-04 -1.48503663e-02]
CPU times: user 2.86 s, sys: 321 ms, total: 3.19 s
Wall time: 1minWe will use the Adam optimizer to adjust the parameters for a number of iterations until we achieve a good performance.
%%time
# hyperparameters
num_steps = 25
learning_rate = 0.05
# initialize adam optimizer with starting parameters
params = params0
optimizer = adam(learning_rate=learning_rate)
opt_state = optimizer.init(params)
# store history
J_history = []
params_history = []
for i in range(num_steps):
# compute gradient and current objective function value
value, gradient = dJ(params)
# outputs
print(f"step = {i + 1}")
print(f"\tJ = {value:.3e}")
print(f"\tgrad_norm = {np.linalg.norm(gradient):.4e}")
# compute and apply updates to the optimizer based on gradient (-1 sign to maximize obj_fn)
updates, opt_state = optimizer.update(-gradient, opt_state, params)
params[:] = apply_updates(params, updates)
# save history
J_history.append(value)
params_history.append(params.copy())08:09:56 UTC Created task 'GC_invdes' with resource_id 'fdve-dff112db-b705-4a4e-bb41-bf77e6896cac' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dff112db-b70 5-4a4e-bb41-bf77e6896cac'.
Task folder: 'default'.
08:09:59 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:10:16 UTC starting up solver
running solver
08:10:20 UTC early shutoff detected at 47%, exiting.
status = postprocess
08:10:27 UTC status = success
08:10:29 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dff112db-b70 5-4a4e-bb41-bf77e6896cac'.
08:10:30 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:10:33 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:11:11 UTC Batch complete.
step = 1
J = 3.390e-01
grad_norm = 8.9328e-0108:11:12 UTC Created task 'GC_invdes' with resource_id 'fdve-3234666f-dfec-4dad-a819-e9508cb01c71' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-3234666f-dfe c-4dad-a819-e9508cb01c71'.
Task folder: 'default'.
08:11:17 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:11:27 UTC status = preprocess
08:11:31 UTC starting up solver
08:11:32 UTC running solver
08:11:35 UTC early shutoff detected at 48%, exiting.
08:11:36 UTC status = postprocess
08:11:43 UTC status = success
08:11:45 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-3234666f-dfe c-4dad-a819-e9508cb01c71'.
08:11:46 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:11:50 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:12:12 UTC Batch complete.
step = 2
J = 3.173e-01
grad_norm = 9.2350e-0108:12:13 UTC Created task 'GC_invdes' with resource_id 'fdve-18acd1d0-2219-473e-8511-1a9bc0cd08d9' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-18acd1d0-221 9-473e-8511-1a9bc0cd08d9'.
Task folder: 'default'.
08:12:19 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:12:36 UTC starting up solver
running solver
08:12:42 UTC early shutoff detected at 47%, exiting.
status = postprocess
08:12:47 UTC status = success
08:12:49 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-18acd1d0-221 9-473e-8511-1a9bc0cd08d9'.
08:12:50 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:12:55 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:13:16 UTC Batch complete.
step = 3
J = 3.793e-01
grad_norm = 5.8213e-0108:13:17 UTC Created task 'GC_invdes' with resource_id 'fdve-e081b902-5d98-485c-b2de-70e9b3165cfe' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-e081b902-5d9 8-485c-b2de-70e9b3165cfe'.
Task folder: 'default'.
08:13:21 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:13:33 UTC status = preprocess
08:13:38 UTC starting up solver
running solver
08:13:42 UTC early shutoff detected at 46%, exiting.
status = postprocess
08:13:49 UTC status = success
08:13:51 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-e081b902-5d9 8-485c-b2de-70e9b3165cfe'.
08:13:52 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:13:56 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:14:18 UTC Batch complete.
step = 4
J = 3.860e-01
grad_norm = 7.2584e-0108:14:20 UTC Created task 'GC_invdes' with resource_id 'fdve-47926689-ad24-434e-922c-4802e892e631' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-47926689-ad2 4-434e-922c-4802e892e631'.
Task folder: 'default'.
08:14:26 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:14:37 UTC status = preprocess
08:14:42 UTC starting up solver
running solver
08:14:47 UTC early shutoff detected at 45%, exiting.
status = postprocess
08:14:55 UTC status = success
08:14:57 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-47926689-ad2 4-434e-922c-4802e892e631'.
08:14:58 UTC Loading results from simulation_data.hdf5
08:14:59 UTC Started working on Batch containing 1 tasks.
08:15:02 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:15:24 UTC Batch complete.
step = 5
J = 3.884e-01
grad_norm = 5.9201e-0108:15:25 UTC Created task 'GC_invdes' with resource_id 'fdve-a965e191-2c2a-4d94-99f5-cf2b2220a1a8' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-a965e191-2c2 a-4d94-99f5-cf2b2220a1a8'.
Task folder: 'default'.
08:15:31 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:15:45 UTC starting up solver
running solver
08:15:46 UTC status = queued
08:16:12 UTC status = running
08:16:19 UTC status = postprocess
08:16:26 UTC status = success
08:16:28 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-a965e191-2c2 a-4d94-99f5-cf2b2220a1a8'.
Loading results from simulation_data.hdf5
08:16:29 UTC Started working on Batch containing 1 tasks.
08:16:32 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:16:52 UTC Batch complete.
step = 6
J = 3.900e-01
grad_norm = 6.7738e-0108:16:53 UTC Created task 'GC_invdes' with resource_id 'fdve-9d87e67b-bb0a-482a-abf3-0ead79fd2e55' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-9d87e67b-bb0 a-482a-abf3-0ead79fd2e55'.
Task folder: 'default'.
08:16:58 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:17:11 UTC starting up solver
running solver
08:17:16 UTC early shutoff detected at 47%, exiting.
status = postprocess
08:17:24 UTC status = success
08:17:26 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-9d87e67b-bb0 a-482a-abf3-0ead79fd2e55'.
Loading results from simulation_data.hdf5
08:17:27 UTC Started working on Batch containing 1 tasks.
08:17:31 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:17:58 UTC Batch complete.
step = 7
J = 4.029e-01
grad_norm = 6.0698e-0108:17:59 UTC Created task 'GC_invdes' with resource_id 'fdve-8bdb1ceb-cad8-4f5e-ae94-2def49657f61' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-8bdb1ceb-cad 8-4f5e-ae94-2def49657f61'.
Task folder: 'default'.
08:18:02 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:18:10 UTC status = preprocess
08:18:14 UTC starting up solver
08:18:15 UTC running solver
08:18:20 UTC early shutoff detected at 47%, exiting.
status = postprocess
08:18:27 UTC status = success
08:18:29 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-8bdb1ceb-cad 8-4f5e-ae94-2def49657f61'.
08:18:30 UTC Loading results from simulation_data.hdf5
08:18:31 UTC Started working on Batch containing 1 tasks.
08:18:34 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:18:56 UTC Batch complete.
step = 8
J = 4.214e-01
grad_norm = 4.2170e-0108:18:57 UTC Created task 'GC_invdes' with resource_id 'fdve-cb454847-4706-4757-bc49-6f6411e16160' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-cb454847-470 6-4757-bc49-6f6411e16160'.
Task folder: 'default'.
08:19:01 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:19:11 UTC status = preprocess
08:19:16 UTC starting up solver
running solver
08:19:20 UTC early shutoff detected at 45%, exiting.
08:19:21 UTC status = postprocess
08:19:26 UTC status = success
08:19:28 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-cb454847-470 6-4757-bc49-6f6411e16160'.
Loading results from simulation_data.hdf5
08:19:29 UTC Started working on Batch containing 1 tasks.
08:19:31 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:19:56 UTC Batch complete.
step = 9
J = 4.287e-01
grad_norm = 3.5623e-0108:19:57 UTC Created task 'GC_invdes' with resource_id 'fdve-51449963-13a8-4d7e-9f69-613a8c95c813' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-51449963-13a 8-4d7e-9f69-613a8c95c813'.
Task folder: 'default'.
08:20:02 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:20:16 UTC starting up solver
running solver
08:20:21 UTC early shutoff detected at 43%, exiting.
08:20:22 UTC status = postprocess
08:20:27 UTC status = success
08:20:29 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-51449963-13a 8-4d7e-9f69-613a8c95c813'.
Loading results from simulation_data.hdf5
08:20:30 UTC Started working on Batch containing 1 tasks.
08:20:33 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:20:55 UTC Batch complete.
step = 10
J = 4.324e-01
grad_norm = 5.1949e-0108:20:56 UTC Created task 'GC_invdes' with resource_id 'fdve-a63c1550-2f22-42d0-af93-8778452ecb02' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-a63c1550-2f2 2-42d0-af93-8778452ecb02'.
Task folder: 'default'.
08:20:59 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:21:21 UTC starting up solver
08:21:22 UTC running solver
08:21:27 UTC early shutoff detected at 43%, exiting.
status = postprocess
08:21:35 UTC status = success
08:21:37 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-a63c1550-2f2 2-42d0-af93-8778452ecb02'.
Loading results from simulation_data.hdf5
08:21:38 UTC Started working on Batch containing 1 tasks.
08:21:41 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:22:03 UTC Batch complete.
step = 11
J = 4.373e-01
grad_norm = 4.7082e-0108:22:04 UTC Created task 'GC_invdes' with resource_id 'fdve-572c6f95-8b63-4a90-8d17-5ecbf514f6d5' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-572c6f95-8b6 3-4a90-8d17-5ecbf514f6d5'.
Task folder: 'default'.
08:22:07 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:22:21 UTC starting up solver
08:22:22 UTC running solver
08:22:27 UTC early shutoff detected at 43%, exiting.
status = postprocess
08:22:34 UTC status = success
08:22:36 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-572c6f95-8b6 3-4a90-8d17-5ecbf514f6d5'.
08:22:37 UTC Loading results from simulation_data.hdf5
08:22:38 UTC Started working on Batch containing 1 tasks.
08:22:40 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:23:08 UTC Batch complete.
step = 12
J = 4.450e-01
grad_norm = 4.3408e-0108:23:09 UTC Created task 'GC_invdes' with resource_id 'fdve-381cd62d-417b-49ce-a4cd-bdac542eb083' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-381cd62d-417 b-49ce-a4cd-bdac542eb083'.
Task folder: 'default'.
08:23:14 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:23:28 UTC starting up solver
running solver
08:23:34 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:23:41 UTC status = success
08:23:43 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-381cd62d-417 b-49ce-a4cd-bdac542eb083'.
08:23:44 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:23:47 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:24:09 UTC Batch complete.
step = 13
J = 4.516e-01
grad_norm = 3.2819e-0108:24:10 UTC Created task 'GC_invdes' with resource_id 'fdve-dd57c35a-b048-474d-8278-a00ae7ada0ff' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dd57c35a-b04 8-474d-8278-a00ae7ada0ff'.
Task folder: 'default'.
08:24:15 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:24:27 UTC status = preprocess
08:24:32 UTC starting up solver
running solver
08:24:36 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:24:43 UTC status = success
08:24:45 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dd57c35a-b04 8-474d-8278-a00ae7ada0ff'.
08:24:46 UTC Loading results from simulation_data.hdf5
08:24:47 UTC Started working on Batch containing 1 tasks.
08:24:49 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:25:09 UTC Batch complete.
step = 14
J = 4.549e-01
grad_norm = 3.5407e-0108:25:10 UTC Created task 'GC_invdes' with resource_id 'fdve-fa00e1f8-b97b-4b50-9a09-23817692d1c5' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-fa00e1f8-b97 b-4b50-9a09-23817692d1c5'.
Task folder: 'default'.
08:25:15 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:25:31 UTC starting up solver
running solver
08:25:37 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:25:42 UTC status = success
08:25:44 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-fa00e1f8-b97 b-4b50-9a09-23817692d1c5'.
08:25:45 UTC Loading results from simulation_data.hdf5
08:25:46 UTC Started working on Batch containing 1 tasks.
08:25:48 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:27:11 UTC Batch complete.
step = 15
J = 4.580e-01
grad_norm = 3.0350e-0108:27:12 UTC Created task 'GC_invdes' with resource_id 'fdve-43a96e60-6067-45d2-a7cc-a748117349f9' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-43a96e60-606 7-45d2-a7cc-a748117349f9'.
Task folder: 'default'.
08:27:15 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:27:25 UTC status = preprocess
08:27:29 UTC starting up solver
08:27:30 UTC running solver
08:27:32 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:27:40 UTC status = success
08:27:42 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-43a96e60-606 7-45d2-a7cc-a748117349f9'.
Loading results from simulation_data.hdf5
08:27:43 UTC Started working on Batch containing 1 tasks.
08:27:45 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:28:07 UTC Batch complete.
step = 16
J = 4.596e-01
grad_norm = 3.4958e-0108:28:08 UTC Created task 'GC_invdes' with resource_id 'fdve-1a2ce4c3-befe-4a41-93f5-3cba7053b837' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-1a2ce4c3-bef e-4a41-93f5-3cba7053b837'.
Task folder: 'default'.
08:28:13 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:28:34 UTC status = preprocess
08:28:38 UTC starting up solver
08:28:39 UTC running solver
08:28:41 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:28:48 UTC status = success
08:28:50 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-1a2ce4c3-bef e-4a41-93f5-3cba7053b837'.
08:28:51 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:28:55 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:29:16 UTC Batch complete.
step = 17
J = 4.650e-01
grad_norm = 2.7680e-0108:29:17 UTC Created task 'GC_invdes' with resource_id 'fdve-9fe21ded-984d-4f7a-b163-2c246577d3cc' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-9fe21ded-984 d-4f7a-b163-2c246577d3cc'.
Task folder: 'default'.
08:29:21 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:29:31 UTC status = preprocess
08:29:36 UTC starting up solver
running solver
08:29:40 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:29:47 UTC status = success
08:29:49 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-9fe21ded-984 d-4f7a-b163-2c246577d3cc'.
08:29:50 UTC Loading results from simulation_data.hdf5
08:29:51 UTC Started working on Batch containing 1 tasks.
08:29:56 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:30:31 UTC Batch complete.
step = 18
J = 4.681e-01
grad_norm = 2.3597e-0108:30:32 UTC Created task 'GC_invdes' with resource_id 'fdve-eedde884-4c61-4394-8057-bd12e69f2cd4' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-eedde884-4c6 1-4394-8057-bd12e69f2cd4'.
Task folder: 'default'.
08:30:37 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:31:31 UTC status = preprocess
08:31:35 UTC starting up solver
running solver
08:31:41 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:31:49 UTC status = success
08:31:51 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-eedde884-4c6 1-4394-8057-bd12e69f2cd4'.
08:31:52 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:31:55 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:32:17 UTC Batch complete.
step = 19
J = 4.691e-01
grad_norm = 2.8568e-0108:32:19 UTC Created task 'GC_invdes' with resource_id 'fdve-41d6ea84-7759-4cf6-99fc-c1b35b82256d' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-41d6ea84-775 9-4cf6-99fc-c1b35b82256d'.
Task folder: 'default'.
08:32:23 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:32:38 UTC starting up solver
running solver
08:32:44 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:32:49 UTC status = success
08:32:51 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-41d6ea84-775 9-4cf6-99fc-c1b35b82256d'.
08:32:52 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:32:55 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:33:16 UTC Batch complete.
step = 20
J = 4.715e-01
grad_norm = 2.6598e-0108:33:17 UTC Created task 'GC_invdes' with resource_id 'fdve-ab0bc3dd-1328-45bf-aa91-59bb7321df20' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-ab0bc3dd-132 8-45bf-aa91-59bb7321df20'.
Task folder: 'default'.
08:33:22 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:33:41 UTC starting up solver
running solver
08:33:47 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:33:54 UTC status = success
08:33:56 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-ab0bc3dd-132 8-45bf-aa91-59bb7321df20'.
08:33:57 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:34:02 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:34:25 UTC Batch complete.
step = 21
J = 4.740e-01
grad_norm = 2.0068e-0108:34:26 UTC Created task 'GC_invdes' with resource_id 'fdve-dd2d544a-f0dd-4bac-abfd-bf645ec1af36' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dd2d544a-f0d d-4bac-abfd-bf645ec1af36'.
Task folder: 'default'.
08:34:29 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:34:41 UTC starting up solver
running solver
08:34:48 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:34:55 UTC status = success
08:34:57 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-dd2d544a-f0d d-4bac-abfd-bf645ec1af36'.
08:34:58 UTC Loading results from simulation_data.hdf5
08:34:59 UTC Started working on Batch containing 1 tasks.
08:35:02 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:35:22 UTC Batch complete.
step = 22
J = 4.742e-01
grad_norm = 2.9551e-0108:35:24 UTC Created task 'GC_invdes' with resource_id 'fdve-d3fb4fbf-4fcf-4f4c-affb-a8777e15ae58' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-d3fb4fbf-4fc f-4f4c-affb-a8777e15ae58'.
Task folder: 'default'.
08:35:28 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:35:40 UTC status = preprocess
08:35:44 UTC starting up solver
running solver
08:35:50 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:35:57 UTC status = success
08:35:59 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-d3fb4fbf-4fc f-4f4c-affb-a8777e15ae58'.
08:36:00 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:36:03 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:36:24 UTC Batch complete.
step = 23
J = 4.774e-01
grad_norm = 2.1671e-0108:36:25 UTC Created task 'GC_invdes' with resource_id 'fdve-646b8467-6092-444d-b59e-c3809d62466b' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-646b8467-609 2-444d-b59e-c3809d62466b'.
Task folder: 'default'.
08:36:29 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:36:41 UTC status = preprocess
08:36:47 UTC starting up solver
08:36:48 UTC running solver
08:36:51 UTC early shutoff detected at 42%, exiting.
08:36:52 UTC status = postprocess
08:37:01 UTC status = success
08:37:03 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-646b8467-609 2-444d-b59e-c3809d62466b'.
08:37:04 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:37:08 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:37:28 UTC Batch complete.
step = 24
J = 4.795e-01
grad_norm = 2.2187e-0108:37:29 UTC Created task 'GC_invdes' with resource_id 'fdve-860cefe7-af39-40fe-ae3c-809be43f559f' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-860cefe7-af3 9-40fe-ae3c-809be43f559f'.
Task folder: 'default'.
08:37:34 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:37:48 UTC starting up solver
running solver
08:37:54 UTC early shutoff detected at 41%, exiting.
status = postprocess
08:38:01 UTC status = success
08:38:03 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-860cefe7-af3 9-40fe-ae3c-809be43f559f'.
08:38:04 UTC Loading results from simulation_data.hdf5
Started working on Batch containing 1 tasks.
08:38:07 UTC Maximum FlexCredit cost: 0.025 for the whole batch.
Use 'Batch.real_cost()' to get the billed FlexCredit cost after completion.
08:38:35 UTC Batch complete.
step = 25
J = 4.801e-01
grad_norm = 1.6251e-01
CPU times: user 1min 5s, sys: 6.18 s, total: 1min 11s
Wall time: 28min 40sPlot the objective function (average coupling efficiency) as a function of the number of iterations. A steady increase of the average coupling efficiency is observed.
plt.plot(J_history, c="black")
plt.xlabel("Iterations")
plt.ylabel("Objective function (average coupling efficiency)")
plt.show()
Design the Apodized GC
For comparison with more established methods, we’ll simulate a device designed using the approach introduced in the reference paper. First of all, the Bragg condition is given by
\[ p_{\text{i}} = \frac{\lambda_\text{c}}{(n_{\text{eff,i}} - n_{\text{c}}\sin \theta_{\text{c}})}, \]
where \(p_{\text{i}}\) is the local period of the grating, \(\lambda_\text{c}\) is the central wavelength, \(n_{\text{c}}\) is the refractive index of the cladding, \(\theta_{\text{c}}\) is the incident angle in the cladding (~ 10 degrees from Snell’s law), and
\[ n_{\text{eff,i}} = f_\text{i} \times n_\text{neff,nonetch} + (1-f_\text{i}) \times n_\text{neff,etch}. \]
Here \(f_\text{i}\) is the local filling fraction, \(n_\text{neff,unetch}\) is the effective index of the slab mode of the unetched waveguide (260 nm thick), and \(n_\text{neff,etch}\) is the effective index of the slab mode of the etched waveguide (100 nm thick). \(n_\text{neff,unetch}\) and \(n_\text{neff,etch}\) are determined to be about 2.9 and 2.2 by mode analysis (not shown). Furthermore, we follow a linear apodization, namely
\[ f_\text{i} = f_\text{0} - R \times x, \]
where \(f_\text{0}\) is the filling fraction of the first period, \(R\) is the linear apodization factor and \(x\) is the position of each tooth from the starting point of the grating. By using an initial local filling fraction \(f_\text{0}\)=0.9 and apodization factor \(R\)=0.0025 µm\(^{-1}\) we can recursively calculate \(p_{\text{i}}\) and \(f_\text{i}\) for \(i=0, 1, ..., N-1\) using the equations above. The result is a grating with apodized pitch length and filling fraction.
For more details on the design principles, refer to the reference paper.
# and then the apodized design parameters to compare to later
R = 0.025 # linear apodization factor
p_list_apodized = [] # list to store all local periodicities
f_list_apodized = [] # list to store all local filling fractions
f0_apodized = 0.86 # filling fraction of the first teeth
f_current = f0_apodized # variable to store the current local filling fraction
# recursively calculate all local periodicities and filling fractions
for i in range(N):
p_current = get_periodicity(f_current)
p_list_apodized.append(p_current)
f_list_apodized.append(f_current)
f_current = f0_apodized - R * sum(p_list_apodized)
widths_apodized = get_widths(p_list_apodized, f_list_apodized)
params_apodized = inverse_project(widths_apodized)
sim_apodized = make_2d_sim(params_apodized)Let’s plot all of the simulations together to see how their features compare.
params_opt = params_history[-1]
sim_opt = make_2d_sim(params_opt)
f, (ax1, ax2, ax3) = plt.subplots(3, 1, tight_layout=True, figsize=(10, 6))
sim0.plot_eps(y=0, freq=freq0, ax=ax1)
sim_apodized.plot_eps(y=0, freq=freq0, ax=ax2)
sim_opt.plot_eps(y=0, freq=freq0, ax=ax3)
ax1.set_title("uniform")
ax2.set_title("apodized")
ax3.set_title("optimized")
plt.show()
Next, we’ll run each of the designs to be able to compare the coupling efficiencies directly.
sim_data_apodized = web.run(simulation=sim_apodized, task_name="apodized design")
sim_data_opt = web.run(simulation=sim_opt, task_name="optimized design")08:38:55 UTC Created task 'apodized design' with resource_id 'fdve-61e70a14-4550-4d5b-874d-4db56cda4d06' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-61e70a14-455 0-4d5b-874d-4db56cda4d06'.
Task folder: 'default'.
08:38:57 UTC Estimated FlexCredit cost: 0.025. This assumes the FDTD solver runs for the full simulation time; if early shutoff is reached, the billed cost can be lower. Use 'web.real_cost(task_id)' to get the billed FlexCredit cost after a simulation run.
08:38:58 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:39:12 UTC status = preprocess
08:39:16 UTC starting up solver
08:39:17 UTC running solver
08:39:20 UTC early shutoff detected at 42%, exiting.
08:39:21 UTC status = postprocess
status = success
08:39:23 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-61e70a14-455 0-4d5b-874d-4db56cda4d06'.
08:39:24 UTC Loading results from simulation_data.hdf5
Created task 'optimized design' with resource_id 'fdve-b6a3e1e8-7c1e-4e22-bff7-90480a190c53' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-b6a3e1e8-7c1 e-4e22-bff7-90480a190c53'.
Task folder: 'default'.
08:39:25 UTC Estimated FlexCredit cost: 0.025. This assumes the FDTD solver runs for the full simulation time; if early shutoff is reached, the billed cost can be lower. Use 'web.real_cost(task_id)' to get the billed FlexCredit cost after a simulation run.
08:39:31 UTC status = queued
To cancel the simulation, use 'web.abort(task_id)' or 'web.delete(task_id)' or abort/delete the task in the web UI. Terminating the Python script will not stop the job running on the cloud.
08:39:43 UTC starting up solver
running solver
08:39:48 UTC early shutoff detected at 42%, exiting.
status = postprocess
08:39:51 UTC status = success
08:39:53 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-b6a3e1e8-7c1 e-4e22-bff7-90480a190c53'.
08:39:54 UTC Loading results from simulation_data.hdf5
ce_opt = compute_ce(sim_data_opt)
ce_apodized = compute_ce(sim_data_apodized)
plt.plot(ldas, dB(ce0), c="red", linewidth=2, label="Initial design")
plt.plot(ldas, dB(ce_apodized), c="orange", linewidth=2, label="Apodized design")
plt.plot(ldas, dB(ce_opt), c="blue", linewidth=2, label="Optimized design")
plt.xlim(min(ldas), max(ldas))
plt.ylim(-10, 0)
plt.xlabel("Wavelength (µm)")
plt.ylabel("Coupling efficiency (dB)")
plt.legend()
plt.grid()
plt.show()
We see that the optimized and apodized designs both exceed that of the uniform with -2dB efficiency. However, the optimized design far outperforms in terms of bandwidth. We can also confirm that the minimal feature size is maintained above 85 nm.
print("Uniform:")
print(f" The 1dB bandwidth is {bandwidth(ldas, ce0):.1f} nm")
print(f" The maximum coupling efficiency is {np.max(dB(ce0)):.2f} dB.")
print(f" The minimal feature size is {1e3 * np.min(project(params0)):.2f} nm.")
print("")
print("Apodized:")
print(f" The 1dB bandwidth is {bandwidth(ldas, ce_apodized):.1f} nm")
print(f" The maximum coupling efficiency is {np.max(dB(ce_apodized)):.2f} dB.")
print(f" The minimal feature size is {1e3 * np.min(project(params_apodized)):.2f} nm.")
print("")
print("Optimized:")
print(f" The 1dB bandwidth is {bandwidth(ldas, ce_opt):.1f} nm")
print(f" The maximum coupling efficiency is {np.max(dB(ce_opt)):.2f} dB.")
print(f" The minimal feature size is {1e3 * np.min(project(params_opt)):.2f} nm.")Uniform:
The 1dB bandwidth is 36.0 nm
The maximum coupling efficiency is -2.74 dB.
The minimal feature size is 86.12 nm.
Apodized:
The 1dB bandwidth is 43.0 nm
The maximum coupling efficiency is -1.99 dB.
The minimal feature size is 85.00 nm.
Optimized:
The 1dB bandwidth is 50.0 nm
The maximum coupling efficiency is -1.99 dB.
The minimal feature size is 86.12 nm.Take the Model Further
In this notebook, we only perform 2D simulations. A good next step would be to convert the 2D design to a 3D linear GC or focusing GC. The coupling efficiency of the 3D GC is likely slightly lower than that in 2D. Inverse design can then be applied to fine-tune it again. To minimize the footprint, one can potentially also apply shape optimization to design a compact and low-loss taper section, as demonstrated in this example.
Further visualizations
In case you would like to gain more insight into how your device has changed over optimization, the following are a few post-processing steps to make animations and visualize the grating teeth locations.
import matplotlib.animation as animation
from IPython.display import HTML
fig, ax = plt.subplots(1, 1, tight_layout=True, figsize=(10, 3))
sim0 = make_2d_sim(params_history[0])
def plot_fn(params, ax=ax):
sim = make_2d_sim(params)
ax = sim.plot_eps(y=0.0, ax=ax, freq=freq0)
ax.set_aspect("equal")
def animate(i):
plot_fn(params_history[i])
# create animation
ani = animation.FuncAnimation(fig, animate, frames=len(params_history))
plt.close()td.config.logging.level = "ERROR"
# display the animation (press "play" to start)
HTML(ani.to_jshtml())<Figure size 672x480 with 0 Axes># uncomment to save to file
# ani.save("/users/twhughes/Desktop/apodized_grating_coupler.gif", fps=10)history_min = [[] for _ in range(len(params_history))]
history_max = [[] for _ in range(len(params_history))]
for i, params in enumerate(params_history):
sim = make_2d_sim(params)
for geometry in sim.structures[2].geometry.geometries:
(xmin, _, _), (xmax, _, _) = geometry.bounds
history_min[i].append(xmin)
history_max[i].append(xmax)
# print(i, j)
history_min = np.array(history_min).T
history_max = np.array(history_max).Tfrom matplotlib.pyplot import fill_between
xs = np.arange(len(params_history))
_, ax = plt.subplots()
for y1, y2 in zip(history_min, history_max):
fill_between(xs, y1, y2)
plt.xlabel("iteration number")
plt.ylabel("x position")
plt.title("tooth position and width over optimization")
plt.show()
TERMS & CONDITIONS
How the process works:
Try this real URL to see an example: