Vibronic and Environmental Effects in Simulations of Optical Spectroscopy.

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy.

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore-environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore-environment configurations with Franck-Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Employing an average Franck-Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores.

In the condensed phase, ions often create heterogeneous local environments around a solute, which may impart chemical reactivity or perturbations to physicochemical properties. Although the former has been the subject of some study, the latter-particularly as is pertains to optical absorption spectroscopy-is much less understood. In this work, the computed UV-vis absorption spectrum is examined for the aqueously solvated chromophore anion of green fluorescent protein for different local ion configurations. The strong ability of water to screen the ions from the chromophore results in little change in excitation energy compared to a purely aqueous environment. However, upon forming a contact ion pair with a sodium ion at either of the two electronegative oxygen sites of the chromophore, there is a spectral shift to either higher or lower energies. Surprisingly, our analysis suggests that the cause of the spectral shift is dominated not by the electrostatic presence of the ion but instead by ion disruption of the hydrogen bond network at the oxygen contact ion pair site.

Nonuniform continuum model for solvatochromism based on frozen-density embedding theory.

Frozen-density embedding theory (FDET) provides the formal framework for multilevel numerical simulations, such that a selected subsystem is described at the quantum mechanical level, whereas its environment is described by means of the electron density (frozen density; ${\rho _{\rm{B}} (\vec r)}$). The frozen density ${\rho _{\rm{B}} (\vec r)}$ is usually obtained from some lower-level quantum mechanical methods applied to the environment, but FDET is not limited to such choices for ${\rho _{\rm{B}} (\vec r)}$. The present work concerns the application of FDET, in which ${\rho _{\rm{B}} (\vec r)}$ is the statistically averaged electron density of the solvent ${\left\langle {\rho _{\rm{B}} (\vec r)} \right\rangle }$. The specific solute-solvent interactions are represented in a statistical manner in ${\left\langle {\rho _{\rm{B}} (\vec r)} \right\rangle }$. A full self-consistent treatment of solvated chromophore, thus involves a single geometry of the chromophore in a given state and the corresponding ${\left\langle {\rho _{\rm{B}} (\vec r)} \right\rangle }$. We show that the coupling between the two descriptors might be made in an approximate manner that is applicable for both absorption and emission. The proposed protocol leads to accurate (error in the range of 0.05 eV) descriptions of the solvatochromic shifts in both absorption and emission.

Non-uniform Continuum Model for Solvated Species Based on Frozen-Density Embedding Theory: The Study Case of Solvatochromism of Coumarin 153.

Recent application of the Frozen-Density Embedding Theory based continuum model of the solvent, which is used for calculating solvatochromic shifts in the UV/Vis range, are reviewed. In this model, the solvent is represented as a non-uniform continuum taking into account both the statistical nature of the solvent and specific solute-solvent interactions. It offers, therefore, a computationally attractive alternative to methods in which the solvent is described at atomistic level. The evaluation of the solvatochromic shift involves only two calculations of excitation energy instead of at least hundreds needed to account for inhomogeneous broadening. The present review provides a detailed graphical analysis of the key quantities of this model: the average charge density of the solvent (<ρB>) and the corresponding Frozen-Density Embedding Theory derived embedding potential for coumarin 153.

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol.

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. The strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

Critical study of the charge transfer parameter for the calculation of interaction energy using the local hard-soft acid-base principle.

Local hard-soft acid-base (HSAB) principle is semiquantitative in nature due to the presence of an ad hoc charge transfer parameter. The accuracy of HSAB principle significantly depends on the definition of this ad hoc parameter. In this paper, for the first time we have introduced the second-order approximation of ΔN (ΔNsecond) as an ad hoc parameter for charge transfer to calculate interaction energies of multiple site based interactions using local hard soft acid base principle. The second-order approximation of ΔN has been derived from Sanderson's electronegativity equalization principle. To validate our approach, we have studied interaction energies of some prototype molecules. The interaction energies obtained from our approach have been further compared with the interaction energies of those obtained using other charge transfer parameters (ΔNfirst and λ) and the conventional methods. We have also discussed the advantages and limitations of the approach.

Molecular polariton electroabsorption.

