Molecular orientations at interfaces by extended polarizable continuum model

This study presents an investigation into orientation of molecular solutes at the interface of liquid water and other media. The calculation of electrostatic free energy of molecular solute is based on an extension of the polarizable continuum model (PCM) to interfacial system. The extended PCM computational scheme is incorporated with the self‐consistent field procedure which is necessary to obtain more accurate electrostatic free energy and charge density distribution. The computation of non‐electrostatic energy for interfacial system is also realized. Applying the numerical procedure to molecular systems, N,N′‐diethyl‐p‐nitroaniline (DEPNA) at air/water interface and p‐nitrophenol (PNP) at cyclohexane/water interface, the average orientational angles are in reasonable agreement with the experimental results. Taking both the electrostatic and the non‐electrostatic energies into account, the analysis on the energy profiles shows that the electrostatic solvation energy is the dominant factor in determining the orientation angle for PNP, whereas for DEPNA, the orientation angle mainly depends on the cavitation energy. This suggests that, in addition to the electrostatic energy, taking the cavitation energy into account may provide a more complete view when we survey the molecular orientation at interface. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Neutral hydrolyses of carbon disulfide: An ab initio study of water catalysis

The water-catalyzed hydrolysis reaction of carbon disulfide (CS(2)) has been investigated at the levels of HF and MP2 with the basis set of 6-311++G(d,p) using the combined supramolecular/continuum models, in which up to six water molecules are involved in the hydrolysis and the effect of water bulk solvent is taken into account according to the polarizable continuum model (PCM). The activation Gibbs free energies in water solution, DeltaG(sol) (not equal) (298 K), for the rate-determining steps of one up to six water hydrolyses are 247.9, 184.2, 152.3, 141.8, 134.4, and 118.9 kJ/mol, respectively. The most favorable hydrolysis path of CS(2) involves a sort of eight-membered ring transition structure formed by six water molecules, among which three water molecules are not involved in the proton transfer, two near to the nonreactive sulfur atom, and one below the parent carbon disulfide. This suggests that the hydrolysis of CS(2) can be mediated with the water molecule(s) and be significantly facilitated by the cooperative effects of the water molecule(s) in the nonreactive region. The catalytic effects of water molecule(s) due to the alleviation of ring strain in the proton transfer process may result from the synergistic effects of rehybridization and charge reorganization from the prereaction complex to the rate-determining transition state structure induced by water molecule(s). PCM solvation models could significantly lower the rate-determining activation Gibbs free energies by 20-38 kJ/mol when two up to six explicit water molecules involved in the neutral hydrolysis of CS(2).     

A study on orientation and absorption spectrum of interfacial molecules by using continuum model

In this work, a numerical procedure based on the continuum model is developed and applied to the solvation energy for ground state and the spectral shift against the position and the orientation of the interfacial molecule. The interface is described as a sharp boundary separating two bulk media. The polarizable continuum model (PCM) allows us to account for both electrostatic and nonelectrostatic solute-solvent interactions when we calculate the solvation energy. In this work we extend PCM to the interfacial system and the information about the position and orientation of the interfacial molecule can be obtained. Based on the developed expression of the electrostatic free energy of a nonequilibrium state, the numerical procedure has been implemented and used to deal with a series of test molecules. The time-dependent density functional theory (TDDFT) associated with PCM is used for the electron structure and the spectroscopy calculations of the test molecules in homogeneous solvents. With the charge distribution of the ground and excited states, the position- and orientation-dependencies of the solvation energy and the spectrum have been investigated for the interfacial systems, taking the electrostatic interaction, the cavitation energy, and the dispersion-repulsion interaction into account. The cavitation energy is paid particular attention, since the interface portion cut off by the occupation of the interfacial molecule contributes an extra part to the stabilization for the interfacial system. The embedding depth, the favorable orientational angle, and the spectral shift for the interfacial molecule have been investigated in detail. From the solvation energy calculations, an explanation has been given on why the interfacial molecule, even if symmetrical in structure, tends to take a tilting manner, rather than perpendicular to the interface.     

Hydrolysis of cyclic phosphates by ribonuclease A: a computational study using a simplified ab initio quantum model

