Fragment distribution of thermal decomposition for PS and PET with QMD calculations by considering the excited and charged model molecules

Simulations by a quantum molecular dynamics (QMD) (MD with MO) method were demonstrated on the thermal decomposition of PS and PET polymers using the model molecules at the ground state including excited and positive charged states. For the excited and positive charged model molecules, we adopted CH₃CHC₆H₅CH₃ and CH₃OCOC₆H₄COOCH₃ of PS and PET monomers, respectively at the singlet and triplet states in single excitation, and at (+2) positive charged state by semiempirical AM1 MO method. Geometry and energy optimized results of the excited and positive charged models by MO calculations were used as the initial MD step of QMD calculations. In the QMD calculations, we controlled the total energy of the system using Nosé-Hoover thermostats in the total energy range of 0.69-0.95 eV, and the sampling position data with a time step of 0.5 fs were carried out up to 5000 steps at 60 different initial conditions. The calculated neutral, positive and negative charged fragment distributions of PS and PET models with 0.82 eV energy control were obtained as (93.5, 2.3, and 4.3%), and (87.8, 5.3, and 6.9%) to the total fragments, respectively. The ratios seem to correspond well to the values observed experimentally in SIMS.     

Simulation of SIMS for monomer and dimer of lignin under the assumption of thermal decomposition using QMD method

The thermal decomposition of the monomer and dimer of lignin has been simulated by a quantum molecular dynamics (DMD) method. In the calculation, we controlled the total energy of the system using Nóse-Hoover thermostats in the total energy range of 0.69-0.95 eV, and the sampling position data with a time step of 0.5 fs were carried out up to 3000 (1.5 ps) or 5000 (2.5 ps) steps in ab initio and semiempirical MO methods, respectively. We obtained the thermally decomposed fragments with positive, neutral and negative charges from SCF MO calculation at each data of the last MD step, and simulated the fragment distribution of the monomer and dimer lignins from the last step in 30-40 runs. Simulated mass numbers of positively and negatively charged fragments for lignin monomer and dimer showed considerably good accordance with the experimental results in TOF-SIMS observed by Saito and co-workers.     

Hydrocarbon film growth by energetic CH3 molecule impact on SiC (0 0 1) surface

Classic molecular dynamics (MD) calculations were performed to investigate the deposition of thin hydrocarbon film. SiC (1 0 0) surfaces were bombarded with energetic CH₃ molecules at impact energies ranging from 50 to 150 eV. The simulated results show that the deposition yield of H atoms decreases with increasing incident energy, which is in good agreement with experiments. During the initial stages, with breaking Si-C bonds in SiC by CH₃ impacting, H atoms preferentially reacts with resulting Si to form Si-H bond. The C/H ratio in the grown films increases with increasing incident energy. In the grown films, CH species are dominant. For 50 eV, H-Csp3 bond is dominant. With increasing energy to 200 eV, the atomic density of H-Csp2 bond increases.     

The adsorption and dissociation of O2 molecular precursors on Cu: the effect of steps

The adsorption and dissociation of dioxygen on Cu steps are studied using periodic self-consistent density functional theory (PW91-GGA) calculations. Cu steps are modeled with a Cu(2 1 1) surface. The results are compared with those on the flat Cu(1 1 1) surface. The adsorption of both atomic and molecular oxygen is enhanced on the stepped surface: the binding energy of atomic oxygen is −4.5 eV at its preferred site on the relaxed Cu(2 1 1) surface vs. −4.3 eV at its preferred site on the relaxed Cu(1 1 1) surface, and the binding energy of the molecular oxygen precursor is increased from ∼−0.6 to ∼−1.0 eV. Several possible O₂ dissociation paths at the edge of the Cu(2 1 1) step have been investigated. The activation energies range from 0.09 to 0.24 eV, comparable to a minimum activation energy of 0.20 eV found on Cu(1 1 1). However, compared to Cu(1 1 1) the paths on Cu(2 1 1) are stabilized in their entirety by the step by ∼0.5 eV in terms of initial state, transition state, and final state energies. The dissociation of O₂ precursors at the foot of the step is close to being barrier-less. Because of the small dissociation barrier on Cu(1 1 1), the effect of steps on O₂ dissociation is nevertheless not expected to be as pronounced as in other gas/metal systems.     

