Nascent Decomposition Pathways of CH4 Pyrolysis in Gas-Phase Metal Halides.

We have performed a combined quantum mechanical and microkinetic modeling study to understand the nascent decomposition pathways of methane pyrolysis, catalyzed by gas-phase ZnCl2, in a constant pressure batch reactor at 1273 K. We find that ZnCl2 catalyzes methane pyrolysis with an apparent activation energy of 227 kJ/mol. We have also performed sensitivity analysis on a reaction network comprising initiation, termination, and primary propagation reactions. The results suggest that the whole reaction network can be simplified to four reactions, which contributes to the initial rate of methane decomposition. Based on these insights, we have also explored the catalyzing effects of gas-phase AlCl3, CoCl2, CuCl2, FeCl2, and NiCl2 for methane decomposition. Our calculations suggest that gas-phase CuCl2 and NiCl2 are the most active catalysts among the metal halides studied in this work.

Rates of adsorption and desorption: Entropic contributions and errors due to mean-field approximations.

We have performed exact classical rate calculations to compute adsorption and desorption rate constants with a model representative of a real system. We compute the desorption rate using transition-state theory by taking the dividing-surface far from the surface of the solid. We find that using a mean-field assumption, i.e., applying potential of mean force to transition state theory, could lead to two orders-of-magnitude errors in the rate constant owing to large fluctuations in the desorption barrier. Furthermore, we compute the adsorption rate by including a dynamical factor which reflects the probability of sticking to the solid surface. We find that the sticking probability is highly sensitive to the coverage. Also, we find that the adsorption rate computed from the mean-field assumption is not very different from the exact adsorption rate. We also compute entropic contribution to desorption rates and compare it to that obtained from two limiting models of adsorption-2D ideal gas and 2D ideal lattice gas. We show that at high temperatures (700 K), the entropic contribution to desorption rates computed from the exact calculations is very close to that obtained from the 2D ideal gas model. However, for lower to intermediate temperatures from 200 K to 500 K, the entropic contributions cover a wide range which lies in between the two limiting models and could lead to over two-orders-of-magnitude errors in the rate coefficient.

Investigation of the active sites of rhodium sulfide for hydrogen evolution/oxidation using carbon monoxide as a probe.

Carbon monoxide (CO) was observed to decrease the activity for hydrogen evolution, hydrogen oxidation, and H2-D2 exchange on rhodium sulfide, platinum, and rhodium metal. The temperature at which the CO was desorbed from the catalyst surface (detected by recovery in the H2-D2 exchange activity of the catalyst) was used as a descriptor for the CO binding energy to the active site. The differences in the CO desorption temperature between the different catalysts showed that the rhodium sulfide active site is not metallic rhodium. Using density functional theory, the binding energy of CO to the Rh sites in rhodium sulfide is found comparable to the binding energy on Pt. Coupled with experiment this supports the proposition that rhodium rather than sulfur atoms in the rhodium sulfide are the active site for the hydrogen reaction. This would indicate the active sites for hydrogen evolution/oxidation as well as oxygen reduction (determined by other groups using X-ray absorption spectroscopy) may be the same.

Catalytic oxidation of methanol to formaldehyde by mass-selected vanadium oxide clusters supported on a TiO2(110) surface.

We report the results of a systematic study of the catalytic activity of mass-selected vanadium oxide clusters deposited on rutile TiO2 surfaces under ultrahigh vacuum (UHV) conditions. Our results show that supported V, VO, and VO2 clusters are not catalytically active for the oxidative dehydrogenation of methanol to formaldehyde but can be made catalytically active by postoxidation. In addition, we found that the postoxidized VO/TiO2 produces the most formaldehyde. Scanning tunneling microscopy (STM) imaging of the postoxidized VO/TiO2 reveals isolated clusters with height and width indicative of VO3 bound to the TiO2 surface. Our results are consistent with previous density functional theory (DFT) calculations that predict that VO3 will be produced by postoxidation of VO and that VO3/TiO2 is an active catalyst.

Transition metal sulfide hydrogen evolution catalysts for hydrobromic acid electrolysis.

Mixed metal sulfides containing combinations of W, Fe, Mo, Ni, and Ru were synthesized and screened for activity and stability for the hydrogen evolution reaction (HER) in aqueous hydrobromic acid (HBr). Co- and Ni-substituted RuS(2) were identified as potentially active HER electrocatalysts by high-throughput screening (HTS), and the specific compositions Co(0.4)Ru(0.6)S(2) and Ni(0.6)Ru(0.4)S(2) were identified by optimization. Hydrogen evolution activity of Co(0.4)Ru(0.6)S(2) in HBr is greater than RuS(2) or CoS(2) and comparable to Pt and commercial Rh(x)S(y). Structural and morphological characterizations of the Co-substituted RuS(2) suggest that the nanoparticulate solids are a homogeneous solid solution with a pyrite crystal structure. No phase separation is detected for Co substitutions below 30% by X-ray diffraction. In 0.5 M HBr electrolyte, the Co-Ru electrode material synthesized with 30% Co rapidly lost approximately 34% of the initial loading of Co; thereafter, it was observed to exhibit stable activity for HER with no further loss of Co. Density functional theory calculations indicate that the S(2)(2-) sites are the most important for HER and the presence of Co influences the S(2)(2-) sites such that the hydrogen binding energy at sufficiently high hydrogen coverage is decreased compared to ruthenium sulfide. Although showing high HER activity in a flow cell, the reverse reaction of hydrogen oxidation is slow on the RuS(2) catalysts tested when compared to platinum and rhodium sulfide, leaving rhodium sulfide as the only suitable tested material for a regenerative HBr cell due its stability compared to platinum.

VO Cluster-Stabilized H2O Adsorption on a TiO2 (110) Surface at Room Temperature.

We probe the adsorption of molecular H2O on a TiO2 (110)-(1 × 1) surface decorated with isolated VO clusters using ultrahigh-vacuum scanning tunneling microscopy (UHV-STM) and temperature-programmed desorption (TPD). Our STM images show that preadsorbed VO clusters on the TiO2 (110)-(1 × 1) surface induce the adsorption of H2O molecules at room temperature (RT). The adsorbed H2O molecules form strings of beads of H2O dimers bound to the 5-fold coordinated Ti atom (5c-Ti) rows and are anchored by VO. This RT adsorption is completely reversible and is unique to the VO-decorated TiO2 surface. TPD spectra reveal two new desorption states for VO stabilized H2O at 395 and 445 K, which is in sharp contrast to the desorption of water due to recombination of hydroxyl groups at 490 K from clean TiO2(110)-(1 × 1) surfaces. Density functional theory (DFT) calculations show that the binding energy of molecular H2O to the VO clusters on the TiO2 (110)-(1 × 1) surface is higher than binding to the bare surface by 0.42 eV, and the resulting H2O-VO-TiO2 (110) complex provides the anchor point for adsorption of the string of beads of H2O dimers.

The properties of small Ag clusters bound to DNA bases.

We study the binding of neutral silver clusters, Ag(n) (n=1-6), to the DNA bases adenine (A), cytosine (C), guanine (G), and thymine (T) and the absorption spectra of the silver cluster-base complexes. Using density functional theory (DFT), we find that the clusters prefer to bind to the doubly bonded ring nitrogens and that binding to T is generally much weaker than to C, G, and A. Ag(3) and Ag(4) make the stronger bonds. Bader charge analysis indicates a mild electron transfer from the base to the clusters for all bases, except T. The donor bases (C, G, and A) bind to the sites on the cluster where the lowest unoccupied molecular orbital has a pronounced protrusion. The site where cluster binds to the base is controlled by the shape of the higher occupied states of the base. Time-dependent DFT calculations show that different base-cluster isomers may have very different absorption spectra. In particular, we find new excitations in base-cluster molecules, at energies well below those of the isolated components, and with strengths that depend strongly on the orientations of planar clusters with respect to the base planes. Our results suggest that geometric constraints on binding, imposed by designed DNA structures, may be a feasible route to engineering the selection of specific cluster-base assemblies.

Inelastic scattering with Chebyshev polynomials and preconditioned conjugate gradient minimization.

We describe and test an implementation, using a basis set of Chebyshev polynomials, of a variational method for solving scattering problems in quantum mechanics. This minimum error method (MEM) determines the wave function Psi by minimizing the least-squares error in the function (H Psi - E Psi), where E is the desired scattering energy. We compare the MEM to an alternative, the Kohn variational principle (KVP), by solving the Secrest-Johnson model of two-dimensional inelastic scattering, which has been studied previously using the KVP and for which other numerical solutions are available. We use a conjugate gradient (CG) method to minimize the error, and by preconditioning the CG search, we are able to greatly reduce the number of iterations necessary; the method is thus faster and more stable than a matrix inversion, as is required in the KVP. Also, we avoid errors due to scattering off of the boundaries, which presents substantial problems for other methods, by matching the wave function in the interaction region to the correct asymptotic states at the specified energy; the use of Chebyshev polynomials allows this boundary condition to be implemented accurately. The use of Chebyshev polynomials allows for a rapid and accurate evaluation of the kinetic energy. This basis set is as efficient as plane waves but does not impose an artificial periodicity on the system. There are problems in surface science and molecular electronics which cannot be solved if periodicity is imposed, and the Chebyshev basis set is a good alternative in such situations.

O2 evolution on a clean partially reduced rutile TiO2(110) surface and on the same surface precovered with Au1 and Au2: the importance of spin conservation.

We have used spin-polarized density functional theory (DFT) to study O(2) evolution on a clean partially reduced rutile TiO(2)(110) surface (i.e., a surface having oxygen vacancies) and its interaction with Au(1) or Au(2) cluster adsorbed on it. We assume that the total spin of the electronic wave function is related to the number of unpaired spins (N(s)) and calculate the binding and the activation energies involved in O(2) evolution for fixed values of N(s). In addition to keeping N(s) constant, we assume that reactions in which the N(s) of the reactants differs from that of the products are very slow. The potential energy surfaces obtained in this way depend strongly on N(s). For example, O(2) dissociation at the vacancy site on a clean partially reduced TiO(2)(110) surface is exothermic by 0.85 eV in the triplet state and the highest activation energy in the chain of reactions leading to the O(2) dissociation is 0.67 eV. In the singlet state, O(2) dissociation is endothermic by 0.11 eV and the activation energy leading to dissociation is 1.30 eV. These observations are in qualitative agreement with scanning tunneling microscopy experiment in which O(2) dissociation on a partially reduced rutile TiO(2)(110) surface is observed at temperature as low as 120 K. In contrast, O(2) dissociation is predicted to be endothermic and is prevented by an activation barrier larger than 1 eV in all the previous DFT calculations, in which the DFT program varies N(s) to get the lowest energy state. We find that on a partially reduced rutile TiO(2)(110) with Au(1) and Au(2) preadsorbed on its surface, O(2) dissociates at the vacancy site: One oxygen atom fills the oxygen vacancy and the other becomes available for oxidation chemistry. This means that Au(1) and Au(2) supported on a partially reduced TiO(2)(110) surface is not an oxidation catalyst since the presence of oxygen turns it into a stoichiometric Au(n)/TiO(2)(110) surface. Finally, we find that the evolution of oxygen on Au(1) and Au(2) in the gas phase is very different from the evolution on the same clusters supported on the partially reduced TiO(2)(110) surface. For example, the molecular adsorption of O(2) is favored in the gas phase (except on Au(1) (-) and Au(2) (-) in the quartet state), while the dissociative adsorption is favored by more than 1 eV when Au(1) and Au(2) are supported on the partially reduced TiO(2)(110). Furthermore, the activation energies associated with O(2) dissociation in the gas phase (DeltaE(act)>2.4 eV) are reduced by at least a factor of 2 when the clusters are supported on TiO(2)(110).

Dynamics of H2O and Na+ in nafion membranes.

We investigate the transport properties of a model of a hydrated Na-Nafion membrane using molecular dynamics simulations. The system consists of several Nafion chains forming a pore with the water and ions inside. At low water content, the hydrophilic domain is not continuous and diffusion is very slow. The diffusion coefficient of both water and Na+ increases with increasing hydration (more strongly so for Na+). The simulations are in qualitative agreement with experimental results for similar systems. The diffusion coefficient is an average over the motion of ions or water molecules located in different environments. To better understand the role of the environment, we calculate the distribution of the residence times of the ion (or water) at different locations in the system. We discuss the transport mechanism in light of this information.

Catalytic molten metals for the direct conversion of methane to hydrogen and separable carbon.

Metals that are active catalysts for methane (Ni, Pt, Pd), when dissolved in inactive low-melting temperature metals (In, Ga, Sn, Pb), produce stable molten metal alloy catalysts for pyrolysis of methane into hydrogen and carbon. All solid catalysts previously used for this reaction have been deactivated by carbon deposition. In the molten alloy system, the insoluble carbon floats to the surface where it can be skimmed off. A 27% Ni-73% Bi alloy achieved 95% methane conversion at 1065°C in a 1.1-meter bubble column and produced pure hydrogen without CO2 or other by-products. Calculations show that the active metals in the molten alloys are atomically dispersed and negatively charged. There is a correlation between the amount of charge on the atoms and their catalytic activity.

Efficient electrocatalyst utilization: electrochemical deposition of Pt nanoparticles using nafion membrane as a template.

We deposit Pt particles electrochemically on an electrode covered with a Nafion membrane. Platinum ions travel through the hydrophilic channels of the membrane, and platinum deposits are formed at the place where the channels make contact with the planar electrode. This procedure deposits the catalyst only at the end of the hydrophilic channels that cross the membrane; no catalyst is placed under the hydrophobic domains, where it would not be in contact with the electrolyte. By performing a series of cyclic voltammograms with this system, we show that deposition of the platinum through the membrane achieves better platinum utilization than deposition of platinum on the naked electrode followed by the placement of the membrane on top.

Pinning mononuclear Au on the surface of titania.

We have deposited Au atoms on the surface of titania without sintering or surface damage. Mass-selected Au+ atoms were deposited from the gas phase at room temperature with kinetic energies from <3 to 190+/-3.5 eV. Scanning tunneling microscopy reveals island formation following deposition at <3 eV, while mainly atomic features are observed for energies between approximately 35 and approximately 190 eV. A mixture of islands and atomic features is observed at a landing energy of 20+/-3.5 eV, suggesting a critical energy above which pinning occurs. Cluster size is also probed as a function of coverage in the deposition of Au+ with 100 eV of energy, revealing that sintering begins at a coverage of only 0.06 ML. These observations suggest a mechanism in which high-energy collision leads to the annealing of any impact-created surface damage and the pinning of Au atoms to the surface. We provide a new method of preparing isolated Au atoms on an oxide surface, which can serve as a platform for catalytic studies.

Electrolithographic investigations of the hydrophilic channels in Nafion membranes.

Nafion membranes are used as semisolid electrolytes in methanol and hydrogen fuel cells. The ion conduction takes place through those hydrophilic channels in the Nafion that can provide continuous pathways through the membrane. There is as yet limited information about the density, the size, and the shape of these channels. We have developed two electrochemical methods of visualizing the pore structure which involve the creation of metal lithographs using the membrane pores as templates. In the experiments, the membrane is supported on a flat solid surface on one side, and is in contact with an electrolyte on the other side. Using hydrogen-terminated n-doped Si(111), we deposited gold from an electrolyte containing a gold salt. The Au ions traverse the membrane through the pores, reach the silicon surface, and are spontaneously reduced. A metallic Au deposit is formed on the silicon surface, at the base of the hydrophilic channel. The Au deposits are imaged after the membrane is dissolved. Another method involves supporting the membrane on a Pt surface and depositing silver wires through the hydrophilic channels of the membrane. The scanning electron microscope pictures of these wires provide an image of the size and the shape of the hydrophilic channels.