The second step in the enzyme-catalyzed hydrolysis of phosphate esters by ribonuclease A (RNase A) was studied using an ab initio quantum-based model of the active site including constrained parts of three critical residues, His-12, His-119, and Lys-41, and a small substrate. The competition between release of the cyclic phosphate intermediate and subsequent hydrolysis following transphosphorylation was explored to determine the electronic factors that contribute to preferential intermediate product release observed experimentally. The structural and energetic results obtained at both the RHF and MP2 levels reveal several contributing factors consistent with experimental observation. Although the intrinsic electronic effects tend to favor hydrolysis slightly with an overall activation free energy of approximately 70 kJ mol(-1), entropic and environmental effects favor release of the cyclic phosphate intermediate over hydrolysis. Exploration of the second, hydrolysis step also revealed interesting similarity with the transphosphorylation step, including the observation of autocatalysis by the substrate. Moreover, both steps of the overall RNase A reaction reveal multiple pathways involving proton transfers to sites of similar proton affinities. The anionic phosphate in both steps can act as a stable proton binding site as protons are moved around the active site throughout the progress of the reaction. These results suggest autocatalysis may be representative of more general behavior in enzymes containing highly charged substrates, especially phosphates.     

The effect of solvent on the internal rotation of formamide: an ab initio study using a polarizable continuum model

The effect of solvent on the barrier to internal rotation of formamide has been studied using an ab initio 4-3IG method coupled with a polarizable continuum model. The trend of increasing rotation barrier with increasing solvent dielectric constant is reproduced. Conformational free energies were calculated by adding dispersion and cavitation terms to the electrostatic energy.     

Solids modeled by ab initio crystal field methods, part 22: the Fock matrix transformed supermolecule model and the structure determination of D-erythronic acid-3,4-carbonate

To eliminate the cutoff effects in the supermolecule model (SM), the Fock matrix transformed supermolecule model (FTSM) is developed. In this model a cyclization of the cluster is performed by means of transformations of elements of the Fock matrix, thereby restoring the translational symmetry. Besides this fundamental enhancement, significant CPU time savings are realized because, in this new procedure, not all Fock matrix elements need to be calculated. The method is applied in a study of the structure of D-erythronic acid-3,4-carbonate in the crystal phase, where the new model is compared to the supermolecule model as well as the XRD experiment. The results are found to be in good agreement with experimental data.     

On the performance of time-dependent density functional theory and polarizable continuum model in the study of aqueous formaldehyde

In this paper Time-Dependent Density Functional Theory (TD-DFT) with Becke's three-parameter hybrid method using the Lee–Yang–Parr correlation functional (B3LYP) in conjunction with different basis sets is applied to calculate the solvatochromic shift on the lowest singlet transition of formaldehyde in vacuum and in solution. Our results, in excellent agreement with solvent excitation energy shift of acetone in water at 298–300 K (no experimental data on formaldehyde in water are available so far) and with previous high level quantum chemical calculations, help in appreciating the potentialities of the combined DFT/PCM approach on aquo-complexes for the study of structural and spectroscopic properties of molecule in solution.     

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.     

Finite‐field method with unbiased polarizable continuum model for evaluation of the second hyperpolarizability of an open‐shell singlet molecule in solvents

The static second hyperpolarizability γ of the complexes composed of open‐shell singlet 1,3‐dipole molecule involving a boron atom and a water molecule in aqueous phase are investigated by the finite‐field (FF) method combined with a standard polarized continuum model (PCM) and with a newly proposed unbiased PCM (UBPCM). On the basis of the comparison with the results calculated by the FF method using the full quantum and the quantum‐mechanical/molecular‐mechanical and molecular‐dynamics (QM/MM‐MD) treatments, the present FF‐UBPCM method is demonstrated to remedy the artificial overestimation of the γ caused by standard FF‐PCM calculations and to well reproduce the FF‐QM/MM‐MD and FF‐full‐QM results with much lower costs. © 2013 Wiley Periodicals, Inc.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

Theoretical study on the hydrolysis mechanism of N,N-dimethyl-N'-(2-oxo-1, 2-dihydro-pyrimidinyl)formamidine: water-assisted mechanism and cluster-continuum model

