Water clusters on graphite: methodology for quantum chemical a priori prediction of reaction rate constants

The properties, interactions, and reactions of cyclic water clusters (H(2)O)(n=1-5) on model systems for a graphite surface have been studied using pure B3LYP, dispersion-augmented density functional tight binding (DFTB-D), and integrated ONIOM(B3LYP:DFTB-D) methods. Coronene C(24)H(12) as well as polycircumcoronenes C(96)H(24) and C(216)H(36) in monolayer, bilayer, and trilayer arrangements were used as model systems to simulate ABA bulk graphite. Structures, binding energies, and vibrational frequencies of water clusters on mono- and bilayer graphite models have been calculated, and structural changes and frequency shifts due to the water cluster-graphite interactions are discussed. ONIOM(B3LYP:DFTB-D) with coronene and water in the high level and C(96)H(24) in the low level mimics the effect of extended graphite pi-conjugation on the water-graphite interaction very reasonably and suggests that water clusters only weakly interact with graphite surfaces, as suggested by the fact that water is an excellent graphite lubricant. We use the ONIOM(B3LYP:DFTB-D) method to predict rate constants for model pathways of water dissociative adsorption on graphite. Quantum chemical molecular dynamics (QM/MD) simulations of water clusters and water addition products on the C(96)H(24) graphite model are presented using the DFTB-D method. A three-stage strategy is devised for a priori investigations of high temperature corrosion processes of graphite surfaces due to interaction with water molecules and fragments.     

Water monomer interaction with gold nanoclusters from van der Waals density functional theory

We investigate the interaction between water molecules and gold nanoclusters Au(n) through a systematic density functional theory study within both the generalized gradient approximation and the nonlocal van der Waals (vdW) density functional theory. Both planar (n = 6-12) and three-dimensional (3D) clusters (n = 17-20) are studied. We find that applying vdW density functional theory leads to an increase in the Au-Au bond length and a decrease in the cohesive energy for all clusters studied. We classify water adsorption on nanoclusters according to the corner, edge, and surface adsorption geometries. In both corner and edge adsorptions, water molecule approaches the cluster through the O atom. For planar clusters, surface adsorption occurs in a O-up/H-down geometry with water plane oriented nearly perpendicular to the cluster. For 3D clusters, water instead favors a near-flat surface adsorption geometry with the water O atom sitting nearly atop a surface Au atom, in agreement with previous study on bulk surfaces. Including vdW interaction increases the adsorption energy for the weak surface adsorption but reduces the adsorption energy for the strong corner adsorption due to increased water-cluster bond length. By analyzing the adsorption induced charge rearrangement through Bader's charge partitioning and electron density difference and the orbital interaction through the projected density of states, we conclude that the bonding between water and gold nanocluster is determined by an interplay between electrostatic interaction and covalent interaction involving both the water lone-pair and in-plane orbitals and the gold 5d and 6s orbitals. Including vdW interaction does not change qualitatively the physical picture but does change quantitatively the adsorption structure due to the fluxionality of gold nanoclusters.     

Ionization of water clusters by collisions with graphite surfaces

We present experimental results on the scattering of neutral water clusters from graphite surfaces. We use cluster beams with an average cluster size up to 3700 molecules and an incident velocity of 1300 m/s, and study the emission of negatively and positively charged cluster fragments from the surface. The ionization probability is found to depend on cluster size and surface temperature, and for a given mean cluster size the emission rate of positive and negative cluster ions follows the Arrhenius equation. In the surface temperature range 950–1450 K, activation energies of 0.52±0.02 and 3.1±0.3 eV are determined for the emission of positive and negative ions, respectively. The emission of negative cluster fragments is attributed to electron transfer from the surface, and we estimate an electron affinity of 1.4±0.3 eV for large water clusters. Positive cluster fragments are proposed to be formed by dissocative ionization inside the cluster, followed by removal of the negative ion during surface contact.     

Analysis of localization sites for an excess electron in neutral methanol clusters using approximate pseudopotential quantum-mechanical calculations

