Solvation Effects on Electronic Transitions: Exploring the Performance of Advanced Solvent Potentials in Polarizable Embedding Calculations

Calculating absorption and fluorescence spectra for chromophores in solution with ensemble Franck-Condon methods

Accurately modeling absorption and fluorescence spectra for molecules in solution poses a challenge due to the need to incorporate both vibronic and environmental effects, as well as the necessity of accurate excited state electronic structure calculations. Nuclear ensemble approaches capture explicit environmental effects, Franck-Condon methods capture vibronic effects, and recently introduced ensemble-Franck-Condon approaches combine the advantages of both methods. In this study, we present and analyze simulated absorption and fluorescence spectra generated with combined ensemble-Franck-Condon approaches for three chromophore-solvent systems and compare them to standard ensemble and Franck-Condon spectra, as well as to the experiment. Employing configurations obtained from ground and excited state ab initio molecular dynamics, three combined ensemble-Franck-Condon approaches are directly compared to each other to assess the accuracy and relative computational time. We find that the approach employing an average finite-temperature Franck-Condon line shape generates spectra nearly identical to the direct summation of an ensemble of Franck-Condon spectra at one-fourth of the computational cost. We analyze how the spectral simulation method, as well as the level of electronic structure theory, affects spectral line shapes and associated Stokes shifts for 7-nitrobenz-2-oxa-1,3-diazol-4-yl and Nile red in dimethyl sulfoxide and 7-methoxy coumarin-4-acetic acid in methanol. For the first time, our studies show the capability of combined ensemble-Franck-Condon methods for both absorption and fluorescence spectroscopy and provide a powerful tool for simulating linear optical spectra.

Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation

The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.

Polarizable Embedding without Artificial Boundary Polarization

We present a fully self-consistent polarizable embedding (PE) model that does not suffer from unphysical boundary polarization. This is achieved through the use of the minimum-image convention (MIC) in the induced electrostatics. It is a simple yet effective approach that includes a more physically accurate description of the polarization throughout the molecular system. Using PE with MIC (PE-MIC), we shed new light on the limitations of commonly employed cutoff models, such as the droplet model, when used in PE calculations. Specifically, we investigate the effects of the unphysical polarization at the outer boundary by comparing induced dipoles and the associated electrostatic potentials, as well as some optical properties of solute-solvent and biomolecular systems. We show that the magnitude of the inaccuracies caused by the unphysical polarization depends on multiple parameters: the nature of the quantum subsystem and of the environment, the cutoff model and distance, and the calculated property.

Trajectory Surface Hopping for a Polarizable Embedding QM/MM Formulation

We present the implementation of trajectory surface-hopping nonadiabatic dynamics for a polarizable embedding QM/MM formulation. Time-dependent density functional theory was used at the quantum mechanical level of theory, whereas the molecular mechanics description involved the polarizable AMOEBA force field. This implementation has been obtained by integrating the surface-hopping program Newton-X NS with an interface between the Gaussian 16 and the Tinker suites of codes to calculate QM/AMOEBA energies and forces. The implementation has been tested on a photoinduced electron-driven proton-transfer reaction involving pyrimidine and a hydrogen-bonded water surrounded by a small cluster of water molecules and within a large water droplet.

Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach

We report on the first formulation of a novel polarizable QM/MM approach, where the density functional tight binding (DFTB) is coupled to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ) is then extended to the linear response within the TD-DFTB framework and challenged to study absorption spectra of large condensed-phase systems.

Harmonic Infrared and Raman Spectra in Molecular Environments Using the Polarizable Embedding Model

We present a fully analytic approach to calculate infrared (IR) and Raman spectra of molecules embedded in complex molecular environments modeled using the fragment-based polarizable embedding (PE) model. We provide the theory for the calculation of analytic second-order geometric derivatives of molecular energies and first-order geometric derivatives of electric dipole moments and dipole-dipole polarizabilities within the PE model. The derivatives are implemented using a general open-ended response theory framework, thus allowing for an extension to higher-order derivatives. The embedding-potential parameters used to describe the environment in the PE model are derived through first-principles calculations, thus allowing a wide variety of systems to be modeled, including solvents, proteins, and other large and complex molecular environments. Here, we present proof-of-principle calculations of IR and Raman spectra of acetone in different solvents. This work is an important step toward calculating accurate vibrational spectra of molecules embedded in realistic environments.

