Nonlinear optics is essential for many recent photonic technologies. Here, a novel multi-scale approach is introduced to simulate the nonlinear optical response of molecular nanomaterials combining ab initio quantum-chemical and classical Maxwell-scattering computations. In this approach, the first hyperpolarizability tensor is computed with time-dependent density-functional theory and incorporated into a multi-scattering formalism that considers the optical interaction between neighboring molecules. Such incorporation is achieved by a novel object: the Hyper-Transition(T)-matrix. With this object at hand, the nonlinear optical response from single molecules and also from entire photonic devices can be computed, including the full tensorial and dispersive nature of the optical response of the molecules, as well as the optical interaction between different molecules as, for example, in the lattice of a molecular crystal. To demonstrate the applicability of the novel approach, the generation of a second-harmonic signal from a thin film of an Urea molecular crystal is computed and compared to more traditional simulations. Furthermore, an optical cavity is designed, which enhances the second-harmonic response of the molecular film up to more than two orders of magnitude. This approach is highly versatile and accurate and can be the working horse for the future exploration of nonlinear photonic molecular materials in structured photonic environments.
Keywords: multi-scale workflow; nonlinear optics; second-harmonic generation; time-dependent density-functional theory; transition matrix formalism.
© 2023 The Authors. Advanced Materials published by Wiley-VCH GmbH.