Search

Structural properties of liquid water at ambient temperature were studied using Car-Parrinello [Phys. Rev. Lett. 55, 2471 (1985)] ab initio molecular dynamics (CPAIMD) simulations combined with the Kohn-Sham (KS) density functional theory and the BLYP exchange-correlation functional for the electronic structure. Unlike other recent work on the same subject, where plane-wave (PW) or hybrid Gaussian/plane-wave basis sets were employed, in the present paper, a discrete variable representation (DVR) basis set is used to expand the KS orbitals, so that with the real-space grid adapted in the present work, the properties of liquid water could be obtained very near the complete basis set limit. Structural properties of liquid water were extracted from a 30 ps CPAIMD-BLYP/DVR trajectory at 300 K. The radial distribution functions (RDFs), spatial distribution functions, and hydrogen bond geometry obtained from the CPAIMD-BLYP/DVR simulation are generally in good agreement with the most up to date experimental measurements. Compared to recent ab initio MD simulations based on PW basis sets, less significant overstructuring was found in the RDFs and the distributions of hydrogen bond angles, suggesting that previous plane-wave and Gaussian basis set calculations have exaggerated the tendency toward overstructuring. [ABSTRACT FROM AUTHOR]

It is a well recognized difficult task to simulate the vibrational dynamics of ices using the density functional theory (DFT), and there has thus been rather limited success in modelling the inelastic neutron scattering (INS) spectra for even the simplest structure of ice, ice Ih, particularly in the translational region below 400 cm-1. The reason is partly due to the complex nature of hydrogen bonding (H-bond) among water-water molecules which require considerable improvement of the quantum mechanical simulation methods, and partly owing to the randomness of protons in ice structures which often requires simulation of large super-lattices. In this report, we present the first series of successful simulation results for ice Ih using DFT methods. On the basis of the recent advancement in the DFT programs, we have achieved for the first time theoretical outcomes that not only reproduce the rotational frequencies between 500 to 1200 cm-1 for ice Ih, but also the two optic peaks at ∼240 and 320 cm-1 in the translational region of the INS spectra [J. C. Li, J. Chem. Phys 105, 6733 (1996)]. Besides, we have also investigated the impact of pairwise configurations of H2O molecules on the H-bond and found that different proton arrangements of pairwise H2O in the ice Ih crystal lattice could not alter the nature of H-bond as significantly as suggested in an early paper [J. C. Li and D. K. Ross, Nature (London) 365, 327 (1993)], i.e., reproducing the two experimental optic peaks do not need to invoke the two H-bonds as proposed in the previous model which led to considerable debates. The results of this work suggest that the observed optic peaks may be attributed to the coupling between the two bands of H-O stretching modes in H2O. The current computational work is expected to shed new light on the nature of the H-bonds in water, and in addition to offer a new approach towards probing the interaction between water and biomaterials for which H-bond is essential. [ABSTRACT FROM AUTHOR]

Hydrogen bond interactions strongly affect vibrational properties and frequencies, the most common consequence being a redshift of the stretching vibration involved; there are, however, few exceptions to this general trend. In previous works, we have proved the effectiveness of ab initio simulations combined with wavelet analysis to investigate these effects and put them into relation to structural environment. In this work, we investigate the hydrogen bond effects on the structural and vibrational properties of 1,3-propanediol in acetonitrile by a combined experimental and computational approach. We explain the appearance of two spectral components in the O-H stretching band on the basis of intra- and intermolecular hydrogen bond interactions. We also elucidate the blueshift of the C≡N stretching band as due to a hydrogen bond interaction between the glycol and acetonitrile that modify the electron density distribution inside the CN group. This effect is well reproduced by ab initio molecular dynamics simulations and density functional calculations reported in this work. [ABSTRACT FROM AUTHOR]