We have used a recently developed electron-methanol molecule pseudopotential in approximate quantum mechanical calculations to evaluate and statistically analyze the physical properties of an excess electron in the field of equilibrated neutral methanol clusters ((CH(3)OH)(n), n=50-500). The methanol clusters were generated in classical molecular dynamics simulations at nominal 100 and 200 K temperatures. Topological analysis of the neutral clusters indicates that methyl groups cover the surface of the clusters almost exclusively, while the associated hydroxyl groups point inside. Since the initial neutral clusters are lacking polarity on the surface and compact inside, the excess electron can barely attach to these structures. Nevertheless, most of the investigated cluster configurations do support weakly stabilized cluster anion states. We find that similarly to water clusters, the pre-existing instantaneous dipole moment of the neutral clusters binds the electron. The localizing electrons occupy diffuse, weakly bound surface states that largely engulf the cluster although their centers are located outside the cluster molecular frame. The initial localization of the excess electron is reflected in its larger radius compared to water due to the lack of free OH hydrogens on the cluster surface. The stabilization of the excess electron increases, while the radius decreases monotonically as the clusters grow in size. Stable, interior bound states of the excess electron are not observed to form neither in finite size methanol clusters nor in the equilibrium bulk.     

Electron solvation in water-ammonia mixed clusters: Structure, energetics, and the nature of localization states of the excess electron

The structure and energetics of water-ammonia mixed clusters with an excess electron, [(H2O)n(NH3)m]- with m=1, n=2-6 and m=2, n=2, and also the corresponding neutral clusters are investigated in detail by means of ab initio quantum chemical calculations. The authors focus on the localization structure of the excess electron with respect to its surface versus interiorlike states, its binding to ammonia versus water molecules, the spatial and orientational arrangement of solvent molecules around the excess electron, the changes of the overall hydrogen-bonded structure of the clusters as compared to those of the neutral ones and associated dipole moment changes, vertical detachment energies of the anionic clusters, and also the vertical attachment energies of the neutral clusters. It is found that the hydrogen-bonded structure of the anionic clusters are very different from those of the neutral clusters unlike the case of water-ammonia dimer anion, and these changes in structural arrangements lead to drastically different dipole moments of the anionic and the neutral clusters. The spatial distribution of the singly occupied molecular orbital holding the excess electron shows only surface states for the smaller clusters. However, for n=5 and 6, both surface and interiorlike binding states are found to exist for the excess electron. For the surface states, the excess electron can be bound to the dangling hydrogens of either an ammonia or a water molecule with different degrees of stability and vertical detachment energies. The interiorlike states, wherever they exist, are found to have a higher vertical detachment energy than any of the surface states of the same cluster. Also, for interiorlike states, the ammonia molecule with its dangling hydrogens is always found to stay on top or on a far side of the charge density of the excess electron without participating in the hydrogen bond network of the cluster; the intermolecular hydrogen bonds are formed by the water molecules only which add to the overall stability of these anionic clusters.     

A scanning tunneling microscope study of single platinum atoms,platinum dimers and trimers on highly-oriented pyrolytic graphite

We report a topographic study of platinum clusters on highly-oriented pyrolytic graphite (HOPG) using a scanning tunneling microscope operating in air. The particles were produced by evaporation of platinum onto the graphite-surface in high vacuum. The simultaneous finding of single platinum atoms, clusters and small particles on an otherwise clean and atomically flat surface shows that the platinum-HOPG surface interaction is strong enough to yield stable images of Pt atoms and yet is not strong enough to annihilate the Pt-Pt interaction. Small flat platinum clusters on HOPG can be imaged with atomic resolution of the cluster and the surrounding graphite lattice. We show the adsorption site distribution for the monomers. The Pt-dimers show a very broad bond length distribution on graphite with an average of 2.46 Å. We found two types of Pt-trimers, one which is almost linear and one of triangular form. The average nearest neighbour distance of the trimers is 2.61 Å.     

Effect of compression of a stearic acid monolayer on interfacial binding of copper ions and cluster formation

The effect of compression of a stearic acid Langmuir monolayer on the interaction of copper ions with the monolayer and on the formation of interface clusters has been studied as functions of pH and ionic content of the water subphase. Copper binding was estimated by direct electron paramagnetic resonance measurements of ion concentrations in the water subphase. A scanning tunnelling microscopy study of monolayer Langmuir-Blodgett films, deposited on graphite substrates, gave evidence for the formation of nanosized clusters on the monolayer surface. The data obtained showed that the interaction of the monolayer with copper ions and the accompaning cluster formation processes were determined by the arrangement, order and mobility of the stearic acid molecules in the monolayer and the electrostatics at the interface.     

