Cavity-enhanced second harmonic generation via nonlinear-overlap optimization

We describe an approach based on topology optimization that enables automatic discovery of wavelength-scale photonic structures for achieving high-efficiency second-harmonic generation (SHG). A key distinction from previous formulation and designs that seek to maximize Purcell factors at individual frequencies is that our method not only aims to achieve frequency matching (across an entire octave) and large radiative lifetimes, but also optimizes the equally important nonlinear--coupling figure of merit $\bar{\beta}$, involving a complicated spatial overlap-integral between modes. We apply this method to the particular problem of optimizing micropost and grating-slab cavities (one-dimensional multilayered structures) and demonstrate that a variety of material platforms can support modes with the requisite frequencies, large lifetimes $Q > 10^4$, small modal volumes $\sim (\lambda/n)^3$, and extremely large $\bar{\beta} \gtrsim 10^{-2}$, leading to orders of magnitude enhancements in SHG efficiency compared to state of the art photonic designs. Such giant $\bar{\beta}$ alleviate the need for ultra-narrow linewidths and thus pave the way for wavelength-scale SHG devices with faster operating timescales and higher tolerance to fabrication imperfections

Efficient construction of maximally localized photonic Wannier functions: locality criterion and initial conditions

Wannier function expansions are well suited for the description of photonic- crystal-based defect structures, but constructing maximally localized Wannier functions by optimizing the phase degree of freedom of the Bloch modes is crucial for the efficiency of the approach. We systematically analyze different locality criteria for maximally localized Wannier functions in two- dimensional square and triangular lattice photonic crystals, employing (local) conjugate-gradient as well as (global) genetic-algorithm-based, stochastic methods. Besides the commonly used second moment (SM) locality measure, we introduce a new locality measure, namely the integrated modulus (IM) of the Wannier function. We show numerically that, in contrast to the SM criterion, the IM criterion leads to an optimization problem with a single extremum, thus allowing for fast and efficient construction of maximally localized Wannier functions using local optimization techniques. We also present an analytical formula for the initial choice of Bloch phases, which under certain conditions represents the global maximum of the IM criterion and, thus, further increases the optimization efficiency in the general case

Topology optimization of freeform large-area metasurfaces

We demonstrate optimization of optical metasurfaces over $10^5$--$10^6$ degrees of freedom in two and three dimensions, 100--1000+ wavelengths ($\lambda$) in diameter, with 100+ parameters per $\lambda^2$. In particular, we show how topology optimization, with one degree of freedom per high-resolution "pixel," can be extended to large areas with the help of a locally periodic approximation that was previously only used for a few parameters per $\lambda^2$. In this way, we can computationally discover completely unexpected metasurface designs for challenging multi-frequency, multi-angle problems, including designs for fully coupled multi-layer structures with arbitrary per-layer patterns. Unlike typical metasurface designs based on subwavelength unit cells, our approach can discover both sub- and supra-wavelength patterns and can obtain both the near and far fields

Coupling techniques between dielectric waveguides and planar photonic crystals

El objetivo de esta tesis es la investigación de estructuras y técnicas de acoplo para minimizar las pérdidas de acoplo entre guías dieléctricas y cristales fotónicos planares. En primer lugar se ha estudiado el modelado del acoplo entre guías dieléctricas y guías en cristal fotónico así como la influencia de los principales parámetros del cristal en la eficiencia de acoplo. Se han obtenido expresiones cerradas para las matrices de reflexión y transmisión que caracterizan totalmente el scattering que ocurre en el interfaz formado entre una guía dieléctrica y una guía en cristal fotónico. A continuación y con el fin de mejorar la eficiencia de acoplo desde guías dieléctrica de anchura arbitraria, se ha propuesto como contribución original una técnica de acoplo basada en la introducción de defectos puntuales en el interior de una estructura de acoplo tipo cuña realizada en el cristal fotónico. Diferentes soluciones, incluida los algoritmos genéticos, han sido propuestas con el objetivo de conseguir el diseño óptimo de la configuración de defectos. Una vez conseguido un acoplo eficiente desde guías dieléctricas a guías en cristal fotónico, se ha investigado el acoplo en guías de cavidades acopladas. Como contribución original se ha propuesto una técnica de acoplo basada en la variación gradual del radio de los defectos situados entre cavidades adyacentes. Además, se ha realizado un riguroso análisis en el dominio del tiempo y la frecuencia de la propagación de pulsos en guías acopladas de longitud finita. Dicho estudio ha tenido como objetivo la caracterización de la influencia de la eficiencia del acoplo en los parámetros del pulso. Finalmente, se han presentado los procesos de fabricación y resultados experimentales de las estructuras de acoplo propuestas.Sanchis Kilders, P. (2005). Coupling techniques between dielectric waveguides and planar photonic crystals [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/1854Palanci

Fast Numerical Methods for the Design of Layered Photonic Structures with Rough Interfaces

Modified boundary conditions (MBC) and a multilayer approach (MA) are proposed as fast and efficient numerical methods for the design of 1D photonic structures with rough interfaces. These methods are applicable for the structures, composed of materials with arbitrary permittivity tensor. MBC and MA are numerically validated on different types of interface roughness and permittivities of the constituent materials. The proposed methods can be combined with the 4x4 scattering matrix method as a field solver and an evolutionary strategy as an optimizer. The resulted optimization procedure is fast, accurate, numerically stable and can be used to design structures for various applications

Fast, Scalable, and Interactive Software for Landau-de Gennes Numerical Modeling of Nematic Topological Defects

Numerical modeling of nematic liquid crystals using the tensorial Landau-de Gennes (LdG) theory provides detailed insights into the structure and energetics of the enormous variety of possible topological defect configurations that may arise when the liquid crystal is in contact with colloidal inclusions or structured boundaries. However, these methods can be computationally expensive, making it challenging to predict (meta)stable configurations involving several colloidal particles, and they are often restricted to system sizes well below the experimental scale. Here we present an open-source software package that exploits the embarrassingly parallel structure of the lattice discretization of the LdG approach. Our implementation, combining CUDA/C++ and OpenMPI, allows users to accelerate simulations using both CPU and GPU resources in either single- or multiple-core configurations. We make use of an efficient minimization algorithm, the Fast Inertial Relaxation Engine (FIRE) method, that is well-suited to large-scale parallelization, requiring little additional memory or computational cost while offering performance competitive with other commonly used methods. In multi-core operation we are able to scale simulations up to supra-micron length scales of experimental relevance, and in single-core operation the simulation package includes a user-friendly GUI environment for rapid prototyping of interfacial features and the multifarious defect states they can promote. To demonstrate this software package, we examine in detail the competition between curvilinear disclinations and point-like hedgehog defects as size scale, material properties, and geometric features are varied. We also study the effects of an interface patterned with an array of topological point-defects.Comment: 16 pages, 6 figures, 1 youtube link. The full catastroph

Inverse Design of Photonic Crystal Devices

This work deals with the inverse design in the field of photonic crystal based devices. Here an inverse method containing a fast and accurate simulation method integrated with a competent optimization method is presented. Two designs yielded from this conjunction of multiple scattering theory with a genetic algorithm is analyzed. The potential of this approach is illustrated by designing a lens that has a very low F-number (F=0.47) and a conversion ratio of 11:1. We have also designed a coupler device that introduces the light from an optical fiber into a PC based wave-guide with a predicted coupling efficiency that exceeds 87%.Comment: 8 pages with figure