Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

Non-linear conductance of disordered quantum wires

The self-consistent electron potential in a current-carrying disordered quantum wire is spatially inhomogeneous due to the formation of resistivity dipoles across scattering centres. In this paper it is argued that these inhomogeneities in the potential result in a suppression of the differential conductance of such a wire at finite applied voltage. A semi-classical argument allows this suppression, quadratic in the voltage, to be related directly to the amount of intrinsic defect scattering in the wire. This result is then tested against numerical calculations.

Methane transformation to carbon and hydrogen on Pd(100): Pathways and energetics from density functional theory calculations

Density functional theory with gradient corrections has been employed to study the reaction pathways and the reaction energetics for the transformations of CH4 to C and H on a Pd(100) surface. On examination of transition state structures identified in each elementary reaction, a clear relationship between the valencies of the CHx fragments and the locations of the transition states emerges. The higher the valency of the CHx fragment, the higher the coordination number of the CHx with the surface atoms. The calculated reaction energetics are in good agreement with the experiments. In addition, calculation results are also used to illustrate an interesting issue concerning the CH3 stability on Pd surfaces. (C) 2002 American Institute of Physics.

Model calculations for copper clusters on Au(111) surfaces

Using the semi-empirical embedded-atom method, the structure of small copper clusters on Au(111) surfaces has been investigated both by static and dynamic calculations. By varying the size of roughly circular clusters, the edge energy per atom is obtained; it agrees quite well with estimates based on experimental results. Small three-dimensional clusters tend to have the shape of a pyramid, whose sides are oriented in the directions of small surface energy. The presence of a cluster is found to distort the underlying lattice of adsorbed copper atoms. (C) 2002 Published by Elsevier Science B.V.

Quantum efficiency in GaAs Schottky photodetectors with enhancement due to surface plasmon excitations

The overall quantum efficiency in surface plasmon (SP) enhanced Schottky barrier photodetectors is examined by considering both the external and internal yield. The external yield is considered through calculations of absorption and transmission of light in a configuration that allows reflectance minimization due to SP excitation. Following a Monte Carlo method, a procedure is presented to estimate the internal yield while taking into account the effect of elastic and inelastic scattering processes on excited carriers subsequent to photon absorption. The relative importance of internal photoemission and band-to-band contributions to the internal yield is highlighted along with the variation of the yield as a function of wavelength, metal thickness and other salient parameters of the detector. (C) 2002 Elsevier Science Ltd. All rights reserved.

Prospects for ultracold carbon via charge exchange reactions and laser cooled carbides

Strategies to produce an ultracold sample of carbon atoms are explored and assessed with the help of quantum chemistry. After a brief discussion of the experimental difficulties using conventional methods, two strategies are investigated. The first attempts to exploit charge exchange reactions between ultracold metal atoms and sympathetically cooled C+ ions. Ab initio calculations including electron correlation have been conducted on the molecular ions [LiC]+ and [BeC]+ to determine whether alkali or alkaline earth metals are a suitable buffer gas for the formation of C atoms but strong spontaneous radiative charge exchange ensure they are not ideal. The second technique involves the stimulated production of ultracold C atoms from a gas of laser cooled carbides. Calculations on LiC suggest that the alkali carbides are not suitable but the CH radical is a possible laser cooling candidate thanks to very favourable Frank-Condon factors. A scheme based on a four pulse STIRAP excitation pathway to a Feshbach resonance is outlined for the production of atomic fragments with near zero centre of mass velocity.

Theory of the Kinetics of Chemical Potentials in Heterogeneous Catalysis

Simple and powerful: The reaction kinetics at surfaces of heterogeneous catalysts is reformulated in terms of the involved chemical potentials. Based on this formulism, an approach of searching for good catalysts is proposed without recourse to extensive calculations of reaction barriers and detailed kinetic analyses. (see picture; R=reactant, I=surface intermediate, P=product, and =standard chemical potential).

A 330-360 GHz spectral survey of G34.3+0.15 .2. Chemical modelling

We describe a detailed depth-and time-dependent model of the molecular cloud associated with the ultracompact H II region G 34.3+0.15. Previous work on observations of NH3 and CS indicates that the molecular cloud has three distinct physical components:- an ultracompact hot core, a compact hot core and an extended halo. We have used the physical parameters derived from these observations as input to our detailed chemical kinetic modelling. The results of the model calculations are discussed with reference to the different chemistries occuring in each component and are compared with abundances derived from our recent spectral line survey of G 34.3+0.15 (Paper I).

Effect of positron-atom interactions on the annihilation gamma spectra of molecules

Calculations of ?-spectra for positron annihilation on a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The annihilation ?-spectra characterise the momentum distribution of the electron-positron pair at the instant of annihilation. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum densities calculated using modern computational quantum chemistry density functional theory tools. The calculation, in its simplest form, effectively treats the low-energy (thermalised, room-temperature) positron as a plane wave and gives annihilation ?-spectra that are about 40% broader than experiment, although the main chemical trends are reproduced. We show that this effective 'narrowing' of the experimental spectra is due to the action of the molecular potential on the positron, chiefly, due to the positron repulsion from the nuclei. It leads to a suppression of the contribution of small positron-nuclear separations where the electron momentum is large. To investigate the effect of the nuclear repulsion, as well as that of short-range electron-positron and positron-molecule correlations, a linear combination of atomic orbital description of the molecular orbitals is employed. It facilitates the incorporation of correction factors which can be calculated from atomic many-body theory and account for the repulsion and correlations. Their inclusion in the calculation gives -spectrum linewidths that are in much better agreement with experiment. Furthermore, it is shown that the effective distortion of the electron momentum density, when it is observed through positron annihilation -spectra, can be approximated by a relatively simple scaling factor. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.

Sonochemical synthesis and measurement of optical properties of zinc sulfide quantum dots

A facile sonochemical method has been developed to prepare very small zinc sulfide nanoparticles (ZnS NPs) of extremely small size about 1. nm in diameter using a set of ionic liquids based on the bis (trifluoromethylsulfonyl) imide anion and different cations of 1-alkyl-3-methyl-imidazolium. The structural features and optical properties of the NPs were determined in depth with X-ray powder diffraction (XRD), transmission electron microscopy (TEM), high-resolution TEM (HRTEM), Fourier transform infrared (FTIR) spectroscopy, dynamic light scattering (DLS) analysis, and UV-vis absorption spectroscopy. The energy band gap measurements of ZnS NPs were calculated by UV-vis absorption spectroscopy. One of the interesting features of the present work is that the wide band gap semiconductor ZnS nanocrystals were prepared which are used in the fabrication of photonic devices.

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Quantum chemistry and in situ FTir spectroscopy studies on potential-dependent properties of CO adsorbed on Pt electrodes

The electronic and vibrational properties of CO adsorbed on Pt electrodes at different potentials have been studied, by using methods of self-consistent-charge discrete variational Xa (SCC-DV-Xa) cluster calculations and in situ FTir spectroscopy. Two new models have been developed and verified to be successful: (1) using a "metallic state cluster" to imitate a metal (electrode) surface; and (2) charging the cluster and shifting its Fermi level (e{lunate}) to simulate, according to the relation of -d e{lunate}e dE, quantitatively the variation of the electrode potential (E). It is shown that the binding of PtCO is dominated by the electric charge transfer of dp ? 2p, while that of s ? Pt is less important in this binding. The electron occupancy of the 2p orbital of CO weakens the CO bond and decreases the v. Variation of E mainly influences the charge transfer process of dp ? 2p, but hardly influences that of s ? Pt. A linear potential-dependence of v has been shown and the calculated dv/dE = 35.0 cm V. All results of calculations coincide with the ir experimental data. © 1993.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. (C) 2014 AIP Publishing LLC.