A polarizable three-layer frozen density embedding/molecular mechanics approach

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding

Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.

A zeroth-order active-space frozen-orbital embedding scheme for multireference calculations

Multireference computations of large-scale chemical systems are typically limited by the computational cost of quantum chemistry methods. In this work, we develop a zeroth-order active space embedding theory [ASET(0)], a simple and automatic approach for embedding any multireference dynamical correlation method based on a frozen-orbital treatment of the environment. ASET(0) is combined with the second-order multireference driven similarity renormalization group and tested on several benchmark problems, including the excitation energy of 1-octene and bond-breaking in ethane and pentyldiazene. Finally, we apply ASET(0) to study the singlet-triplet gap of p-benzyne and 9,10-anthracyne diradicals adsorbed on a NaCl surface. Our results show that despite its simplicity, ASET(0) is a powerful and sufficiently accurate embedding scheme applicable when the coupling between the fragment and the environment is in the weak to medium regime.

Coupled Cluster Theory with Induced Dipole Polarizable Embedding for Ground and Excited States

In this work, we present the theory and implementation of the coupled cluster single and double excitations (CCSD) method combined with a classical polarizable molecular mechanics force field (MMPol) based on the induced dipole model. The method is developed to compute electronic excitation energies within the state specific (SS) and linear response (LR) formalisms for the interaction of the quantum mechanical and classical regions. Furthermore, we consider an approximate expression of the correlation energy, originally developed for CCSD with implicit solvation models, where the interaction term is linear in the coupled cluster density. This approximation allows us to include the explicit contribution of the environment to the CC equations without increasing the computational effort. The test calculations on microsolvated systems, where the CCSD/MMPol method is compared to full CCSD calculations, demonstrates the reliability of this computational protocol for all interaction schemes (errors < 2%). We also show that it is important to include induced dipoles on all atom centers of the classical region and that too diffuse functions in the basis set may be problematic due to too strong interaction with the environment.

Discrete Solvent Reaction Field Calculations for One- and Two-Photon Absorptions of Solution-Phase Dimethylaminonitrostilbene Molecule

Based on the configurations generated by molecular dynamics (MD) simulations using the on-the-fly density-functional tight-bonding (DFTB) force field, we investigated performance of the discrete solvent reaction field (DRF) model coupled to time-dependent density functional theory (TD-DFT) for solvatochromic effect of one- and two-photon absorption phenomena. Dimethylaminonitrostilbene (DANS) molecule solvated in chloroform, dichloromethane, and dimethyl sulfoxide solvents was selected as a model system for our research purpose. For every selected MD/DFTB configuration, within the context of the DRF, solute molecule is represented by TD-DFT and solvent molecules are described by atomic charges and polarizabilities. The calculated one-photon absorption energies reproduce well the positive solvatochromic behavior of solvated DANS and are in good agreement with available experimental data. For the two-photon absorption cross section, even though our approach overshot the experimental data by about 20% in absolute magnitude, experimentally observed solvatochromic change was captured qualitatively in this work. At last, we examined the contributions of atomic charges and polarizabilities of solvent molecules to the solvatochromic shifts of properties of interest.

Mutually polarizable QM/MM model with in situ optimized localized basis functions