The molecular approach to heterogeneous nucleation

A molecular approach to heterogeneous nucleation has been developed. The expressions for the equilibrium cluster distribution, the reversible work of the cluster formation, and the nucleation rate have been derived. Two separate statements for the work of formation were formulated. If the equilibrium cluster distribution is normalized on the monomer concentration near the substrate surface, the reversible work of formation is expressed by DeltaG(het) (I) = (F(n) (het)-F(n) (hom))-(F(1) (het)-F(1) (hom)) + DeltaG(hom) where F(n) (het) and F(n) (hom) are the Helmholtz free energies of a cluster interacting with a substrate and a cluster not interacting with the substrate, respectively. If the equilibrium cluster distribution is normalized on the monomer concentration far from the substrate surface, the work of cluster formation is given by DeltaG(het) (II) = (F(n) (het)-F(n) (hom)) + DeltaG(hom). The former expression corresponds to the approach of the classical heterogeneous nucleation theory. The cluster partition function appears to be dependent on the location of a virtual plane, which separates the volume, where the interaction of the clusters with the substrate is effective from the one where interaction is negligible. Our Monte Carlo simulations have shown that the dependence is rather weak and thus the location of the plane is not very important. According to the simulations the variation of the plane position in the range from 20 to 50 Angstroms does not lead to a considerable change of the heterogeneous nucleation rate.     

Growth pattern of Ag(n) (n = 1-8) clusters on the α-Al2O3(0001) surface: a first principles study

We report an extensive first-principles study of the structure and electronic properties of Ag(n) (n = 1-8) clusters isolated in gas phase and deposited on the α-Al(2)O(3) surface. We have used the plane wave based pseudopotential method within the framework of density functional theory. The electron ion interaction has been described using projector augmented wave (PAW), and the spin-polarized GGA scheme was used for the exchange correlation energy. The results reveal that, albeit interacting with support alumina, the Ag atoms prefers to remain bonded together suggesting an island growth motif is preferred over wetting the surface. When compared the equilibrium structures of Ag clusters between free and on alumina substrate, a significant difference was observed starting from n = 7 onward. While Ag(7) forms a three-dimensional (3D) pentagonal bipyramid in the isolated gas phase, on alumina support it forms a planar hexagonal structure parallel to the surface plane. Moreover, the spin moment of the Ag(7) cluster was found to be fully quenched. This has been attributed to higher delocalization of electron density as the size of the cluster increases. Furthermore, a comparison of chemical bonding analysis through electronic density of state (EDOS) shows that the EDOS of the deposited Ag(n) cluster is significantly broader, which has been ascribed to the enhanced spd hybridization. On the basis of the energetics, it is found that the adsorption energy of Ag clusters on the α-Al(2)O(3) surface decreases with cluster size.     

Pinning of size-selected Pd nanoclusters on graphite

The production of stable cluster arrays on smooth surfaces has several potential technological applications. We report a study of the pinning of size-selected palladium nanoclusters on the graphite surface. The clusters formed during gas aggregation in vacuum are projected with sufficient kinetic energy to create a defect in the graphite surface. The energy necessary to create such an immobilizing defect is investigated as a function of the palladium cluster size. The palladium pinning energy is found to deviate from the simple binary collision model as appropriate to previously reported silver and gold results. This finding is in agreement with the deviation of nickel clusters and points to the influence of the interatomic cluster bonding on the mechanics of the collision.     

Excess electron localization sites in neutral water clusters

We present approximate pseudopotential quantum-mechanical calculations of the excess electron states of equilibrated neutral water clusters sampled by classical molecular dynamics simulations. The internal energy of the clusters are representative of those present at temperatures of 200 and 300 K. Correlated electronic structure calculations are used to validate the pseudopotential for this purpose. We find that the neutral clusters support localized, bound excess electron ground states in about 50% of the configurations for the smallest cluster size studied (n = 20), and in almost all configurations for larger clusters (n > 66). The state is always exterior to the molecular frame, forming typically a diffuse surface state. Both cluster size and temperature dependence of energetic and structural properties of the clusters and the electron distribution are explored. We show that the stabilization of the electron is strongly correlated with the preexisting instantaneous dipole moment of the neutral clusters, and its ground state energy is reflected in the electronic radius. The findings are consistent with electron attachment via an initial surface state. The hypothetical spectral dynamics following such attachment is also discussed.     

