Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Density functional/Monte Carlo study of ring-opening polymerization

Density functional calculations of the structure, potential energy surface and reactivity for organic systems closely related to bisphenol-A-polycarbonate (BPA-PC) provide the basis for a model describing the ring-opening polymerization of its cyclic oligomers by nucleophilic molecules. Monte Carlo simulations using this model show a strong tendency to polymerize that is increased by increasing density and temperature, and is greater in 3D than in 2D. Entropy in the distribution of inter-particle bonds is the driving force for chain formation. (C) 2002 Elsevier Science B.V. All rights reserved.

Reactions of polycarbonate with cyclohexene oxide and phosphites: A density functional study

Epoxides and phosphites are often used as additives to stabilize the properties of polymers, including bisphenol A polycarbonate (BPA-PC). We describe density functional (DF) calculations of the reactions of cyclohexene oxide (CHO, cyclohexane epoxide) and phosphites with chain segments of BPA-PC, with the aim of identifying possible reaction paths and energy barriers. The reactions of CHO with the OH-terminated PC chains and with the carbonate group are exothermic, although there is an energy barrier in each case of more than 10 kcal/mol. A comparison of results for different CHO isomers demonstrates the importance of steric effects. The reactions between the same groups of the PC chain and the phosphites 2-[2,4-bis(tert-butyl)phenoxy]-5,5-dimethyl-1,3,2-dioxaphosphorinane] (BPDD) and trimethyl phosphite (TMP), and their phosphonate isomers are characterized by large energy barriers.

Hydrogen bonding and collective proton modes in clusters and periodic layers of squaric acid: A density functional study

Hydrogen bonding in clusters and extended layers of squaric acid molecules has been investigated by density functional computations. Equilibrium geometries, harmonic vibrational frequencies, and energy barriers for proton transfer along hydrogen bonds have been determined using the Car-Parrinello method. The results provide crucial parameters for a first principles modeling of the potential energy surface, and highlight the role of collective modes in the low-energy proton dynamics. The importance of quantum effects in condensed squaric acid systems has been investigated, and shown to be negligible for the lowest-energy collective proton modes. This information provides a quantitative basis for improved atomistic models of the order-disorder and displacive transitions undergone by squaric acid crystals as a function of temperature and pressure. (C) 2001 American Institute of Physics.

Density functional study of reactions of phenoxides with polycarbonate

Density functional calculations with simulated annealing have been used to study the reactions of chains of bisphenol A polycarbonate (BPA-PC) with sodium phenoxide (NaOPh), diphenyl carbonate (DPC), and tetraphenylphosphonium phenoxide (PPh4OPh). These calculations extend our work on the reactions of LiOPh, NaOPh, and phenol with the cyclic tetramer of BPA-PC. We study, in particular, chain growth catalyzed by NaOPh and PPh4OH. The energy barriers for reactions with PPh4OPh are somewhat larger than those involving LiOPh and NaOPh, but they are significantly lower than those involving phenol (HOPh), due in part to the collective rearrangement of phenyl groups in the reacting molecules. We discuss in the Appendix the bonds between alkali metal atoms (Na in the present calculations) and other atoms (here oxygen) that are analogous to the more familiar "hydrogen bonds".

Fast heating of ultrahigh-density plasma as a step towards laser fusion ignition

ate studies(2) and fusion energy research(3,4). Laser-driven implosions of spherical polymer shells have, for example, achieved an increase in density of 1,000 times relative to the solid state(5). These densities are large enough to enable controlled fusion, but to achieve energy gain a small volume of compressed fuel (known as the 'spark') must be heated to temperatures of about 10(8) K (corresponding to thermal energies in excess of 10 keV). In the conventional approach to controlled fusion, the spark is both produced and heated by accurately timed shock waves(4), but this process requires both precise implosion symmetry and a very large drive energy. In principle, these requirements can be significantly relaxed by performing the compression and fast heating separately(6-10); however, this 'fast ignitor' approach(7) also suffers drawbacks, such as propagation losses and deflection of the ultra-intense laser pulse by the plasma surrounding the compressed fuel. Here we employ a new compression geometry that eliminates these problems; we combine production of compressed matter in a laser-driven implosion with picosecond-fast heating by a laser pulse timed to coincide with the peak compression. Our approach therefore permits efficient compression and heating to be carried out simultaneously, providing a route to efficient fusion energy production.

Understanding the dehydrogenation mechanism of tetrahydrocarbazole over palladium using a combined experimental and density functional theory approach

The mechanism of the dehydrogenation of tetrahydrocarbazole to carbazole over palladium has been examined for the first time. By use of a combination of deuterium exchange experiments and density functional theory calculations, a detailed reaction profile for the aromatization of tetrahydrocarbazole has been identified and validated by experiment. As with many dehydrogenation reactions, the initial hydrogen abstraction is found to have the highest reaction barrier. Tetrahydrocarbazole has four hydrogens which can, in principle, be cleaved initially; however, the theory and experiment show that the reaction is dominated by the cleavage of the carbon hydrogens at the carbon atoms in positions 1 and 4. The two pathways originating from these two C-H bond cleavage processes are found to have similar reaction energy profiles and both contribute to the overall reaction.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