We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with onetep linear-scaling density functional theory and the classical subsystem - with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully mutually polarizable. A total energy minimization scheme is employed for the Hamiltonian of the coupled QM/MM system. We demonstrate that, compared to simpler models using fixed basis sets, the additional flexibility offered by in situ optimized basis functions improves the accuracy of the QM/MM interface, but also poses new challenges, making the QM subsystem more prone to overpolarization and unphysical charge transfer due to increased charge penetration. We show how these issues can be efficiently solved by replacing the classical repulsive van der Waals term for QM/MM interactions with an interaction of the electronic density with a fixed, repulsive MM potential that mimics Pauli repulsion, together with a modest increase in the damping of QM/MM polarization. We validate our method, with particular attention paid to the hydrogen bond, in tests on water-ion pairs, the water dimer, first solvation shells of neutral and charged species, and solute-solvent interaction energies. As a proof of principle, we determine suitable repulsive potential parameters for water, K⁺, and Cl^-. The mechanisms we employed to counteract the unphysical overpolarization of the QM subsystem are demonstrated to be adequate, and our approach is robust. We find that the inclusion of explicit polarization in the MM part of QM/MM improves agreement with fully QM calculations. Our model permits the use of minimal size QM regions and, remarkably, yields good energetics across the well-balanced QM/MM interface.

Solvent effects on excitation energies obtained using the state-specific TD-DFT method with a polarizable continuum model based on constrained equilibrium thermodynamics

Nonequilibrium solvation effects need to be treated properly in the study of electronic absorption processes of solutes since solvent polarization is not in equilibrium with the excited-state charge density of the solute. In this work, we developed a state specific (SS) method based on the novel nonequilibrium solvation model with constrained equilibrium manipulation to account for solvation effects in electronic absorption processes. Time-dependent density functional theory (TD-DFT) is adopted to calculate electronic excitation energies and a polarizable continuum model is employed in the treatment of bulk solvent effects on both the ground and excited electronic states. The equations based on this novel nonequilibrium solvation model in the framework of TDDFT to calculate vertical excitation energy are presented and implemented in the Q-Chem package. The implementation is validated by comparing reorganization energies for charge transfer excitations between two atoms obtained from Q-Chem and those obtained using a two-sphere model. Solvent effects on electronic transitions of coumarin 153 (C153), acetone, pyridine, (2E)-3-(3,4-dimethoxyphenyl)-1-(2-hydroxyphenyl)prop-2-en-1-one (DMHP), and uracil in different solvents are investigated using the newly developed code. Our results show that the obtained vertical excitation energies as well as spectral shifts generally agree better with the available experimental values than those obtained using the traditional nonequlibrium solvation model. This new model is thus appropriate to study nonequilibrium excitation processes in solution.

Dispersion Interactions in QM/EFP

The dispersion energy term between quantum-mechanical (QM) and classical (represented by effective fragment potentials, EFP) subsystems is developed and implemented. A new formulation is based on long-range perturbation theory and uses dynamic polarizability tensors of the effective fragments and electric field integrals and orbital energies of the quantum-mechanical subsystem. No parametrization is involved. The accuracy of the QM-EFP dispersion energy is tested on a number of model systems; the average mean unsigned error is 0.8 kcal/mol or 13% with respect to the symmetry adapted perturbation theory on the S22 data set of noncovalent interactions. The computational cost of the dispersion energy computation is low compared to the self-consistent field calculation of the QM subsystem. The dispersion energy is sensitive to the level of theory employed for the QM part and to the electrostatic interactions in the system. The latter means that the dispersion interactions in the QM/EFP method are not purely two-body but have more complex many-body behavior.

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis

The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH₃, NH₄⁺, CN^-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.

Electronic absorption spectra of pyridine and nicotine in aqueous solution with a combined molecular dynamics and polarizable QM/MM approach

The electronic absorption spectra of pyridine and nicotine in aqueous solution have been computed using a multistep approach. The computational protocol consists in studying the solute solvation with accurate molecular dynamics simulations, characterizing the hydrogen bond interactions, and calculating electronic transitions for a series of configurations extracted from the molecular dynamics trajectories with a polarizable QM/MM scheme based on the fluctuating charge model. Molecular dynamics simulations and electronic transition calculations have been performed on both pyridine and nicotine. Furthermore, the contributions of solute vibrational effect on electronic absorption spectra have been taken into account in the so called vertical gradient approximation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Averaged Condensed Phase Model for Simulating Molecules in Complex Environments

