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.

Calculation of gamma spectra for positron annihilation on molecules

Calculations of gamma spectra for positron annihilation for a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum distributions calculated using the density functional theory (DFT) based B3LYP/TZVP model. For positrons thermalised to room temperature, the calculation, in its simplest form, effectively treats the positron as a plane wave and gives positron annihilation ?-spectra linewidths that are broader (30-40%) than experiment, although the main chemical trends are reproduced. The main physical reason for this is the neglect of positron repulsion from the nuclei. We show that this effect can be incorporated through momentum-dependent correction factors, determined from positron-atom calculations, e.g., many-body perturbation theory. Inclusion of these factors in the calculation gives linewidths that are in improved agreement with experiment.

Radiative charge transfer in cold and ultracold sulphur atoms colliding with protons

Radiative decay processes at cold and ultra cold temperatures for sulfur atoms colliding with protons are investigated. The MOLPRO quantum chemistry suite of codes was used to obtain accurate potential energies and transition dipole moments, as a function of internuclear distance, between low-lying states of the SH+ molecular cation. A multi-reference configuration-interaction approximation together with the Davidson correction is used to determine the potential energy curves and transition dipole moments, between the states of interest, where the molecular orbitals are obtained from state-averaged multi-configuration-self-consistent field calculations. The collision problem is solved approximately using an optical potential method to obtain radiative loss, and a fully two-channel quantum approach for radiative charge transfer. Cross sections and rate coefficients are determined for the first time for temperatures ranging from 10 μK up to 10 000 K. Results are obtained for all isotopes of sulfur, colliding with H+ and D+ ions and comparison is made to a number of other collision systems.

A density functional theory study of CO and atomic oxygen chemisorption on Pt(111)

Ab initio total energy calculations within a density functional theory framework have been performed for CO and atomic oxygen chemisorbed on the Pt(111) surface. Optimised geometries and chemisorption energies for CO and O on four high-symmetry sites, namely the top, bridge, fee hollow and hcp hollow sites, are presented, the coverage in all cases being 0.25 ML. The differences in CO adsorption energies between these sites are found to be small, suggesting that the potential energy surface for CO diffusion across Pt(111) is relatively flat. The 5 sigma and 2 pi molecular orbitals of CO are found to contribute to bonding with the metal. Some mixing of the 4 sigma and 1 pi molecular orbitals with metal states is also observed. For atomic oxygen, the most stable adsorption site is found to be the fee hollow site, followed in decreasing order of stability by the hcp hollow and bridge sites, with the top site being the least stable. The differences in chemisorption energies between sites for oxygen are larger than in the case of CO, suggesting a higher barrier to diffusion for atomic oxygen. The co-adsorption of CO and O has also been investigated. Calculated chemisorption energies for CO on an O/fcc-precovered surface show that of the available chemisorption sites, the top site at the oxygen atom's next-nearest neighbour surface metal atom is the most stable, with the other four sites calculated bring at least 0.29 eV less stable. The trend of CO site stability in the coadsorption system is explained in terms of a 'bonding competition' model. (C) 2000 Elsevier Science B.V. All rights reserved.

Rydberg structure associated with core-excited autoionizing states of the LiH radical

We have investigated inner-shell excitation of the LiH + molecular ion by electron impact within several different collision models to delineate Rydberg autoionizing resonance structure associated with the LiH + (1σ2σ 2 2 Σ + ) core-excited threshold. The minimal representation requires only the retention of the 1σ and 2σ molecular orbitals, in which the core-excited state involves the promotion of a single electron into the 2σ orbital. This model is extended to include two further representations, in which both the 3σ and 4σ orbitals obtained from a self-consistent field calculation improve target representation, correlation and support additional autoionization channels. This affects the autoionization widths and to a lesser degree the positions of the LiH (1σ2σ 2 n s, n p 1,3 Σ + ) resonance series. Comparing our work with calculations on the counterpart atomic Be system assists in the assignment of the core-excited molecular resonance states. The results from our investigation provide helpful insights into the study of inner-shell transitions produced by electron or photon impact in more complex diatomic molecules.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.