A Monte Carlo–quantum mechanics study of the spectroscopic properties of molecules in solution

Solvent effects for vertical absorption and emission processes in solution using a self-consistent state specific method based on constrained equilibrium thermodynamics

A self-consistent state specific (SS) method in the framework of TDDFT is presented to account for solvent effects on absorption and emission processes for molecules in solution. In these processes, the initial state is an equilibrium state, while the polarization of the solvent is in nonequilibrium with the electron density of the solute in the final state. Nonequilibrium solvation free energy is calculated based on a novel nonequilibrium solvation model with constrained equilibrium manipulation. The bulk solvent effects are considered using the polarizable continuum method (PCM), where the solvent-solute interaction is described with a reaction field. Molecular orbitals and orbital energies in the presence of the reaction field corresponding to the excited state are employed and the response of the solvent is not included in the TDDFT calculations. A self-consistent procedure is designed to obtain the excited state reaction field. The equations based on this new nonequilibrium solvation model in the framework of the self-consistent SS-PCM/TDDFT method for calculation of vertical absorption and emission energies are presented and implemented in the Q-Chem package. Vertical absorption and emission energies for several small molecules in solution using the newly developed code are calculated and compared with available experimental data and the results of other theoretical studies. Solvent shifts of absorption and emission energies are reasonably reproduced with this approach. The new model is a promising approach to study nonequilibrium absorption and emission processes in solution.

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.

A DFT study of infrared spectra and Monte Carlo predictions of the solvation shell of Praziquantel and β-cyclodextrin inclusion complex in liquid water

In this paper, we report a theoretical study of the inclusion complexes of Praziquantel (PZQ) and β-cyclodextrin (β-CD) in liquid water. The starting geometry has been carried out by molecular mechanics simulations, and afterwards optimized in B3LYP level with a 6-311G(d) basis set. Monte Carlo simulations have been used to calculate the solvation shell of the PZQ/β-CD inclusion complexes. Moreover, the vibrational frequencies and the infrared intensities for the PZQ/β-CD complex were computed using the B3LYP method. It is demonstrated that this combined model can yield well-converged thermodynamic data even for a modest number of sample configurations, which makes the methodology particularly adequate for understanding the solute-solvent interaction used for generating the liquid structures of one solute surrounded by solvent molecules. The complex solvation shell showed an increase of the water molecule level in relation to the isolated PZQ molecule because of the hydrophilic effect of the CD molecule. The infrared spectra showed that the contribution that originated in the PZQ molecule was not predominant in the upper-wave number region in the drug/β-CD. The movement that purely originated in the PZQ molecule was localized in the absorption band, ranging from 1328 to 1688cm(-1).

A Multiscale Treatment of Angeli's Salt Decomposition

Sodium trioxodinitrate's (Na2N2O3, Angeli's salt) unique cardiovascular effects have been associated with its ability to yield HNO upon dissociation under physiological conditions. Due to its potential applications in new therapies for heart failure, the dissociation of Angeli's salt has recently received increased attention. The decomposition mechanism has been previously studied by quantum mechanical methods using a continuum approximation (PCM) for the solvent effects. In this work we use our recently developed interface of the Amber and Gaussian packages via the PUPIL package to study Angeli's salt dissociation in a hybrid QM/MM scheme where the water solvent molecules are treated explicitly with classical mechanics while the solute is treated with full quantum mechanics (UB3LYP/6-31+G(d) and UMP2/6-31+G(d)) level. Multiple steered molecular dynamics was used with the Jarzynski relationship to extract the free energy profile for the process. We obtain 4.8 kcal mol(-1) and 6.4 kcal mol(-1) free energy barriers for the N-N bond breaking for UB3LYP and UMP2, respectively. The geometries and Mulliken charges for reactant, transition state, and products have been characterized through a number of hybrid QM/MM molecular dynamics runs with the N-N distance restrained to representative values of each species. The results highlight the role of individual solvent molecules for the reaction energetics and provide a comparison point against implicit solvation methods.

Polarizable empirical force field for nitrogen-containing heteroaromatic compounds based on the classical Drude oscillator

