Doubly resonant photonic crystal cavities for efficient second-harmonic generation in III-V semiconductors

Second-order nonlinear effects, such as second-harmonic generation, can be strongly enhanced in nanofabricated photonic materials when both fundamental and harmonic frequencies are spatially and temporally confined. Practically designing low-volume and doubly resonant nanoresonators in conventional semiconductor compounds is challenging owing to their intrinsic refractive index dispersion. In this work we review a recently developed strategy to design doubly resonant nanocavities with low mode volume and large quality factor by localized defects in a photonic crystal structure. We build on this approach by applying an evolutionary optimisation algorithm in connection with Maxwell equations solvers, showing that the proposed design recipe can be applied to any material platform. We explicitly calculate the second-harmonic generation efficiency for doubly resonant photonic crystal cavity designs in typical III-V semiconductor materials, such as GaN and AlGaAs, targeting a fundamental harmonic at telecom wavelengths, and fully accounting for the tensor nature of the respective nonlinear susceptibilities. These results may stimulate the realisation of small footprint photonic nanostructures in leading semiconductor material platforms to achieve unprecedented nonlinear efficiencies.

Doubly resonant second-harmonic generation of a vortex beam from a bound state in the continuum

Second harmonic generation in nonlinear materials can be greatly enhanced by realizing doubly-resonant cavities with high quality factors. However, fulfilling such doubly resonant condition in photonic crystal (PhC) cavities is a long-standing challenge, because of the difficulty in engineering photonic bandgaps around both frequencies. Here, by implementing a second-harmonic bound state in the continuum (BIC) and confining it with a heterostructure design, we show the first doubly-resonant PhC slab cavity with $2.4times10^{-2}$ W$^{-1}$ conversion efficiency under continuous wave excitation. We also report the confirmation of highly normal-direction concentrated far-field emission pattern with radial polarization at the second harmonic frequency. These results represent a solid verification of previous theoretical predictions and a cornerstone achievement, not only for nonlinear frequency conversion but also for vortex beam generation and prospective nonclassical sources of radiation.

Inverse design of photonic crystals through automatic differentiation

Gradient-based inverse design in photonics has already achieved remarkable results in designing small-footprint, high-performance optical devices. The adjoint variable method, which allows for the efficient computation of gradients, has played a major role in this success. However, gradient-based optimization has not yet been applied to the mode-expansion methods that are the most common approach to studying periodic optical structures like photonic crystals. This is because, in such simulations, the adjoint variable method cannot be defined as explicitly as in standard finite-difference or finite-element time- or frequency-domain methods. Here, we overcome this through the use of automatic differentiation, which is a generalization of the adjoint variable method to arbitrary computational graphs. We implement the plane-wave expansion and the guided-mode expansion methods using an automatic differentiation library, and show that the gradient of any simulation output can be computed efficiently and in parallel with respect to all input parameters. We then use this implementation to optimize the dispersion of a photonic crystal waveguide, and the quality factor of an ultra-small cavity in a lithium niobate slab. This extends photonic inverse design to a whole new class of simulations, and more broadly highlights the importance that automatic differentiation could play in the future for tracking and optimizing complicated physical models.

Higher-order topological insulators in synthetic dimensions

Conventional topological insulators support boundary states that have one dimension lower than the bulk system that hosts them, and these states are topologically protected due to quantized bulk dipole moments. Recently, higher-order topological insulators have been proposed as a way of realizing topological states that are two or more dimensions lower than the bulk, due to the quantization of bulk quadrupole or octupole moments. However, all these proposals as well as experimental realizations have been restricted to real-space dimensions. Here we construct photonic higher-order topological insulators (PHOTI) in synthetic dimensions. We show the emergence of a quadrupole PHOTI supporting topologically protected corner modes in an array of modulated photonic molecules with a synthetic frequency dimension, where each photonic molecule comprises two coupled rings. By changing the phase difference of the modulation between adjacently coupled photonic molecules, we predict a dynamical topological phase transition in the PHOTI. Furthermore, we show that the concept of synthetic dimensions can be exploited to realize even higher-order multipole moments such as a 4th order hexadecapole (16-pole) insulator, supporting 0D corner modes in a 4D hypercubic synthetic lattice that cannot be realized in real-space lattices.

