DFT:B3LYP ab initio molecular dynamics study of the Zundel and Eigen proton complexes, H5O2+ and H9O4+, in the triplet state in gas phase and solution

DFT:B3LYP ab initio molecular dynamics (MD) approach is used to elucidate the properties of the Zundel and Eigen, H5O2+ and H9O4+, proton complexes in the triplet state. The simulation considers the complexes in the gas phase (isolated complexes) and inside the clusters composed of 32, 64, and 128 water molecules, mimicking the behavior of aqueous solutions. MD simulations reveal three distinct periods. For the complex in solutions, the periods are smoothed out. The H5O2+ and H9O4+ complexes in the triplet state undergo structural rearrangements, which eventually result in hydrogen elimination. For the H5O2+, the hydrogen is eliminated from the center of the water cluster, whereas for the H9O4+ it is removed from a near-surface water molecule. The rate of hydrogen elimination decreases with increasing the number of water molecules surrounding the complex.     

Dynamics and structural changes of small water clusters on ionization

Despite utmost importance in understanding water ionization process, reliable theoretical results of structural changes and molecular dynamics (MD) of water clusters on ionization have hardly been reported yet. Here, we investigate the water cations [(H₂O)_{n = 2–6}⁺] with density functional theory (DFT), Möller–Plesset second‐order perturbation theory (MP2), and coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The complete basis set limits of interaction energies at the CCSD(T) level are reported, and the geometrical structures, electronic properties, and infrared spectra are investigated. The characteristics of structures and spectra of the water cluster cations reflect the formation of the hydronium cation moiety (H₃O⁺) and the hydroxyl radical. Although most density functionals fail to predict reasonable energetics of the water cations, some functionals are found to be reliable, in reasonable agreement with high‐level ab initio results. To understand the ionization process of water clusters, DFT‐ and MP2‐based Born‐Oppenheimer MD (BOMD) simulations are performed on ionization. On ionization, the water clusters tend to have an Eigen‐like form with the hydronium cation instead of a Zundel‐like form, based on reliable BOMD simulations. For the vertically ionized water hexamer, the relatively stable (H₂O)₅⁺ (5sL4A) cluster tends to form with a detached water molecule (H₂O). © 2013 Wiley Periodicals, Inc.     

Theoretical studies on the hydration of formic acid by ab initio and ABEEMσπ fluctuating charge model

The interaction between formic acid (FA) and water was systemically investigated by atom-bond electronegativity equalization method fused into molecular mechanics (ABEEMσπ/MM) and ab initio methods. The geometries of 20 formic acid–water complexes (FA–water) were obtained using B3LYP/aug-cc-pVTZ level optimizations, and the energies were determined at the MP2/aug-cc-pVTZ level with basis set superposition error (BSSE) and zero-point vibrational energy (ZPVE) corrections. The ABEEMσπ potential model gives reasonable properties of these clusters when compared with the present ab initio data. For interaction energies, the root mean square deviation is 0.74 kcal/mol, and the linear coefficient reaches 0.993. Next, FA in aqueous solution was also studied. The hydrogen-bonding pattern due to the interactions with water has been analyzed in detail. Furthermore, the ABEEMσπ charges changed when H₂O interacted with the FA molecule, especially at the sites where the hydrogen bonds form. These results show that the ABEEMσπ fluctuating charge model is fine giving the overall characteristic hydration properties of FA–water systems in good agreement with the high-level ab initio calculations.     

Water trimer cation

