Solvent Induced Shifts in the UV Spectrum of Amides

Solvent effects on the electronic spectra of formamide and trans-N-methylacetamide are studied using four different levels of theory: singly excited configuration interaction (CIS), equations of motion coupled-cluster theory with singles and doubles (EOM-CCSD), completely renormalized coupled-cluster theory with singles and doubles with perturbative triple excitations (CR-EOM-CCSD(T)), and time-dependent density functional theory (TDDFT), employing small clusters of water molecules. The simulated electronic spectrum is obtained via molecular dynamics simulations with 100 waters modeled with the effective fragment potential method and exhibits a blue-shift and red-shift, respectively, for the n → π* and πnb → π* vertical excitation energies, in good agreement with the experimental electronic spectra of amides

Ionization energies in solution with the QM : QM approach

We discuss a fragment-based QM:QM scheme as a practical way to access the energetics of vertical electronic processes in the condensed phase. In the QM:QM scheme, we decompose the large molecular system into small fragments, which interact solely electrostatically. The energies of the fragments are calculated in a self-consistent field generated by the other fragments and the total energy of the system is calculated as a sum of the fragment energies. We show on two test cases (cytosine and a sodium cation) that the method allows one to accurately simulate the shift of vertical ionization energies (VIE) while going from the gas phase to the bulk. For both examples, the predicted solvent shifts and peak widths estimated at the DFT level agree well with the experimental observations. We argue that the QM:QM approach is more suitable than either an electrostatic embedding based QM/MM approach, a full quantum description at the DFT level with a generally used functional or a combination of both. We also discuss the potential scope of the applicability for other electronic processes such as Auger decay.Peer reviewe

Polarizable Density Embedding for Large Biomolecular Systems

This document is the Accepted Manuscript version of a Published Work that appeared in final form in Journal of Chemical Theory and Computation, Copyright © 2020 American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see https://doi.org/10.1021/acs.jctc.0c00763.We present an efficient and robust fragment-based quantum–classical embedding model capable of accurately capturing effects from complex environments such as proteins and nucleic acids. This is realized by combining the molecular fractionation with conjugate caps (MFCC) procedure with the polarizable density embedding (PDE) model at the level of Fock matrix construction. The PDE contributions to the Fock matrix of the core region are constructed using the local molecular basis of the individual fragments rather than the supermolecular basis of the entire system. Thereby, we avoid complications associated with the application of the MFCC procedure on environment quantities such as electronic densities and molecular-orbital energies. Moreover, the computational cost associated with solving self-consistent field (SCF) equations of the core region remains unchanged from that of purely classical polarized embedding models. We analyze the performance of the resulting model in terms of the reproduction of the electrostatic potential of an insulin monomer protein and further in the context of solving problems related to electron spill-out. Finally, we showcase the model for the calculation of one- and two-photon properties of the Nile red molecule in a protein environment. Based on our analyses, we find that the combination of the MFCC approach with the PDE model is an efficient, yet accurate approach for calculating molecular properties of molecules embedded in structured biomolecular environments

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid phase environment are significantly renormalized by the surrounding solvent molecules. Herein, we employ the many-body Green's function approach (in the GW approximation) to investigate the solvation effects on the ionization energy of phenol in various solvents with distinct polarizability. The many-body effects among the investigated solvents differ by up to 0.4 eV, and this difference is not simply owed to the macroscopic solvent polarizability. Utilizing orbital localization and projection, we define an electronic subspace for a fragment, i.e., a solvation shell, in the solvent environment. The resulting fragment correlation self-energy is shown to decay with respect to the intermolecular distance and vanish at ~9 angstroms. This decaying pattern is independent of the ionization state and the solvent type. The 9-angstrom cut-off distance defines an effective interacting volume, within which we find the quasiparticle energy shift per solvent molecule is directly related to the polarizability of the solvent molecules. Finally, we propose a simple solvation model for computing the quasiparticle energies of solvated systems

Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions

We present the extension of the quantum/classical polarizable fluctuating charge model to the calculation of single residues of quadratic response functions, as required for the computational modeling of two-photon absorption cross-sections. By virtue of a variational formulation of the quantum/classical polarizable coupling, we are able to exploit an atomic orbital-based quasienergy formalism to derive the additional coupling terms in the response equations. Our formalism can be extended to the calculation of arbitrary order response functions and their residues. The approach has been applied to the challenging problem of one- and two-photon spectra of rhodamine 6G (R6G) in aqueous solution. Solvent effects on one- and two-photon spectra of R6G in aqueous solution have been analyzed by considering three different approaches, from a continuum (QM/PCM) to two QM/MM models (non-polarizable QM/TIP3P and polarizable QM/FQ). Both QM/TIP3P and QM/FQ simulated OPA and TPA spectra show that the inclusion of discrete water solvent molecules is essential to increase the agreement between theory and experiment. QM/FQ has been shown to give the best agreement with experiments

Multiscale Models for Light-Driven Processes

Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes