Structure of the water/platinum interface--a first principles simulation under bias potential

Ab initio molecular dynamics simulations have been performed on the water/Pt interface. When the surface is neutral, water is found to form a contact layer directing its O atom toward the surface, i.e., O-down configuration. When the surface is negatively biased, the contact layer shows a significant structural change. The O-down configuration is converted mostly to the H-down configuration. As the surface is biased more strongly, we find that a hydrophobic double layer is formed in the contact layer.     

A first principles simulation study of fluctuations of hydrogen bonds and vibrational frequencies of water at liquid-vapor interface

We present a theoretical study of the structure and dynamics of water-vapor interface by means of ab initio molecular dynamics simulations. The inhomogeneous density, hydrogen bond and orientational profiles, voids and vibrational frequency distributions are investigated. We have also studied various dynamical properties of the interface such as diffusion, orientational relaxation, hydrogen bond dynamics and vibrational frequency fluctuations. The diffusion and orientational relaxation of water molecules are found to be faster at the interface which can be correlated with the voids present in the system. The hydrogen bond dynamics, however, is found to be slightly slower at the interface than that in bulk water. The correlations of hydrogen bond relaxation with the dynamics of vibrational frequency fluctuations are also discussed.     

Simulating fluid-phase equilibria of water from first principles

Efficient Monte Carlo algorithms and a mixed-basis set electronic structure program were used to compute from first principles the vapor-liquid coexistence curve of water. A water representation based on the Becke-Lee-Yang-Parr exchange and correlation functionals yields a saturated liquid density of 900 kg/m3 at 323 K and normal boiling and critical temperatures of 350 and 550 K, respectively. An analysis of the structural and electronic properties of the saturated liquid phase shows an increase of the asymmetry of the local hydrogen-bonded structure despite the persistence of a 4-fold coordination and decreases of the molecular dipole moment and of the spread of the lowest unoccupied molecular orbital with increasing temperature.     

Properties of liquid crystal molecules from first principles computer simulation

We explore the valence charge distribution, equilibrium geometry and harmonic force fields of the 4-pentyl-4-cyanobiphenyl (5CB) molecule and the benzene (C6H6) molecule, which provides an important mesogenic fragment, using first principles techniques adapted from large scale electronic structure calculations of periodic solids. We present for the first time accurate structural data for the isolated 5CB molecule and observe subtle broken symmetries relative to the constituent mesogenic fragments. The dynamic properties of these molecules are determined by diagonalization of dynamical matrices, the elements of which are obtained directly from quantum mechanical Hellmann-Feynman forces. Results for both molecules are in excellent agreement with available spectroscopic data, and for benzene are comparable to the most accurate quantum chemistry calculations to date. For 5CB we also present values for the molecular dipole and quadrupole moments.     

Reversible redox reaction and water configuration on a positively charged platinum surface: first principles molecular dynamics simulation

The water dissociation reaction and water molecule configuration on a positively charged platinum (111) surface were investigated by means of first principles molecular dynamics under periodic boundary conditions. Water molecules on the Pt surface were mostly in the O-down orientation but some H-down structures were also found. OH(-) ion, generated by removing H from H(2)O in the bulk region, moved to the Pt surface, on which a positive charge is induced, by a Grotthuss-like proton-relay mechanism and adsorbed on it as OH(Pt). Hydrogen atom exchange between OH(Pt) and a near-by water molecule frequently occurred on the Pt surface and had a low activation energy of the same order as room temperature energy. When a positive charge (7 μC cm(-2)) was added to the Pt surface, H(3)O(+) and OH(Pt) were generated from 2H(2)O on the Pt. This may be coupled with an electron transfer to the Pt electrode [2H(2)O → H(3)O(+) + OH(Pt) + e(-)]. The opposite reaction was also observed on the same charged surface during a simulation of duration about 10 ps; it is a reversible redox reaction. When further positive charge (14 μC cm(-2)) was added, the reaction shifted to the right hand side completely. Thus, this one-electron transfer reaction, which is a part of the oxygen electrode reaction in fuel cells and water electrolysis, was confirmed to be a low activation energy process.     

A first principles approach to optimal control

This article shows that by using ab initio or first principle calculations it is possible to obtain reliable ingredients needed to simulate pump-probe and optimal control experiments. Our experimental challenge is to elucidate the reaction mechanism behind an optimal pulse tailored to maximize ionization in the system CpMn(CO)₃, while avoiding CO dissociation. Starting from MRCI/CASSCF potential energy curves calculated along the relevant CO fragmentation channel, we use the resulting MRCI wave function to estimate non-adiabatic couplings, as well as neutral-to-neutral and neutral-to-ionic dipole couplings. The state-of-the-art potentials and couplings serve to perform wave packet propagations which simulate the femtosecond pump-probe spectra that explain the features shown in the experimental optimal pulse.     

A first principles study of water oxidation catalyzed by a tetraruthenium-oxo core embedded in polyoxometalate ligands

