Ab initio molecular dynamics with a continuum solvation model

Pyrrolizidine alkaloids from Jacobaea vulgaris Gaertn and theoretical studies on intramolecular interactions

Two undescribed pyrrolizidine alkaloids, 13-dehydrosenkirkine (1) and chloromethylretrorsine (2), along with three known analogues, onetine (3), retrorsine (4), and usaramine N-oxide (5), were isolated from Jacobaea vulgaris Gaertn. The structures of the undescribed compounds were elucidated by extensive spectrometric and spectroscopic techniques, including HRESIMS, NMR, calculated 13C-NMR DP4+ analysis and comparison with experimental and calculated ECD spectra. The undescribed compounds were evaluated for their antitumour activity against HT29, HeLa, and HepG2 cells. In addition, the intramolecular interactions and quadrupolar couplings were revealed by investigating the geometrical and electronic properties of three typical otonecine-type PAs in DFT theory.

Toward theoretical terahertz spectroscopy of glassy aqueous solutions: partially frozen solute-solvent couplings of glycine in water

The molecular-level understanding of THz spectra of aqueous solutions under ambient conditions has been greatly advanced in recent years. Here, we go beyond previous analyses by performing ab initio molecular dynamics simulations of glycine in water with artificially frozen solute or solvent molecules, respectively, while computing the total THz response as well as its decomposition into mode-specific resonances based on the "supermolecular solvation complex" technique. Clamping the water molecules and keeping glycine moving breaks the coupling of glycine to the structural dynamics of the solvent, however, the polarization and dielectric solvation effects in the static solvation cage are still at work since the full electronic structure of the quenched solvent is taken into account. The complementary approach of fixing glycine reveals both the dynamical and electronic response of the solvation cage at the level of its THz response. Moreover, to quantitatively account for the electronic contribution solely due to solvent embedding, the solute species is "vertically desolvated", thus preserving the fully coupled solute-solvent motion in terms of the solute's structural dynamics in solution, while its electronic structure is no longer subject to solute-solvent polarization and charge transfer effects. When referenced to the free simulation of Gly(aq), this three-fold approach allows us to decompose the THz spectral contributions due to the correlated solute-solvent dynamics into entirely structural and purely electronic effects. Beyond providing hitherto unknown insights, the observed systematic changes of THz spectra in terms of peak shifts and lineshape modulations due to conformational freezing and frozen solvation cages might be useful to investigate the solvation of molecules in highly viscous H-bonding solvents such as ionic liquids and even in cryogenic ices as relevant to polar stratospheric and dark interstellar clouds.

Selective Solvent-Induced Stabilization of Polar Oxide Surfaces in an Electrochemical Environment

The impact of an electrochemical environment on the thermodynamic stability of polar oxide surfaces is investigated for the example of ZnO(0001) surfaces immersed in water using density functional theory calculations. We show that solvation effects are highly selective: They have little effect on surfaces showing a metallic character, but largely stabilize semiconducting structures, particularly those that have a high electrostatic penalty in vacuum. The high selectivity is shown to have direct consequences for the surface phase diagram and explains, e.g., why certain surface structures could be observed only in an electrochemical environment.

Continuous and smooth potential energy surface for conductorlike screening solvation model using fixed points with variable areas

Rigorously continuous and smooth potential energy surfaces, as well as exact analytic gradients, are obtained for a conductorlike screening solvation model (CPCM, a variant of the general COSMO) with Hartree-Fock (RHF, ROHF, UHF, and MCSCF) and density functional theory (R-DFT, RO-DFT, and U-DFT) methods using a new tessellation scheme, fixed points with variable areas (FIXPVA). In FIXPVA, spheres centered at atoms are used to define the molecular cavity and surface. The surface of each sphere is divided into 60, 240, or 960 tesserae, which have positions fixed relative to the sphere center and areas scaled by switching functions of their distances to neighboring spheres. Analytic derivatives of the positions and areas of the surface tesserae with respect to atomic coordinates can be obtained and used to evaluate the solvation energy gradients. Due to the accurate analytic gradients and smooth potential energy surface, geometry optimization processes using these methods are stable and convergent.

A unified electrostatic and cavitation model for first-principles molecular dynamics in solution

The electrostatic continuum solvent model developed by [Fattebert and Gygi J. Comput. Chem. 23, 662 (2002); Int. J. Quantum Chem. 93, 139 (2003)] is combined with a first-principles formulation of the cavitation energy based on a natural quantum-mechanical definition for the surface of a solute. Despite its simplicity, the cavitation contribution calculated by this approach is found to be in remarkable agreement with that obtained by more complex algorithms relying on a large set of parameters. Our model allows for very efficient Car-Parrinello simulations of finite or extended systems in solution and demonstrates a level of accuracy as good as that of established quantum-chemistry continuum solvent methods. We apply this approach to the study of tetracyanoethylene dimers in dichloromethane, providing valuable structural and dynamical insights on the dimerization phenomenon.

Hydration of Cyclohexylamines: CPCM Calculation of Hydration Gibbs Energy of the Conformers