Identifying an O-2 supply pathway in CO oxidation on Au/TiO2(110): A density functional theory study on the intrinsic role of water

Au catalysis has been one of the hottest topics in chemistry in the last 10 years or so. How O-2 is supplied and what role water plays in CO oxidation are the two challenging issues in the field at the moment. In this study, using density functional theory we show that these two issues are in fact related to each other. The following observations are revealed: (i) water that can dissociate readily into OH groups can facilitate O-2 adsorption on TiO2; (ii) the effect of OH group on the O-2 adsorption is surprisingly long-ranged; and (iii) O-2 can also diffuse along the channel of Ti (5c) atoms on TiO2(1 10), and this may well be the rate-limiting step for the CO oxidation. We provide direct evidence that O-2 is supplied by O-2 adsorption on TiO2 in the presence of OH and can diffuse to the interface of Au/TiO2 to participate in CO oxidation. Furthermore, the physical origin of the water effects on Au catalysis has been identified by electronic structure analyses: There is a charge transfer from TiO2 in the presence of OH to O-2, and the O-2 adsorption energy depends linearly on the 02 charge. These results are of importance to understand water effects in general in heterogeneous catalysis.

The importance of hydrogen's potential-energy surface and the strength of the forming R-H bond in surface hydrogenation reactions

An understanding of surface hydrogenation reactivity is a prevailing issue in chemistry and vital to the rational design of future catalysts. In this density-functional theory study, we address hydrogenation reactivity by examining the reaction pathways for N+H -> NH and NH+H -> NH2 over the close-packed surfaces of the 4d transition metals from Zr-Pd. It is found that the minimum-energy reaction pathway is dictated by the ease with which H can relocate between hollow-site and top-site adsorption geometries. A transition state where H is close to a top site reduces the instability associated with bond sharing of metal atoms by H and N (NH) (bonding competition). However, if the energy difference between hollow-site and top-site adsorption energies (Delta E-H) is large this type of transition state is unfavorable. Thus we have determined that hydrogenation reactivity is primarily controlled by the potential-energy surface of H on the metal, which is approximated by Delta E-H, and that the strength of N (NH) chemisorption energy is of less importance. Delta E-H has also enabled us to make predictions regarding the structure sensitivity of these reactions. Furthermore, we have found that the degree of bonding competition at the transition state is responsible for the trend in reaction barriers (E-a) across the transition series. When this effect is quantified a very good linear correlation is found with E-a. In addition, we find that when considering a particular type of reaction pathway, a good linear correlation is found between the destabilizing effects of bonding competition at the transition state and the strength of the forming N-H (HN-H) bond. (c) 2006 American Institute of Physics.

CHx hydrogenation on Co(0001): A density functional theory study

Hydrogenation is an important process in the Fischer-Tropsch synthesis. In this work, all the elementary steps of the hydrogenation from C to CH4 are studied on both flat and stepped Co(0001) using density functional theory (DFT). We found that (i) CH3 hydrogenation (CH3+H-->CH4) is the most difficult one among all the elementary reactions on both surfaces, possessing barriers of around 1.0 eV; (ii) the other elementary reactions have the barriers below 0.9 eV on the flat and stepped surfaces; (iii) CH2 is the least stable species among all the CHx(x=1-3) species on both surfaces; and (iv) surface restructuring may have little effect on the CHx(x=0-3) hydrogenation. The barriers of each elementary step on both flat and stepped surfaces are similar and energy profiles are also similar. The reason as to why CHx hydrogenation is not structure-sensitive is also discussed. (C) 2005 American Institute of Physics.

General insight into CO oxidation: A density functional theory study of the reaction mechanism on platinum oxides

CO oxidation on PtO2(110) has been studied using density functional theory calculations. Four possible reaction mechanisms were investigated and the most feasible one is the following: (i) the O at the bridge site of PtO2(110) reacts with CO on the coordinatively unsaturated site (CUS) with a negligible barrier; (ii) O-2 adsorbs on the bridge site and then interacts with CO on the CUS to form an OO-CO complex; (iii) the bond of O-OCO breaks to produce CO2 with a small barrier (0.01 eV). The CO oxidation mechanisms on metals and metal oxides are rationalized by a simple model: The O-surface bonding determines the reactivity on surfaces; it also determines whether the atomic or molecular mechanism is preferred. The reactivity on metal oxides is further found to be related to the 3rd ionization energy of the metal atom.

The mechanism of N2O formation via the (NO)(2) dimer: A density functional theory study

