Molecular simulation of alkyl monolayers on the Si(111) surface YUAN

The preparation of monolayers on silicon surface is of growing interest for potential applica-tions in biosensor or semiconductor technology[1—5]. The alkyl modified Si(111) surfaces[6—10] can be obtained using the thermal, catalyzed, or photochemical reaction of hydrogen-terminated sili-con with alkenes, Grignard reagents, and so on. At the same time, the monolayer properties on Si(111) surface have been studied by a variety of experimental methods[8—10] such as X-ray photo-electron spect…     

Disulfide passivation of the Ge(100)-2 × 1 surface

Understanding the bonding of sulfur at the germanium surface is important to developing good passivation routes for germanium-based electronic devices. The adsorption behavior of ethyl disulfide (EDS) and 1,8-naphthalene disulfide (NDS) at the Ge(100)-2 × 1 surface has been studied under ultrahigh vacuum conditions to investigate both their fundamental reactivity and their effectiveness as passivants of this surface. X-ray photoelectron spectroscopy, multiple internal reflection-infrared spectroscopy, and density functional theory results indicate that both molecules adsorb via S-S dissociation at room temperature. Upon exposure to ambient air, the thiolate adlayer remains intact for both EDS- and NDS-functionalized surfaces, indicating the stability of this surface attachment. Although both systems resist oxidation compared to the bare Ge(100)-2 × 1 surface, the Ge substrate is significantly oxidized in all cases (17-57% relative to the control), with the NDS-passivated surface undergoing up to two times more oxidation than the EDS-passivated surface at the longest air exposure times studied. The difference in passivation capability is attributed to the difference in surface coverage on Ge(100)-2 × 1, where EDS adsorption leads to a saturation coverage 17% higher than that for NDS/Ge(100)-2 × 1.     

Structure and stability of one-dimensional o-phthalaldehyde lines on the Si(100)-2 × 1:H surface

We present the structural models for the o-phthalaldehyde (OP) molecular lines on the H-terminated Si(100) surface which were recently observed by scanning tunnelling microscopy. Our first-principles density-functional theory calculations show that the formation of OP lines is not only kinetically more facile but also thermodynamically more stable than those of previously reported alkene lines.     

Mechanism of•OH Formation in the thermal decomposition of (≡Si–O)2Si(C2H5)(O–O•) radicals

ESR andIR spectroscopy and quantum chemical calculations are used to obtain the mechanistic data on the reaction . The IR bands are assigned by simulating the vibrational spectra of model low-molecular compounds. Quantum chemical calculations provided the data on the shapes of potential energy surface for the systems under study and transition states. These data are used to interpret the experimental data. The title reaction occurs much more readily in the case of organosilicon peroxy radicals than in the case of their hydrocarbon analogs. Surface silanone groups of silica react with ethylene molecules to form the siloxacyclobutane group.     

Investigation of the Ti/MgCl2 interface on a Si(111) 7 × 7 substrate

Photoelectron spectroscopy with synchrotron radiation, low energy electron diffraction, and ion-scattering spectroscopy were used in order to study the Ti/MgCl(2) interface grown on an atomically clean Si(111) 7 × 7 substrate. A series of high resolution spectra after deposition of a thick MgCl(2) layer, step by step deposition of Ti and gradual annealing, indicated a very reactive interface even at room temperature. Strong interaction between the incoming Ti atoms and the MgCl(2) layer, leads to the formation of Ti(2+) and Ti(4+) oxidation states. The interfacial interaction continues even at multilayer Ti coverage mainly by the partial disruption of Mg-Cl bonds and the formation of Ti-Cl sites, rendering this interface a very promising UHV-compatible model of a pre-catalyst for olefin polymerization. After the final annealing, the MgCl(2) multilayers desorb while Ti remains on the surface forming a silicide layer on which Cl and Mg atoms are attached.     

Stereo-isomerism controls surface reactivity: 1-chloropentane-pairs on Si(100)-2×1

Chloropentane forms asymmetric ('A') and symmetric ('S') pairs on Si(100)-2×1, differing in the direction of curvature of one pentane tail. Surprisingly this renders the rate of thermal reaction of 'A' fifteen times greater than 'S' in chlorinating room-temperature silicon. Correspondingly, for electron-induced reaction the energy threshold for A is 1 eV less than for S.     

Atomic structure and formation mechanism of identically sized Au clusters grown on Si(111)-(7×7) surface