By using density functional theory (DFT) and high-level ab initio theory, we have investigated the structure, interaction energy, electronic property, and IR spectra of the water trimer cation [(H₂O) ₃ ⁺ ]. Two structures of the water trimer cation [the H₃O⁺ containing linear (3Lp) structure versus the ring (3OO) structure] are compared. For the complete basis set (CBS) limit of coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)], the 3Lp structure is 11.9?kcal/mol more stable than the 3OO structure. This indicates that the ionization of water clusters produce the hydronium cation moiety (H₃O⁺) and the hydroxyl radical. It is interesting to note that the calculation results of the water trimer cation vary seriously depending on the calculation level. At the level of M?ller?CPlesset second-order perturbation (MP2) theory, the stability of 3OO is underestimated due to the underestimated O??O hemibonding energy. This stability is also underestimated even for the CCSD(T) single point calculations on the MP2-optimized geometry. For the 3OO structure, the MP2 and CCSD(T) calculations give closed-ring structures with a hemi-bond between two O atoms, while the DFT calculations show open-ring structures due to the overestimated O??O hemibonding energy. Thus, in order to obtain reliable stabilities and frequencies of the water trimer cation, the CCSD(T) geometry optimizations and frequency calculations are necessary. In this regard, the DFT functionals need to be improved to take into account the proper O??O hemibonding energy.     

Density functional complete study of hydrogen bonding between the water molecule and the hydroxyl radical (H2O · HO)

The hydrogen bonding complexes formed between the H₂O and OH radical have been completely investigated for the first time in this study using density functional theory (DFT). A larger basis set 6‐311++G(2d,2p) has been employed in conjunction with a hybrid density functional method, namely, UB3LYP/6‐311++G(2d,2p). The two degenerate components of the OH radical ²Π ground electronic state give rise to independent states upon interaction with the water molecule, with hydrogen bonding occurring between the oxygen atom of H₂O and the hydrogen atom of the OH radical. Another hydrogen bond occurs between one of the H atoms of H₂O and the O atom of the OH radical. The extensive calculation reveals that there is still more hydrogen bonding form found first in this investigation, in which two or three hydrogen bonds occur at the same time. The optimized geometry parameter and interaction energy for various isomers at the present level of theory was estimated. The infrared (IR) spectrum frequencies, IR intensities, and vibrational frequency shifts are reported. The estimates of the H₂O · OH complex's vibrational modes and predicted IR spectra for these structures are also made. It should be noted that a total of 10 stationary points have been confirmed to be genuine minima and transition states on the potential energy hypersurface of the H₂O · HO system. Among them, four genuine minima were located. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Density functional theory and post-Hartree-Fock studies on molecular structure and harmonic vibrational spectrum of formaldehyde

Density functional theory (DFT) with the Becke's three-parameter exchange correlation functional and the functional of Lee, Yang and Parr, gradient-corrected functionals of Perdew, and Perdew and Wang [the DFT(B3LYP), DFT(B3P86) and DFT(B3PW91) methods, respectively], and several levels of conventional ab initio post-Hartree-Fock theory (second- and fourth-order perturbation theory M?ller-Plesset MP2 and MP4(SDTQ), coupled cluster with the single and double excitations (CCSD), and CCSD with perturbative triple excitation [CCSD(T)], configuration interaction with the single and double excitations [CISD], and quadratic configuration interaction method [QCISD(T)], using several basis sets [ranging from a simple 6-31G(d,p) basis set to a 6-311+ +G(3df, 2pd) one], were applied to study of the molecular structure (geometrical parameters, rotational constants, dipole moment) and harmonized infrared (IR) spectrum of formaldehyde (CH₂O). High-level ab initio methods CCSD(T) and QCISD(T) with the 6-311+ +G(3df, 2pd) predict correctly molecular parameters, vibrational harmonic wavenumbers and the shifts of the harmonic IR spectrum of ¹²CH₂ ¹⁶O upon isotopic substitution. Received: 30 January 1997 / Accepted: 7 May 1997     

Direct ab initio MD study on the hydrogen abstraction reaction of triplet state acetone from methanol molecule

