Abstract:CONSPECTUS: Organic glasses containing chromophores with large first hyperpolarizabilities (β) are promising for compact, high-bandwidth, and energy-efficient electro-optic devices. Systematic optimization of device performance requires development of materials with high acentric order and enhanced hyperpolarizability at operating wavelengths. One essential component of the design process is the accurate calculation of optical transition frequencies and hyperpolarizability. These properties can be computed with a wide range of electronic structure methods implemented within commercial and open-source software packages. A wide variety of methods, especially hybrid density-functional theory (DFT) variants have been used for this purpose. However, in order to provide predictions useful to chromophore designers, a method must be able to consistently predict the relative ordering of standard and novel materials. Moreover, it is important to distinguish between the resonant and nonresonant contribution to the hyperpolarizabiliy and be able to estimate the trade-off between improved β and unwanted absorbance (optical loss) at the target device's operating wavelength. Therefore, we have surveyed a large variety of common methods for computing the properties of modern high-performance chromophores and compared these results with prior experimental hyper-Rayleigh scattering (HRS) and absorbance data. We focused on hybrid DFT methods, supplemented by more computationally intensive Møller-Plesset (MP2) calculations, to determine the relative accuracy of these methods. Our work compares computed hyperpolarizabilities in chloroform relative to standard chromophore EZ-FTC against HRS data versus the same reference. We categorized DFT methods used by the amount of Hartree-Fock (HF) exchange energy incorporated into each functional. Our results suggest that the relationship between percentage of long-range HF exchange and both βHRS and λmax is nearly linear, decreasing as the fraction of long-range HF exchange increases. Mild hybrid DFT methods are satisfactory for prediction of λmax. However, mild hybrid methods provided qualitatively incorrect predictions of the relative hyperpolarizabilities of three high-performance chromophores. DFT methods with approximately 50% HF exchange, and especially the Truhlar M062X functional, provide superior predictions of relative βHRS values but poorer predictions of λmax. The observed trends for these functionals, as well as range-separated hybrids, are similar to MP2, though predicting smaller absolute magnitudes for βHRS. Frequency dependence for βHRS can be calculated using time-dependent DFT and HF methods. However, calculation quality is sensitive not only to a method's ability to predict static hyperpolarizability but also to its prediction of optical resonances. Due to the apparent trade-off in accuracy of prediction of these two properties and the need to use static finite-field methods for MP2 and higher-level hyperpolarizability calculations in most codes, we suggest that composite methods could greatly improve the accuracy of calculations of β and λmax.

Application to nonlinear optical properties of the RSX‐QIDH double‐hybrid range‐separated functional

Non-empirical Double-Hybrid Density Functionals as Reliable Tools for Electronic Structure Calculations

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

RPA, an accurate and fast method for the computation of static non-linear optical properties

Range-separated double-hybrid density-functional theory applied to periodic systems

Assessment of range-separated exchange functionals and nonempirical functional tuning for calculating the static second hyperpolarizabilities of streptocyanines.

A general range-separated double-hybrid density-functional theory

Local approach to coupled cluster evaluation of polarizabilities for long conjugated molecules.

Inverse molecular design and parameter optimization with Hückel theory using automatic differentiation

Range-separated density-functional theory with random phase approximation applied to noncovalent intermolecular interactions

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

Optimizing calculations of electronic excitations and relative hyperpolarizabilities of electrooptic chromophores

Transferable screened range-separated hybrid functionals for electronic and optical properties of van der Waals materials

Reaching a Uniform Accuracy for Complex Molecular Systems: Long-Range-Corrected XYG3 Doubly Hybrid Density Functional

Accurate Electronic and Optical Properties of Organic Doublet Radicals Using Machine Learned Range-Separated Functionals

Which options exist for NISQ-friendly linear response formulations?

Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS.

Towards polaritonic molecular orbitals for large molecular systems

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Symmetry Dilemma of Doubly Hybrid Density Functionals for Equilibrium Molecular Property Calculations.