Identically sized Au clusters are grown on the Si(111)-(7×7) surface by room temperature deposition of Au atoms and subsequent annealing at low-temperature. The topographical images investigated by in situ scanning tunneling microscopy show a bias-dependent feature. The current-voltage properties measured by scanning tunneling spectroscopy indicate some semiconducting characteristics of the Au adsorbed surface, which is attributable to the saturation of Si dangling bonds. These experimental results, combined with the simulated scanning tunneling microscopy images and the first-principles adsorption energy calculations, show that the Au cluster is most likely to have a Au(6)Si(3) structure. In the Au(6)Si(3) cluster, three adsorbed Au atoms replace the three Si center adatoms, forming a hollow triangle, while the replaced Si atoms and other three Au atoms connect into a hexagon locating within the triangle. The formation mechanism of this atomic configuration is intimately associated with the complicated chemical valences of Au and the specific annealing conditions.     

Thermal decomposition mechanisms of methylamine, ethylamine, and 1-propylamine on Si(100)-2 × 1 surface

The thermal decomposition reactions of methylamine, ethylamine, and 1-propylamine absorbed on Si(100)-2 × 1 surface were theoretically investigated. Eight decomposition channels were found leading to desorption products of imine, H(2), alkyl cyanide, ammonia, aziridine, alkene, azetidine, and cyclopropane, which supports the experimental assignments. Our mechanistic studies strongly suggest that the alkyl cyanide (hydrogen cyanide in the case of methylamine) channel is coupled with the hydrogen desorption step. The β-hydrogen of ethylamine and 1-propylamine was found to undergo additional decomposition reactions producing aziridine and alkene, which were classified as γ- and β-eliminations, respectively. It was also found that the γ-hydrogen of 1-propylamine undergoes azetidine and cyclopropane producing decompositions, which were classified as δ- and γ-eliminations. In general, γ- and δ-hydrogen involved decomposition reactions are kinetically less favorable than β-hydrogen involved ones. Consequently, it is expected that the thermal decompositions of the primary alkyl amines with longer alkyl chains would not add additional favorable decomposition channels. Except alkyl cyanide and ammonia desorption channels, the decompositions occur in a concerted fashion.     

Dynamics of H2 dissociation on the 1/2 ML c(2 × 2)-Ti/Al(100) surface

The dissociation of H(2) on Ti-covered Al surfaces is relevant to the rehydrogenation and dehydrogenation of the NaAlH(4) hydrogen storage material. The energetically most stable structure for a 1/2 monolayer of Ti deposited on the Al(100) surface has the Ti atoms in the second layer with a c(2 × 2) structure, as has been confirmed by both low-energy electron diffraction and low-energy ion scattering experiments and density functional theory studies. In this work, we investigate the dynamics of H(2) dissociation on a slab model of this Ti/Al(100) surface. Two six-dimensional potential energy surfaces (PESs) have been built for this H(2) + Ti/Al(100) system, based on the density functional theory PW91 and RPBE exchange-correlation functionals. In the PW91 (RPBE) PES, the lowest H(2) dissociation barrier is found to be 0.65 (0.84) eV, with the minimum energy path occurring for H(2) dissociating above the bridge to top sites. Using both PESs, H(2) dissociation probabilities are calculated using the classical trajectory (CT), the quasi-classical trajectory (QCT), and the time-dependent wave-packet methods. We find that the QCT H(2) dissociation probabilities are in good agreement with the quantum dynamics results in the collision energy range studied up to 1.0 eV. We have also performed molecular beam simulations and present predictions for molecular beam experiments. Our molecular beam simulations show that H(2) dissociation on the 1/2 ML Ti/Al(100) surface is an activated process, and the reaction probability is found to be 6.9% for the PW91 functional and 1.8% for the RPBE at a nozzle temperature of 1700 K. Finally, we have also calculated H(2) dissociation rate constants by applying transition state theory and the QCT method, which could be relevant to modeling Ti-catalyzed rehydrogenation and dehydrogenation of NaAlH(4).     

甲亚胺1,2—脱氢过程中化学键变化的量子拓扑研究

在IRC解析的基础上,用电子密度的量子拓扑分析方法研究了甲亚胺1,2-脱氢过程中化学键的变化.反应途径中各点化学键的拓扑性质可以清楚地反映出化学键断裂及生成过程,计算结果进一步证明该反应为协同非同步反应.本工作为研究化学反应过程提供了一种新的方法。     

Topological Analysis of Electron Densities: Is the Presence of an Atomic Interaction Line in an Equilibrium Geometry a Sufficient Condition for the Existence of a Chemical Bond?

The structure, energetics, and electron density in the inclusion complex of He in adamantane, C10H16, have been studied by density functional theory calculations at the B3LYP6-311++G(2p,2d) level. Topological analysis of the electron density shows that the He atom is connected to the four tertiary tC atoms in the cage by atomic interaction lines with (3,-1) critical points. The calculated dissociation energy of the complex He@adamantane(g)=adamantane(g) + He(g) of DeltaE=-645 kJ mol(-1) nevertheless shows that the He-tC interactions are antibonding.     

Topology and Chemistry