We investigate electroabsorption (EA) in organic semiconductor microcavities to understand whether strong light-matter coupling non-trivially alters their nonlinear optical [[Formula: see text]] response. Focusing on strongly-absorbing squaraine (SQ) molecules dispersed in a wide-gap host matrix, we find that classical transfer matrix modeling accurately captures the EA response of low concentration SQ microcavities with a vacuum Rabi splitting of [Formula: see text] meV, but fails for high concentration cavities with [Formula: see text] meV. Rather than new physics in the ultrastrong coupling regime, however, we attribute the discrepancy at high SQ concentration to a nearly dark H-aggregate state below the SQ exciton transition, which goes undetected in the optical constant dispersion on which the transfer matrix model is based, but nonetheless interacts with and enhances the EA response of the lower polariton mode. These results indicate that strong coupling can be used to manipulate EA (and presumably other optical nonlinearities) from organic microcavities by controlling the energy of polariton modes relative to other states in the system, but it does not alter the intrinsic optical nonlinearity of the organic semiconductor inside the cavity.

Comparison of the auxiliary density perturbation theory and the noniterative approximation to the coupled perturbed Kohn-Sham method: case study of the polarizabilities of disubstituted azoarene molecules.

We present a theoretical study of the polarizabilities of free and disubstituted azoarenes employing auxiliary density perturbation theory (ADPT) and the noniterative approximation to the coupled perturbed Kohn-Sham (NIA-CPKS) method. Both methods are noniterative but use different approaches to obtain the perturbed density matrix. NIA-CPKS is different from the conventional CPKS approach in that the perturbed Kohn-Sham matrix is obtained numerically, thereby yielding a single-step solution to CPKS. ADPT is an alternative approach to the analytical CPKS method in the framework of the auxiliary density functional theory. It is shown that the polarizabilities obtained using these two methods are in good agreement with each other. Comparisons are made for disubstituted azoarenes, which give support to the push-pull mechanism. Both methods reproduce the same trend for polarizabilities because of the substitution pattern of the azoarene moiety. Our results are consistent with the standard organic chemistry "activating/deactivating" sequence. We present the polarizabilities of the above molecules calculated with three different exchange-correlation functionals and two different auxiliary function sets. The computational advantages of both methods are also discussed.

Excitonic Hamiltonians for Calculating Optical Absorption Spectra and Optoelectronic Properties of Molecular Aggregates and Solids.

Rational design of disordered molecular aggregates and solids for optoelectronic applications relies on our ability to predict the properties of such materials using theoretical and computational methods. However, large molecular systems where disorder is too significant to be considered in the perturbative limit cannot be described using either first principles quantum chemistry or band theory. Multiscale modeling is a promising approach to understanding and optimizing the optoelectronic properties of such systems. It uses first-principles quantum chemical methods to calculate the properties of individual molecules, then constructs model Hamiltonians of molecular aggregates or bulk materials based on these calculations. In this paper, we present a protocol for constructing a tight-binding Hamiltonian that represents the excited states of a molecular material in the basis of Frenckel excitons: electron-hole pairs that are localized on individual molecules that make up the material. The Hamiltonian parametrization proposed here accounts for excitonic couplings between molecules, as well as for electrostatic polarization of the electron density on a molecule by the charge distribution on surrounding molecules. Such model Hamiltonians can be used to calculate optical absorption spectra and other optoelectronic properties of molecular aggregates and solids.

Noniterative density functional response approach: application to nonlinear optical properties of p-nitroaniline and its methyl-substituted derivatives.

We report the effect of substitution, position of the substituent, and the symmetry on the nonlinear optical properties of p-nitroanline (PNA) and its derivatives using our implementation of the noniterative approximation of couple-perturbed Kohn-Sham (CPKS) equation in the deMon2k. Dipole moment, static polarizability, and first hyperpolarizability of these pi-conjugated donor-acceptor organic derivatives of PNA and its methyl-substituted analogs are calculated using our method at different exchange correlation functionals, namely, BP86, BPW91, and BLYP, using 6-31++G** basis set. A comparison of results obtained by our method with those obtained by MP2 (finite-field perturbation) method is presented in this paper. The effect of optical gap on charge transfer and subsequently on polarizabilities has been illustrated.