We derive the coarse-grained interactions between DNA nucleotides from ab initio total-energy calculations based on density functional theory (DFT). The interactions take into account base and sequence specificity, and are decomposed into physically distinct contributions that include hydrogen bonding, stacking interactions, backbone, and backbone-base interactions. The interaction energies of each contribution are calculated from DFT for a wide range of configurations and are fitted by simple analytical expressions for use in the coarse-grained model, which reduces each nucleotide into two sites. This model is not derived from experimental data, yet it successfully reproduces the stable B-DNA structure and gives good predictions for the persistence length. It may be used to realistically probe dynamics of DNA strands in various environments at the μs time scale and the μm length scale. [ABSTRACT FROM AUTHOR]

We present a simplified implementation of the non-local van der Waals correlation functional introduced by Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] and reformulated by Román-Pérez et al. [Phys. Rev. Lett. 103, 096102 (2009)]. The proposed numerical approach removes the logarithmic singularity of the kernel function. Complete expressions of the self-consistent correlation potential and of the stress tensor are given. Combined with various choices of exchange functionals, five versions of van der Waals density functionals are implemented. Applications to the computation of the interaction energy of the benzene-water complex and to the computation of the equilibrium cell parameters of the benzene crystal are presented. As an example of crystal structure calculation involving a mixture of hydrogen bonding and dispersion interactions, we compute the equilibrium structure of two polymorphs of aspirin (2-acetoxybenzoic acid, C9H8O4) in the P21/c monoclinic structure. [ABSTRACT FROM AUTHOR]

Proton jump processes in the hydration layer on the iso-structural TiO2 rutile (110) and SnO2 cassiterite (110) surfaces were studied with density functional theory molecular dynamics. We find that the proton jump rate is more than three times faster on cassiterite compared with rutile. A local analysis based on the correlation between the stretching band of the O-H vibrations and the strength of H-bonds indicates that the faster proton jump activity on cassiterite is produced by a stronger H-bond formation between the surface and the hydration layer above the surface. The origin of the increased H-bond strength on cassiterite is a combined effect of stronger covalent bonding and stronger electrostatic interactions due to differences of its electronic structure. The bridging oxygens form the strongest H-bonds between the surface and the hydration layer. This higher proton jump rate is likely to affect reactivity and catalytic activity on the surface. A better understanding of its origins will enable methods to control these rates. [ABSTRACT FROM AUTHOR]

We have studied substituent effects on the properties of the intramolecular hydrogen bond of some ortho-hydroxy Schiff bases using density functional theory (DFT) based first-principle molecular dynamics (FPMD) and path integral molecular dynamics. The studied compounds possess a strong intramolecular hydrogen bond (r(O···N) ≤ 2.6 Å), which can be tuned by substitution to either (i) enhance the basicity of the acceptor moiety by induction effects or (ii) decrease the hydrogen bond length through steric repulsion. DFT calculations and FPMD were employed to investigate structural and dynamical properties of the selected molecules, while quantum effects on the structural properties were assessed using path integral FPMD. The simulations were performed in vacuo and in the solid state to study the influence of the environment on the hydrogen bond and spectroscopic properties. We give computational support to the suggestion that induction effects are less effective to tune the intramolecular hydrogen bond properties of the discussed ortho-hydroxy Schiff bases than the steric or the environmental effects. [ABSTRACT FROM AUTHOR]

Density functional theory (DFT) in the commonly used local density or generalized gradient approximation fails to describe van der Waals (vdW) interactions that are vital to organic, biological, and other molecular systems. Here, we propose a simple, efficient, yet accurate local atomic potential (LAP) approach, named DFT+LAP, for including vdW interactions in the framework of DFT. The LAPs for H, C, N, and O are generated by fitting the DFT+LAP potential energy curves of small molecule dimers to those obtained from coupled cluster calculations with single, double, and perturbatively treated triple excitations, CCSD(T). Excellent transferability of the LAPs is demonstrated by remarkable agreement with the JSCH–2005 benchmark database [P. Jurecˇka et al. Phys. Chem. Chem. Phys. 8, 1985 (2006)], which provides the interaction energies of CCSD(T) quality for 165 vdW and hydrogen-bonded complexes. For over 100 vdW dominant complexes in this database, our DFT+LAP calculations give a mean absolute deviation from the benchmark results less than 0.5 kcal/mol. The DFT+LAP approach involves no extra computational cost other than standard DFT calculations and no modification of existing DFT codes, which enables straightforward quantum simulations, such as ab initio molecular dynamics, on biomolecular systems, as well as on other organic systems. [ABSTRACT FROM AUTHOR]