Structure of hydrated Na-Nafion polymer membranes.

We use molecular dynamics simulations to investigate the structure of the hydrated Na-Nafion membranes. The membrane is "prepared" by starting with the Nafion chains placed on a cylinder having the water inside it. Minimizing the energy of the system leads to a filamentary hydrophilic domain whose structure depends on the degree of hydration. At 5 wt % water the system does not have enough water molecules to solvate all the ions that could be formed by the dissociation of the -SO3Na groups. As a result, the -SO3Na groups aggregate with the water to form very small droplets that do not join into a continuous phase. The size of the droplets is between 5 and 8 A. As the amount of water present in the membrane is increased, the membrane swells, and SO3Na has an increasing tendency to dissociate into ions. Furthermore, a transition to a percolating hydrophilic network is observed. In the percolating structure, the water forms irregular curvilinear channels branching in all directions. The typical dimension of the cross section of these channels is about 10-20 A. Calculated neutron scattering from the simulated system is in qualitative agreement with experiment. In all simulations, the pendant sulfonated perfluorovinyl side chains of the Nafion hug the walls of the hydrophilic channel, while the sulfonate groups point toward the center of the hydrophilic phase. The expulsion of the side chains from the hydrophilic domain is favored because it allows better interaction between the water molecules. We have also examined the probability of finding water molecules around the Na+ and the -SO3(-) ions as well as the probability of finding other water molecules next to a given water molecule. These probabilities are much broader than those found in bulk water or for one ion in bulk water (calculated with the potentials used in the present simulation). This is due to the highly inhomogeneous nature of the material contained in the small hydrophilic pores.

Excess proton solvation and delocalization in a hydrophilic pocket of the proton conducting polymer membrane nafion.

Solvation properties of the hydrated excess proton are studied in a hydrophilic pocket of Nafion 117 through a series of molecular dynamics simulations. The multistate empirical valence bond (MS-EVB) methodology, which enables the delocalization of the excess proton through the Grotthuss hopping mechanism, was employed for one of the excess protons in the simulation cell. Simulations were performed such that "classical" nondissociable hydronium cations and a single excess proton treated with the MS-EVB methodology were at a concentration ratio of 39:1. Two degrees of hydration of the Nafion polymer electrolyte membrane were simulated, each displaying the same marked difference between the solvation structures of the classical versus MS-EVB treated (Grotthuss shuttling) excess proton species. These differences are attributed to the solvent dynamics needed to transfer the cation between the solvent separated and contact pair positions about the sulfonic acid counterion. The results demonstrate that it is generally impossible to describe the low pH conditions in the hydrophilic domains of Nafion without the explicit treatment of Grotthuss delocalization in the underlying molecular dynamics model for the excess protons.

Catechol and HCl Adsorption on TiO2(110) in Vacuum and at the Water-TiO2 Interface.

Coadsorbed water is often unavoidable in electrochemistry and low-temperature catalysis. In addition, water influences the adsorption of biomolecules on surfaces. We use ab initio DFT molecular dynamics and ground-state calculations to study the adsorption of HCl and catechol on the rutile TiO2(110) surface and at a water-rutile interface. We find that a coadsorbed water film reduces the adsorption energy of both catechol and HCl significantly because water molecules must be displaced from the surface before catechol or HCl can adsorb. The adsorption energy of catechol (or HCl) at the water-rutile interface can be estimated as the adsorption energy in vacuum minus the energy to remove two water molecules (respectively, one water molecule) from the rutile surface in vacuum and place them in liquid water. This estimate predicts the effect of a surface water film on adsorption without the need of molecular dynamics.

