Prediction of geometrical parameters given a dispersion curve.
We aim for advanced dispersion engineering design of microresonators to achieve broad (ideally octave-spanning) frequency combs.
Frequency comb is an array of equidistant spikes in frequency domain. Recently, FC were discovered in microresonators, where they can be generated due to the optical parametric oscillations, i.e four-wave mixing.
However, this nonlinear process requires high optical confinement of the light, so that non-linear processes can play a significant role, and an appropriate dispersion of the light modes in microresonator. The former can be obtained by increasing the Q-factor of the microresonator, while the latter one is a problem of dispersion engineering, where the geometry of the microresonator is designed in a special way to guide the light.
It appears that one of the key conditions for comb generation - is an anomalous dispersion (where the distance between two consecutive resonances of the resonator increases with increasing optical frequency, D2>0 (on the figures) ).
The dispersion profile (Dint vs frequency) determines the spectrum of the comb one can achieve in microresonator.
Generally, we want Dint (see figure above) to be slightly > 0 and very broad, thus the comb is flat and broad. If |Dint| is significantly nore than 0, which means that the mode is far from the corresponding microresonator resonance, the comb substantially decays.
Conventional way to engineer dispersion is to play with geometrical parameters of the single waveguide.
However, this is already well established technique, which has already met its limit. New approaches has been proposed: combining several resonators to locally induce anomalous/normal dispersion thus extending overall Dint profile (see for example: [Dispersion engineering and frequency comb generation in thin silicon nitride concentric microresonators]). This happens due to the mode hybridization between two resonators at a given range of frequincie, where one can effect the dispersion.
On the figure above one can observe extension of the Dint (Mode deviation) profile, due to locally induced anomalous dispersion (around 220THz). For the red curve comb would vast until until 230 THz on the low frequency side, while for hybridized mode it would go 190 THz. For high frequencies red and blue curves are almost the same.
Generally, to predict dispersion profile, one usually simulates distribution of eigenmodes in microresonator, which is done by solving system of Maxwell equations in commercial softwear like COMSOL or ANSYS (via Finite Elements Method). Combining two resonators already becomes hard problem as the number of geometrical parametters increases. Now one has: widths of 2 resonators, heights, gap distance and radius. One round of simulations takes about 10 min for quite sparse array of frrequencies taken to calculate modes at (we calculate modes in range from 180THz up to 380THz each 10Thz), (8 CPUs and many RAM).
Here one might think to apply some fancy algorithms to be able to predict parameters of the system, given the known dispersion curve. Our idea is to feed a Neural Network with simulated batch of data, so it is later able to map dispersion curve that we draw to real parameters of the microresonator.
Acquiring the data was the toughest challenge. One can estimate, that 10 mins/sample turns to be around 140 samples per day. Due to the license restrictions it is possible to use the softwear on one computer only. Also there is a tradeoff between quality of the data and quantity of the samples. Trying to get about 6000 samples we simulated quite sparse dispersion array which results in low precesion.
One might think about optimization problem, where it is possible to assign kind off score of each dispersion curve based on its width and slight positive deviations from zero, or even multiple zero-crossings ([which leads to dispersive waves]) and try to optimize it (where generative algorithms might be useful). But here we focus on a Neural Network implementation, eventhough it might be not optimal way to solve this sort of problem.
As an input for the network we use dispersion data and for the output - geometrical parameters, corresponding to this curve.
import pickle
## images data
with open('data/rede/raw/5070-62-111-images.pkl', 'rb') as f:
img = pickle.load(f)
## vectors of Dispersion and frequency grid
with open('data/rede/raw/5070-frequencies_modes.pkl', 'rb') as f:
Disp = pickle.load(f)
## Geometrical parameters
with open('data/rede/raw/5070-5-parameters.pkl', 'rb') as f:
parameters_list = pickle.load(f)
To get processed data you can visit [prepare_dataset_keras.ipynb]. (These 2 pictures correspond to differnt parameters).
NN implementation, visit [prepare_dataset_keras.ipynb]
- First, we feed our network with one-dimensional vectors of Dint corresponding to given frequencies (same for the whole vectors)
- Second, we implemented convolutional NN, where it is fed with images.
- Third, we combined 2 networks.
- Anton Lukashchuk <anton.lukashchuk@epfl.ch>
- Anton Karazeev <anton.karazeev@phystech.edu>
- Maxim Karpov <maxim.karpov@epfl.ch>