Dynamical properties of liquid water were studied using Car-Parrinello [Phys. Rev. Lett. 55, 2471 (1985)] ab initio molecular dynamics (AIMD) simulations within the Kohn-Sham (KS) density functional theory employing the Becke-Lee-Yang-Parr exchange-correlation functional for the electronic structure. The KS orbitals were expanded in a discrete variable representation basis set, wherein the complete basis set limit can be easily reached and which, therefore, provides complete convergence of ionic forces. In order to minimize possible nonergodic behavior of the simulated water system in a constant energy (NVE) ensemble, a long equilibration run (30 ps) preceded a 60 ps long production run. The temperature drift during the entire 60 ps trajectory was found to be minimal. The diffusion coefficient [0.055 Å2/ps] obtained from the present work for 32 D2O molecules is a factor of 4 smaller than the most up to date experimental value, but significantly larger than those of other recent AIMD studies. Adjusting the experimental result so as to match the finite-sized system used in the present study brings the comparison between theory and experiment to within a factor of 3. More importantly, the system is not observed to become “glassy” as has been reported in previous AIMD studies. The computed infrared spectrum is in good agreement with experimental data, especially in the low frequency regime where the translational and librational motions of water are manifested. The long simulation length also made it possible to perform detailed studies of hydrogen bond dynamics. The relaxation dynamics of hydrogen bonds observed in the present AIMD simulation is slower than those of popular force fields, such as the TIP4P potential, but comparable to that of the TIP5P potential. [ABSTRACT FROM AUTHOR]

The properties of aqueous sulfuric acid have been studied employing density functional theory–based molecular dynamics simulations in conjunction with norm-conserving pseudopotentials. The simulations were carried out for two different concentrations whose molar concentrations were fixed at 0.84 and 10.2 mol/l. The structural features of aqueous sulfuric acid solutions show a strong dependency on the concentration. The Grötthuss-type proton transfer mechanism is not effectively operative at the higher concentration because of the broken hydrogen bond network of water induced by ions generated by the dissociation of sulfuric acid. In addition, to evaluate electrical properties, we carried out a simulation that takes an electric field into account. Results are compared with those of the simulation undertaken with no external electric field. [ABSTRACT FROM AUTHOR]

We have carried out a series of molecular dynamics simulations to investigate the hydrogen bond and residence dynamics of X-–water (X=F, Cl, and I) and pairs in aqueous solutions at a temperature of 673 K. The calculations are done at six different water densities ranging from 1.0 to 0.15 g cm-3. The hydrogen bonds are defined by using a set of configurational criteria with respect to the anion(oxygen)–oxygen and anion(oxygen)–hydrogen distances and the anion(oxygen)–oxygen–hydrogen angle for an anion(water)–water pair. The F-–water hydrogen bonds are found to have a longer lifetime than all other hydrogen bonds considered in the present study. The lifetime of Cl-–water hydrogen bonds is shorter than that of F-–water hydrogen bonds but longer than the lifetime of water–water hydrogen bonds. The lifetimes of I-–water and water–water hydrogen bonds are found to be very similar. Generally, the lifetimes of both anion–water and water–water hydrogen bonds are found to be significantly shorter than those found under ambient conditions. In addition to hydrogen bond lifetimes, we have also calculated the residence times and the orientational relaxation times of water molecules in ion(water) hydration shells and have discussed the correlations of these dynamical quantities with the observed dynamics of anion(water)-water hydrogen bonds as functions of the ion size and density of the supercritical solutions. [ABSTRACT FROM AUTHOR]