The hydrolysis reaction of N,N-dimethyl-N'-(2-oxo-1, 2-dihydro-pyrimidinyl)formamidine (DMPFA), a model compound of the antivirus drug amidine-3TC (3TC = 2', 3'-dideoxy-3'-thiacytidine), is investigated by the hybrid density functional theory B3LYP/6-31+G (d,p) method. The hydrolysis reaction of the title compound is predicted to undergo via two pathways, each of which is a stepwise process. Path A is the addition of H2O to the C=N double bond in the amidine group to form a tetrahedral structure in its first step, and then the transfer of the H atom of hydroxyl leads to the corresponding products via four possible channels. Path B simultaneously involves the nucleophilic attack of H2O to the C atom of the C=N bond and the proton transfer to the N atom of amino group leading to the cleavage of the C-N single bond in the amidine group. The results indicate that path A is more favorable than path B in the gas phase. Moreover, to simulate the title reaction in aqueous solution, water-assisted mechanism and the cluster-continuum model, based on the SCRF/CPCM model, are taken into account in our work. The results indicate that it is rational for two water molecules served as a bridge to assist in the first step of path A and that cytosine rather than the cytosine-substituted formamide should be released from the tetrahedral intermediate via s six-membered cycle transition state (channel 2). Our calculations exhibit that the process toward the tetrahedral intermediate is the rate-determining step both in the gas phase and in aqueous solution.     

Hydrogen bonding and dimeric self-association of 2-pyrrolidinone: An ab initio study

Ab initio calculations on the different associated structures of 2-pyrrolidinone with water and with itself were carried out using 3-21G and 6-31G^* basis sets at the Hartree–Fock level, including electron correlation using second-order Møller–Plesset perturbation theory. The calculated free energy changes for the intermolecular hydrogen bonded dimer and hydrated species indicated that the molecular systems with cyclic dimerization and association with two water molecules are dominant. The results are compared to the available experimental data in the literature.     

Building cavities in a fluid of spherical or rod-like particles: a contribution to the solvation free energy in isotropic and anisotropic polarizable continuum model

A general formalism for the calculation of cavitation energies in the framework of the scaled particle theory has been implemented in the Polarizable Continuum Model (PCM), contributing to the nonelectrostatic part of the molecular free energy in solution. The solute cavity and the solvent molecules are described as hard spherocylinders, whose radius and length are related to the actual molecular shape, while the solvent density is estimated from experimental data, or from the solvent molecular volume, suitably scaled. The present model can describe isotropic solutions of spherical and rod-like molecules in spherical or rod-like solvents, and also anisotropic solutions in which the solvent molecules are oriented in space: in this case, the cavitation energy also depends on the relative orientation of solute and solvent molecules. Test calculations have been performed on simple systems to evaluate the accuracy of the present approach, in comparison with other methods and with the available experimental estimates of the cavitation energy, giving encouraging results.     

A new extension of the polarizable continuum model: Toward a quantum chemical description of chemical reactions at extreme high pressure

A quantum chemical method for studying potential energy surfaces of reactive molecular systems at extreme high pressures is presented. The method is an extension of the standard Polarizable Continuum Model that is usually used for Quantum Chemical study of chemical reactions at a standard condition of pressure. The physical basis of the method and the corresponding computational protocol are described in necessary detail, and an application of the method to the dimerization of cyclopentadiene (up to 20 GPa) is reported. © 2015 Wiley Periodicals, Inc.     

Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: new energy gradients and molecular surface tessellation

New equations are derived and implemented for efficient and accurate computation of solvation energy derivatives for the conductor-like polarizable continuum model (C-PCM) and the isotropic integral equation formalism polarizable continuum model (IEF-PCM). Two new molecular surface tessellation procedures GEPOL-RT and GEPOL-AS that generate near continuous potential energy surfaces are proposed for PCM geometry optimization. The combined use of these new techniques leads to efficient and convergent geometry optimizations with the PCMs.     

Exchange coupling and magnetic anisotropy in a family of bipyrimidyl radical‐bridged dilanthanide complexes: Density functional theory and ab initio calculations

The origin of the magnetic anisotropy energy barriers in a series of bpym^? (bpym = 2,2′‐bipyrimidine) radical‐bridged dilanthanide complexes [(Cp*₂Ln)₂(μ‐bpym)]⁺ [Cp* = pentamethylcyclopentadienyl; Ln = Gd^III ( 1 ), Tb^III ( 2 ), Dy^III ( 3 ), Ho^III ( 4 ), Er^III ( 5 )] has been explored using density functional theory (DFT) and ab initio methods. DFT calculations show that the exchange coupling between the two lanthanide ions for each complex is very weak, but the antiferromagnetic Ln‐bpym^? couplings are strong. Ab initio calculations show that the effective energy barrier of 2 or 3 mainly comes from the contribution of a single Tb^III or Dy^III fragment, which is only about one third of a single Ln energy barrier. For 4 or 5 , however, both of the two Ho^III or Er^III fragments contribute to the total energy barrier. Thus, it is insufficient to only increase the magnetic anisotropy energy barrier of a single Ln ion, while enhancing the Ln‐bpym^? couplings is also very important. © 2014 Wiley Periodicals, Inc.     

