Search

In this work, we consider a new combination of vibrational analysis and normal-like mode decomposition of Debye-Waller factors of solvated ions entirely based on molecular dynamics data. Such a novel time-dependent analysis procedure provides a direct link between x-ray absorption fine structure parameters and normal mode contributions for an ion-solvent system. The potentialities of such a methodology rely on two fundamental aspects which distinguish it from already available tools. First, a general vibrational analysis that does not require any Gaussian or harmonic model for describing atomic fluctuations in liquids. Second, a very accurate sampling of the short range motions around the structural probe via the recently developed atom centered density matrix propagation/general liquid optimized boundary method. This novel molecular dynamics methodology is based on an integrated ab initio/classical potential using localized basis functions and nonperiodic boundary conditions. As a case study we have chosen the Zn(II) ion in aqueous solution. The consistency of our results and the observed good agreement with experiments show how the key support to advanced structural techniques from molecular dynamics can be further expanded and investigated. [ABSTRACT FROM AUTHOR]

In this work, we present a hybrid explicit/implicit solvation model, well suited for first-principles molecular dynamics simulations of solute-solvent systems. An effective procedure is presented that allows to reliably model a solute with a few explicit solvation shells, ensuring solvent bulk behavior at the boundary with the continuum. Such an approach is integrated with high-level ab initio methods using localized basis functions to perform first-principles or mixed quantum mechanics/molecular mechanics simulations within the extended-Lagrangian formalism. A careful validation of the model along with illustrative applications to solutions of acetone and glycine radical are presented, considering two solvents of different polarity, namely, water and chloroform. Results show that the present model describes dynamical and solvent effects with an accuracy at least comparable to that of conventional approaches based on periodic boundary conditions. [ABSTRACT FROM AUTHOR]

A novel molecular dynamics methodology recently proposed by our group [Rega et al., Chem. Phys. Lett. 422, 367 (2006)], which is based on an integrated hybrid potential rooted in high level quantum mechanical methods using localized basis functions and nonperiodic boundary conditions, has been applied to study acrolein in aqueous solution. The solute structural rearrangement and its hydrogen-bonding pattern due to the interactions with water have been analyzed in some detail. Moreover, the solvent effects on the UV n→π* vertical transition and on the NMR 13C and 17O shielding constants of acrolein have been investigated theoretically by performing a posteriori quantum mechanical calculations on a statistically significant number of snapshots extracted from both gas-phase and aqueous solution simulations. Results show that such effective computational strategy can be successfully used to improve our understanding, at atomic level, of important spectroscopic observables. [ABSTRACT FROM AUTHOR]

In this work, the mean field (MF) method, a continuum-based model designed for treating complex molecular systems, such as liquids and solutions, recently presented by Brancato et al. [J. Chem. Phys. 122, 154109 (2005)], has been further developed and improved especially in the treatment of the electrostatics. The revised model has been used to investigate the size effects on several physical properties of various solute-solvent systems by increasing the number of explicitly included solvent molecules from few tens up to thousands. Results on simple ions, such as sodium and chloride ions, and on a small peptide, such as alanine dipeptide analog (AcAlaNHMe), have shown that solvation structures and dynamics, as well as solvent-induced changes in the solute conformation, can be correctly reproduced by the MF model, providing that only two or three solvent layers are treated explicitly. [ABSTRACT FROM AUTHOR]

We present a computational method, exploiting some features of the polarizable continuum model (PCM) to describe heterogeneous media; it belongs to the family of electrostatic embedding mixed methods, such as the more common quantum-mechanical (QM)/molecular mechanics approaches, with the electrostatic long range effects accounted for by a polarized continuum instead of atomic point charges. Provided effective dielectric constants are determined for the various parts of the system, the method is much faster than its atomistic counterpart, and allows for high-level QM calculations on the fragment of interest, using all the highly efficient computational tools developed for homogeneous PCM. Two case studies (the calculation of the pKa of solvent exposed acidic residues in a model protein, and the calculation of the electron spin resonance spectrum of a typical spin probe partially embedded in a membrane) are analyzed in some detail, to illustrate the application of the method to complex systems. [ABSTRACT FROM AUTHOR]

The polarizable continuum model (PCM), used for the calculation of molecular energies, structures, and properties in liquid solution has been deeply revised, in order to extend its range of applications and to improve its accuracy. The main changes effect the definition of solute cavities, of solvation charges and of the PCM operator added to the molecular Hamiltonian, as well as the calculation of energy gradients, to be used in geometry optimizations. The procedure can be equally applied to quantum mechanical and to classical calculations; as shown also with a number of numerical tests, this PCM formulation is very efficient and reliable. It can also be applied to very large solutes, since all the bottlenecks have been eliminated to obtain a procedure whose time and memory requirements scale linearly with solute size. The present procedure can be used to compute solvent effects at a number of different levels of theory on almost all the chemical systems which can be studied in vacuo. [ABSTRACT FROM AUTHOR]

The hyperfine parameters of a number of representative free radicals have been computed by post-Hartree-Fock and density functional approaches including averaging effects from large amplitude vibrations and solvent effects through a recent implementation of the polarizable continuum model. Our results show that fully ab initio hyperfine splittings are accurate enough to back the interpretation of experimental data and to allow an unbiased judgement of the role played by electronic, vibrational, and environmental effects in determining the observed value. The very good results obtained by a density functional approach including some Hartree-Fock exchange both for intrinsic values and for solvent shifts pave the route for the investigation of large biologically significant radicals in their natural aqueous medium. [ABSTRACT FROM AUTHOR]

An approximate method, recently proposed to include in continuum solvation models the effects of electronic charge lying outside the solute cavity, has been adapted and implemented in the framework of the polarizable continuum model (PCM). This formulation exploits all the features already developed for the other PCM versions; it provides molecular free energies, gradients and second derivatives with respect to nuclear coordinates. The performances of this method have been tested in comparison with other PCM versions, in particular, we examined the reliability of this technique to reproduce actual volume charge distribution effects, compared to traditional procedures based on Gauss' Law. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]