Undissociated versus dissociated structures for water clusters and ammonia-water clusters: (H2O)n and NH3(H2O)n-1 (n = 5, 8, 9, 21). Theoretical study

To understand the autoionization of pure water and the solvation of ammonia in water, we investigated the undissociated and dissociated (ion-pair) structures of (H2O) n and NH3(H2O)n-1 (n = 5, 8, 9, 21) using density functional theory (DFT) and second order Moller-Plesset perturbation theory (MP2). The stability, thermodynamic properties, and infrared spectra were also studied. The dissociated (ion-pair) form of the clusters tends to favor the solvent-separated ion-pair of H3O+/NH4+ and OH-. As for the NH3(H2O)20 cluster, the undissociated structure has the internal conformation, in contrast to the surface conformation for the (H2O)21 cluster, whereas the dissociated structure of NH3(H2O)20 has the surface conformation. As the cluster size of (H2O)n/NH3(H2O)n-1 increases, the difference in standard free energy between undissociated and dissociated (ion-pair) clusters is asymptotically well corroborated with the experimental free energy change at infinite dilution of H3O+/NH4+ and OH-. The predicted NH and OH stretching frequencies of the undissociated and dissociated (ion-pair) clusters are discussed.     

Effect of formic acid addition on water cluster stability and structure

Computational chemistry simulations were performed to determine the effect that the addition of a single formic acid molecule has on the structure and stability of protonated water clusters. Previous experimental studies showed that addition of formic acid to protonated pure water results in higher intensities of large-sized clusters when compared to pure water and methanol-water mixed clusters. For larger, protonated clusters, molecular dynamics simulations were performed on H(+)(H(2)O)(n), H(+)(H(2)O)(n)CH(3)OH, and H(+)(H(2)O)(n)CHOOH clusters, 19-28 molecules in size, using a reactive force field (ReaxFF). Based on these computations, formic acid-water clusters were found to have significantly higher binding energies per molecule. Addition of formic acid to a water cluster was found to alter the structure of the hydrogen-bonding network, creating selective sites within the cluster, enabling the formation of new hydrogen bonds, and increasing both the stability of the cluster and its rate of growth.     

Mid-infrared characterization of the NH4 +(H2O)n clusters in the neighborhood of the n=20 "magic" number

Vibrational predissociation spectra are reported for size-selected NH4+ (H2O)n clusters (n=5-22) in the 2500-3900 cm(-1) region. We concentrate on the sharp free OH stretching bands to deduce the local H-bonding configurations of water molecules on the cluster surface. As in the spectra of the protonated water clusters, the free OH bands in NH4+ (H2O)n evolve from a quartet at small sizes (n<7), to a doublet around n=9, and then to a single peak at the n=20 magic number cluster, before the doublet re-emerges at larger sizes. This spectral simplification at the magic number cluster mirrors that found earlier in the H+(H2O)n clusters. We characterize the likely structures at play for the n=19 and 20 clusters with electronic structure calculations. The most stable form of the n=20 cluster is predicted to have a surface-solvated NH4+ ion that lies considerably lower in energy than isomers with the NH4+ in the interior.     

DFT studies of the interactions of a graphene layer with small water aggregates

We have investigated the structure, adsorption, electronic states, and charge transfer of small water aggregates on the surface of a graphene layer using density functional theory. Our calculations were focused on water adsorbates containing up to five water molecules interacting with one and both sides of a perfect freestanding sheet. Different orientations of the aggregates with respect to the graphene sites were considered. The results show that the adsorption energy of one water molecule is primarily determined by its orientation, although it is also strongly dependent on the implemented functional scheme. Despite its intrinsic difficulties with dispersion interactions, the Perdew and Wang's exchange-correlation functional may be a viable alternative to investigate the adsorption of large molecular aggregates on a graphene surface. Although water physisorption is expected to occur in the regime of droplets, we found no induced impurity states close to the Fermi level of graphene interacting with small water clusters. In order to investigate the donor/acceptor tendency of the water clusters on graphene, we have performed a Bader charge analysis. Considering the charge transfer mechanism, we have noticed that it should preferentially occur from water to graphene only when the oxygen atom is pointing toward the surface. Otherwise, and in the case of larger adsorbed clusters, charge transfers systematically occur from graphene to water.     