This work presents a theoretical study on the hydration of cyclohexylamine and isomers of cyclohexyldiamine. All possible conformers were fully optimized in solution using the conductor-like polarizable continuum model (CPCM) and density functional theory. Values of the Gibbs energy of solvation, its respective contributions (electrostatic, nonelectrostatic and conformational change), and the relative Gibbs energy of the conformers in aqueous solution and gas phase are reported. From these values and the Boltzmann populations of the conformers in both phases, the weighted mean values of DeltaG(solv) for the compounds are calculated. Three structural features were found to be important for the hydration of these compounds: the distance between the two NH2 groups (proximity disfavors hydration), their position relative to the ring (equatorial is preferred over axial), and the orientation of the nitrogen lone-pairs (gauche is more favorable to hydration than trans). In the particular case of vicinal cyclohexyldiamines, in addition to these two factors, the relative orientation of one group to the other should also be taken into account.

Ab initio molecular dynamics study of glycine intramolecular proton transfer in water

We use ab initio molecular-dynamics simulations to quantify structural and thermodynamic properties of a model proton transfer reaction that converts a neutral glycine molecule, stable in the gas phase, to the zwitterion that predominates in aqueous solution. We compute the potential of mean force associated with the direct intramolecular proton transfer event in glycine. Structural analyses show that the average hydration number (N(w)) of glycine is not constant along the reaction coordinate, but rather progresses from N(w) = 5 in the neutral molecule to N(w) = 8 for the zwitterion. We report the free-energy difference between the neutral and charged glycine molecules, and the free-energy barrier to proton transfer. Finally, we identify the approximations inherent in our method and estimate the corresponding corrections to our reported thermodynamic predictions.

[CpRu((R)-Binop-F)(H2O)][SbF6], a new fluxional chiral Lewis acid catalyst: synthesis, dynamic nmr, asymmetric catalysis, and theoretical studies

The C(2)-symmetric electron-poor ligand (R)-BINOP-F (4) was prepared by reaction of (R)-BINOL with bis(pentafluorophenyl)-phosphorus bromide in the presence of triethylamine. The iodo complex [CpRu((R)-BINOP-F)(I)] ((R)-6) was obtained by substitution of two carbonyl ligands by (R)-4 in the in situ-prepared [CpRu(CO)(2)H] complex followed by reaction with iodoform. Complex 6 was reacted with [Ag(SbF(6))] in acetone to yield [CpRu((R)-BINOP-F)(acetone)][SbF(6)] ((R)-7). X-ray structures were obtained for both (R)-6 and (R)-7. The chiral one-point binding Lewis acid [CpRu((R)-BINOP-F)][SbF(6)] derived from either (R)-7 or the corresponding aquo complex (R)-8 activates methacrolein and catalyzes the Diels-Alder reaction with cyclopentadiene to give the [4 + 2] cycloadduct with an exo/endo ratio of 99:1 and an ee of 92% of the exo product. Addition occurs predominantly to the methacrolein C(alpha)-Re face. In solution, water in (R)-8 exchanges readily. Moreover, a second exchange process renders the diastereotopic BINOP-F phosphorus atoms equivalent. These processes were studied by the application of variable-temperature (1)H, (31)P, and (17)O NMR spectroscopy, variable-pressure (31)P and(17)O NMR spectroscopy, and, using a simpler model complex, density functional theory (DFT) calculations. The results point to a dissociative mechanism of the aquo ligand and a pendular motion of the BINOP-F ligand. NMR experiments show an energy barrier of 50.7 kJ mol(-1) (12.2 kcal mol(-1)) for the inversion of the pseudo-chirality at the ruthenium center.

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.

Solvent effects on glycine. I. A supermolecule modeling of tautomerization via intramolecular proton transfer

The relative stabilities of glycine tautomers involved in the intramolecular proton transfer are investigated computationally by considering glycine-water complexes containing up to five water molecules. The supermolecule results are compared with continuum calculations. Specific solute-solvent interactions and solvent induced changes in the solute wave function are considered using the natural bond orbitals (NBO) method. The stabilization of the zwitterion upon solvation is explained by the changes in the wave functions localized on the forming and breaking bonds as well as by the different interaction energies in the zwitterionic and neutral clusters. Only the neutral species exist in mono- and dihydrated clusters and in the gas phase. In the smaller clusters, zwitterions are mainly stabilized by conformational effects, whereas in larger clusters, in particular when glycine is solvated on both sides of its heavy atom backbone, polarization effects dominate the stability of a given tautomer. Generally, the strength of the solute-solvent interactions is governed by the intermolecular charge transfer interactions. As the solvation progresses, the hypothetical gaseous zwitterion is better solvated than the gaseous neutral, making zwitterion to neutral tautomerization progressively less exothermic for clusters containing up to three water molecules, and endothermic for larger clusters. The neutral isomer does not exist for some solvent arrangements with five water molecules. Only solvent arrangements in which water molecules do not interact with the reactive proton are considered. Hence, the experimentally observed double well potential energy surface may be due to such an interaction or to a different reaction mechanism.