Forward-Mode Differentiation of Maxwells Equations

We present a previously unexplored forward-mode differentiation method for Maxwells equations, with applications in the field of sensitivity analysis. This approach yields exact gradients and is similar to the popular adjoint variable method, but provides a significant improvement in both memory and speed scaling for problems involving several output parameters, as we analyze in the context of finite-difference time-domain (FDTD) simulations. Furthermore, it provides an exact alternative to numerical derivative methods, based on finite-difference approximations. To demonstrate the usefulness of the method, we perform sensitivity analysis of two problems. First we compute how the spatial near-field intensity distribution of a scatterer changes with respect to its dielectric constant. Then, we compute how the spectral power and coupling efficiency of a surface grating coupler changes with respect to its fill factor.

Experimental demonstration of dynamical input isolation in nonadiabatically modulated photonic cavities

Modulated optical cavities have been proposed and demonstrated for applications in communications, laser frequency stabilization, microwave-to-optical conversion and frequency comb generation. However, most studies are restricted to the adiabatic regime, where either the maximum excursion of the modulation or the modulation frequency itself is below the linewidth of the cavity. Here, using a fiber ring resonator with an embedded electro-optic phase modulator, we investigate the nonadiabatic regime. By strongly driving the modulator at frequencies that are significantly smaller than the free-spectral range of the ring resonator, but well beyond the linewidth of the resonator, we experimentally observe counterintuitive behavior predicted in a recent theoretical study by Minkov et al. [APL Photonics 2, 076101 (2017)], such as the complete suppression of drop-port transmission even when the input laser wavelength is on resonance with the optical cavity. This can be understood as dynamical isolation of the cavity from the input light. We also show qualitative differences in the steady-state responses of the system between the adiabatic and nonadiabatic limits. Our experiments probe a seldom explored regime of operation that is promising for applications in integrated photonic systems with current state-of-the-art technology.

Adjoint method and inverse design for nonlinear nanophotonic devices

The development of inverse design, where computational optimization techniques are used to design devices based on certain specifications, has led to the discovery of many compact, non-intuitive structures with superior performance. Among various methods, large-scale, gradient-based optimization techniques have been one of the most important ways to design a structure containing a vast number of degrees of freedom. These techniques are made possible by the adjoint method, in which the gradient of an objective function with respect to all design degrees of freedom can be computed using only two full-field simulations. However, this approach has so far mostly been applied to linear photonic devices. Here, we present an extension of this method to modeling nonlinear devices in the frequency domain, with the nonlinear response directly included in the gradient computation. As illustrations, we use the method to devise compact photonic switches in a Kerr nonlinear material, in which low-power and high-power pulses are routed in different directions. Our technique may lead to the development of novel compact nonlinear photonic devices.

Training of photonic neural networks through in situ backpropagation

Recently, integrated optics has gained interest as a hardware platform for implementing machine learning algorithms. Of particular interest are artificial neural networks, since matrix-vector multi- plications, which are used heavily in artificial neural networks, can be done efficiently in photonic circuits. The training of an artificial neural network is a crucial step in its application. However, currently on the integrated photonics platform there is no efficient protocol for the training of these networks. In this work, we introduce a method that enables highly efficient, in situ training of a photonic neural network. We use adjoint variable methods to derive the photonic analogue of the backpropagation algorithm, which is the standard method for computing gradients of conventional neural networks. We further show how these gradients may be obtained exactly by performing intensity measurements within the device. As an application, we demonstrate the training of a numerically simulated photonic artificial neural network. Beyond the training of photonic machine learning implementations, our method may also be of broad interest to experimental sensitivity analysis of photonic systems and the optimization of reconfigurable optics platforms.

p-q-superstrings in Anti-de-Sitter space-time

We study a field theoretical model for p-q-superstrings in a fixed Anti-de-Sitter background. We find that the presence of the negative cosmological constant tends to decrease the core radius of the strings. Moreover, the binding energy decreases with the increase of the absolute value of the cosmological constant. Studying the effect of the p-q-strings on Anti-de-Sitter space, we observe that the presence of the negative cosmological constant tends to decrease the deficit angle as compared to asymptotically flat space-time.