Catalytic formation of N2O via a (NO)(2) intermediate was studied employing density functional theory with generalized gradient approximations. Dimer formation was not favored on Pt(111), in agreement with previous reports. On Pt(211) a variety of dimer structures were studied, including trans-(NO)(2) and cis-(NO)(2) configurations. A possible pathway involving (NO)(2) formation at the terrace near to a Pt step is identified as the possible mechanism for low-temperature N2O formation. The dimer is stabilized by bond formation between one O atom of the dimer and two Pt step atoms. The overall mechanism has a low barrier of approximately 0.32 eV. The mechanism is also put into the context of the overall NO+H-2 reaction. A consideration of the step-wise hydrogenation of O-(ads) from the step is also presented. Removal of O-(ads) from the step is significantly different from O-(ads) hydrogenation on Pt(111). The energetically favored structure at the transition state for OH(ads) formation has an activation energy of 0.63 eV. Further hydrogenation of OH(ads) has an activation energy of 0.80 eV. (C) 2004 American Institute of Physics.

Insight into association reactions on metal surfaces: Density-functional theory studies of hydrogenation reactions on Rh(111)

Hydrogenation reaction, as one of the simplest association reactions on surfaces, is of great importance both scientifically and technologically. They are essential steps in many industrial processes in heterogeneous catalysis, such as ammonia synthesis (N-2+3H(2)-->2NH(3)). Many issues in hydrogenation reactions remain largely elusive. In this work, the NHx (x=0,1,2) hydrogenation reactions (N+H-->NH, NH+H-->NH2 and NH2+H-->NH3) on Rh(111) are used as a model system to study the hydrogenation reactions on metal surfaces in general using density-functional theory. In addition, C and O hydrogenation (C+H-->CH and O+H-->OH) and several oxygenation reactions, i.e., C+O, N+O, O+O reactions, are also calculated in order to provide a further understanding of the barrier of association reactions. The reaction pathways and the barriers of all these reactions are determined and reported. For the C, N, NH, and O hydrogenation reactions, it is found that there is a linear relationship between the barrier and the valency of R (R=C, N, NH, and O). Detailed analyses are carried out to rationalize the barriers of the reactions, which shows that: (i) The interaction energy between two reactants in the transition state plays an important role in determining the trend in the barriers; (ii) there are two major components in the interaction energy: The bonding competition and the direct Pauli repulsion; and (iii) the Pauli repulsion effect is responsible for the linear valency-barrier trend in the C, N, NH, and O hydrogenation reactions. For the NH2+H reaction, which is different from other hydrogenation reactions studied, the energy cost of the NH2 activation from the IS to the TS is the main part of the barrier. The potential energy surface of the NH2 on metal surfaces is thus crucial to the barrier of NH2+H reaction. Three important factors that can affect the barrier of association reactions are generalized: (i) The bonding competition effect; (ii) the local charge densities of the reactants along the reaction direction; and (iii) the potential energy surface of the reactants on the surface. The lowest energy pathway for a surface association reaction should correspond to the one with the best compromise of these three factors. (C) 2003 American Institute of Physics.

Catalytic role of gold in gold-based catalysts: A density functional theory study on the CO oxidation on gold

Gold-based catalysts have been of intense interests in recent years, being regarded as a new generation of catalysts due to their unusually high catalytic performance. For example, CO oxidation on Au/TiO2 has been found to occur at a temperature as low as 200 K. Despite extensive studies in the field, the microscopic mechanism of CO oxidation on Au-based catalysts remains controversial. Aiming to provide insight into the catalytic roles of Au, we have performed extensive density functional theory calculations for the elementary steps in CO oxidation on Au surfaces. O atom adsorption, CO adsorption, O-2 dissociation, and CO oxidation on a series of Au surfaces, including flat surfaces, defects and small clusters, have been investigated in detail. Many transition states involved are located, and the lowest energy pathways are determined. We find the following: (i) the most stable site for O atom on Au is the bridge site of step edge, not a kink site; (ii) O-2 dissociation on Au (O-2-->20(ad)) is hindered by high barriers with the lowest barrier being 0.93 eV on a step edge; (iii) CO can react with atomic O with a substantially lower barrier, 0.25 eV, on Au steps where CO can adsorb; (iv) CO can react with molecular O-2 on Au steps with a low barrier of 0.46 eV, which features an unsymmetrical four-center intermediate state (O-O-CO); and (v) O-2 can adsorb on the interface of Au/TiO2 with a reasonable chemisorption energy. On the basis of our calculations, we suggest that (i) O-2 dissociation on Au surfaces including particles cannot occur at low temperatures; (ii) CO oxidation on Au/inactive-materials occurs on Au steps via a two-step mechanism: CO+O-2-->CO2+O, and CO+O-->CO2; and (iii) CO oxidation on Au/active-materials also follows the two-step mechanism with reactions occurring at the interface.