Analog photonic solutions offer unique opportunities to address complex computational tasks with unprecedented performance in terms of energy dissipation and speeds, overcoming current limitations of modern computing architectures based on electron flows and digital approaches. The lack of modularization and lumped element reconfigurability in photonics has prevented the transition to an alloptical analog computing platform. Here, we explore, using numerical simulation, a nanophotonic platform based on epsilonnearzero materials capable of solving in the analog domain partial differential equations (PDE). Wavelength stretching in zeroindex media enables highly nonlocal interactions within the board based on the conduction of electric displacement, which can be monitored to extract the solution of a broad class of PDE problems. By exploiting the experimentally achieved control of deposition technique through process parameters, used in our simulations, we demonstrate the possibility of implementing the proposed nanooptic processor using CMOScompatible indiumtinoxide, whose optical properties can be tuned by carrier injection to obtain programmability at high speeds and low energy requirements. Our nanooptical analog processor can be integrated at chipscale, processing arbitrary inputs at the speed of light.
 Award ID(s):
 1748294
 NSFPAR ID:
 10295749
 Date Published:
 Journal Name:
 Nanophotonics
 Volume:
 10
 Issue:
 6
 ISSN:
 21928606
 Page Range / eLocation ID:
 1711 to 1721
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Abstract 
Finch, a domain specific language and code generation framework for partial differential equations (PDEs), is demonstrated here to solve two classical problems: steadystate advection diffusion equation (single PDE) and the phonon Boltzmann transport equation (coupled PDEs). Both finite volume and finite element methods are explored. In addition to work presented at the 2022 International Conference on Computational Science (Heisler et al., 2022), we include recent developments for solving nonlinear equations using both automatic and symbolic differentiation, and demonstrate the capability for the Bratu (nonlinear Poisson) equation.more » « less

Abstract Three control problems for the system of two coupled differential equations governing the dynamics of an energy harvesting model are studied. The system consists of the equation of an Euler–Bernoulli beam model and the equation representing the Kirchhoff's electric circuit law. Both equations contain coupling terms representing the inverse and direct piezoelectric effects. The system is reformulated as a single evolution equation in the state space of 3‐component functions. The control is introduced as a separable forcing term on the right‐hand side of the operator equation. The first control problem deals with an explicit construction of that steers an initial state to zero on a time interval [0,
T ]. The second control problem deals with the construction of such that the voltage output is equal to some given function (with being given as well). The third control problem deals with an explicit construction of both the force profile, , and the control, , which generate the desired voltage output . Interpolation theory in the Hardy space of analytic functions is used in the solution of the second and third problems. 
null ; null ; null ; null (Ed.)Many imaging problems can be formulated as inverse problems expressed as finitedimensional optimization problems. These optimization problems generally consist of minimizing the sum of a data fidelity and regularization terms. In Darbon (SIAM J. Imag. Sci. 8:2268–2293, 2015), Darbon and Meng, (On decomposition models in imaging sciences and multitime HamiltonJacobi partial differential equations, arXiv preprint arXiv:1906.09502, 2019), connections between these optimization problems and (multitime) HamiltonJacobi partial differential equations have been proposed under the convexity assumptions of both the data fidelity and regularization terms. In particular, under these convexity assumptions, some representation formulas for a minimizer can be obtained. From a Bayesian perspective, such a minimizer can be seen as a maximum a posteriori estimator. In this chapter, we consider a certain class of nonconvex regularizations and show that similar representation formulas for the minimizer can also be obtained. This is achieved by leveraging minplus algebra techniques that have been originally developed for solving certain HamiltonJacobi partial differential equations arising in optimal control. Note that connections between viscous HamiltonJacobi partial differential equations and Bayesian posterior mean estimators with Gaussian data fidelity terms and logconcave priors have been highlighted in Darbon and Langlois, (On Bayesian posterior mean estimators in imaging sciences and HamiltonJacobi partial differential equations, arXiv preprint arXiv:2003.05572, 2020). We also present similar results for certain Bayesian posterior mean estimators with Gaussian data fidelity and certain nonlogconcave priors using an analogue of minplus algebra techniques.more » « less

Abstract. Geoscientific models are facing increasing challenges to exploit growing datasets coming from remote sensing. Universal differential equations (UDEs), aided by differentiable programming, provide a new scientific modelling paradigm enabling both complex functional inversions to potentially discover new physical laws and data assimilation from heterogeneous and sparse observations. We demonstrate an application of UDEs as a proof of concept to learn the creep component of ice flow, i.e. a nonlinear diffusivity differential equation, of a glacier evolution model. By combining a mechanistic model based on a twodimensional shallowice approximation partial differential equation with an embedded neural network, i.e. a UDE, we can learn parts of an equation as nonlinear functions that then can be translated into mathematical expressions. We implemented this modelling framework as ODINN.jl, a package in the Julia programming language, providing high performance, sourcetosource automatic differentiation (AD) and seamless integration with tools and global datasets from the Open Global Glacier Model in Python. We demonstrate this concept for 17 different glaciers around the world, for which we successfully recover a prescribed artificial law describing ice creep variability by solving ∼ 500 000 ordinary differential equations in parallel. Furthermore, we investigate which are the best tools in the scientific machine learning ecosystem in Julia to differentiate and optimize large nonlinear diffusivity UDEs. This study represents a proof of concept for a new modelling framework aiming at discovering empirical laws for largescale glacier processes, such as the variability in ice creep and basal sliding for ice flow, and new hybrid surface mass balance models.