Solvent re-orientation process of triplet acetone/methanol complex and intermolecular hydrogen atom abstraction reaction on the triplet state energy surface, (CH₃)₂C=O (T₁) + CH₃OH → (CH₃)₂C–OH + CH₂OH in gas phase, have been investigated by means of density functional theory (DFT) and direct ab initio molecular dynamics (MD) methods. The static DFT calculation of hydrogen abstraction reaction at the T₁ state showed that the transition state is 16.4 and 30.9 kcal/mol lower than the energy levels of S₁ and S₂ states, respectively, and 9.2 kcal/mol higher than the bottom of T₁ state. The product state, (CH₃)₂C–OH⋯CH₂OH, is 8.4 kcal/mol lower in energy than the level of T₁ state. The direct ab initio MD calculation showed that the product is rapidly formed within 150 fs and the separated products (CH₃)₂C–OH + CH₂OH were formed. The mechanism of reaction dynamics of the triplet acetone/methanol complex was discussed on the basis of theoretical results.     

Opposing Electronic and Nuclear Quantum Effects on Hydrogen Bonds in H2O and D2O

The effect of extending the O−H bond length(s) in water on the hydrogen-bonding strength has been investigated using static ab initio molecular orbital calculations. The “polar flattening” effect that causes a slight σ-hole to form on hydrogen atoms is strengthened when the bond is stretched, so that the σ-hole becomes more positive and hydrogen bonding stronger. In opposition to this electronic effect, path-integral ab initio molecular-dynamics simulations show that the nuclear quantum effect weakens the hydrogen bond in the water dimer. Thus, static electronic effects strengthen the hydrogen bond in H₂O relative to D₂O, whereas nuclear quantum effects weaken it. These quantum fluctuations are stronger for the water dimer than in bulk water.     

Theoretical study on structures and vibrational spectra of M+(H2O)Ar (M = Cu,Ag, Au)

A theoretical study on the structures and vibrational spectra of M⁺(H₂O)Ar_0‐1 (M = Cu, Ag, Au) complexes was performed using ab initio method. Geometrical structures, binding energies (BEs), OH stretching vibrational frequencies, and infrared (IR) absorption intensities are investigated in detail for various isomers with Ar atom bound to different binding sites of M⁺(H₂O). CCSD(T) calculations predict that BEs are 14.5, 7.5, and 14.4 kcal/mol for Ar atom bound to the noble metal ion in M⁺(H₂O)Ar (M = Cu, Ag, Au) complexes, respectively, and the corresponding values have been computed to be 1.5, 1.3, and 2.1 kcal/mol when Ar atom attaches to a H atom of water molecule. The former structure is predicted to be more stable than the latter structure. Moreover, when compared with the M⁺(H₂O) species, tagging Ar atom to metal cation yields a minor perturbation on the IR spectra, whereas binding Ar atom to an OH site leads to a large redshift in OH stretching vibrations. The relationships between isomers and vibrational spectra are discussed. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Vibrational analysis of free and hydrogen bonded complexes of nicotinamide and picolinamide

Harmonic and anharmonic vibrations of free nicotinamide (NIA) and picolinamide (PIA) molecules together with their hydrogen bonded complexes H₂O–NIA and H₂O–PIA have been studied by means of density functional method. The calculation results of the vibrational spectra of free molecules have been investigated and are compared to the available experimental spectra. The vibrational wavenumbers of both molecules have also been calculated by polarizable continuum model (PCM) that represents the solvent as a polarizable continuum and places the solute in a cavity within the solvent (water is chosen as the solvent in this study). The results of PCM calculations and the H₂O–NIA, H₂O–PIA complexes, are used to investigate the H-bonding interactions of both molecules with the water molecule. The harmonic wavenumbers have been scaled by proper factors obtained by comparing the observed versus calculated wavenumbers and it is shown that anharmonic corrections on the vibrational spectra provided a better agreement between the observed and calculated wavenumbers compared to the results obtained by scaling factor method.     

Theoretical characterization of H2 adsorption on AuPt clusters