The determinants of chemical bonding are the chemical properties of the atoms and the constraints of three-dimensional (3-D) space into which the atoms must fit, but topology provides a convenient way of describing the resultant structure. This paper explores the topologies of various scalar fields associated with atoms in molecules and crystals and what they can tell us about chemical bonding. The scalar fields examined are the electron density, the electrostatic potential, and two simplified electrostatic potentials in which the contributions of the electron cores have been removed, namely the Madelung and the covalent field. Not all of the information contained in these fields is present in the topology but, since the topology is insensitive to the details of the field, it can often be determined using simplified calculations. Although the same topological model is used to explore all four fields, each field has its own distinctive topology and each provides different information about the nature of chemical bonding and structure. The analysis of these topologies, when combined with simple electrostatic theory and a few empirical observations, leads to a quantitative model of localized chemical bonding. In the process, the analysis provides insights into the nature of chemical bonding.     

Structures and bonding patterns of nanoannular carbon clusters (C4–C20) through AIM analyses

Structures, binding energies, harmonic frequencies, dipole moments, HOMO–LUMO energy gaps and particularly atoms in molecules (AIM) analyses of some nanoannular carbon clusters (C₄–C₂₀) are investigated at B3LYP/6-31+G(d) level of theory. No correlation is found by plotting the calculated binding energies as a functional number of carbon atoms of carbon clusters. The calculated binding energies sharply increase from C₄ to C₁₀ while slowly from C₁₀ to C₂₀. The binding energies of C_4n+2 clusters including C₆, C₁₀, C₁₄, and C₁₈ have a clear increase when compared with others indicating their aromatic characters which is confirmed by results of HOMO–LUMO energy gaps and geometrical parameters. AIM analyses show that most of our carbon clusters are topologically normal (non-conflict) with stable structures. Nevertheless, the topological networks of small antiaromatic rings, C₄ and C₈, at their equilibrium geometries may change via molecular vibrations. The existence of straight bond paths in 3D molecular graphs of carbon clusters with n > 10 implies that ring strains are decreased as the ring sizes grow. Except for C₄ and C₈, the ellipticity values for the remaining carbon clusters are small indicating that the C–C bond is conserved in these clusters. Dipole moments of even-numbered structures are negligible, whereas odd-numbered ones have μ values of 0.09−0.73 D.     

亚烷基卡宾与亚烷基锂氟类卡宾化学键结构的量子拓扑研究

用ab initio(3-21G)方法对亚烷基卡宾H_2C—C的单线态及三线态结构进行了电子密度拓扑分析,说明了它们的亲电、亲核反应方向,讨论了亚烷基锂氟类卡宾H_2C—CLiF的4种构型,论证了该分子中不存在四元环结构、Li—C键以静电作用为主的特性,并预测了加成反应机理。     

HNO与(HF)1≤n≤3分子间的蓝移与红移氢键

运用量子化学从头算方法研究了HNO与分子簇(HF)1≤n≤移氢键, 重极化与重杂化和分子内超共轭导致了氢键的蓝移; 所有的X…H—F(X=O, N, F)氢键都是红移的, 分子间超共轭导致了氢键的红移. 在多分子体系形成的氢键链中, 分子间超共轭作用呈现规律性递变, 它导致了氢键强度与频率位移的规律性变化, 电子密度拓扑分析结果反映和支持了这种规律性变化.     

Manganese,Iron, Cobalt,and Nickel Oxo‐, Peroxo‐, and Superoxoclusters: A Density Functional Theory Study

The 3d‐transition‐metal dioxo‐, peroxo‐, and superoxoclusters with the general composition MO₂, M(O₂), and MOO (M=Mn, Fe, Co, and Ni) were studied by DFT by the B1LYP functional. The dioxides in their ground states represent the global minima for the M+O₂ system. Both ground‐state dioxides and the lowest‐energy peroxides are in their (d‐only) highest spin states. The ⁶A₁ state of Co(O₂) exceeds the d‐only spin‐multiplicity value (quartet), being nearly isoenergetic with the ⁴A₁ state of Co(O₂). The energy gain on transforming the peroxides to the corresponding dioxides decreases in the order Mn(O₂)>Fe(O₂)>Co(O₂)>Ni(O₂) and varies in the range 0.27–1.8 eV. The dissociation energy to M+O₂ for all studied peroxides is less than 1 eV being the lowest (0.47 eV) for Mn(O₂). The Mn and Fe peroxides need less than 0.3 eV to rupture one of the MO bonds to form the corresponding superoxide. Mn and Fe superoxides are less stable than the corresponding peroxides; the superoxide of Co is more stable than its peroxide, while Ni superoxide is unstable—its energy is above the limit of dissociation to Ni+O₂. According to the electrostatic potential maps, the oxygen atoms in the peroxides are more nucleophilic than those in the dioxides and superoxides, in which the terminal oxygen atom is more nucleophilic than the M‐bonded oxygen atom. This result differs from the expectations based on charge‐distribution analysis.