水在石墨(0001)面簇模型桥位上吸附的量子化学研究

用从头计算方法对水在石墨(0001)面桥位上的吸附进行了研究.用C₆H₈原子簇模拟石墨表面,在6-31G^*水平上计算了水在不同方向和位置上的吸附能量.研究表明:水在石墨面上的吸附很弱,属于物理吸附;在中性或带负电荷的石墨表面,当水分子中的氢原子靠近石墨面时,体系存在能量最小值,而在带正电荷的表面,当氧原子靠近石墨面时存在稳定的吸附点;不论表面带正电荷还是带负电荷,均对水分子的吸附起增强作用.     

Microsolvation of the sodium and iodide ions and their ion pair in acetonitrile clusters: a theoretical study

The structural and thermodynamic properties of Na+(CH3CN)n, I-(CH3CN)n, and NaI(CH3CN)n clusters have been investigated by means of room-temperature Monte Carlo simulations with model potentials developed to reproduce the properties of small clusters predicted by quantum chemistry. Ions are found to adopt an interior solvation shell structure, with a first solvation shell containing approximately 6 and approximately 8 acetonitrile molecules for large Na+(CH3CN)n and I-(CH3CN)n clusters, respectively. Structural features of Na+(CH3CN)n are found to be similar to those of Na+(H2O)n clusters, but those of I-(CH3CN)n contrast with those of I-(H2O)n, for which "surface" solvation structures were observed. The potential of mean force calculations demonstrates that the NaI ion pair is thermodynamically stable with respect to ground-state ionic dissociation in acetonitrile clusters. The properties of NaI(CH3CN)n clusters exhibit some similarities with NaI(H2O)n clusters, with the existence of contact ion pair and solvent-separated ion pair structures, but, in contrast to water clusters, both types of ion pairs adopt a well-defined interior ionic solvation shell structure in acetonitrile clusters. Whereas contact ion pair species are thermodynamically favored in small clusters, solvent-separated ion pairs tend to become thermodynamically more stable above a cluster size of approximately 26. Hence, ground-state charge separation appears to occur at larger cluster sizes for acetonitrile clusters than for water clusters. We propose that the lack of a large Na+(CH3CN)n product signal in NaI(CH3CN)n multiphoton ionization experiments could arise from extensive stabilization of the ground ionic state by the solvent and possible inhibition of the photoexcitation mechanism, which may be less pronounced for NaI(H2O)n clusters because of surface solvation structures. Alternatively, increased solvent evaporation resulting from larger excess energies upon photoexcitation or major solvent reorganization on the ionized state could account for the observed solvent-selectivity in NaI cluster multiphoton ionization.     

Loss of chlorine in mass spectra of DCl picked up by water clusters in a beam

Water clusters (H(2)O)(n), n相似文献   

Ab initio study of structural stability of Mo-S clusters and size specific stoichiometries of magic clusters

Using first-principles calculations with ultrasoft pseudopotential formalism and the generalized gradient approximation for the exchange-correlation functional, we study the stability of MonSm (n =1-6 and m ranging from n to 3n) clusters and obtain the optimal stoichiometry for each n corresponding to the magic cluster. It is found that in this size range, the lowest-energy structures favor a core of metal atoms, which is covered by sulfur. In particular, we observe that for Mo6S14 isolated clusters, a 3D structure is significantly lower in energy as compared to platelet structures found recently on Au (111) surface. The composition ratio between S and Mo in the magic clusters is less than 2 for n=3 and greater than 2 for n<3. The structural stability of the magic clusters arises from the optimization of the Mo-Mo and Mo- S bonding as well as the symmetry of the cluster. Addition of a terminal sulfur in a magic cluster generally lowers its binding energy. The presence of partially occupied d-orbitals in Mo atoms contributes to Mo-Mo bonding and for higher S concentration it leads to sulfur-sulfur bond formation. The variation in energy due to a change in the sulfur composition suggests that sulfurization of the magic clusters is generally more favorable than desulfurization.