Abstract
Topology opens many new horizons for photonics, from integrated optics to lasers. The complexity of largescale devices asks for an effective solution of the inverse problem: how best to engineer the topology for a specific application? We introduce a machinelearning approach applicable in general to numerous topological problems. As a toy model, we train a neural network with the Aubry–Andre–Harper band structure model and then adopt the network for solving the inverse problem. Our application is able to identify the parameters of a complex topological insulator in order to obtain protected edge states at target frequencies. One challenging aspect is handling the multivalued branches of the direct problem and discarding unphysical solutions. We overcome this problem by adopting a selfconsistent method to only select physically relevant solutions. We demonstrate our technique in a realistic design and by resorting to the widely available opensource TensorFlow library.
Introduction
The rapidly growing interest in topological photonics^{1,2} is leading to the design of complex structures for many applications of optical topological insulators^{3}. One leading goal of topological photonics is photon transport protected from unwanted random scattering. This is achieved by realizing analogs of the quantum Hall effect^{4,5,6} through magneticlike Hamiltonians in photonic systems^{7}. In the optical domain, topological insulators^{8} have been implemented in modulated honeycomb lattices^{7}, in arrays of coupled opticalring resonators^{9}, and optical quantum walks^{10}. Geometryindependent topological structures have been proposed to obtain nonreciprocal single mode lasing^{11,12,13,14} as well as systems with balanced gain and loss for paritytime symmetric structures with topological order^{15,16}. Emulations of fourdimensional (4D) physics have also been reported^{17,18}. By using onedimensional (1D) Harper modulations, it is possible to simulate twodimensional (2D) topological systems. Similarly, by 2D topological systems, one can simulate 4D ones, as recently investigated in refs.^{17,18}.
One challenge in this field is to find an effective methodology for the inverse problem in which the target optical properties result from topological characteristics. Although various computational techniques are available, these require specific implementations tailored to the task at hand. Machine learning (ML)^{19,20,21} has recently been proposed as an encompassing technology for dealing with greatly differing problems through a unified approach. ML techniques have shown a remarkable growth in sophistication and application scope in multiple fields^{22,23,24}; ML offers exciting perspectives in topological photonics. ML is applied in two main classes of problems: (i) classification for categorizing information and (ii) regression to predict continuous values, both typically performed by supervised training. Unlike parametric regression—in which a best fit of the data is determined on the basis of a specific function—ML regression employs a neural network (NN) emulating the behavior of the data on which it has been trained: “the NN learns the model”.
In this paper, we employ ML regression for solving the inverse problem in topological photonics. We apply advanced ML techniques to design photonic topological insulators enabling innovative applications through custom tailoring of desired optical parameters. In our approach, we introduce a twist in order to ensure that only physically possible solutions are found. This twist is based on a selfconsistent cycle in which a tentative solution obtained from the inverse problem NN is run through the direct problem NN in order to ensure that the solution obtained is indeed viable. This has the added benefit of checking that multivalued degeneracy has been effectively removed.
Results
We consider one of the simplest structures that support nontrivial topological properties. In 1D systems, synthetic magnetic fields occur by lattice modulation^{25} of the optical structure. In the Aubry–Andre–Harper (AAH) model^{26,27}, identical sites—resonators, twolevel atoms, waveguides, etc.—are centered at positions \(z_n = d_{\rm o}\left( {n + \eta \delta _n^H} \right)\), with n an integer label, d_{o} the primary lattice period, η the modulation strength, and \(\delta _n^H = {\mathrm{cos}}\left( {2\pi \beta n + \phi } \right)\) the Harper modulation^{27}. The parameter β is the frequency of the Harper modulation. Together, β and the phase shift ϕ furnish the topological properties by a “2D ancestor” mapping^{28}. The 2D ancestor is characterized by the dependence of the dielectric function on the coordinate z and on the parameter ϕ, which acts as a periodic artificial coordinate. Hence, the phase ϕ can be treated as a wave vector in a fictitious auxiliary direction^{28}. For β = p/q with p > 0 and q > 0 integers, the lattice displays two commensurate periods with q sites z_{n} in the unitcell. Properly chosen parameters give rise to nontrivial topological phases with protected states at the border of the structure. These “edgestates” are hallmarks of topological insulators. The phase ϕ tunes edgestate eigenfrequency in the photonic bandgaps.
Our photonic topological insulator is an array of layers A of normalized thickness ξ = L_{A}/d_{o}, centered in z_{n}, in an homogeneous bulk of material B. This kind of structure can be effectively modeled by the transfer matrix technique^{16,29}, as reported in Fig. 1a. In this figure A_{0} and A_{n} are the initial and final amplitudes of the righttraveling waves; while B_{0} and B_{n} are their equivalent for the lefttraveling wave amplitudes. As detailed in Methods, we obtain the transfer matrix for the single period T^{(1)}(ω, ϕ, ξ) with elements \(T_{11}^{(1)}\), \(T_{12}^{(1)}\), \(T_{21}^{(1)}\), and \(T_{22}^{(1)}\). Figure 1a shows the final wave amplitudes A_{n}, B_{n} by the nfold repeated action of T^{(1)}(ω, ϕ, ξ) on A_{0}, B_{0}. The dielectric constant profile  for the case β = 1/3 is schematically illustrated in Fig. 1b.
For η = 0, we have a periodical unmodulated structure with frequency bandgaps labeled by an integer i. For η ≠ 0, each gap of the unmodulated structure splits into q gaps, each one labeled by indices (i, j) (j = 1, …, q)^{30}. This splitting is shown in Fig. 1c for β = 1/3 with respect to the variable χ = ϕ + π(2β − 1)/2.
As detailed in Methods and illustrated in Fig. 1d, enforcing boundary conditions at the left edge^{31,32} and defining the function Q(ω, ϕ, ξ) enables one to establish the presence of edge states corresponding to poles ω_{t} of the reflection coefficient. However, the function ω_{t} = ω(χ, ξ) cannot be analytically inverted to express the geometrical parameters χ and ξ in terms of the variable ω_{t}. Exploiting ML techniques, we solve this inverse problem and design topological insulators with target edge modes. The inverse problem in artificial NN theory—and therefore in ML—is widely discussed in numerical modeling, engineering, and other fields^{33,34,35,36,37}. Regression in ML optimizes an NN so that a given vector input (\({\Bbb R}^n\)) will result in a scalar (\({\Bbb R}\)) output, emulating the behavior of the training data. A regressive NN is a configuration of computational layers such that a specific set of input nodes \(\underline I\) is connected to a single output node, through a configurable set of N_{h} hidden layers each containing n_{i} nodes h_{ij}, where i = 1, …. N_{h} and j = 1, …. n_{i}. Examples of such regressive NNs are shown in Fig. 2a, b. A generic node k + 1, j, shown in Fig. 2c, receiving as inputs h_{ki}, with i = 1, …. n_{k}, yields on output \(h_{k \,+\, 1j} = g\left( {\mathop {\sum}\nolimits_l {\kern 1pt} w_{k \,+ \, 1jkl}h_{kl} + b_{k \,+\, 1j}} \right)\), with g(x) being a nonlinear activation function, w_{k + 1jkl} the weight of h_{kl} on h_{k+1j} with a bias term b_{k + 1j}. Following accepted practice, our activation function is g(x) = tanh(x).
Optimization of the NN is performed by minimizing a cost function by a gradient descent method that updates weights and biases. In the initial state, weights w_{ijkl} are selected from a truncated normal and biases are set to zero. Training applies this procedure to a dataset randomly split into two separate classes: (i) an actual training set and (ii) a validation set. The network is iteratively updated until the error on the validating dataset converges to a given rate.
The inverse topological problem at hand is to obtain the desired optical behavior: a target edgestate at frequency ω_{t}, which is an input to the design (Fig. 2a). ML techniques achieve this result by modeling the multidimensional nonlinear relationships among all the structure parameters ω_{t}, χ, β, \(\epsilon _{\rm A}\), \(\epsilon _{\rm B}\), and ξ. In our specific case, the dataset fixes \(\epsilon _{\rm A}\), \(\epsilon _{\rm B}\), β at the values \(\epsilon _{\rm A} = 9\), \(\epsilon _{\rm B} = 4\) and β = 1/3.
First, we generate a dataset to train our NNs by numerically computing the complex roots of \(T_{12}^{(1)}(\omega ,\chi ,\xi )\) covering the region of interest for parameters χ and ξ. The real part of these roots, shown in Fig. 3a, represents the edge states dispersion. Interestingly the same dataset can be used both for the inverse and direct NN training phase, by suitably selecting the features and target fields. The inverse problem NN (Fig. 2a) targets a value χ = χ_{o}, a topological parameter on the basis of features including ω_{t}. For a direct problem (Fig. 2b), the mode frequency ω_{t} would be the target of a network whose features include the topological parameters (χ, ξ).
The dataset contains various branches since there exist an edge state for each band gap (i, j) with j ≠ 3, as results by Eq. (2) in Methods. Due to the folding of the Brillouin zones, the edge state frequency ω(χ, ξ) is then a multimode function, which we unfold by introducing a label \(m_{ij}^ \pm\) for each mode; here i = 1, …∞ and j = 1, …q, while the sign ± indicates modes in the positive/negative χ domain. In Fig. 3a, data points with different ij values are identified with different colors, and solving the inverse problem is a matter of determining when these surfaces intercept a specific target value of the ω axis. Three outcomes are possible: a single value for χ and ξ when a monotonic mode surface is intercepted, no solution for values of ω laying between surfaces, and multiple solutions in other cases. This implies that the feature set (χ, ξ, \(m_{ij}^ \pm\)) is insufficient. To tackle this problem, we take into account the trend s_{±} = sgn(dω_{t}/dχ) as an additional variable. The NNs with this enlarged feature set are illustrated in Fig. 2a, b.
In the terminology used in ML, the mode index \(m_{ij}^ \pm\) and trend s_{±} labels are "categorical features'' and lead to two possible courses of action for the actual implementation of the NNs used in our problem. One in which a single NN is constructed in a hybrid feature space with both continuous variables (real valued ξ's and χ's) and categorical features, as illustrated in Fig. 2b. Another course is to adopt multiple independent NNs, one NN for each mode and each trend.
The single NN approach is hindered by the presence of discontinuities in the features domain: with respect to the ω variable they are a consequence of the fact that edge states fall within the bulk energy gaps; with respect to the χ variable these arise from considering only the leftedge states. Figure 3a clarifies this aspect. Due to these discontinuities, we have chosen to use multiple independent NNs.
Moreover, when considering the solution provided by the inverse NNs, we identify a specific problem in the use of ML as they may furnish solutions that are not physical. An example of this issue is given in Fig. 3b where—for a fixed band and a fixed ξ—the curve representing ω as a function of χ is shown together with its inverse (Fig. 3c). Inverting the function ω(χ), we consider an interval of values for ω spanning from its minimum ω_{min} to the maximum ω_{max}, but for the two branches of the inverse function χ(ω)—identified by colors in Fig. 3c—the range of ω is different. For example, for the red branch, the maximal value of ω is \({\omega{\prime}} _{{\rm max}} < \omega _{{\rm max}}\). When the target frequency is outside of this range, the NN produces an output outside of the physically acceptable range for χ (see details in Supplementary Information: Supplementary Figs. 1–3). The inverse NN can furnish spurious nonphysical solutions.
Our approach tackles this issue by a twostep selfconsistent cycle, detailed in the Supplementary Information (Supplementary Fig. 4): (i) in the first stage, a desired input ω_{t} forms part of the feature set \(\left( {\omega _t,m_{ij}^ \pm ,s_ \pm } \right)\) resulting in the output χ_{o} of the inverse NN; this set is used as input \(\left( {\chi _{\rm o},m_{ij}^ \pm ,s_ \pm } \right)\) to a direct problem network; (ii) in the second stage, the target of this direct network ω_{sc} is compared with the input value ω_{t} and χ_{o} is retained as a solution of the inverse model if \(\left {\omega _{sc}  \omega _t} \right < \delta\) with δ a userdefined small positive quantity. The value of δ affects the model accuracy (see Supplementary Fig. 5 and related comments). A reasonable choice can be \(\delta \sim E_j^{{\rm max}}\) (with j = I,D), i.e., the maximum value of the squared error functions for the inverse (I) and the direct (D) networks.
The training dataset was generated with 11 sets of ξ ranging from 0.10 to 0.20 in steps of 0.01 and for each set χ spans −π to π with 997 equally spaced values. Results based on using an array of NNs each composed of 5 hidden layer of 131 nodes are shown in Fig. 4 together with its training set (colored lines). The model was developed using 80% of the dataset randomly chosen, the rest being used for validation and comprising of 250,000 steps. Training each model takes about 8 min on our hardware using a single Nvidia GPGPU Tesla K20c. The purple dots in Fig. 4 are based on 100 values of ξ while exploring the ω domain with a resolution of 10^{−5}. Each array element is trained for a specific value of the categorical features and pertains to either the positive or the negative χ domain.
The results of applying the direct and inverse NNs, portrayed in Fig. 4a, b, respectively, show that the proposed method gives accurate solutions matching the original data in the whole range of interest. Figure 4 clearly shows that our ML strategy solves the inverse topological design problem.
Discussion
The inverse problem in topological design is solved by a supervised ML regression technique. We employ a selfconsistent procedure to rule out unphysical solutions enabling tailored engineering of protected edgestates. We successfully tackle multivalued functions introducing categorical features, as the trend, which tags training data according to their gradient’s sign. Discontinuous domains are effectively treated by adopting multiple independent NNs each one specific to its domain. Our general method can be extensively applied—well beyond the example considered in this work—and may also be exploited for other physical systems in topological science, as polaritonics^{38,39}, quantum technologies and ultracold atoms^{40,41}. The method is scalable to very complex structures involving hundreds of topological devices, as those recently considered for large scale synchronization^{42}, and frequency comb generation^{43}, eventually including nonHermitian systems^{44,45}. Further applications include 2D and 3D topological systems^{11} and quantum sources and simulations^{17,18}.
Methods
TensorFlow
Tensorflow is Google’s versatile opensource multiplatform dataflow library capable of efficiently performing ML tasks such as implementing NNs (https://tensorflow.org). Multidimensional data arrays, referred to as “tensors” are executed on the basis of stateful dataflow graphs, hence the name TensorFlow. For our final code implementation, Tensorflow version 1.3 with python API bindings was used.
The nature of our problem is such that there is a discontinuity in ξ = 0 which cannot be correctly handled by a single NN bridging this point; this is relevant to both the inverse and direct cases. Breaking up the dataset into two parts to be used for two separate NNs is the simplest solution to this problem.
Another interesting aspect is related to the fact that the feature set in our inverse and direct NNs contain both continuous and discrete variables. The discrete variables can either be treated as such or handled by constructing multiple NNs each relative to a specific value of the discrete variable. The trend variable which has two possible values is one such case as is the mode number. In our code, we have implemented a flexible system which allows one to decide which discrete variables are to be included in each NN, the others being broken up into arrays of NNs one for each value of the variable. Once the bookkeeping issues have been tackled, this generalized approach allows one to tailor the problem to the given dataset.
Transfer matrix
Given the stepped and periodic dielectric function of period D = qd_{o}:
in each layer, the electric field can be represented as the superposition of a left and a righttraveling wave. Applying the boundary conditions, the matrices
with α, γ = A or B and \(r_{\alpha \gamma } = {\textstyle{{q_\gamma  q_\alpha } \over {q_\gamma + q_\alpha }}}\), describe the light propagation through the interfaces, having introduced \(q_\alpha = (\omega {\mathrm{/}}c)\sqrt {\epsilon _\alpha }\), while the propagation within each layer A and B is given by:
where s_{n} = [z_{n+1} − z_{n} − L_{A}]/d_{o} are the normalized thicknesses of the B layers.
From these, we obtain the transfer matrix for the single period T^{(1)}(ω), the matrix connecting the fields in the left side of the elementary cell to the ones in the right side:
with M = M_{AB}T_{A}M_{BA}. The quantity \(\rho =  \frac{1}{2}TrT^{(1)}(\omega ,\phi ,\xi )\) allows one to locate bulk bands in the regions where \(\rho ^2\leqslant 1\), and gaps where ρ^{2} > 1. Alternatively, the amplitude \(\left {r_\infty (\omega ,\phi ,\xi )} \right^2\) of the reflection coefficient of the structure^{28}
where e^{ik(ω)D} is an eigenvalue of the matrix T^{(1)}(ω, ϕ, ξ), can also be used to locate the gaps of the system.
Band structure of the unmodulated system
The unmodulated structure (η = 0) features stopbands at \(\tilde \omega _0\) = \(\omega _0d_0{\mathrm{/}}c\) = \(\pi {\mathrm{/}}\left( {\sqrt {\varepsilon _{\rm A}} + (1  \xi )\sqrt {\varepsilon _{\rm B}} } \right)\), where ξ = L_{A}/d_{o} is the characteristic size ratio.
Q(ω, ϕ, ξ) function
To determine the existence of the edge states, one needs to specify the boundary conditions on each edge of the structure. For the left edge, this condition is given by:
where A_{1} and B_{1} are the amplitudes of the right and lefttraveling waves in the first layer of the structure. This condition can be reformulated as
with b_{1} = ((q_{a} − q_{b}), (q_{a} + q_{b}))^{T} and a_{1} = (A_{1}, B_{1})^{T}, and together with the eigenvalues λ_{±} and eigenvectors \(v_ \pm = \left( {T_{12}^{(1)},\lambda _ \pm  T_{11}^{(1)}} \right)\) of the transfer matrix T^{(1)}, it is possible to determine the existence and dispersion of edge states.
Following refs.^{31,32}, it can be in fact shown that a proportionality relation exists between the boundary vector b_{1} and the eigenvectors v_{±} of the transfer matrix. So the condition for the existence of the edge states is given by det(b_{1}, v_{±}) = 0 in a gap where \(\left {\lambda _ \pm } \right < 1\). This entails searching for the zeros of the function F_{l,±} = \(\left( {q_{\rm A} q_{\rm B}} \right)\left( {\lambda _ \pm  T_{11}^{(1)}} \right)\) − \(T_{12}^{(1)}\left( {q_{\rm A} + q_{\rm B}} \right)\).
Specifically, the real part of F_{l,±} = 0 yields the function Q(ω, ϕ, ξ) = \({\rm Re}\left\{ {T_{12}^{(1)}\left( {q_{\rm A} + q_{\rm B}} \right)} \right.\) − \(\left. {\left( {q_{\rm A}  q_{\rm B}} \right)\left( {T_{22}^{(1)}  T_{11}^{(1)}} \right){\mathrm{/}}2} \right\}\) and, as shown in Fig. 1c, this implies that edge states exist only in the gaps where ρ > 1 and Q(ω, ϕ, ξ) · ρ > 0. At the same time, edge states cannot exist in gaps where Q(ω, ϕ, ξ) does not change sign. Moreover, due to a bulkboundary correspondence^{46}, the number of these edge modes is equal to the modulus of the associated topological invariant ν_{ij}, given by the winding number of the reflection coefficient:
i.e., the extra phase (divided by 2π) of r_{∞} (ω, χ) when χ varies in the range (−π, π) with ω in the stop band^{47}.
By relying on the transfer matrix method, our approach can be applied to a general class of problems and thus makes it suitable for a wide range of systems beyond our baseline AAH model. Specifically, it can be extended to many physical systems whose behavior is described by a gapped unitary operator, e.g., photonic Floquet topological insulators^{7,48} and photonic topological quantum walks^{10}. Analogously to the AAH model, the edge states of these systems can be defined with an equivalent F_{l,±}(ω, p_{1}, ..p_{n}) function, where (p_{1}, ..p_{n}) are relevant parameters describing the structure. The imaginary component of F_{l,±}(ω, p_{1}, ..p_{n}) = 0 furnishes the dispersion relations of the edge modes and hence the training dataset of our ML inverse problem.
Code availability
The code developed for the present study is available from the corresponding author on reasonable request.
Data availability
The datasets generated during the current study are available from the corresponding author on reasonable request.
References
 1.
Lu, L., Joannopoulos, J. D. & Soljai, M. Topological photonics. Nat. Photonics 8, 821 (2014).
 2.
Wu, Y. et al. Applications of topological photonics in integrated photonic devices. Adv. Opt. Mater. 5, 1700357 (2017).
 3.
Ozawa, T. et al. Topological photonics. Preprint at ArXiv (2018).
 4.
Haldane, F. D. M. & Raghu, S. Possible realization of directional optical waveguides in photonic crystals with broken timereversal symmetry. Phys. Rev. Lett. 100, 013904 (2008).
 5.
Raghu, S. & Haldane, F. D. M. Analogs of quantumhalleffect edge states in photonic crystals. Phys. Rev. A. 78, 033834 (2008).
 6.
Wang, Z., Chong, Y. D., Joannopoulos, J. D. & Soljacic, M. Reflectionfree oneway edge modes in a gyromagnetic photonic crystal. Phys. Rev. Lett. 100, 013905 (2008).
 7.
Rechtsman, M. C. et al. Photonic floquet topological insulators. Nature 496, 196 (2013).
 8.
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045 (2010).
 9.
Hafezi, M., Mittal, S., Fan, J., Migdall, A. & Taylor, J. Imaging topological edge states in silicon photonics. Nat. Photonics 7, 1001 (2013).
 10.
Kitagawa, T. et al. Observation of topologically protected bound states in photonic quantum walks. Nat. Commun. 3, 882 (2012).
 11.
Bahari, B. et al. Nonreciprocal lasing in topological cavities of arbitrary geometries. Science 358, 636 (2017).
 12.
Bandres, M. A. et al. Topological insulator laser: experiments. Science 359, eaar4005 (2018).
 13.
Harari, G. et al. Topological insulator laser: theory. Science, http://science.sciencemag.org/content/early/2018/01/31/science.aar4003 (2018.)
 14.
StJean, P. et al. Lasing in topological edge states of a onedimensional lattice. Nat. Photonics 11, 651–656 (2017).
 15.
Rivolta, N. X. A., Benisty, H. & Maes, B. Topological edge modes with ptsymmetry in a quasiperiodic structure. Phys. Rev. A 96, 023864 (2017).
 16.
Pilozzi, L. & Conti, C. Topological lasing in resonant photonic structures. Phys. Rev. B 93, 195317 (2016).
 17.
Zilberberg, O. et al. Photonic topological boundary pumping as a probe of 4d quantum Hall physics. Nature 553, 59 (2018).
 18.
Lohse, M., Schweizer, C., Price, H. M., Zilberberg, O. & Bloch, I. Exploring 4d quantum Hall physics with a 2d topological charge pump. Nature 553, 55 (2018).
 19.
Bishop, C. Pattern Recognition and Machine Learning (Cambridge, Springer, 2006).
 20.
Duda, R. Pattern Classification. (Wiley, New York, 2001).
 21.
Murphy, K. Machine Learning: A Probabilistic Perspective. (The MIT Press, Cambridge, 2012).
 22.
Zdeborova, L. Machine learning: new tool in the box. Nat. Phys. 13, 420 (2017).
 23.
Carrasquilla, J. & Melko, R. G. Machine learning phases of matter. Nat. Phys. 13, 431 (2017).
 24.
Zhang, Y. & Kim, E.A. Quantum loop topography for machine learning. Phys. Rev. Lett. 118, 216401 (2017).
 25.
Kraus, Y. E. & Zilberberg, O. Topological equivalence between the fibonacci quasicrystal and the harper model. Phys. Rev. Lett. 109, 116404 (2012).
 26.
Aubry, S. & André, G. Analyticity breaking and anderson localization in incommensurate lattices. Ann. Isr. Phys. Soc. 3, 133 (1980).
 27.
Harper, P. G. Single band motion of conduction electrons in a uniform magnetic field. Proc. Phys. Soc. Lond. A 68, 874 (1955).
 28.
Poshakinskiy, A. V., Poddubny, A. N., Pilozzi, L. & Ivchenko, E. L. Radiative topological states in resonant photonic crystals. Phys. Rev. Lett. 112, 107403 (2014).
 29.
Chew, W. C. Waves and Fields in Inhomogeneous Media. (WileyIEEE Press, New York, 1999).
 30.
Hofstadter, D. R. Energy levels and wave functions of bloch electrons in rational and irrational magnetic fields. Phys. Rev. B 14, 2239–2249 (1976).
 31.
Hatsugai, Y. Edge states in the integer quantum Hall effect and the riemann surface of the Bloch function. Phys. Rev. B 48, 11851 (1993).
 32.
Tauber, C. & Delplace, P. Topological edge states in twogap unitary systems: a transfer matrix approach. New J. Phys. 17, 115008 (2015).
 33.
Kabir, H., Wang, Y., Yu, M. & Zhang, Q. Neural network inverse modeling and applications to microwave filter design. IEEE 56, 867 (2008).
 34.
Gosal, G., Almajali, E., McNamara, D. & Yagoub, M. Transmitarray antenna design using forward and inverse neural network modeling. IEEE, Antennas Wirel. Propag. Lett. 15, 1483 (2016).
 35.
Aoad, A., Simsek, M. & Aydin, Z. Knowledge based response correction method for design of reconfigurable nshaped microstrip patch antenna using inverse anns. Int. J. Numer. Model. Electron. Netw. Devices Fields 30, e2129–e2129 (2017).
 36.
Liu, D., Tan, Y., Khoram, E. & Yu, Z. Training deep neural networks for the inverse design of nanophotonic structures. ACS Photonics 5, 1365–1369 (2018).
 37.
Adler, J. & Öktem, O. Solving illposed inverse problems using iterative deep neural networks. Inverse Probl. 33, 124007 (2017).
 38.
Kartashov, Y. V. & Skryabin, D. V. Bistable topological insulator with excitonpolaritons. Phys. Rev. Lett. 119, 253904 (2017).
 39.
Mihalache, D. et al. Stable topological modes in twodimensional ginzburglandau models with trapping potentials. Phys. Rev. A. 82, 023813 (2010).
 40.
Jünemann, J. et al. Exploring interacting topological insulators with ultracold atoms: The synthetic creutzhubbard model. Phys. Rev. X 7, 031057 (2017).
 41.
Mancini, M. et al. Observation of chiral edge states with neutral fermions in synthetic hall ribbons. Science 349, 1510–1513 (2015).
 42.
Parto, M. et al. Complex edgestate phase transitions in 1d topological laser arrays. Preprint at arXiv:1709.00523 (2017).
 43.
Pilozzi, L. & Conti, C. Topological cascade laser for frequency comb generation in PTsymmetric structures. Opt. Lett. 42, 5174 (2017).
 44.
Longhi, S. Paritytime symmetry meets photonics: a new twist in nonHermitian optics. Preprint at ArXiv (2018).
 45.
Zeuner, J. M. et al. Observation of a topological transition in the bulk of a nonhermitian system. Phys. Rev. Lett. 115, 040402 (2015).
 46.
Graf, G. & Porta, M. Bulkedge correspondence for twodimensional topological insulators. Commun. Math. Phys. 324, 851 (2013).
 47.
Poshakinskiy, A. V., Poddubny, A. N. & Hafezi, M. Phase spectroscopy of topological invariants in photonic crystals. Phys. Rev. A 91, 043830 (2015).
 48.
Graf, G. & Tauber, C. BulkEdge correspondence for twodimensional Floquet topological insulators. Preprint at ArXiv (2017).
Acknowledgements
We acknowledge support from the Templeton foundation (grant number 58277), the PRIN2015 NEMO project (2015KEZNYM grant), the H2020 QuantERA project QUOMPLEX (grant number 731473), and the Italian MAE project NECST. We thank Dr. Alexander Poshakinskiy for the fruitful comments regarding the training dataset generation.
Author information
Affiliations
Contributions
C.C. conceived the initial idea and supervised the project. F.F. expanded the concept and developed the code. L.P., G.M. and C.C. developed the theoretical part. F.F. and L.P. carried out the simulations. F.F., L.P. and G.M. contributed to data analysis and figure preparation. All the authors contributed to manuscript writing.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pilozzi, L., Farrelly, F.A., Marcucci, G. et al. Machine learning inverse problem for topological photonics. Commun Phys 1, 57 (2018). https://doi.org/10.1038/s4200501800588
Received:
Accepted:
Published:
Further reading

Highlighting photonics: looking into the next decade
eLight (2021)

Broadband vectorial ultrathin optics with experimental efficiency up to 99% in the visible region via universal approximators
Light: Science & Applications (2021)

Deep learning for the design of photonic structures
Nature Photonics (2021)

Inferring the dynamics of “blackbox” systems using a learning machine
Science China Physics, Mechanics & Astronomy (2021)

Recent advances in 2D, 3D and higherorder topological photonics
Light: Science & Applications (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.