The polarizable empirical CHARMM force field based on the classical Drude oscillator has been extended to the nitrogen-containing heteroaromatic compounds pyridine, pyrimidine, pyrrole, imidazole, indole, and purine. Initial parameters for the six-membered rings were based on benzene with nonbond parameter optimization focused on the nitrogen atoms and adjacent carbons and attached hydrogens. In the case of five-member rings, parameters were first developed for imidazole and transferred to pyrrole. Optimization of all parameters was performed against an extensive set of quantum mechanical and experimental data. Ab initio data were used for the determination of initial electrostatic parameters, the vibrational analysis, and in the optimization of the relative magnitudes of the Lennard-Jones (LJ) parameters, through computations of the interactions of dimers of model compounds, model compound-water interactions, and interactions of rare gases with model compounds. The absolute values of the LJ parameters were determined targeting experimental heats of vaporization, molecular volumes, heats of sublimation, crystal lattice parameters, and free energies of hydration. Final scaling of the polarizabilities from the gas-phase values by 0.85 was determined by reproduction of the dielectric constants of pyridine and pyrrole. The developed parameter set was extensively validated against additional experimental data such as diffusion constants, heat capacities, and isothermal compressibilities, including data as a function of temperature.

Statistical mechanically averaged molecular properties of liquid water calculated using the combined coupled cluster/molecular dynamics method

Liquid water is investigated theoretically using combined molecular dynamics (MD) simulations and accurate electronic structure methods. The statistical mechanically averaged molecular properties of liquid water are calculated using the combined coupled cluster/molecular mechanics (CC/MM) method for a large number of configurations generated from MD simulations. The method includes electron correlation effects at the coupled cluster singles and doubles level and the use of a large correlation consistent basis set. A polarizable force field has been used for the molecular dynamics part in both the CC/MM method and in the MD simulation. We describe how the methodology can be optimized with respect to computational costs while maintaining the quality of the results. Using the optimized method we study the energetic properties including the heat of vaporization and electronic excitation energies as well as electric dipole and quadrupole moments, the frequency dependent electric (dipole) polarizability, and electric-field-induced second harmonic generation first and second hyperpolarizabilities. Comparisons with experiments are performed where reliable data are available. Furthermore, we discuss the important issue on how to compare the calculated microscopic nonlocal properties to the experimental macroscopic measurements.

Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde

A combined quantum mechanics/molecular mechanics (QM/MM) method is described, where the polarization between the solvent and solute is accounted for using a self-consistent scheme linear in the solvent polarization. The QM/MM method is implemented for calculation of energies and molecular response properties including the calculation of linear and quadratic response functions using the density-functional theory (DFT) and the Hartree-Fock (HF) theory. Sample calculations presented for ground-state energies, first-order ground-state properties, excitation energies, first-order excited state properties, polarizabilities, first-hyperpolarizabilities, and two-photon absorptions strengths of formaldehyde suggests that DFT may in some cases be a sufficiently reliable alternative to high-level theory, such as coupled-cluster (CC) theory, in modeling solvent shifts, whereas results obtained with the HF wave function deviate significantly from the CC results. Calculations carried out on water gives results that also are comparable with CC calculations in accuracy for ground-state and first-order properties. However, to obtain such accuracy an exchange-correlation functional capable of describing the diffuse Rydberg states must be chosen.

Solvent effects on the n→π[sup ∗] electronic transition in formaldehyde: A combined coupled cluster/molecular dynamics study

We present a study of the blueshift of the n-->pi* electronic transition in formaldehyde in aqueous solution using a combined coupled cluster/molecular mechanics model including mutual polarization effects in the Hamiltonian. In addition, we report ground and excited state dipole moments. Configurations are generated from molecular dynamics simulations with two different force fields, one with and one without an explicit polarization contribution. A statistical analysis using 1200 configurations is presented. Effects of explicit polarization contributions are found to be significant. It is found that the main difference in the effects on the excitation energies arises from the fact that the two force fields result in different liquid structures, and thus a different set of configurations is generated for the coupled cluster/molecular mechanics calculations.

QM/MM Car-Parrinello Molecular Dynamics Study of the Solvent Effects on the Ground State and on the First Excited Singlet State of Acetone in Water

We present a hybrid Car-Parrinello quantum mechanical/molecular mechanical (QM/MM) approach that is capable of treating the dynamics of molecular systems in electronically excited states in complex environments. The potential energy surface in the excited state is described either within the restricted open-shell Kohn-Sham (ROKS) formalism or within time-dependent density functional theory (TDDFT). As a test case, we apply this technique to the study of the solvent effects on the ground state and on the first excited singlet state of acetone in water. Our results demonstrate that for this system a purely classical description of the solvent is sufficient, since inclusion of the first solvent shell of 12 water molecules into the quantum system does not show a significant effect on this transition. The excited-state energies calculated with ROKS are red shifted by a constant value compared to the TDDFT results, while the relative variations of the excitation energy for different configurations are in very good agreement. The experimentally observed blue shift of the excitation energy in going from gas phase to condensed phase is well reproduced. Excited-state dynamics carried out with ROKS yield the relaxation of the solute and the rearrangement of the solvent structure on a picosecond timescale. The calculated Stokes shift is in reasonable agreement with experimental data.