Structural and electronic properties of H-passivated graphene

The atomic and electronic structures of graphane (hydrogen-passivated graphene) are theoretically investigated using the local density approximation (LDA) of the density functional theory (DFT) and the pseudopotential method. Our total energy calculations suggest that the chairlike configuration for graphane is more energetically stable than the boatlike and tablelike configurations by approximately 0.129 eV/cell and 0.655 eV/cell, respectively. Our calculations suggest that the LDA band gap of the chairlike structure is approximately 3.9 eV. The equilibrium geometry and the band structure of the chairlike conformer are investigated and compared with the available experimental and theoretical data. We further present total and partial charge density to reveal the orbital nature of the highest occupied and the lowest unoccupied states.     

Ab-initio calculations of electronic structure and optical properties of TiAl alloy

The electronic structures and optical properties of TiAl intermetallic alloy system are studied by the first-principle orthogonalized linear combination of atomic orbitals method. Results on the band structure, total and partial density of states, localization index, effective atomic charges, and optical conductivity are presented and discussed in detail. Total density of states spectra reveal that (near the Fermi level) the majority of the contribution is from Ti-3d states. The effective charge calculations show an average charge transfer of 0.52 electrons from Ti to Al in primitive cell calculations of TiAl alloy. On the other hand, calculations using supercell approach reveal an average charge transfer of 0.48 electrons from Ti to Al. The localization index calculations, of primitive cell as well as of supercell, show the presence of relatively localized states even above the Fermi level for this alloy. The calculated optical conductivity spectra of TiAl alloy are rich in structures, showing the highest peak at 5.73 eV for supercell calculations. Calculations of the imaginary part of the linear dielectric function show a prominent peak at 5.71 eV and a plateau in the range 1.1-3.5 eV.     

Coverage effects on phenol adsorption on Si(1 0 0)2 × 1 as: A first principle calculation

We investigated the adsorption of a 6-dimers Si(1 0 0)2 × 1 surface as a function of coverage and adsorption type (molecular/dissociative) by first principle calculations. In particular, we performed calculations on models with 2, 3, 4 and 6 phenol molecules, corresponding to coverage Θ = 0.34, 0.5, 0.67 and 1. We found that total adsorption energy, when at least one phenol is in a molecular state is lower than the sum of the corresponding singly adsorbed molecules. The dissociative adsorption of multiple molecules, both in parallel and switched configuration is most favoured for a coverage Θ = 0.34 (2.6 eV per adsorbed molecule). This values decreases to 2.0 eV and remains constant till the coverage 1 is reached.The energy barrier for the molecular-to-dissociated transition of a phenol molecule, in presence of another dissociatively adsorbed molecule is ∼0.008 eV and it is similar to the value in case of single adsorption. Possible hydrogen displacements were also considered.     

Half-metallic ferromagnetism in rocksalt and zinc-blende MS (M=Li, Na and K): A first-principles study

First-principles full-potential linearized augmented plane-wave method is used to investigate the electronic structure and magnetic properties of hypothetical zinc-blende and rocksalt LiS, NaS and KS. We find that all the compounds except rocksalt LiS exhibit half-metallic ferromagnetism with an integer magnetic moment of 1.00 μ_B per formula unit. The ferromagnetism results from the spin-polarization of p states of anion S. Total energies calculations indicate the rocksalt phase is lower in energy than the zinc-blende one. The total energy differences are about 0.38, 0.36 and 0.32 eV per formula unit for LiS, NaS and KS, respectively. Meanwhile, it is shown that rocksalt NaS and KS have the half-metallic gaps of 0.22 and 0.41 eV, respectively, and the half-metallic gaps are 0.03, 0.46 and 0.65 eV for zinc-blende LiS, NaS and KS, respectively. We also find the half-metallicity is robust against the lattice contraction up to 7% and 13% for rocksalt NaS and KS, respectively. Although rocksalt LiS is nonmagnetic and metallic at the equilibrium lattice constant, it shows half-metallic ferromagnetism when the lattice constant is larger than 5.40 Å.     

Pd diffusion on MgO(1 0 0): The role of defects and small cluster mobility

Density functional theory is used to explore the energy landscape of Pd atoms adsorbed on the terrace of MgO(1 0 0) and at oxygen vacancy sites. Saddle point finding methods reveal that small Pd clusters diffuse on the terrace in interesting ways. The monomer and dimer diffuse via single atom hops between oxygen sites with barriers of 0.34 eV and 0.43 eV respectively. The trimer and tetramer, however, form 3D clusters by overcoming a 2D-3D transition barrier of less than 60 meV. The trimer diffuses along the surface either by a walking or flipping motion, with comparable barriers of ca. 0.5 eV. The tetramer rolls along the terrace with a lower barrier of 0.42 eV. Soft rotational modes at the saddle point lead to an anomalously high prefactor of 1.3 × 10¹⁴ s⁻¹ for tetramer diffusion. This prefactor is two order of magnitude higher than for monomer diffusion, making the tetramer the fastest diffusing species on the terrace at all temperatures for which diffusion is active (above 200 K). Neutral oxygen vacancy sites are found to bind Pd monomers with a 2.63 eV stronger binding energy than the terrace. A second Pd atom, however, binds to this trapped monomer with a smaller energy of 0.56 eV, so that dimers at defects dissociate on a time scale of milliseconds at room temperature. Larger clusters bind more strongly at defects. Trimers and tetramers dissociate from monomer-bound-defects at elevated temperatures of ca. 600 K. These species are also mobile on the terrace, suggesting they are important for the ripening observed at ?600 K during Pd vapor deposition on MgO(1 0 0) by Haas et al. [G. Haas, A. Menck, H. Brune, J.V. Barth, J.A. Venables, K. Kern, Phys. Rev. B 61 (2000) 11105].     

Aluminum and nitrogen impurities in Wurtzite ZnO: first-principles studies

First-principles calculations have been performed to investigate the doping behaviors of Al and N dopant impurities in ZnO. According to the results, in the Al mono-doping case, the impurity states are quite delocalized, the corresponding effective masses are small, and the formation energy is as low as −9.71 eV. In the N mono-doping case, the impurity states are localized, the effective masses are large, and the formation energy is high (4.55 eV in the most favorable extreme O-rich conditions). In the Al-N codoping case, the corresponding effective masses are marked decreased compared to the N mono-doping situations, and the formation energy of the N-Al-N system is as low as −2.54 eV in the O-rich condition. The above results can explain the electrical behaviors of the doped or codoped ZnO systems observed in experiments.     

KLL Auger resonant Raman transition in metallic Cu and Ni

High energy resolution KL₂₃L₂₃ Auger spectra of polycrystalline Cu and Ni were measured using photon energies up to about 50 eV above the K-absorption edge and down to 5 eV (Cu KLL) and 4 eV (Ni KLL) below threshold. The spectra show strong satellite structures varying considerably as a function of the photon energy. In the sub-threshold region the linear dispersion of the diagram line energy positions and a distortion of the line shape as a function of photon energy, attributable to the Auger resonant Raman process, is clearly observed, indicating the one-step nature of the Auger emission. These changes in the resonant spectra are interpreted using a simple model based on resonant scattering theory in combination with partial density of states obtained from cluster molecular orbital (DV-Xα) calculations.     

Reaction mechanism for CO oxidation on Cu(3 1 1): A density functional theory study