Linear and triangular AuPt_n cluster reactions with H₂ are studied theoretically, using ab initio multiconfiguration self-consistent field (MC-SCF) calculations, followed by extensive multireference configuration interaction (MR-CI) variational and perturbative. Both the linear dimer and the triangular trimers capture the hydrogen molecule by the Pt cluster side and by the Au cluster side. Gold has an electronic effect on the Pt activity, more important than a geometrical one, poisoning the Pt activity to dissociate H₂ and lowering the adsorption heats. This effect is stronger in the AuPt dimer, where the H₂ capture occurs only at the molecular level, without showing hydrogen bond dissociation. The trimers look more active, relaxing the H(SINGLE BOND)H bond until breakage. No activation barriers are observed in all the cases considered. The hydrogen molecule is not able to cross the clusters due to the large barriers present. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 399–409, 1998     

ChemInform Abstract: The Infrared Spectrum of F-·H2O.

The geometry and vibrational frequencies of the F-·H₂O system are determined by ab initio calculations (SCF and CISD levels of theory).     

Ring-polymer Molecular Dynamics Simulation for the Adsorption of H2 on Ice Clusters (H2O)n (n=8, 10, and 12)

In the interstellar medium, the H₂ adsorption and desorption on the solid water ice are crucial for chemical and physical processes. We have recently investigated the probabilities of H₂ sticking on the (H₂O)₈ ice, which has quadrilateral surfaces. We have extended the previous work using classical MD and ring-polymer molecular dynamics (RPMD) simulations to the larger ice clusters, (H₂O)₁₀ and (H₂O)₁₂, which have pentagonal and hexagonal surfaces, respectively. The H₂ sticking probabilities decreased as the temperature increased for both cluster cases, whereas the cluster-size-independent profiles were observed. It is thought that the size independence of the probabilities is qualitatively understood from the similar binding energies for all the three cluster systems. Furthermore, the RPMD sticking probabilities are smaller than the classical ones because of the reduction in the binding energies owing to nuclear quantum effects, such as vibrational quantization.     

Molecular analysis of water clusters: Calculation of the cluster structures and vibrational spectrum using density functional theory

Density functional theory (DFT) is applied to obtain absorption spectra at THz frequencies for molecular clusters of H₂O. The vibrational modes of the clusters are calculated. Coupling among molecular vibrational modes explains their spectral features associated with THz excitation. THz excitation is associated with vibrational frequencies which are here calculated within the DFT approximation of electronic states. This is done for both isolated molecules and collections of molecules in a cluster. The principal result of the paper is that a crystal-like cluster of 38 water molecules together with a continuum solvent background is sufficient to replicate well the experimental vibrational frequencies.     

Electron density of H2O in a crystalline environment

The influence of the surroundings on the water molecule m LiOH·H₂O has been studied by ab initio MO LCAO SCF techniques. The main features are an enhancement of the polarity of the molecule and a decrease in the lone-pair density.     

Reaction mechanism of the preferential oxidation of the CO reaction in an H<Subscript>2</Subscript> stream over Cu–Ni bimetallic catalysts: A computational study

The preferential oxidation (PROX, CO + H₂ + O₂ → CO₂ + H₂O) of the CO reaction in an H₂ stream is the simplest and most cost-effective method to remove CO gas to less than 10 ppm in reformed fuel gas. We study the mechanism of PROX of the CO reaction in the H₂ stream catalyzed by Cu _n Ni (n = 3-12) clusters using a density functional theory (DFT) calculation to investigate bimetallic effects on the catalytic activation. Our results indicate that the Cu₁₂Ni cluster is the most efficient catalyst for H₂ dissociation and the Cu₆Ni cluster is the most efficient catalyst for CO-PROX in excess hydrogen among Cu _n Ni (n = 3-12) clusters. To gain insight into the adsorption and dissociation of the H₂ molecule effect in the catalytic activity over the Cu₁₂Ni cluster and the potential energy surfaces about PROX of CO oxidation on the Cu₆Ni cluster, the nature of the interaction between the adsorbate and substrate is analyzed by detailed electron local densities of states (LDOS) as well as molecular structures.     

Study of water clusters in the n = 2–34 size regime, based on the ABEEM/MM model