The need for configurational sampling dramatically increases the cost of combined quantum mechanics/molecular mechanics (QM/MM) simulations of chemical processes in solution. We developed an averaged condensed phase environment (ACPE) model that constructs an effective polarizable environment directly from explicitly sampled molecular dynamics configurations via the K-means++ algorithm and a mathematically rigorous translation of the molecular mechanics parameters. The model captures detailed heterogeneous features in the environment that may be difficult to describe using a conventional polarizable continuum model. Instead of performing repeated QM/MM calculations for each new configuration of the environment, the ACPE approach allows one to perform a single QM calculation on an averaged configuration. Here, we demonstrate the model by computing electronic excitation energies for several small molecules in solution. The ACPE model predicts the excitation energies in excellent agreement with conventional configurational averaging yet with orders of magnitude of reduction in the computational cost.

An averaged polarizable potential for multiscale modeling in phospholipid membranes

A set of average atom-centered charges and polarizabilities has been developed for three types of phospholipids for use in polarizable embedding calculations. The lipids investigated are 1,2-dimyristoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, and 1-palmitoyl-2-oleoyl-sn-glycerol-3-phospho-L-serine given their common use both in experimental and computational studies. The charges, and to a lesser extent the polarizabilities, are found to depend strongly on the molecular conformation of the lipids. Furthermore, the importance of explicit polarization is underlined for the description of larger assemblies of lipids, that is, membranes. In conclusion, we find that specially developed polarizable parameters are needed for embedding calculations in membranes, while common non-polarizable point-charge force fields usually perform well enough for structural and dynamical studies. © 2017 Wiley Periodicals, Inc.

Automated Fragmentation Polarizable Embedding Density Functional Theory (PE-DFT) Calculations of Nuclear Magnetic Resonance (NMR) Shielding Constants of Proteins with Application to Chemical Shift Predictions

Full-protein nuclear magnetic resonance (NMR) shielding constants based on ab initio calculations are desirable, because they can assist in elucidating protein structures from NMR experiments. In this work, we present NMR shielding constants computed using a new automated fragmentation (J. Phys. Chem. B 2009, 113, 10380-10388) approach in the framework of polarizable embedding density functional theory. We extend our previous work to give both basis set recommendations and comment on how large the quantum mechanical region should be to successfully compute ¹³C NMR shielding constants that are comparable with experiment. The introduction of a probabilistic linear regression model allows us to substantially reduce the number of snapshots that are needed to make comparisons with experiment. This approach is further improved by augmenting snapshot selection with chemical shift predictions by which we can obtain a representative subset of snapshots that gives the smallest predicted error, compared to experiment. Finally, we use this subset of snapshots to calculate the NMR shielding constants at the PE-KT3/pcSseg-2 level of theory for all atoms in the protein GB3.

Polarizable Embedding Density Matrix Renormalization Group

The polarizable embedding (PE) approach is a flexible embedding model where a preselected region out of a larger system is described quantum mechanically, while the interaction with the surrounding environment is modeled through an effective operator. This effective operator represents the environment by atom-centered multipoles and polarizabilities derived from quantum mechanical calculations on (fragments of) the environment. Thereby, the polarization of the environment is explicitly accounted for. Here, we present the coupling of the PE approach with the density matrix renormalization group (DMRG). This PE-DMRG method is particularly suitable for embedded subsystems that feature a dense manifold of frontier orbitals which requires large active spaces. Recovering such static electron-correlation effects in multiconfigurational electronic structure problems, while accounting for both electrostatics and polarization of a surrounding environment, allows us to describe strongly correlated electronic structures in complex molecular environments. We investigate various embedding potentials for the well-studied first excited state of water with active spaces that correspond to a full configuration-interaction treatment. Moreover, we study the environment effect on the first excited state of a retinylidene Schiff base within a channelrhodopsin protein. For this system, we also investigate the effect of dynamical correlation included through short-range density functional theory.

Efficient determination of accurate atomic polarizabilities for polarizeable embedding calculations

We evaluate embedding potentials, obtained via various methods, used for polarizable embedding computations of excitation energies of para-nitroaniline in water and organic solvents as well as of the green fluorescent protein. We found that isotropic polarizabilities derived from DFTD3 dispersion coefficients correlate well with those obtained via the LoProp method. We show that these polarizabilities in conjunction with appropriately derived point charges are in good agreement with calculations employing static multipole moments up to quadrupoles and anisotropic polarizabilities for both computed systems. The (partial) use of these easily-accessible parameters drastically reduces the computational effort to obtain accurate embedding potentials especially for proteins. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Revealing Nucleic Acid Mutations Using Förster Resonance Energy Transfer-Based Probes