Dispersion and repulsion contributions to the solvation free energy: comparison of quantum mechanical and classical approaches in the polarizable continuum model

We report a systematic comparison of the dispersion and repulsion contributions to the free energy of solvation determined using quantum mechanical self-consistent reaction field (QM-SCRF) and classical methods. In particular, QM-SCRF computations have been performed using the dispersion and repulsion expressions developed in the framework of the integral equation formalism of the polarizable continuum model, whereas classical methods involve both empirical pairwise potential and surface-dependent approaches. Calculations have been performed for a series of aliphatic and aromatic compounds containing prototypical functional groups in four solvents: water, octanol, chloroform, and carbon tetrachloride. The analysis is focused on the dependence of the dispersion and repulsion components on the level of theory used in QM-SCRF computations, the contribution of those terms in different solvents, and the magnitude of the coupling between electrostatic and dispersion-repulsion components. Finally, comparison is made between the dispersion-repulsion contributions obtained from QM-SCRF calculations and the results determined from classical approaches.     

A comparative molecular field analysis (CoMFA) study using semiempirical, density functional, ab initio methods and pharmacophore derivation using DISCOtech on sigma 1 ligands

The Comparative Molecular Field Analysis (CoMFA) was developed to investigate a three-dimensional quantitative structure activity relationship (3D-QSAR) model of ligands for the sigma 1 receptor. The starting geometry of sigma-1 receptor ligands was obtained from the Tripos force field minimizations and conformations were decided from DISCOtech using the SYBYL 6.8. program. The structures of 48 molecules were fully optimized at the ab initio HF/3-21G* and semiempirical AM1 calculations using GAUSSIAN 98. The electrostatic charges were calculated using several methods such as semiempirical AM1, density functional B3LYP/3-21G*, and ab initio HF/3-21G*, MP2/3-21G* calculations within GAUSSIAN 98. Using the optimized geometries, the CoMFA results derived from the HF/3-21G method were better than those from AM1. The best CoMFA was obtained from HF/3-21G* optimized geometry and charges (R2 = 0.977). Using the optimized geometries, the CoMFA results derived from the HF/3-21G methods were better than those from AM1 calculations. The training set of 43 molecules gave higher R2 (0.989-0.977) from HF/3-21G* optimized geometries than R2 (0.966-0.911) values from AM1 optimized geometries. The test set of five molecules also suggested that HF/3-21G* optimized geometries produced good CoMFA models to predict bioactivity of sigma 1 receptor ligands but AM1 optimized geometries failed to predict reasonable bioactivity of sigma 1 receptor ligands using different calculations for atomic charges.     

Second-order Mø øller–Plesset second derivatives for the polarizable continuum model: theoretical bases and application to solvent effects in electrophilic bromination of ethylene

An analytical theory for the energy second derivatives at second-order Mø øller–Plesset level for solvated molecules described within the polarizable continuum model is presented. The method, which is based on the differentiation of relaxed density equations for the first derivatives, is firstly presented in its formal aspects and is then applied to the study of the rate-determining step of the electrophilic bromination of ethylene in aqueous solution.Contribution to the Jacopo Tomasi Honorary Issue     

Ab initio and DFT modeling of stereoselective deamination of aziridines by nitrosyl chloride

The stereochemical course of the deamination of cis‐2,3‐dimethylaziridine by nitrosyl chloride was investigated at the QCISD/6‐31G(d) level. Calculations reveal that the reaction takes place in two steps. In the first step, the reactants form a pre‐reactive complex, followed by a spiro‐type bicyclic transition state, which on dissociative cycloelimination gives the N‐nitrosoaziridine intermediate. In the second step, this intermediate undergoes cycloreversion through a slightly asynchronous concerted transition state to form an alkene with the same stereochemistry, which is in total agreement with experiment. In the whole reaction, the denitrosation step is found to be rate‐determining. For comparison, geometry optimizations and energies were also obtained at the B3LYP/6‐31G(d) level. It was found that the B3LYP energy results differed significantly from the QCISD ones. To analyze the reason for this difference, B3LYP calculations were repeated by varying the contribution of exact exchange in the Becke functional. With respect to the QCISD results, it has been shown that the functional with 0% exact exchange yields the best activation barriers, whereas the functional with 30% exact exchange is the most suitable one to carry out the complexation and reaction energy calculations. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005