Various properties of water clusters in the n = 2–34 size regime with the change of cluster size have been systemically explored based on the newly developed flexible-body and charge-fluctuating ABEEM/MM water potential model. The ABEEM/MM water model is to take ABEEM charges of all atoms, bonds, and lone-pairs of water molecules into the intermolecular electrostatic interaction term in molecular mechanics. The computed correlating properties characterizing water clusters (H₂O) _n (n = 2–34) include optimal structures, structural parameters, ABEEM charge distributions, binding energies, hydrogen bonds, dipole moments, and so on. The study of optimal structures shows that the ABEEM/MM model can correctly predict the following important structural features, such as the transition from two-dimensional (from dimer to pentamer) to three-dimensional (for clusters larger than the hexamer) structures at hexamer region, the transition from cubes to cages at dodecamer (H₂O)₁₂, the transition from all-surface (all water molecules on the surface of the cluster) to one water-centered (one water molecule at the center of the cluster, fully solvated) structures at (H₂O)₁₇, the transition from one to two internal molecules in the cage at (H₂O)₃₃, and so on. The first three structural transitions are in good agreement with those obtained from previous work, while the fourth transition is different from that identified by Hartke. Subsequently, a systematic investigation of structural parameters, ABEEM charges, energetic properties, and dipole moments of water clusters with increasing cluster size can provide important reference for describing the objective trait of hydrogen bonds in water cluster system, and also provide a strong impetus toward understanding how the water clusters approach the bulk limit.     

Ab initio study of the internal rotation barrier of formamide and the formamide–H2O complex

The barriers to internal rotation about the N? C bond of formamide and the formamide–H₂O complex have been examined by ab initio quantum chemistry methods. Both self-consistent field and correlated approaches have been utilized to determine the geometries, energies, and local harmonic vibrational frequencies of the minimum-energy and transition-state structures of these two systems. We find that formamide's rotation barrier of 14–15 kcal/mol is increased to 16–18 kcal/mol when a single H₂O molecule is attached. This result contrasts with the effect of a single H₂O solvent molecule on the barrier to tautomerization of formamide (to form formamidic acid HN? CHOH) for which a barrier lowering of ca. 20 kcal/mol has been observed. The rotation barriers obtained for formamide and for its H₂O complex are compared with barriers obtained experimentally in various solvents. © 1993 John Wiley & Sons, Inc.     

MNDO and INDO geometrical dependences of the hydrogen diamagnetic shielding tensor of some small molecules

The structural dependence of the hydrogen diamagnetic shielding tensor in H₂, CH₄, NH₄⁺, H₂O and HCl has been calculated using MNDO and INDO semiempirical methods. The influence of inner shells as well as two-center integrals are taken into account. Results are found to be in good agreement with ab initio predictions, and in the HCl molecule, with experimental data.     

Molecular Dynamics Simulations for Ice Modifications II and IX

Molecular vibrations of water (H₂O and D₂O) in crystals of ice II and ice IX are studied by molecular dynamics in a rigid bond approximation with a fixed bond angle. Using an atom-atomic potential PM for describing the interactions between water molecules in ice II (N = 576 molecules) and ice IX (N = 768) in an NVE ensemble leads to reproduction of the structure of both types of ice. For all water molecules and separately for each system of crystallographically equivalent water molecules in ice crystals, we defined the time dependence of the mean-square displacement of the center of mass of the molecule, the autocorrelation function of velocity for the center of mass, and the autocorrelation function of velocity for hydrogen (deuterium) atoms. The densities of vibrational states are calculated as Fourier integrals of the corresponding autocorrelation functions. In the case of ice II, the densities of states agree well with the experimental incoherent inelastic neutron scattering spectra. In the case of ice IX, agreement is worse. For both polymorphs, the mean-square displacement and the densities of vibrational states of the center of mass of the molecule and the hydrogen (deuterium) atom differ slightly between molecules belonging to different systems of crystallographic positions. This is explained by the difference in their environments.