Enhanced adsorption energy of Au1 and O2 on the stoichiometric TiO2(110) surface by coadsorption with other molecules.

During heterogeneous catalysis the surface is simultaneously covered by several adsorbed molecules. The manner in which the presence of one kind of molecule affects the adsorption of a molecule of another kind has been of interest for a long time. In most cases the presence of one adsorbate does not change substantially the binding energy of another adsorbate. The calculations presented here show that the stoichiometric rutile TiO(2)(110) surface, on which one of the compounds -OH, Au(3), Au(5), Au(7), Na, K, or Cs or two different gold strips was preadsorbed, behaves differently: the binding energy of Au(1) or O(2) to such a surface is much stronger than the binding to the clean stoichiometric TiO(2)(110) surface. Moreover, the binding energy of Au(1) or O(2) and the amount of charge they take from the surface when they adsorb are the same, regardless of which of the above species is preadsorbed. The preadsorbed species donate electrons to the conduction band of the oxide, and these electrons are used by Au(1) or O(2) to make stronger bonds with the surface. This suggests that adding an electron to the conduction band of the clean stoichiometric TiO(2)(110) slab used in the calculation will affect similarly the adsorption energy of Au(1) or O(2). Our calculations show that it does. We have also studied how the preadsorption of Au(4) or Au(6) affects the binding of Au(1) or O(2) to the surface. These two gold clusters do not donate electrons to the surface when they bind to it and therefore should not influence substantially the binding energy of Au(1) or O(2) to the surface. However, adsorbing O(2) or Au(1) on the surface forces the clusters to change their structure into that of isomers that donate charge to the oxide. This charge is used by Au(1) or O(2) to bind to the surface and the energy of this bond exceeds the isomerization energy. As a result the surface with the isomerized cluster is the lowest energy state of the system. We believe that these results can be generalized as follows. The molecules that we coadsorbed with Au(1) or O(2) donate electrons to the oxide and are Lewis bases. By giving the surface high energy electrons, they turn it into a Lewis base and this increases its ability to bind strong Lewis acids such as Au(1) and O(2). We speculate that this kind of interaction is general and may be observed for other oxides and for other coadsorbed Lewis base-Lewis acid pairs.

Density functional study of the interaction between small Au clusters, Au(n) (n=1-7) and the rutile TiO(2) surface. I. Adsorption on the stoichiometric surface.

This is the first paper in a series of four dealing with the adsorption site, electronic structure, and chemistry of small Au clusters, Au(n) (n=1-7), supported on stoichiometric, partially reduced, or partially hydroxylated rutile TiO(2)(110) surfaces. Analysis of the electronic structure reveals that the main contribution to the binding energy is the overlap between the highest occupied molecular orbitals of Au clusters and the Kohn-Sham orbitals localized on the bridging and the in-plane oxygen of the rutile TiO(2)(110) surface. The structure of adsorbed Au(n) differs from that in the gas phase mostly because the cluster wants to maximize this orbital overlap and to increase the number of Au-O bonds. For example, the equilibrium structures of Au(5) and Au(7) are planar in the gas phase, while the adsorbed Au(5) has a distorted two-dimensional structure and the adsorbed Au(7) is three-dimensional. The dissociation of an adsorbed cluster into two adsorbed fragments is endothermic, for all clusters, by at least 0.8 eV. This does not mean that the gas-phase clusters hitting the surface with kinetic energy greater than 0.8 eV will fragment. To place enough energy in the reaction coordinate for fragmentation, the impact kinetic energy needs to be substantially higher than 0.8 eV. We have also calculated the interaction energy between all pairs of Au clusters. These interactions are small except when a Au monomer is coadsorbed with a Au(n) with odd n. In this case the interaction energy is of the order of 0.7 eV and the two clusters interact through the support even when they are fairly far apart. This happens because the adsorption of a Au(n) cluster places electrons in the states of the bottom of the conduction band and these electrons help the Au monomer to bind to the five-coordinated Ti atoms on the surface.