Nucleic acid mutations are of tremendous importance in modern clinical work, biotechnology and in fundamental studies of nucleic acids. Therefore, rapid, cost-effective and reliable detection of mutations is an object of extensive research. Today, Förster resonance energy transfer (FRET) probes are among the most often used tools for the detection of nucleic acids and in particular, for the detection of mutations. However, multiple parameters must be taken into account in order to create efficient FRET probes that are sensitive to nucleic acid mutations. In this review; we focus on the design principles for such probes and available computational methods that allow for their rational design. Applications of advanced, rationally designed FRET probes range from new insights into cellular heterogeneity to gaining new knowledge of nucleic acid structures directly in living cells.

LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields

We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER-HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single-point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo-bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.

Averaged Solvent Embedding Potential Parameters for Multiscale Modeling of Molecular Properties

We derive and validate averaged solvent parameters for embedding potentials to be used in polarizable embedding quantum mechanics/molecular mechanics (QM/MM) molecular property calculations of solutes in organic solvents. The parameters are solvent-specific atom-centered partial charges and isotropic polarizabilities averaged over a large number of geometries of solvent molecules. The use of averaged parameters reduces the computational cost to obtain the embedding potential, which can otherwise be a rate-limiting step in calculations involving large environments. The parameters are evaluated by analyzing the quality of the resulting molecular electrostatic potentials with respect to full QM potentials. We show that a combination of geometry-specific parameters for solvent molecules close to the QM region and averaged parameters for solvent molecules further away allows for efficient polarizable embedding multiscale modeling without compromising the accuracy. The results are promising for the development of general embedding parameters for biomolecules, where the reduction in computational cost can be considerable.

Open-ended response theory with polarizable embedding: multiphoton absorption in biomolecular systems

We present the theory and implementation of an open-ended framework for electric response properties that includes effects from the molecular environment modeled by the polarizable embedding model.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

Analytical gradients for MP2, double hybrid functionals, and TD-DFT with polarizable embedding described by fluctuating charges

A polarizable quantum mechanics (QM)/ molecular mechanics (MM) approach recently developed for Hartree-Fock (HF) and Kohn-Sham (KS) methods has been extended to energies and analytical gradients for MP2, double hybrid functionals, and TD-DFT models, thus allowing the computation of equilibrium structures for excited electronic states together with more accurate results for ground electronic states. After a detailed presentation of the theoretical background and of some implementation details, a number of test cases are analyzed to show that the polarizable embedding model based on fluctuating charges (FQ) is remarkably more accurate than the corresponding electronic embedding based on a fixed charge (FX) description. In particular, a set of electronegativities and hardnesses has been optimized for interactions between QM and FQ regions together with new repulsion-dispersion parameters. After validation of both the numerical implementation and of the new parameters, absorption electronic spectra have been computed for representative model systems including vibronic effects. The results show remarkable agreement with full QM computations and significant improvement with respect to the corresponding FX results. The last part of the article provides some hints about computation of solvatochromic effects on absorption spectra in aqueous solution as a function of the number of FQ water molecules and on the use of FX external shells to improve the convergence of the results. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

Assessing Molecular Dynamics Simulations with Solvatochromism Modeling

For the modeling of solvatochromism with an explicit representation of the solvent molecules, the quality of preceding molecular dynamics simulations is crucial. Therefore, the possibility to apply force fields which are derived with as little empiricism as possible seems desirable. Such an approach is tested here by exploiting the sensitive solvatochromism of p-nitroaniline, and the use of reliable excitation energies based on approximate second-order coupled cluster results within a polarizable embedding scheme. The quality of the various MD settings for four different solvents, water, methanol, ethanol, and dichloromethane, is assessed. In general, good agreement with the experiment is observed when polarizable force fields and special treatment of hydrogen bonding are applied.