We present a computational study addressing the catalytic cycle of a recently-synthesized all-inorganic homogeneous catalyst capable to promote water oxidation with low overpotential and high turnover frequency [Sartorel et al., J. Am. Chem. Soc., 2008, 130, 5006; Geletii et al., Angew. Chem., Int. Ed., 2008, 47, 3896]. This catalyst consists of a tetraruthenium-oxo core [Ru(4)O(4)(OH)(2)·(H(2)O)(4)](6+)capped by two polyoxometalate [SiW(10)O(36)](8-) units. The reaction mechanism underpinning its efficiency is currently under debate. We study a reaction cycle involving four consecutive proton-coupled electron transfer (PCET) processes that successively oxidize the four Ru(IV)-H(2)O units of the initial state (S(0)) to the four Ru(V)-OH centers of the activated intermediate (S(4)). The energetics of these electrochemical processes as well as the structural and electronic properties of the reaction intermediates are studied with ab initio Density Functional Theory (DFT) calculations. After characterizing these reaction intermediates in the gas phase, we show that the solvated tetraruthenate core undergoes a solvent-induced structural distortion that brings the predicted molecular geometry to excellent agreement with the experimental X-ray diffraction data. The calculated electronic properties of the catalyst are instead weakly dependent on the presence of the solvent. The frontier orbitals of the initial state as well as the electronic states involved in the PCET steps are shown to be localized on the tetraruthenium-oxo core. The reaction thermodynamics predicted for the intermediate reaction steps is in good agreement with the available cyclic voltammetry measurements up to S(3), but the calculated free energy difference between the initial and the activated state (S(0)/S(4)) turns out to be significantly lower than the thermodynamic limit for water oxidation. Since the oxidizing power of the S(0)/S(4) couple is not sufficient to split water, we suggest that promoting this reaction would require cycling between higher oxidation states.     

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.     

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.     

Optical properties of silicon clusters in the presence of water: a first principles theoretical analysis

We investigate the impact of water on the optical absorption of prototypical silicon clusters. Our clusters contain 5 silicon atoms, tetrahedrally coordinated and passivated with either hydrogen or oxygen. We approach this complex problem by assessing the contributions of three factors: chemical reactivity, thermal equilibration, and dielectric screening. We find that the silanone (Si=O) functional group is not chemically stable in the presence of water and exclude this as a source of significant red shift in absorption in aqueous environments. We perform first principles molecular dynamics simulations of the solvation of a chemically stable, oxygenated silicon cluster with explicit water molecules at 300 K. We find a systematic 0.7 eV red shift in the absorption gap of this cluster, which we attribute to thermally induced fluctuations in the molecular structure. Surprisingly, we find no observable screening impact of the solvent, in contrast with consistent blue shifts observed for similarly sized organic molecules in polar solvents. The predicted red shift is expected to be significantly smaller for larger Si quantum dots produced experimentally, guaranteeing that their vacuum optical properties are preserved even in aqueous environments.     

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.     

A first principles study on organic molecule encapsulated boron nitride nanotubes

The electronic structures of boron nitride nanotubes (BNNTs) doped with organic molecules are investigated using density functional theory. An electrophilic molecule introduces acceptor states in the wide gap of BNNT close to the valence band edge, which makes the doped system a p-type semiconductor. However, with typical nucleophilic organic molecules encapsulation, only deep occupied molecular states but no shallow donor states are observed. There is a significant electron transfer from a BNNT to an electrophilic molecule, while the charge transfer between a nucleophilic molecule and a BNNT is negligible. When both electrophilic and nucleophilic molecules are encapsulated in the same BNNT, a large charge transfer between the two kinds of molecules occurs. The resulting small energy gap can strongly modify the transport and optical properties of the system.     

A first principles theory of nuclear magnetic resonance J-coupling in solid-state systems

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has also been applied to the silicophosphate, Si(5)O(PO(4))(6), giving (31)P-(29)Si-couplings which are in excellent agreement with experiment.     

A study of gas phase and surface formaldehyde polymerisation from first principles

First principles computational studies have been carried out on the gas phase and surface polymerisation of formaldehyde. Any gas phase polymerisation in dark conditions is shown to occur by way of neutral species, as opposed to earlier suggestions that a charged species may be involved. The resulting paraformaldehyde chains are shown to be stable only when fully bonded at both chain ends, as opposed to some of the earlier suggested mechanisms. Polymerisation is shown to be thermodynamically favourable in the gas phase, but kinetically limited. Derived rate constants are insufficient to account for experimental results, in agreement with prior studies suggesting that gas phase polymerisation actually occurs on surfaces rather than directly in the gas phase. The TiO(2) (110) surface is employed as a model surface to study formaldehyde adsorption and possible polymerisation mechanisms. Formaldehyde monomers are shown to weakly bind by way of coordinate bonding through the carbonyl group to surface Ti. A particularly strong dimerisation configuration is also found where the two ends of paraformaldehyde are terminated by a surface Ti and surface O, consistent with the earlier observation that paraformaldehyde chains are only stable when saturated at both ends.     

Guanine crystals: a first principles study

We report first principles density functional theory studies on the basic ground state characteristics, dynamic properties, and the electronic structure of guanine crystals. The effect of water molecules within the crystal is studied in detail, and we discuss their influence on the structural, vibrational, and electronic properties. The geometries calculated for various crystal structures are compared with gas-phase calculations and available experimental data. Phonon frequencies and eigenvectors are predicted for intermolecular and intramolecular lattice vibrations. Vibrational and electronic density-of-states are presented and analyzed. The electronic band structure near the fundamental gap is calculated from the Kohn-Sham approach. We find that the former molecular HOMO states form a dispersive band in the pi-pi stacking direction upon condensation resulting in a large bandwidth of 0.83 eV. Consequences for the charge transport in layered van der Waals bonded organic molecular crystals are discussed.