The microscopic reaction mechanism for CO oxidation on Cu(3 1 1) surface has been investigated by means of comprehensive density functional theory (DFT) calculations. The elementary steps studied include O₂ adsorption and dissociation, dissociated O atom adsorption and diffusion, as well as CO adsorption and oxidation on the metal. Our results reveal that O₂ is considerably reactive on the Cu(3 1 1) surface and will spontaneously dissociate at several adsorption states, which process are highly dependent on the orientation and site of the adsorbed oxygen molecule. The dissociated O atom may likely diffuse via inner terrace sites or from a terrace site to a step site due to the low barriers. Furthermore, we find that the energetically most favorable site for CO molecule on Cu(3 1 1) is the step edge site. According to our calculations, the reaction barrier of CO + O → CO₂ is about 0.3 eV lower in energy than that of CO + O₂ → CO₂ + O, suggesting the former mechanism play a main role in CO oxidation on the Cu(3 1 1) surface.     

A density-functional study on the atomic geometry and adsorption of the Cu(1 0 0) c(2 × 2)/N Surface

In this work we have performed total-energy calculations on the geometric structure and adsorption properties of Cu(1 0 0) c(2 × 2)/N surface by using the density-functional theory and the projector-augmented wave method. It is concluded that nitrogen atom was adsorbed on a FFH site with a vertical distance of 0.2 Å towards from surface Cu layer. The bond length of the shortest Cu-N bonding is calculated to be 1.83 Å. Geometry optimization calculations exclude out the possibilities of adsorbate induced reconstruction mode suggested by Driver and Woodruff and the atop structural model. The calculated workfunction for this absorbate-adsorbent system is 4.63 eV which is quite close to that of a clean Cu(1 0 0) surface. The total-energy calculations showed that the average adsorption energy per nitrogen in the case of Cu(1 0 0) c(2 × 2)-N is about 4.88 eV with respect to an isolated N atom. The absorption of nitrogen on Cu(1 0 0) surface yields the hybridization between surface Cu atoms and N, and generates the localized surface states at −1.0 eV relative to Fermi energy E_F. The stretch mode of the adsorbed nitrogen at FFH site is about 30.8 meV. The present study provides a strong criterion to account for the local surface geometry in Cu(1 0 0) c(2 × 2)/N surface.     

Local structure and electronic properties of BaTaO2N with perovskite-type structure

First-principles calculation based on density-functional theory in the pseudo-potential approach have been performed for the total energy and crystal structure of BaTaO₂N. The calculations indicate a random occupation of the anionic positions by O and N in a cubic structure, in agreement with neutron diffraction measurements and infrared spectra. The local symmetry in the crystal is broken, maintaining a space group Pm3?m, as used in structure refinement, which represents only the statistically averaged result. The calculations also indicate displacive disordering in the crystal. The average Ta-N distance is smaller (2.003 Å), while the average Ta-O distance becomes larger (2.089 Å). The local relaxation of the atoms has an influence on the electronic structure, especially on the energy gap. BaTaO₂N is calculated to be a semiconductor with an energy gap of about 0.5 eV. The upper part of the valence band is dominated by N 2p states, while O 2p states are mainly in the lower part. The conduction band is dominated by Ta 5d states.     

A density functional study of atomic oxygen and water molecule adsorption on Ni(1 1 1) and chromium-substituted Ni(1 1 1) surfaces

Density functional theory (DFT) for generalized gradient approximation calculations has been used to study the adsorption of atomic oxygen and water molecules on Ni(1 1 1) and different kind of Ni-Cr(1 1 1) surfaces. The fcc hollow site is energetically the most favorable for atomic oxygen adsorption and on top site is favorable for water adsorption. The Ni-Cr surface has the highest absorption energy for oxygen at 6.86 eV, followed by the hcp site, whereas the absorption energy is 5.56 eV for the Ni surface. The Ni-O bond distance is 1.85 Å for the Ni surface. On the other hand, the result concerning the Ni-Cr surface implies that the bond distances are 1.93-1.95 Å and 1.75 Å for Ni-O and Cr-O, respectively. The surface adsorption energy for water on top site for two Cr atom substituted Ni-Cr surface is 0.85 eV. Oxygen atoms prefer to bond with Cr rather than Ni atoms. Atomic charge analysis demonstrates that charge transfer increases due to the addition of Cr. Moreover, a local density of states (LDOS) study examines the hybridization occurring between the metal d orbital and the oxygen p orbital; the bonding is mainly ionic, and water bonds weakly in both cases.     

A study of core and valence levels in β-PbO2 by hard X-ray photoemission

The core and valence levels of β-PbO₂ have been studied using hard X-ray photoemission spectroscopy (hν = 6000 eV and 7700 eV). The Pb 4f core levels display an asymmetric lineshape which may be fitted with components associated with screened and unscreened final states. It is found that intrinsic final state screening is suppressed in the near-surface region. A shift in the O 1s binding energy due to recoil effects is observed under excitation at 7700 eV. It is shown that conduction band states have substantial 6s character and are selectively enhanced in hard X-ray photoemission spectra. However, the maximum amplitude in the Pb 6s partial density of states is found at the bottom of the valence band and the associated photoemission peak shows the most pronounced enhancement in intensity at high photon energy.     

Diffusion of Pb adatom and ad-dimer on Si(1 0 0) from ab initio studies

The diffusion pathways of Pb adatoms and ad-dimers on Si(1 0 0) are investigated by first-principles calculations. Pb adatoms are found to diffuse on top of the Si(1 0 0) dimer row with an energy barrier of 0.31 eV. However, Pb dimers are energetically more stable. Pb dimers on top of the dimer row have a high energy barrier (0.95 eV) to rotate from the lowest energy configuration to the orientation parallel to the underlying Si(1 0 0) dimer row. Once the ad-dimer is oriented parallel to Si(1 0 0) dimer row, they can diffuse along the dimer row with an energy barrier of only 0.32 eV.     

First-principles study of the optical properties of BeO in its ambient and high-pressure phases

Optical properties such as the dynamic dielectric function, reflectance, and energy-loss function of beryllium oxide (BeO) in its ambient and high-pressure phases are reported for a wide energy range of 0-50 eV. The calculations of optical properties employ first-principles methods based on all-electron density functional theory together with sum over states and finite-field methods. Our results show subtle differences in the calculated optical properties of the wurtzite, zincblende, rocksalt and CsCl phases of BeO, which may be attributed to the higher symmetry and packing density of these phases. For the wurtzite phase, the calculated band gap of 10.4 eV corresponds well with the experimental value of 10.6 eV and the calculated (average) index of refraction of 1.70 shows excellent agreement with the experimental value of 1.72.     

Theoretical treatment of excited electronic states of adsorbates on metals: Electron attachment to CO2 adsorbed on Pt(1 1 1)

Photochemistry involving adsorbates on metals often proceeds by photoexcitation of the metal followed by transient attachment of photoemitted electrons to the adsorbate. First principles theoretical methods suitable for describing electronic states embedded in a near continuum of metal to metal excitations are described and an application to electron attachment to CO₂ adsorbed on Pt(1 1 1) is reported. Wavefunctions are constructed by ab initio configuration interaction methods which allow a rigorous resolution of states and differentiation between competing pathways of molecular desorption and dissociation. An embedding theory is used to achieve high accuracy in the adsorbate-surface region. The energy required to form the electron attached state is 5.2 eV for excitation to bent CO₂ and 6.8 eV for excitation to linear CO₂, hence both energies are near the work function of the metal (5.7 eV). The process also involves localization of the metal hole and attraction of the charged adsorbate to the metal. Optimum geometries are calculated and pathways that results in desorption, dissociation by bond rupture directly in the excited electronic state, or dissociation after return to the ground state potential energy surface via vibrational processes are explored.