# The Bayesian optimizer uses the bayesian-optimization external package version 1.5.1.
# Uncomment the following line to install the package
# pip install bayesian-optimization==1.5.1
import gdstk
import matplotlib.pyplot as plt
import numpy as np
import tidy3d as td
import tidy3d.plugins.design as tdd
import tidy3d.web as web
from scipy.interpolate import make_interp_spline
TIDY3D
LEARNING CENTER
Example Library
Bayesian optimization of a Y branch¶
Note: the cost of running the entire notebook is larger than 10 FlexCredits.
Bayesian optimization is a popular technique for searching and optimizing a design space. It works by building a probabilistic model, such as a Gaussian process, to predict the outcome of a more complex objective function. Unlike other optimizers like genetic algorithms or particle swarm optimization, Bayesian optimization uses an acquisition function to select new potential solutions, focusing on areas of the design space that the Gaussian process predicts will have high objective values and high uncertainty. Over the course of the optimization, the model’s uncertainty reduces, improving its accuracy as a surrogate for the true objective function. This makes Bayesian optimization a direct approach to optimization, particularly well-suited to problems with small (< 30 dimensions) and expensive-to-evaluate objective functions.
However, as Bayesian optimization is an iterative process with many individual components, it can be challenging to design and parallelize, leading to long run times and high computational costs. Fortunately, the Design plugin in Tidy3D includes the Bayesian optimization method MethodBayOpt which eliminates much of this complexity. Users can quickly and easily run a Bayesian optimization for complex FDTD simulations and analyze the results. Further details on Bayesian optimization and the methods used in Tidy3D are available in the open-source Python library bayesian_optimization.
This notebook details how to perform a Bayesian optimization with the Tidy3D Design plugin for the development of a Y-junction. It is based on the work of Zhengqi Gao, Zhengxing Zhang, and Duane S. Boning, "Automatic Synthesis of Broadband Silicon Photonic Devices via Bayesian Optimization" Journal of Lightwave Technology 40, 7879-7892 (2022) DOI:10.1109/JLT.2022.3207052. A detailed description of how to build a Y-junction can be found in the Waveguide Y junction.
If you are curious about other Design plugin features, see these notebooks:
Simulation Setup
The simulation is defined across the 1.5 \(\mu m\) to 1.6 \(\mu m\) wavelength range with 100 sampling points. The base of the model is silicon, the top of the model is silicon dioxide; we use the material constants included in the Tidy3D Tidy3D’s material library.
lda0 = 1.55 # central wavelength
freq0 = td.C_0 / lda0 # central frequency
n_wav = 100 # Number of wavelengths to sample in the range
ldas = np.linspace(1.5, 1.6, n_wav) # wavelength range
freqs = td.C_0 / ldas # frequency range
fwidth = 0.5 * (np.max(freqs) - np.min(freqs)) # width of the source frequency range
si = td.material_library["cSi"]["Palik_LowLoss"]
sio2 = td.material_library["SiO2"]["Palik_LowLoss"]The following parameters describe the geometry of the waveguide. The optimizable design space is the junction which is split into 13 discrete segments.
t = 0.22 # thickness of the silicon layer
num_d = 13 # dimensional space of the design region
l_in = 1 # input waveguide length
l_junction = 2 # length of the junction
l_bend = 6 # horizontal length of the waveguide bend
h_bend = 2 # vertical offset of the waveguide bend
l_out = 1 # output waveguide length
branch_width = 0.5 # width of one Y branch
branch_sep = 0.2 # distance between y branches at the junction
inf_eff = 100 # effective infinityThe most effective way to use the Design plugin is to split the workflow into “pre” and “post” functions which can surround a call to the Tidy3D cloud that carries out the computation. This takes advantage of automated simulation batching, allowing for parallelization which saves a considerable amount of time. The pre-function returns a Simulation; the post-function then analyzes the corresponding SimulationData. Together they can be considered the fitness function (or objective function / figure of merit) which the Bayesian optimization process is working to predict.
The function fn_pre needs to take the parameters that are being optimized: these are the 13 width segments of the junction. They always input as a dictionary and can be unpacked within the function (as below) or included as keyword arguments. These parameters are used to build a Y-junction Simulation object.
The function fn_post then takes the SimulationData object output by the simulation and computes the objective function. The value of objective is then fed back into the Bayesian optimization to inform the probabilistic model. In this case, we extract the power passing through the Y-junction and the power reflected back to the source. These values are evaluated in a custom loss function described by Gao et al. to determine the effectiveness of the junction design:
\[ -\frac{N_\lambda}{3} \sum_\lambda \left( R^2 + 2(T-0.5)^2 \right) \]
where the summation is performed over the simulated wavelengths, \(N_\lambda\) is the number of wavelength points, \(R\) is the reflected power and \(T\) is the transmitted power in one branch.
The aim of this function is to achieve a transmitted power of 0.5 through each branch, whilst driving the power reflected towards the source to zero. Note there is a minus sign in front of the output value as this loss function is a minimizing function, whilst all the optimizers in the Design plugin are built to maximize the objective function.
def fn_pre(**w_params: dict) -> td.Simulation:
"""Create a Simulation of a Y splitter from a series of junction widths.
Includes mode monitors to measure the power transmitted and reflected to source.
"""
w_start = 0.5
w_end = branch_width * 2 + branch_sep
widths = [w_start] # Ensures input waveguide is included in spline for first point of junction
widths.extend(list(w_params.values()))
widths.append(w_end) # Ensures final point of junction smoothly converts to the branches
x_junction = np.linspace(
l_in, l_in + l_junction, num_d + 2
) # x coordinates of the top edge vertices
y_junction = np.array(widths) # y coordinates of the top edge vertices
# pass vertices through spline and increase sampling to smooth the geometry
new_x_junction = np.linspace(
l_in, l_in + l_junction, 100
) # x coordinates of the top edge vertices
spline = make_interp_spline(x_junction, y_junction, k=2)
spline_yjunction = spline(new_x_junction)
# using concatenate to include bottom edge vertices
x_junction = np.concatenate((new_x_junction, np.flipud(new_x_junction)))
y_junction = np.concatenate((spline_yjunction / 2, -np.flipud(spline_yjunction / 2)))
# stacking x and y coordinates to form vertices pairs
vertices = np.transpose(np.vstack((x_junction, y_junction)))
junction = td.Structure(
geometry=td.PolySlab(vertices=vertices, axis=2, slab_bounds=(0, t)), medium=si
)
x_start = l_in + l_junction # x coordinate of the starting point of the waveguide bends
x_bend = np.linspace(x_start, x_start + l_bend, 100) # x coordinates of the top edge vertices
y_bend = (
(x_bend - x_start) * h_bend / l_bend
- h_bend * np.sin(2 * np.pi * (x_bend - x_start) / l_bend) / (np.pi * 2)
+ w_end / 2
- w_start / 2
) # y coordinates of the top edge vertices
# adding the last point to include the straight waveguide at the output
x_bend = np.append(x_bend, inf_eff)
y_bend = np.append(y_bend, y_bend[-1])
# add path to the cell
cell = gdstk.Cell("bends")
cell.add(
gdstk.FlexPath(x_bend + 1j * y_bend, branch_width, layer=1, datatype=0)
) # top waveguide bend
cell.add(
gdstk.FlexPath(x_bend - 1j * y_bend, branch_width, layer=1, datatype=0)
) # bottom waveguide bend
# define top waveguide bend structure
wg_bend_1 = td.Structure(
geometry=td.PolySlab.from_gds(
cell,
gds_layer=1,
axis=2,
slab_bounds=(0, t),
)[0],
medium=si,
)
# define bottom waveguide bend structure
wg_bend_2 = td.Structure(
geometry=td.PolySlab.from_gds(
cell,
gds_layer=1,
axis=2,
slab_bounds=(0, t),
)[1],
medium=si,
)
# straight input waveguide
wg_in = td.Structure(
geometry=td.Box.from_bounds(rmin=(-inf_eff, -w_start / 2, 0), rmax=(l_in, w_start / 2, t)),
medium=si,
)
# the entire model is the collection of all structures defined so far
model_structure = [wg_in, junction, wg_bend_1, wg_bend_2]
Lx = l_in + l_junction + l_out + l_bend # simulation domain size in x direction
Ly = w_end + 2 * h_bend + 1.5 * lda0 # simulation domain size in y direction
Lz = 10 * t # simulation domain size in z direction
sim_size = (Lx, Ly, Lz)
# add a mode source as excitation
mode_spec = td.ModeSpec(num_modes=1, target_neff=3.5)
mode_source = td.ModeSource(
center=(l_in / 2, 0, t / 2),
size=(0, 4 * w_start, 6 * t),
source_time=td.GaussianPulse(freq0=freq0, fwidth=fwidth),
direction="+",
mode_spec=mode_spec,
mode_index=0,
)
# add a mode monitor to measure transmission at the output waveguide
mode_monitor_11 = td.ModeMonitor(
center=(l_in / 3, 0, t / 2),
size=(0, 4 * w_start, 6 * t),
freqs=freqs,
mode_spec=mode_spec,
name="mode_11",
)
mode_monitor_12 = td.ModeMonitor(
center=(l_in + l_junction + l_bend + l_out / 2, w_end / 2 - w_start / 2 + h_bend, t / 2),
size=(0, 4 * w_start, 6 * t),
freqs=freqs,
mode_spec=mode_spec,
name="mode_12",
)
# add a field monitor to visualize field distribution at z=t/2
field_monitor = td.FieldMonitor(
center=(0, 0, t / 2), size=(td.inf, td.inf, 0), freqs=[freq0], name="field"
)
run_time = 5e-13 # simulation run time
# construct simulation
sim = td.Simulation(
center=(Lx / 2, 0, 0),
size=sim_size,
grid_spec=td.GridSpec.auto(min_steps_per_wvl=20, wavelength=lda0),
structures=model_structure,
sources=[mode_source],
monitors=[mode_monitor_11, mode_monitor_12, field_monitor],
run_time=run_time,
boundary_spec=td.BoundarySpec.all_sides(boundary=td.PML()),
medium=sio2,
)
return sim
def fn_post(sim_data: td.SimulationData) -> float:
"""Calculate the loss function from the power at the mode monitors in the SimulationData."""
# Calculate the power reflected back to source and transmitted to one branch
power_reflected = np.squeeze(
np.abs(sim_data["mode_11"].amps.sel(direction="-", mode_index=0)) ** 2
)
power_transmitted = np.squeeze(
np.abs(sim_data["mode_12"].amps.sel(direction="+", mode_index=0)) ** 2
)
# Loss function proposed by Gao et al. which takes advantage of branch symmetry
loss_fn = 1 / 3 * n_wav * np.sum(power_reflected**2 + 2 * (power_transmitted - 0.5) ** 2)
output = -float(loss_fn.values) # Negative value as this is a minimizing loss function
return outputWe can quickly check that fn_pre is working correctly by passing a set of potential test_params and plotting the result. Note that the gap between the Polyslab junction and the branches is a plotting artifact and doesn’t exist in the Simulation.
test_params = {
"w1": 0.5,
"w2": 0.5,
"w3": 0.6,
"w4": 0.7,
"w5": 0.9,
"w6": 1.26,
"w7": 1.4,
"w8": 1.4,
"w9": 1.4,
"w10": 1.4,
"w11": 1.31,
"w12": 0.5,
"w13": 0.5,
}
sim = fn_pre(**test_params)
sim.plot(z=0)
plt.show()
Next, we setup the Bayesian optimization method. This is done with the MethodBayOpt object. We don’t have the lcb (lower confidence bound) acquisition function used in the paper available to us, but by making our loss function negative and using the ucb (upper confidence bound) acquisition we achieve the same optimizer design. The initial_iter and n_iter options control the number of initial random samples and subsequent optimization iterations respectively. We can also set the random seed to ensure reproducible results (optional).
We also create a list of ParameterFloat objects corresponding to the 13 segment widths. The span option defines the bounds of each parameter.
The method and parameters are then passed to a DesignSpace object which contains all we need to run the Bayesian optimization.
method = tdd.MethodBayOpt(
initial_iter=30,
n_iter=70,
acq_func="ucb",
kappa=0.3,
seed=1,
)
parameters = [tdd.ParameterFloat(name=f"w_{i}", span=(0.5, 1.6)) for i in range(num_d)]
design_space = tdd.DesignSpace(
method=method, parameters=parameters, task_name="bay_opt_notebook", path_dir="./data"
)It is then very easy to the launch this optimization with design_space.run(). This launches an initial random batch of 30 simulations, as specified by initial_iter, followed by sequential computation of 70 simulations, as specified by n_iter. In the latter phase, the Bayesian optimizer chooses potential candidates based on simultaneously maximizing the objective value and efficiently exploring the design space. The total of 100 simulations takes around 90 minutes to compute. Once complete, the results are returned in a pandas dataframe for analysis.
results = design_space.run(fn_pre, fn_post, verbose=True)
df = results.to_dataframe()10:14:59 UTC Running 30 Simulations
| 1 | -180.6539 | 0.9587242 | 1.2923569 | 0.5001258 | 0.8325658 | 0.6614314 | 0.6015724 | 0.7048862 | 0.8801167 | 0.9364442 | 1.0926984 | 0.9611139 | 1.2537414 | 0.7248974 | | 2 | -704.4466 | 1.4659291 | 0.5301263 | 1.2375142 | 0.9590352 | 1.1145588 | 0.6544256 | 0.7179116 | 1.3808190 | 1.5650877 | 0.8447665 | 1.2615548 | 1.4640280 | 1.4840673 | | 3 | -576.6328 | 0.5935486 | 0.5429602 | 0.6868134 | 1.4659567 | 0.6081815 | 0.9632183 | 1.5536784 | 1.0864818 | 1.2610648 | 0.8470671 | 1.2551510 | 1.4180882 | 0.5201171 | | 4 | -1225.541 | 1.3251587 | 1.5877471 | 1.3229822 | 0.8084883 | 1.3682072 | 0.6135486 | 0.9926828 | 1.4994550 | 0.8229755 | 0.8165528 | 0.6430314 | 0.5213036 | 1.2467190 | | 5 | -639.9370 | 0.7327909 | 0.7921013 | 1.0407304 | 0.5586987 | 1.1315293 | 0.6614014 | 1.1482360 | 1.2697341 | 0.6125678 | 0.9554615 | 1.2638401 | 0.9555971 | 0.5549488 | | 6 | -1134.121 | 1.0894860 | 1.2301741 | 1.0663780 | 1.5390542 | 1.1452105 | 1.4937421 | 0.6512221 | 0.6532039 | 1.3881304 | 0.9374445 | 0.6818896 | 1.5202594 | 0.8825424 | | 7 | -1053.332 | 1.3258933 | 1.2985977 | 1.4716367 | 1.1860394 | 1.3260366 | 0.8837881 | 0.7969206 | 1.4854748 | 0.9709003 | 1.5613240 | 1.2297856 | 1.1838652 | 0.6262205 | | 8 | -616.3112 | 1.5444381 | 0.9949033 | 1.1362285 | 0.9489504 | 0.7607296 | 1.4937174 | 1.1310474 | 0.5031573 | 1.1788594 | 0.8593093 | 1.0797639 | 1.4745363 | 0.8929967 | | 9 | -1020.071 | 1.4993886 | 1.1856961 | 0.5174033 | 1.5223809 | 1.2599866 | 1.5970551 | 0.6895745 | 0.6508493 | 1.5258550 | 1.2664999 | 0.5726001 | 1.3310093 | 1.3292638 | | 10 | -426.7613 | 1.5153269 | 1.2826772 | 0.6366980 | 0.5218681 | 0.5288320 | 0.5311371 | 0.7708321 | 1.4460307 | 1.0927141 | 1.1081041 | 1.4262339 | 0.6365906 | 0.8071020 | | 11 | -1135.771 | 1.1443351 | 1.5665553 | 1.1171332 | 0.5205120 | 1.3806959 | 0.7562717 | 1.3878157 | 0.9266467 | 1.4498960 | 1.3218338 | 1.1118642 | 0.6501007 | 0.5659094 | | 12 | -872.5194 | 0.6334778 | 0.5490070 | 0.6182435 | 0.7482802 | 1.2842878 | 1.1156886 | 0.5138115 | 0.5791717 | 1.5640039 | 1.1249105 | 0.7236225 | 0.7775583 | 1.3182084 | | 13 | -1325.809 | 0.7149724 | 1.1394948 | 1.5670219 | 1.4315116 | 0.7638325 | 1.0431466 | 1.1819512 | 1.4118789 | 0.6724705 | 0.5204338 | 0.5770243 | 1.0349796 | 1.1669624 | | 14 | -968.4757 | 1.1257365 | 0.8490986 | 1.5874777 | 1.1377197 | 0.9181552 | 1.1060430 | 1.3198678 | 1.2361561 | 0.7914115 | 0.5729683 | 0.9070926 | 1.1926892 | 0.7311914 | | 15 | -665.7242 | 1.3280311 | 0.5731901 | 0.7863466 | 1.3852300 | 0.7127777 | 1.2034069 | 1.0771373 | 1.5172887 | 0.7896264 | 0.5725571 | 1.3085725 | 1.3493958 | 1.4985974 | | 16 | -407.8049 | 1.5251692 | 0.5153467 | 0.7577982 | 1.1784561 | 1.5439179 | 1.5451937 | 1.1123185 | 1.5071669 | 1.2057228 | 0.9290084 | 1.0345897 | 1.1647415 | 1.1045027 | | 17 | -727.7176 | 1.5187995 | 1.5106067 | 0.9343631 | 1.5595887 | 0.6913512 | 0.6389624 | 0.6485870 | 1.0562283 | 0.5236772 | 1.5427672 | 1.4098270 | 0.5165208 | 0.6938158 | | 18 | -352.4688 | 0.8652699 | 0.6440965 | 1.3904397 | 0.8792103 | 1.5341182 | 1.1402155 | 1.4667151 | 1.4292078 | 1.4959315 | 1.0058682 | 1.1009814 | 1.3784639 | 0.8142907 | | 19 | -826.0877 | 1.0392788 | 1.1590213 | 0.5170866 | 1.1528295 | 0.9770439 | 1.3880965 | 0.8467692 | 1.4821775 | 1.1356429 | 0.7024112 | 1.3667221 | 1.1732342 | 0.5593001 | | 20 | -804.4139 | 0.9622130 | 1.2469757 | 1.5104619 | 0.5004422 | 1.5744350 | 0.9142383 | 1.5711618 | 1.1651877 | 1.4117303 | 1.1321826 | 1.1908838 | 0.8141339 | 1.1455166 | | 21 | -1467.759 | 1.3250239 | 1.4441452 | 1.3305904 | 1.2678629 | 1.4509273 | 0.8549490 | 1.2378676 | 0.9959613 | 0.9203130 | 0.9518924 | 0.9416275 | 0.8491223 | 1.1841113 | | 22 | -769.7994 | 0.9732719 | 1.5711822 | 1.2455809 | 0.7184268 | 0.9693711 | 0.8776808 | 1.3774026 | 1.4679981 | 1.4942261 | 1.2289917 | 0.7972290 | 0.7776033 | 1.4403877 | | 23 | -1501.670 | 1.0804861 | 1.3823771 | 1.1297373 | 1.3064567 | 1.0709127 | 1.3479723 | 1.1257437 | 1.0122808 | 0.8769577 | 0.5750302 | 0.9157165 | 0.5875886 | 1.5810988 | | 24 | -1013.943 | 0.6997741 | 1.3930445 | 1.4624578 | 1.2572545 | 1.1264438 | 0.6770685 | 1.0135680 | 0.8796892 | 0.7475439 | 1.1517630 | 0.8434968 | 1.5079361 | 1.5005990 | | 25 | -955.6246 | 0.7828301 | 0.6219804 | 0.7122590 | 1.0495425 | 1.3014442 | 0.7290138 | 0.7728369 | 1.4368390 | 0.9574335 | 1.1783535 | 0.7570327 | 0.6121639 | 1.0674427 | | 26 | -851.0590 | 1.0248550 | 0.6679388 | 1.1839868 | 1.0984111 | 1.2195510 | 0.6590000 | 1.3266805 | 0.7442540 | 1.0712870 | 1.3638256 | 0.5245634 | 0.8567987 | 1.4602146 | | 27 | -944.8509 | 1.4291805 | 1.0922846 | 1.4532691 | 1.5447865 | 1.4090476 | 1.4395269 | 0.6086177 | 1.2164347 | 1.2738686 | 1.1712648 | 1.3795767 | 0.5380283 | 1.3472626 | | 28 | -744.9340 | 1.3049014 | 0.7856682 | 0.7827762 | 1.1955336 | 0.8798272 | 1.3762475 | 0.9907608 | 1.3610243 | 1.5895189 | 0.8302731 | 0.6573064 | 1.4914392 | 1.0957153 | | 29 | -903.7163 | 1.5722144 | 1.2002648 | 1.5933043 | 1.1006778 | 1.0790685 | 0.6489706 | 0.8912756 | 0.5288404 | 0.6764346 | 1.3202009 | 0.5334396 | 0.9031974 | 1.4485808 | | 30 | -785.8112 | 1.2619454 | 1.2600363 | 0.7075004 | 0.9860947 | 1.1397351 | 1.5887268 | 0.7242968 | 0.7725061 | 0.7883903 | 1.3251896 | 1.0026728 | 0.5626223 | 1.0593678 |
10:17:43 UTC Best Fit from Initial Solutions: -180.654
10:17:44 UTC Running 1 Simulations
| 31 | -190.0621 | 1.0522683 | 1.1760757 | 0.5 | 0.6626245 | 0.5 | 0.5432270 | 0.7165272 | 0.8747580 | 0.8993268 | 1.1667382 | 1.0827083 | 1.3027878 | 0.7196552 |
10:18:22 UTC Running 1 Simulations
| 32 | -232.3610 | 0.9338372 | 1.3319437 | 0.5 | 0.7773070 | 0.5031274 | 0.5139294 | 0.6574052 | 0.8545373 | 0.9144821 | 1.1169142 | 0.9583798 | 1.2707503 | 0.6817543 |
10:20:36 UTC Running 1 Simulations
| 33 | -150.2680 | 1.0677360 | 1.1090421 | 0.5 | 0.7270105 | 0.7019969 | 0.6587111 | 0.8049887 | 0.8681509 | 0.9182693 | 1.1306363 | 1.0459074 | 1.3692381 | 0.7706255 |
10:21:40 UTC Latest Best Fit on Iter 2: -150.268
Running 1 Simulations
| 34 | -134.7010 | 1.1017868 | 1.1594999 | 0.5 | 0.7326107 | 0.6980480 | 0.6564651 | 0.7027951 | 0.8535271 | 0.8751159 | 1.1578098 | 1.0796527 | 1.2857289 | 0.7910355 |
10:22:15 UTC Latest Best Fit on Iter 3: -134.701
10:22:16 UTC Running 1 Simulations
| 35 | -122.1144 | 1.0926640 | 1.2030089 | 0.5 | 0.6711740 | 0.7161665 | 0.6823690 | 0.6843437 | 0.8785732 | 0.8405495 | 1.1999396 | 1.0666363 | 1.3098749 | 0.8441808 |
10:23:27 UTC Latest Best Fit on Iter 4: -122.114
10:23:28 UTC Running 1 Simulations
| 36 | -157.0558 | 1.1543354 | 1.2041388 | 0.5 | 0.6124372 | 0.7226210 | 0.7050298 | 0.6679719 | 0.8953800 | 0.7858015 | 1.2018185 | 0.9963484 | 1.3166285 | 0.8354802 |
10:25:11 UTC Running 1 Simulations
| 37 | -80.47734 | 0.9582419 | 1.2301734 | 0.5 | 0.6773326 | 0.7216479 | 0.6689466 | 0.7065752 | 0.8517428 | 0.8915236 | 1.2929234 | 1.2622954 | 1.3626252 | 0.9695925 |
10:25:59 UTC Latest Best Fit on Iter 6: -80.477
10:26:00 UTC Running 1 Simulations
| 38 | -89.34242 | 0.9805040 | 1.2104089 | 0.5 | 0.6967425 | 0.6724455 | 0.7139549 | 0.7431892 | 0.8552651 | 0.8887100 | 1.3662391 | 1.2515370 | 1.3332569 | 0.9972348 |
10:26:56 UTC Running 1 Simulations
| 39 | -92.97644 | 0.9979875 | 1.2509106 | 0.5 | 0.6867444 | 0.7439457 | 0.5848105 | 0.6984067 | 0.8677333 | 0.9052262 | 1.2938998 | 1.2515848 | 1.3731695 | 1.0098429 |
10:27:35 UTC Running 1 Simulations
| 40 | -97.82107 | 0.8909229 | 1.2092209 | 0.5 | 0.6733788 | 0.7713976 | 0.6948590 | 0.6624629 | 0.8054665 | 0.8426538 | 1.3202261 | 1.2669906 | 1.3867445 | 0.8938276 |
10:28:36 UTC Running 1 Simulations
| 41 | -107.5207 | 0.9663897 | 1.2456090 | 0.5 | 0.5677728 | 0.6675065 | 0.7541698 | 0.7230917 | 0.8274898 | 0.9117789 | 1.1839259 | 1.3331353 | 1.3405449 | 1.0382272 |
10:29:37 UTC Running 1 Simulations
| 42 | -74.34442 | 0.9391022 | 1.2475747 | 0.5 | 0.7154334 | 0.6920238 | 0.7053409 | 0.7041178 | 0.9852475 | 0.9201795 | 1.2659107 | 1.3040579 | 1.3952222 | 0.9313445 |
10:30:27 UTC Latest Best Fit on Iter 11: -74.344
10:30:29 UTC Running 1 Simulations
| 43 | -68.42979 | 1.1607788 | 0.5 | 1.0843618 | 0.5063023 | 1.5289877 | 1.5503850 | 1.6 | 1.6 | 1.4759911 | 1.2081043 | 1.3267847 | 1.6 | 1.0916108 |
10:31:51 UTC Latest Best Fit on Iter 12: -68.43
10:31:52 UTC Running 1 Simulations
| 44 | -63.62362 | 1.3714582 | 0.5 | 0.9013382 | 0.5386194 | 1.5035745 | 1.6 | 1.6 | 1.6 | 1.3508395 | 1.2482270 | 1.3576625 | 1.6 | 1.1359782 |
10:32:39 UTC Latest Best Fit on Iter 13: -63.624
10:32:40 UTC Running 1 Simulations
| 45 | -85.05320 | 1.4453884 | 0.5 | 1.0564988 | 0.5393353 | 1.6 | 1.6 | 1.6 | 1.6 | 1.4686618 | 1.1611994 | 1.3454274 | 1.6 | 1.0969831 |
10:34:04 UTC Running 1 Simulations
| 46 | -97.95916 | 1.2362564 | 0.5 | 0.9295396 | 0.5 | 1.2925307 | 1.5185495 | 1.6 | 1.6 | 1.3808034 | 1.3894596 | 1.4229902 | 1.6 | 1.1928763 |
10:35:40 UTC Running 1 Simulations
| 47 | -48.36476 | 1.1898082 | 0.5 | 0.7519351 | 0.5 | 1.6 | 1.6 | 1.6 | 1.6 | 1.2780639 | 1.1684355 | 1.2521264 | 1.6 | 1.1016447 |
10:36:59 UTC Latest Best Fit on Iter 16: -48.365
10:37:01 UTC Running 1 Simulations
| 48 | -54.01582 | 1.2188875 | 0.5 | 0.9027782 | 0.5 | 1.6 | 1.6 | 1.6 | 1.6 | 1.2074730 | 1.2610002 | 1.1841865 | 1.6 | 1.1143966 |
10:37:53 UTC Running 1 Simulations
| 49 | -63.57638 | 1.1762434 | 0.5 | 0.8516728 | 0.5 | 1.6 | 1.6 | 1.6 | 1.6 | 1.1978430 | 1.1600914 | 1.4094348 | 1.6 | 1.1271603 |
10:38:47 UTC Running 1 Simulations
| 50 | -48.50647 | 1.2100349 | 0.5 | 0.7693759 | 0.5 | 1.6 | 1.6 | 1.6 | 1.6 | 1.3515782 | 1.3061366 | 1.2489936 | 1.6 | 0.9930843 |
10:39:40 UTC Running 1 Simulations
| 51 | -43.50288 | 1.2029518 | 0.5 | 0.7543247 | 0.5477109 | 1.6 | 1.6 | 1.6 | 1.6 | 1.3822642 | 1.3182152 | 1.2423992 | 1.6 | 1.1872446 |
10:40:51 UTC Latest Best Fit on Iter 20: -43.503
10:40:52 UTC Running 1 Simulations
| 52 | -52.05282 | 1.3659270 | 0.5 | 0.5974036 | 0.5 | 1.6 | 1.6 | 1.6 | 1.6 | 1.2667933 | 1.3228177 | 1.2378488 | 1.6 | 1.1850347 |
10:41:43 UTC Running 1 Simulations
| 53 | -37.89904 | 1.2279095 | 0.5 | 0.7487606 | 0.6572372 | 1.6 | 1.6 | 1.6 | 1.6 | 1.3368683 | 1.2829257 | 1.2466565 | 1.6 | 1.1085687 |
10:44:00 UTC Latest Best Fit on Iter 22: -37.899
10:44:01 UTC Running 1 Simulations
| 54 | -38.96237 | 1.2143083 | 0.5 | 0.7456376 | 0.6365094 | 1.6 | 1.4982227 | 1.6 | 1.6 | 1.3458085 | 1.2750961 | 1.2490804 | 1.6 | 1.1058533 |
10:44:33 UTC Running 1 Simulations
| 55 | -43.64383 | 1.2593036 | 0.5 | 0.7509741 | 0.6500018 | 1.6 | 1.5762487 | 1.6 | 1.6 | 1.3775599 | 1.2617909 | 1.1994043 | 1.6 | 1.1026213 |
10:45:27 UTC Running 1 Simulations
| 56 | -43.31062 | 1.1709196 | 0.5 | 0.6267335 | 0.6807331 | 1.6 | 1.5580424 | 1.6 | 1.6 | 1.2603170 | 1.3529759 | 1.3250584 | 1.6 | 1.1229842 |
10:46:24 UTC Running 1 Simulations
| 57 | -35.31533 | 1.1932773 | 0.5 | 0.7782393 | 0.6450345 | 1.6 | 1.5676684 | 1.6 | 1.5102760 | 1.3124962 | 1.3048318 | 1.3062498 | 1.6 | 1.1078715 |
10:47:04 UTC Latest Best Fit on Iter 26: -35.315
10:47:06 UTC Running 1 Simulations
| 58 | -31.26061 | 1.1179085 | 0.5 | 0.8340106 | 0.6695061 | 1.6 | 1.5574909 | 1.6 | 1.5649421 | 1.3573147 | 1.3114341 | 1.3145082 | 1.6 | 1.0969044 |
10:47:36 UTC Latest Best Fit on Iter 27: -31.261
10:47:37 UTC Running 1 Simulations
| 59 | -20.36655 | 0.9379617 | 0.5 | 0.9076883 | 0.6829084 | 1.6 | 1.5161154 | 1.6 | 1.5749746 | 1.4409073 | 1.3372378 | 1.3232654 | 1.6 | 1.0765260 |
10:48:08 UTC Latest Best Fit on Iter 28: -20.367
10:48:09 UTC Running 1 Simulations
| 60 | -12.95504 | 0.5772258 | 0.5 | 1.0395362 | 0.6300330 | 1.6 | 1.4315432 | 1.6 | 1.6 | 1.6 | 1.4140505 | 1.3491333 | 1.6 | 1.0763863 |
10:49:02 UTC Latest Best Fit on Iter 29: -12.955
10:49:03 UTC Running 1 Simulations
| 61 | -11.86397 | 0.6839974 | 0.5 | 1.0648263 | 0.7190414 | 1.6 | 1.5441775 | 1.6 | 1.6 | 1.6 | 1.4357676 | 1.3597450 | 1.6 | 1.0802142 |
10:49:44 UTC Latest Best Fit on Iter 30: -11.864
10:49:45 UTC Running 1 Simulations
| 62 | -15.25406 | 0.6941479 | 0.5 | 0.9799304 | 0.7098980 | 1.6 | 1.4477466 | 1.6 | 1.6 | 1.6 | 1.3823878 | 1.4224568 | 1.6 | 1.0735731 |
10:51:36 UTC Running 1 Simulations
| 63 | -12.66334 | 0.6392819 | 0.5 | 1.0083912 | 0.6571356 | 1.6 | 1.5375297 | 1.6 | 1.6 | 1.6 | 1.3603541 | 1.2899690 | 1.6 | 1.0325465 |
10:52:09 UTC Running 1 Simulations
| 64 | -13.28191 | 0.6944038 | 0.5 | 1.0255076 | 0.6730651 | 1.6 | 1.4555691 | 1.6 | 1.6 | 1.6 | 1.4367655 | 1.2704071 | 1.6 | 1.0783252 |
10:52:56 UTC Running 1 Simulations
| 65 | -10.95505 | 0.6485171 | 0.5 | 1.0478013 | 0.6763221 | 1.6 | 1.5112866 | 1.6 | 1.6 | 1.5574590 | 1.3707955 | 1.3151897 | 1.6 | 1.1418639 |
10:53:32 UTC Latest Best Fit on Iter 34: -10.955
10:53:33 UTC Running 1 Simulations
| 66 | -9.825168 | 0.6574209 | 0.5 | 1.0329899 | 0.6571028 | 1.6 | 1.5119578 | 1.6 | 1.5043804 | 1.5614290 | 1.3995650 | 1.3324339 | 1.6 | 1.0918041 |
10:55:23 UTC Latest Best Fit on Iter 35: -9.825
10:55:24 UTC Running 1 Simulations
| 67 | -8.158044 | 0.6499813 | 0.5 | 1.0341994 | 0.6846829 | 1.6 | 1.4993326 | 1.6 | 1.5548008 | 1.5084516 | 1.4020800 | 1.3204288 | 1.6 | 1.0614925 |
10:56:16 UTC Latest Best Fit on Iter 36: -8.158
10:56:17 UTC Running 1 Simulations
| 68 | -8.022129 | 0.6262763 | 0.5 | 0.9793887 | 0.6862959 | 1.6 | 1.5214481 | 1.6 | 1.5571684 | 1.5106994 | 1.4405897 | 1.3215476 | 1.6 | 1.0878739 |
10:56:52 UTC Latest Best Fit on Iter 37: -8.022
10:56:53 UTC Running 1 Simulations
| 69 | -5.884046 | 0.6019348 | 0.5 | 1.0083738 | 0.7255650 | 1.6 | 1.5008131 | 1.6 | 1.5390440 | 1.5348409 | 1.4184816 | 1.3175751 | 1.6 | 1.0788607 |
10:57:24 UTC Latest Best Fit on Iter 38: -5.884
10:57:26 UTC Running 1 Simulations
| 70 | -4.947835 | 0.5105898 | 0.5 | 1.0554180 | 0.7493046 | 1.6 | 1.5134779 | 1.6 | 1.5309070 | 1.5647630 | 1.4487672 | 1.3164432 | 1.6 | 1.0829239 |
10:58:17 UTC Latest Best Fit on Iter 39: -4.948
10:58:18 UTC Running 1 Simulations
| 71 | -5.179303 | 0.5104270 | 0.5 | 1.0634768 | 0.7342946 | 1.6 | 1.5158677 | 1.6 | 1.5286429 | 1.5532518 | 1.4545660 | 1.2794042 | 1.6 | 1.0838874 |
10:59:12 UTC Running 1 Simulations
| 72 | -4.963722 | 0.5 | 0.5 | 1.0345226 | 0.7326799 | 1.6 | 1.4868094 | 1.6 | 1.5082944 | 1.5519470 | 1.4432093 | 1.3085739 | 1.6 | 1.0764912 |
11:00:04 UTC Running 1 Simulations
| 73 | -4.696639 | 0.5 | 0.5 | 1.0554169 | 0.7488442 | 1.6 | 1.4839414 | 1.6 | 1.5169158 | 1.5607067 | 1.4732963 | 1.3140500 | 1.6 | 1.1024051 |
11:01:28 UTC Latest Best Fit on Iter 42: -4.697
11:01:29 UTC Running 1 Simulations
| 74 | -4.437693 | 0.5082366 | 0.5 | 1.0609732 | 0.7533314 | 1.6 | 1.5073188 | 1.6 | 1.4894295 | 1.5458410 | 1.4807656 | 1.3230833 | 1.6 | 1.0923043 |
11:02:09 UTC Latest Best Fit on Iter 43: -4.438
11:02:10 UTC Running 1 Simulations
| 75 | -4.961593 | 0.5 | 0.5 | 1.0556458 | 0.7578543 | 1.6 | 1.5020699 | 1.6 | 1.4881595 | 1.5654444 | 1.5016678 | 1.3196584 | 1.6 | 1.0704190 |
11:03:56 UTC Running 1 Simulations
| 76 | -4.321766 | 0.5 | 0.5 | 1.0495750 | 0.7582289 | 1.6 | 1.5223067 | 1.6 | 1.4830800 | 1.5582939 | 1.4825048 | 1.3208318 | 1.6 | 1.1154044 |
11:05:17 UTC Latest Best Fit on Iter 45: -4.322
11:05:18 UTC Running 1 Simulations
| 77 | -4.339920 | 0.5 | 0.5 | 1.0808272 | 0.7552558 | 1.6 | 1.5149165 | 1.6 | 1.4739244 | 1.5737245 | 1.4744957 | 1.3238965 | 1.6 | 1.1124266 |
11:06:23 UTC Running 1 Simulations
| 78 | -5.075434 | 0.5 | 0.5 | 1.0692195 | 0.7373111 | 1.6 | 1.5217397 | 1.6 | 1.4836556 | 1.5631093 | 1.4861041 | 1.3282006 | 1.6 | 1.1144092 |
11:07:19 UTC Running 1 Simulations
| 79 | -2.838741 | 0.5 | 0.5 | 1.0826137 | 0.8206927 | 1.6 | 1.5293956 | 1.6 | 1.4729495 | 1.5655795 | 1.4827757 | 1.3203595 | 1.6 | 1.1150366 |
11:08:01 UTC Latest Best Fit on Iter 48: -2.839
11:08:02 UTC Running 1 Simulations
| 80 | -5.021302 | 0.5 | 0.5 | 1.1780496 | 0.9193464 | 1.6 | 1.6 | 1.6 | 1.4804456 | 1.6 | 1.5545033 | 1.3650320 | 1.6 | 1.1610276 |
11:08:44 UTC Running 1 Simulations
| 81 | -3.021010 | 0.5 | 0.5 | 1.0900474 | 0.8700573 | 1.6 | 1.5522617 | 1.6 | 1.4544093 | 1.5761775 | 1.4908380 | 1.3553925 | 1.6 | 1.1194879 |
11:09:21 UTC Running 1 Simulations
| 82 | -2.655108 | 0.5 | 0.5 | 1.1153931 | 0.8558697 | 1.6 | 1.5762891 | 1.6 | 1.4643217 | 1.5608396 | 1.4910057 | 1.3198305 | 1.6 | 1.1132735 |
11:10:17 UTC Latest Best Fit on Iter 51: -2.655
11:10:18 UTC Running 1 Simulations
| 83 | -2.977021 | 0.5 | 0.5 | 1.1034215 | 0.8628779 | 1.6 | 1.5744785 | 1.6 | 1.4712531 | 1.5856606 | 1.4958391 | 1.3133095 | 1.6 | 1.1196191 |
11:11:06 UTC Running 1 Simulations
| 84 | -2.518340 | 0.5 | 0.5 | 1.1015773 | 0.8581018 | 1.6 | 1.5603784 | 1.6 | 1.4738437 | 1.5429247 | 1.4976516 | 1.3241458 | 1.6 | 1.1277677 |
11:12:42 UTC Latest Best Fit on Iter 53: -2.518
11:12:44 UTC Running 1 Simulations
| 85 | -2.548534 | 0.5 | 0.5 | 1.1071853 | 0.8580468 | 1.6 | 1.5597170 | 1.6 | 1.4469135 | 1.5451412 | 1.5110344 | 1.3200677 | 1.6 | 1.1277982 |
11:13:28 UTC Running 1 Simulations
| 86 | -2.446199 | 0.5385742 | 0.5 | 1.1026333 | 0.8565676 | 1.6 | 1.5702467 | 1.6 | 1.4623088 | 1.5404111 | 1.4962952 | 1.3243839 | 1.6 | 1.1297109 |
11:15:18 UTC Latest Best Fit on Iter 55: -2.446
11:15:20 UTC Running 1 Simulations
| 87 | -2.629518 | 0.5300255 | 0.5 | 1.0933592 | 0.8614966 | 1.6 | 1.5898284 | 1.6 | 1.4618496 | 1.5361581 | 1.5130843 | 1.3296225 | 1.6 | 1.1252615 |
11:16:00 UTC Running 1 Simulations
| 88 | -2.462754 | 0.5197740 | 0.5 | 1.1090171 | 0.8562094 | 1.6 | 1.5792717 | 1.6 | 1.4550578 | 1.5443719 | 1.4915863 | 1.3269635 | 1.6 | 1.1469812 |
11:17:03 UTC Running 1 Simulations
| 89 | -2.372081 | 0.5305973 | 0.5 | 1.0981983 | 0.8570101 | 1.6 | 1.5657280 | 1.6 | 1.4578734 | 1.5369034 | 1.4845761 | 1.3196305 | 1.6 | 1.1301083 |
11:17:54 UTC Latest Best Fit on Iter 58: -2.372
11:17:56 UTC Running 1 Simulations
| 90 | -2.418177 | 0.5328262 | 0.5 | 1.1134120 | 0.8553410 | 1.6 | 1.5656784 | 1.6 | 1.4591374 | 1.5422412 | 1.4893737 | 1.3251099 | 1.6 | 1.1301498 |
11:18:31 UTC Running 1 Simulations
| 91 | -2.407531 | 0.5291998 | 0.5 | 1.0973336 | 0.8606927 | 1.6 | 1.5630820 | 1.6 | 1.4651320 | 1.5418543 | 1.4930209 | 1.3182749 | 1.6 | 1.1396851 |
11:19:03 UTC Running 1 Simulations
| 92 | -2.360315 | 0.5284394 | 0.5 | 1.0943213 | 0.8547125 | 1.6 | 1.5690066 | 1.6 | 1.4532833 | 1.5339190 | 1.4826115 | 1.3196519 | 1.6 | 1.1281383 |
11:19:35 UTC Latest Best Fit on Iter 61: -2.36
11:19:37 UTC Running 1 Simulations
| 93 | -2.386844 | 0.5276519 | 0.5 | 1.0984220 | 0.8593770 | 1.6 | 1.5723469 | 1.6 | 1.4601598 | 1.5333727 | 1.4815933 | 1.3216155 | 1.6 | 1.1295624 |
11:20:14 UTC Running 1 Simulations
| 94 | -2.364264 | 0.5304796 | 0.5 | 1.0941835 | 0.8525050 | 1.6 | 1.5666038 | 1.6 | 1.4480961 | 1.5375164 | 1.4868082 | 1.3187537 | 1.6 | 1.1296253 |
11:21:26 UTC Running 1 Simulations
| 95 | -2.359467 | 0.5284423 | 0.5 | 1.0903633 | 0.8543516 | 1.6 | 1.5634501 | 1.6 | 1.4534515 | 1.5368020 | 1.4875853 | 1.3170847 | 1.6 | 1.1319118 |
11:23:13 UTC Latest Best Fit on Iter 64: -2.359
11:23:15 UTC Running 1 Simulations
| 96 | -2.400466 | 0.5303865 | 0.5 | 1.1006098 | 0.8525986 | 1.6 | 1.5701997 | 1.6 | 1.4477451 | 1.5374539 | 1.4851559 | 1.3229099 | 1.6 | 1.1281442 |
11:24:33 UTC Running 1 Simulations
| 97 | -2.344397 | 0.5326708 | 0.5 | 1.0980564 | 0.8530274 | 1.6 | 1.5709758 | 1.6 | 1.4609866 | 1.5337871 | 1.4869915 | 1.3108677 | 1.6 | 1.1317702 |
11:25:24 UTC Latest Best Fit on Iter 66: -2.344
11:25:27 UTC Running 1 Simulations
| 98 | -2.367201 | 0.5241184 | 0.5 | 1.0926104 | 0.8583426 | 1.6 | 1.5677291 | 1.6 | 1.4592777 | 1.5365898 | 1.4876420 | 1.3160742 | 1.6 | 1.1361725 |
11:27:18 UTC Running 1 Simulations
| 99 | -2.338875 | 0.5356358 | 0.5 | 1.1007470 | 0.8524496 | 1.6 | 1.5734875 | 1.6 | 1.4604786 | 1.5329858 | 1.4868023 | 1.3104059 | 1.6 | 1.1307632 |
11:28:20 UTC Latest Best Fit on Iter 68: -2.339
11:28:22 UTC Running 1 Simulations
| 100 | -2.340785 | 0.5345387 | 0.5 | 1.0990839 | 0.8594800 | 1.6 | 1.5661724 | 1.6 | 1.4608347 | 1.5353546 | 1.4869680 | 1.3138533 | 1.6 | 1.1269174 |
11:28:52 UTC Best Result: -2.338875812804492 Best Parameters: w_0: 0.5356358395811599 w_1: 0.5 w_2: 1.1007470922035667 w_3: 0.8524496873280988 w_4: 1.6 w_5: 1.5734875546193445 w_6: 1.6 w_7: 1.4604786696611824 w_8: 1.5329858004284218 w_9: 1.4868023890954158 w_10: 1.310405955631697 w_11: 1.6 w_12: 1.130763290988251
Results
The best result can be extracted directly from the optimizer object. Plotting this, we see what the optimizer has returned as the optimized structure for this design of Y-junction.
best_params = results.optimizer.max["params"]
print(f"Best fitness: {results.optimizer.max['target']}")
sim = fn_pre(**best_params)
sim.plot(z=0)
plt.show()Best fitness: -2.338875812804492
We can then create a plot to evaluate if the Bayesian optimization has converged on a fitness value. The first 30 iterations are from the random initialization, so the fitness values are expected to be more varied. The fitness has converged before the end of the remaining 70 iterations; an early stop criteria could have been included to finish the optimization sooner.
ax = df["output"].plot(xlabel="Simulation Number", ylabel="Fitness")
The mean power at the reflected and transmitted monitors can be calculated and compared to the paper from Gao et al.
idx_best_result = df["output"].idxmax()
best_sim_filename = results.task_paths[idx_best_result]
best_sim = td.SimulationData.from_file(best_sim_filename)
power_reflected = np.array(
np.squeeze(np.abs(best_sim["mode_11"].amps.sel(direction="-", mode_index=0)) ** 2)
).mean()
power_transmitted = np.array(
np.squeeze(np.abs(best_sim["mode_12"].amps.sel(direction="+", mode_index=0)) ** 2)
).mean()
print(f"Mean Reflected Power: {round(power_reflected, 3)} (Paper: 0.004)")
print(f"Mean Transmitted Power: {round(power_transmitted, 3)} (Paper: 0.460)")Mean Reflected Power: 0.001 (Paper: 0.004)
Mean Transmitted Power: 0.482 (Paper: 0.460)Final Optimized Design
Finally, we can simulate the optimal solution with an additional FieldMonitor to visualise the field distribution throughout the splitter.
final_sim = fn_pre(**best_params)
# Define a field monitor to help visualize the field distribution
field_monitor = td.FieldMonitor(
center=(0, 0, 0), size=(td.inf, td.inf, 0), freqs=[freq0], name="field"
)
mode_12 = final_sim.monitors[1]
final_sim = final_sim.copy(update={"monitors": (field_monitor, mode_12)})
final_sim_data = web.run(final_sim, task_name="BO_notebook_final_sim")Created task 'BO_notebook_final_sim' with resource_id 'fdve-c9a0b168-cf11-4d44-8c09-8b21fc780959' and task_type 'FDTD'.
View task using web UI at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-c9a0b168-cf1 1-4d44-8c09-8b21fc780959'.
Task folder: 'default'.
11:28:54 UTC Estimated FlexCredit cost: 0.116. 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.
11:28:55 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.
11:29:07 UTC starting up solver
running solver
11:29:18 UTC early shutoff detected at 76%, exiting.
status = postprocess
11:29:21 UTC status = success
11:29:23 UTC View simulation result at 'https://tidy3d.simulation.cloud/workbench?taskId=fdve-c9a0b168-cf1 1-4d44-8c09-8b21fc780959'.
11:29:24 UTC Loading results from simulation_data.hdf5
Plotting the field shows how it splits evenly between each branch.
final_sim_data.plot_field("field", "E", "abs^2")
plt.show()
And we can visualise how the power transmitted varies over the frequency range.
power_transmitted = np.squeeze(
np.abs(final_sim_data["mode_12"].amps.sel(direction="+", mode_index=0)) ** 2
)
plt.plot(freqs, power_transmitted)
plt.title("Power transmitted across frequency range")
plt.xlabel("Frequency / Hz")
plt.ylabel("Power")
plt.show()
Conclusion
Through comparison of the junction geometry and the transmitted and reflected power, we can show that these results closely follow the results published by Gao et al. for the design of a Y-junction. This notebook demonstrates how to carry out Bayesian optimization with the Design plugin, and can be readily adapted to other use cases.
Learn more about the Design plugin:
TERMS & CONDITIONS
How the process works:
Try this real URL to see an example: