Hydration properties of magnesium and calcium ions from constrained first principles molecular dynamics

We studied the solvation structures of the divalent metal cations Mg(2+) and Ca(2+) in ambient water by applying a Car-Parrinello-based constrained molecular dynamics method. By employing the metal-water oxygen coordination number as a reaction coordinate, we could identify distinct aqua complexes characterized by structural variations of the first coordination shell. In particular, our estimated free-energy profile clearly shows that the global minimum for Mg(2+) is represented by a rather stable sixfold coordination in the octahedral arrangement, in agreement with experiments. Conversely, for Ca(2+) the free-energy curve shows several shallow local minima, suggesting that the hydration structure of Ca(2+) is highly variable. Implications for water exchange reactions are also discussed.     

Communication: Hydration structure and polarization of heavy alkali ions: A first principles molecular dynamics study of Rb(+) and Cs(+)

Hydration structure and polarization of Rb(+) and Cs(+) in liquid water at ambient conditions were studied by first principles molecular dynamics. Our systematic analysis of the relevant electronic structures, based on maximally localized Wannier functions, revealed that the dipole moment of H(2)O molecules in the first solvation shell of the ions slightly increases with increasing the atomic number. We also found that the polarization of heavy alkali ions, particularly Cs(+), tends to stabilize a peculiar asymmetric hydration structure with relevant consequences in the extraction of the harmful (137)Cs resulting from nuclear wastes.     

High pressure reactivity of propene by first principles molecular dynamics calculations

The reactivity of propene under high pressure has been investigated in the framework of Car-Parrinello molecular dynamics. Changes in structural and electronic properties due to pressure have been analyzed in systems with a density ranging from 0.855 to 2.151 g/cm(3). A ionic collective mechanism which leads to the formation of oligomers has been found by both spin restricted and spin polarized formalism. The maximally localized Wannier centers analysis has allowed us to characterize the addition scheme and to identify a Wannier center with a high spread value involved in the formation of the principal reaction products.     

A recent development in computational chemistry: chemical reactions from first principles molecular dynamics simulations

First-principles molecular dynamics simulations have recently been found an effective tool to study a large variety of chemical problems. Finite temperature simulations reveal unique information, including explicit dynamical effects and the evaluation of proper free energy differences. Moreover, dynamics simulations reveal information on the flexibility of molecular systems, and elucidate, often otherwise inaccessible, mechanistic details of chemical reactions. In addition this methodology allows the study of larger, periodic, systems, revealing computationally unique information which may be directly compared to experiments on realistic chemical systems. A variety of examples will be given, although most focus on the important field of catalysis.     

Hydration structures of U(III) and U(IV) ions from ab initio molecular dynamics simulations

We apply DFT+U-based ab initio molecular dynamics simulations to study the hydration structures of U(III) and U(IV) ions, pertinent to redox reactions associated with uranium salts in aqueous media. U(III) is predicted to be coordinated to 8 water molecules, while U(IV) has a hydration number between 7 and 8. At least one of the innershell water molecules of the hydrated U(IV) complex becomes spontaneously deprotonated. As a result, the U(IV)-O pair correlation function exhibits a satellite peak at 2.15 A? associated with the shorter U(IV)-(OH(-)) bond. This feature is not accounted for in analysis of extended x-ray absorption fine structure and x-ray adsorption near edge structure measurements, which yield higher estimates of U(IV) hydration numbers. This suggests that it may be useful to include the effect of possible hydrolysis in future interpretation of experiments, especially when the experimental pH is close to the reported hydrolysis equilibrium constant value.     

Hydration of superoxide studied by molecular dynamics simulation

Molecular dynamics was used to study the hydration of superoxide (O). The Helmholtz free energy of hydration of O was estimated by the thermodynamic integration method. The diffusion of O and the water structure around O were also studied. Two water models were used in the calculations and the results were compared to experiments.     

The mixed alkali effect in ionically conducting glasses revisited: a study by molecular dynamics simulation

When more than two kinds of mobile ions are mixed in ionic conducting glasses and crystals, there is a non-linear decrease of the transport coefficients of either type of ion. This phenomenon is known as the mixed mobile ion effect or Mixed Alkali Effect (MAE), and remains an unsolved problem. We use molecular dynamics simulation to study the complex ion dynamics in ionically conducting glasses including the MAE. In the mixed alkali lithium-potassium silicate glasses and related systems, a distinct part of the van Hove functions reveals that jumps from one kind of site to another are suppressed. Although, consensus for the existence of preferential jump paths for each kind of mobile ions seems to have been reached amongst researchers, the role of network formers and the number of unoccupied ion sites remain controversial in explaining the MAE. In principle, these factors when incorporated into a theory can generate the MAE, but in reality they are not essential for a viable explanation of the ion dynamics and the MAE. Instead, dynamical heterogeneity and "cooperativity blockage" originating from ion-ion interaction and correlation are fundamental for the observed ion dynamics and the MAE. Suppression of long range motion with increased back-correlated motions is shown to be a cause of the large decrease of the diffusivity especially in dilute foreign alkali regions. Support for our conclusion also comes from the fact that these features of ion dynamics are common to other ionic conductors, which have no glassy networks, and yet they all exhibit the MAE.     

Mobility of a single alkali metal atom on fullerene C60: first principles molecular dynamical study

The dynamical behavior of a single K or Na alkali metal atom on the surface of a C(60) molecule is investigated via Car-Parrinello molecular dynamics simulations in a temperature range up to 300 K. These provide direct evidence for the heteroatom motion, postulated earlier in pioneering experiments, and show that an alkali metal atom can move both on the surface and radially outward from the surface, resulting in a dynamics ranging from diffusive to free orbital motion, on time-scales of up to a few picoseconds.     

Thermal conductivity of molten alkali halides from equilibrium molecular dynamics simulations

The thermal conductivity of molten sodium chloride and potassium chloride has been computed through equilibrium molecular dynamics Green-Kubo simulations in the microcanonical ensemble (N,V,E). In order to access the temperature dependence of the thermal conductivity coefficient of these materials, the simulations were performed at five different state points. The form of the microscopic energy flux for ionic systems whose Coulombic interactions are calculated through the Ewald method is discussed in detail and an efficient formula is used by analogy with the methods used to evaluate the stress tensor in Coulombic systems. The results show that the Born-Mayer-Huggins-Tosi-Fumi potential predicts a weak negative temperature dependence for the thermal conductivity of NaCl and KCl. The simulation results are in agreement with part of the experimental data available in the literature with simulation values generally overpredicting the thermal conductivity by 10%-20%.     

Hydration shell structure and dynamics of curium(III) in aqueous solution: first principles and empirical studies

Results of ab initio molecular dynamics (AIMD), quantum mechanics/molecular mechanics (QM/MM), and classical molecular dynamics (CMD) simulations of Cm(3+) in liquid water at a temperature of 300 K are reported. The AIMD simulation was based on the Car-Parrinello MD scheme and GGA-PBE formulation of density functional theory. Two QM/MM simulations were performed by treating Cm(3+) and the water molecules in the first shell quantum mechanically using the PBE (QM/MM-PBE) and the hybrid PBE0 density functionals (QM/MM-PBE0). Two CMD simulations were carried out using ab initio derived pair plus three-body potentials (CMD-3B) and empirical Lennard-Jones pair potential (CMD-LJ). The AIMD and QM/MM-PBE simulations predict average first shell hydration numbers of 8, both of which disagree with recent experimental EXAFS and TRLFS value of 9. On the other hand, the average first shell hydration numbers obtained in the QM/MM-PBE0 and CMD simulations was 9, which agrees with experiment. All the simulations predicted an average first shell and second shell Cm-O bond distance of 2.49-2.53 ? and 4.67-4.75 ? respectively, both of which are in fair agreement with corresponding experimental values of 2.45-2.48 and 4.65 ?. The geometric arrangement of the 8-fold and 9-fold coordinated first shell structures corresponded to the square antiprism and tricapped trigonal prisms, respectively. The second shell hydration number for AIMD QM/MM-PBE, QM/MM-PBE0, CMD-3B, and CMD-LJ, were 15.8, 17.2, 17.7, 17.4, and 16.4 respectively, which indicates second hydration shell overcoordination compared to a recent EXAFS experimental value of 13. Save the EXAFS spectra CMD-LJ simulation, all the computed EXAFS spectra agree fairly well with experiment and a clear distinction could not be made between configurations with 8-fold and 9-fold coordinated first shells. The mechanisms responsible for the first shell associative and dissociative ligand exchange in the classical simulations have been analyzed. The first shell mean residence time was predicted to be on the nanosecond time scale. The computed diffusion constants of Cm(3+) and water are in good agreement with experimental data.     

First principles molecular dynamics of molten NaCl

First principles Hellmann-Feynman molecular dynamics (HFMD) results for molten NaCl at a single state point are reported. The effect of induction forces on the structure and dynamics of the system is studied by comparison of the partial radial distribution functions and the velocity and force autocorrelation functions with those calculated from classical MD based on rigid-ion and shell-model potentials. The first principles results reproduce the main structural features of the molten salt observed experimentally, whereas they are incorrectly described by both rigid-ion and shell-model potentials. Moreover, HFMD Green-Kubo self-diffusion coefficients are in closer agreement with experimental data than those predicted by classical MD. A comprehensive discussion of MD results for molten NaCl based on different ab initio parametrized polarizable interionic potentials is also given.     

Metastable dissociation dynamics of molecular cluster ions

Metastable uni-cluster dissociation for several hydrogen-bonded and van der Waals cluster ions are observed via resonance-enhanced two-photon ionization reflectron time-of-flight (TOF) mass spectrometry. All of the cluster ions studied show evaporation of a single molecule from the respective parent cluster ions as dominant metastable decay processes. Furthermore, the averaged metastable evaporation rate constants (k _evap) of these cluster ions in a fixed time domain of 0.2–50 µs are obtained by analyzing the relative intensity of metastable ion peaks due to evaporation in the acceleration and the field-free drift regions of the TOF mass spectrometer. An intensity anomaly in some of the observed metastable ion peaks, indicative of magic number stability of the cluster ion, is also presented.     

First principles molecular dynamics simulation of liquid rubidium

We present a first principles molecular dynamics simulation of liquid rubidium. The atomic forces are obtained from a quantum mechanical calculation of its electronic structure within the local density approximation of the density functional formalism and using the pseudopotential plane-wave method. We compare our results with the structure and dynamics predicted by classical molecular dynamics simulations.     

Structural fluctuation and dynamics of ribose puckering in aqueous solution from first principles

Using the method of ab initio molecular dynamics, we examine the structural fluctuation and the low-frequency dynamics of beta-ribofuranose puckering in aqueous solution. Our analysis suggests that the distance between the anomeric and hydroxymethyl oxygens is a simple relevant geometrical parameter that dynamically correlates with the phase angle in the north region. The time-frequency analysis using the Hilbert-Huang transform also confirms the correlation, and most of the instantaneous frequencies for the phase angle and the above distance are found to be concentrated on the region below about 100 cm(-1). Our analysis of ab initio molecular dynamics trajectories suggests that the molecular origin of the hydration effects on the low-frequency dynamics of beta-ribofuranose puckering is closely related to this correlation and thus primarily attributed to the relatively local interactions among the anomeric and hydroxymethyl oxygens and the surrounding water molecules near them. Additionally, we discuss the difference in the low-frequency dynamics of beta-ribofuranose puckering between two hydroxymethyl rotamers.     

Tryptophan–water interaction in Monellin: Hydration patterns from molecular dynamics simulation

Femtosecond spectroscopy carried out earlier on Monellin and some other systems has given insights into the hydration dynamics of the proteins. In the present work, molecular dynamics simulations have been performed on Monellin to study the hydration dynamics. A method has been described to follow up the molecular events of the protein–water interactions in detail. The time constants of the survival correlation function match well with the reported experimental values. This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general.     

Understanding the solubility of triamino-trinitrobenzene in hydrous tetramethylammonium fluoride: a first principles molecular dynamics simulation study

With the aim to understand the relatively high solubility of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), an important energetic material with a high degree of inter- and intra-molecular hydrogen bonding, in fluoride anion containing ionic liquids (ILs), first principles molecular dynamics simulations in the isobaric-isothermal ensemble were carried out for a system using hydrous tetramethylammonium fluoride as the prototypical solvent. Simulations initiated from both molecular TATB and its Meisenheimer complex (i.e., a σ-complex of the fluoride and the electrophilic ring of TATB) yield a Zundel-type complex where a proton is shared between an amino group and an F(-) ion, whereas the Meisenheimer complex is found to be only transiently stable. An analysis of the electronic structure probing the Wannier function centers supports the finding of a proton-sharing complex with a three-center four-electron like bond. The Zundel-type complex also yields an electronic absorption spectrum consistent with the experimentally observed color change. This study provides evidence that the remarkable solubility of otherwise hard-to-dissolve molecular crystals in ILs can be aided by chemical modification of the solute.     

Electrochemical potentials from first principles

The electrochemical potential is the fundamental parameter in the theory of electrochemistry. Not only does it determine the position of electrochemical equilibria but also it acts as the driving force for electron transfer reactions, diffusion-migration phenomena, and phase transformations of all kinds. In the present work, the electrochemical potential is defined as the total work done in transferring a single particle of a substance from a universal reference state to a specified location, at constant temperature and pressure. It is the sum of two scalar fields: the electrostatic potential energy and the chemical potential energy. The electrochemical potential is widely underutilized within the fields of solid-state science and electrochemical engineering. For historical reasons, many authors prefer to analyze driving forces in terms of electrode potentials, concentration gradients, or Gibbs free energies. In this paper, the author provides a short introduction to the electrochemical potential and then shows how some of the major branches of electrochemistry can benefit from using it. Topics examined include the Volta potential difference, the membrane potential difference, the scanning Kelvin probe microscope, the electromotive force, the proton motive force, and the activation of electron transfer.

     

Tafel slopes from first principles

Tafel slopes for multistep electrochemical reactions are derived from first principles. The derivation takes place in two stages. First, Dirac’s perturbation theory is used to solve the Schrödinger equation. Second, current–voltage curves are obtained by integrating the single-state results over the full density of states in electrolyte solutions. Thermal equilibrium is assumed throughout. Somewhat surprisingly, it is found that the symmetry factor that appears in the Butler–Volmer equation is different from the symmetry factor that appears in electron transfer theory, and a conversion formula is given. Finally, the Tafel slopes are compiled in a convenient look-up table.     

Beta-sheet preferences from first principles

The natural amino acids have different preferences of occurring in specific types of secondary protein structure. Simulations are performed on periodic model beta-sheets of 14 different amino acids, at the level of density functional theory, employing the generalized gradient approximation. We find that the statistically observed beta-sheet propensities correlate very well with the calculated binding energies. Analysis of the calculations shows that the beta-sheet propensities are determined by the local flexibility of the individual polypeptide strands.     

Vibrational spectral diffusion and hydrogen bond dynamics in heavy water from first principles

We present a first-principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in heavy water without using any empirical model potentials. The calculations are based on ab initio molecular dynamics simulations for trajectory generation and a time series analysis using the wavelet method for frequency calculations. It is found that, in deuterated water, although a one-to-one relation does not exist between the instantaneous frequency of an OD bond and the distance of its associated hydrogen bond, such a relation does hold on average. The dynamics of spectral diffusion is investigated by means of frequency-time correlation and spectral hole dynamics calculations. Both of these functions are found to have a short-time decay with a time scale of approximately 100 fs corresponding to dynamics of intact hydrogen bonds and a slower long-time decay with a time constant of approximately 2 ps corresponding to lifetimes of hydrogen bonds. The connection of the slower time scale to the dynamics of local structural relaxation is also discussed. The dynamics of hydrogen bond making is shown to have a rather fast time scale of approximately 100 fs; hence, it can also contribute to the short-time dynamics of spectral diffusion. A damped oscillation is also found at around 150-200 fs, which is shown to have come from underdamped intermolecular vibrations of a hydrogen-bonded water pair. Such assignments are confirmed by independent calculations of power spectra of intermolecular motion and hydrogen bond kinetics using the population correlation function formalism. The details of the time constants of frequency correlations and spectral shifts are found to depend on the frequencies of chosen OD bonds and are analyzed in terms of the dynamics of hydrogen bonds of varying strengths and also of free non-hydrogen-bonded OD groups.