Exchange-correlation energy and potential as approximate functionals of occupied and virtual Kohn-Sham orbitals: Application to dissociating H-2

The standard local density approximation and generalized gradient approximations fail to properly describe the dissociation of an electron pair bond, yielding large errors (on the order of 50 kcal/mol) at long bond distances. To remedy this failure, a self-consistent Kohn-Sham (KS) method is proposed with the exchange-correlation (xc) energy and potential depending on both occupied and virtual KS orbitals. The xc energy functional of Buijse and Baerends [Mol. Phys. 100, 401 (2002); Phys. Rev. Lett. 87, 133004 (2001)] is employed, which, based on an ansatz for the xc-hole amplitude, is able to reproduce the important dynamical and nondynamical effects of Coulomb correlation through the efficient use of virtual orbitals. Self-consistent calculations require the corresponding xc potential to be obtained, to which end the optimized effective potential (OEP) method is used within the common energy denominator approximation for the static orbital Green's function. The problem of the asymptotic divergence of the xc potential of the OEP when a finite number of virtual orbitals is used is addressed. The self-consistent calculations reproduce very well the entire H-2 potential curve, describing correctly the gradual buildup of strong left-right correlation in stretched H-2. (C) 2003 American Institute of Physics.

Positronium scattering by lithium

We report cross sections for Ps(1s)-Li(2s) scattering in the energy range up to 30 eV. The calculations have been carried out in a coupled state approximation. The Ps states consist of both eigenstates and pseudostates. the latter to allow for ionization of the Ps. The atom is treated as a frozen core represented by it model potential which supports the valence orbitals. The coupled state expansion includes only the 2s and 2p states of the atom as well as in unphysical Is state which exists in the model potential. The inclusion of this Is state is necessary in order to avoid pronounced false pseudostructure. Results are presented for excitation and ionization of the Ps as well as collisions in which the Ps(1s) remains unchanged. These results also differentiate between the case where the Li(2s) remains unexcited and where it is excited to the 2p level. (c) 2005 Published by Elsevier B.V.

Energy levels and transition probabilities for nitrogen-like Fe xx

Energies of the 700 lowest levels in Fe XX have been obtained using the multiconfiguration Dirac-Fock method. Configuration interaction method on the basis set of transformed radial orbitals with variable parameters taking into account relativistic corrections in the Breit-Pauli approximation was used to crosscheck our presented results. Transition probabilities, oscillator and line strengths are presented for electric dipole (E1), electric quadrupole (E2) and magnetic dipole (M1) transitions among these levels. The total radiative transition probabilities from each level are also provided. Results are compared with data compiled by NIST and with other theoretical work.

Electron recombination with multicharged ions via chaotic many-electron states

We show that a dense spectrum of chaotic multiply excited eigenstates can play a major role in collision processes involving many-electron multicharged ions. A statistical theory based on chaotic properties of the eigenstates enables one to obtain relevant energy-averaged cross sections in terms of sums over single-electron orbitals. Our calculation of low-energy electron recombination of Au25+ shows that the resonant process is 200 times more intense than direct radiative recombination, which explains the recent experimental results of Hoffknecht [J. Phys. B 31, 2415 (1998)].

Virtual state scattering with cold electrons: para-xylene and para-difluorobenzene

The scattering of electrons with kinetic energies down to a few meV by para-xylene and para-difluorobenzene has been observed experimentally with an electron beam energy resolution of 0.95 to 1.5 meV (full width half maximum). At low electron energies the collisions can be considered as cold scattering events because the de Broglie wavelength of the electron is considerably larger than the target dimensions. The scattering cross sections measured rise rapidly at low energy due to virtual state scattering. The nature of this scattering process is discussed using s- and p-wave phase shifts derived from the experimental data. Scattering lengths are derived of, respectively, -9.5+/-0.5 and -8.0+/-0.5 a.u. for para-xylene and para-difluorobenzene. The virtual state effect is interpreted in terms of nuclear diabatic and partially adiabatic models, involving the electronic and vibronic symmetries of the unoccupied orbitals in the target species. The concept of direct and indirect virtual state scattering is introduced, through which the present species, in common with carbon dioxide and benzene, scatter through an indirect virtual state process, whereas other species, such as perfluorobenzene, scatter through a direct process. (C) 2005 American Institute of Physics.

Dissociative electron attachment to HCCCN

Dissociative electron attachment to cyanoacetylene (propiolonitrile) HCCCN has been observed in the electron energy range 0-12 eV. Negative ions are formed in two main bands with maxima at similar to 1.6 eV (CCCN-) and similar to 5.3 eV (CCCN-, CN-, HCC- and CC-). There are also weaker resonances which lead to dissociative electron attachment to form CN-, HCC- and CC- with a maximum intensity at similar to 8.1 eV and CCCN-, CN- and CC- at similar to 11.2 eV. A trace of CCN- is observed at similar to 9.1 eV. The positions of the main dissociative attachment bands observed are close to positions of pi* resonances recently calculated by Sommerfeld and Knecht. Calculations have also been performed in this work, which confirm the position of the p* orbitals. The electron affinity of the CCCN radical is determined as 4.59 +/- 0.25 eV from the threshold for CCCN- formation at 1.32 +/- 0.15 eV. Dissociative electron attachment to this molecule will act as a source of negative ions in extraterrestrial environments where electrons are present with more than 1.3 eV energy.

Tools for analysing configuration interaction wavefunctions

The configuration interaction (CI) approach to quantum chemical calculations is a well-established means of calculating accurately the solution to the Schrodinger equation for many-electron systems. It represents the many-body electron wavefunction as a sum of spin-projected Slater determinants of orthogonal one-body spin-orbitals. The CI wavefunction becomes the exact solution of the Schrodinger equation as the length of the expansion becomes infinite, however, it is a difficult quantity to visualise and analyse for many-electron problems. We describe a method for efficiently calculating the spin-averaged one- and two-body reduced density matrices rho(psi)((r) over bar; (r) over bar' ) and Gamma(psi)((r) over bar (1), (r) over bar (2); (r) over bar'(1), (r) over bar'(2)) of an arbitrary CI wavefunction Psi. These low-dimensional functions are helpful tools for analysing many-body wavefunctions; we illustrate this for the case of the electron-electron cusp. From rho and Gamma one can calculate the matrix elements of any one- or two-body spin-free operator (O) over cap. For example, if (O) over cap is an applied electric field, this field can be included into the CI Hamiltonian and polarisation or gating effects may be studied for finite electron systems. (C) 2003 Elsevier B.V. All rights reserved.

Optimal basis set for electronic structure calculations in periodic systems

An efficient method for calculating the electronic structure of systems that need a very fine sampling of the Brillouin zone is presented. The method is based on the variational optimization of a single (i.e., common to all points in the Brillouin zone) basis set for the expansion of the electronic orbitals. Considerations from k.p-approximation theory help to understand the efficiency of the method. The accuracy and the convergence properties of the method as a function of the optimal basis set size are analyzed for a test calculation on a 16-atom Na supercell.

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.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Spontaneous spin polarization and electron localization in constrained geometries: The Wigner transition in nanowires

The Wigner transition in a jellium model of cylindrical nanowires has been investigated by density-functional computations using the local spin-density approximation. A wide range of background densities rho(b) has been explored from the nearly ideal metallic regime (r(s)=[3/4 pi rho(b)](1/3)=1) to the high correlation limit (r(s)=100). Computations have been performed using an unconstrained plane wave expansion for the Kohn-Sham orbitals and a large simulation cell with up to 480 electrons. The electron and spin distributions retain the cylindrical symmetry of the Hamiltonian at high density, while electron localization and spin polarization arise nearly simultaneously in low-density wires (r(s)similar to 30). At sufficiently low density (r(s)>= 40), the ground-state electron distribution is the superposition of well defined and nearly disjoint droplets, whose charge and spin densities integrate almost exactly to one electron and 1/2 mu(B), respectively. Droplets are arranged on radial shells and define a distorted lattice whose structure is intermediate between bcc and fcc. Dislocations and grain boundaries are apparent in the droplets' configuration found by our simulations. Our computations aim at modeling the behavior of experimental low-carried density systems made of lightly doped semiconductor nanostructures or conducting polymers.

Choline saccharinate and choline acesulfamate: Ionic liquids with low toxicities

Choline saccharinate and choline acesulfamate are two examples of hydrophilic ionic liquids, which can be prepared from easily available starting materials (choline chloride and a non-nutritive sweetener). The (eco)toxicity of these ionic liquids in aqueous solution is very low in comparison to other types of ionic liquids. A general method for the synthesis and purification of hydrophilic ionic liquids is presented. The method consists of a silver-free metathesis reaction, followed by purification of the ionic liquid by ion-exchange chromatography. The crystal structures show a marked difference in hydrogen bonding between the two ionic liquids, although the saccharinate and the acesulfamate anions show structural similarities. The optimized structures, the energetics, and the charge distribution of cation-anion pairs in the ionic liquids were studied by density functional theory (DFT) and second-order (Moller-Plesset) perturbation theory calculations. The occupation of the non-Lewis orbitals was considered to obtain a qualitative picture of the Lewis structures. The calculated interaction energies and the dipole moments for the ion pairs in the gas phase were discussed.

Convergence of partial-wave expansions for energies, scattering amplitudes and positron annihilation rates

We use many-body theory to find the asymptotic behaviour of second-order correlation corrections to the energies and positron annihilation rates in many- electron systems with respect to the angular momenta l of the single-particle orbitals included. The energy corrections decrease as 1/(l+1/2)4, in agreement with the result of Schwartz, whereas the positron annihilation rate has a slower 1/(l+1/2)2 convergence rate. We illustrate these results by numerical calculations of the energies of Ne and Kr and by examining results from extensive con?guration-interaction calculations of PsH binding and annihilation.

Long-range electronic coupling in bis(cyclometalated) ruthenium complexes

Symmetrical and unsymmetrical ligands containing terpyridyl coordinating units (N, N, N) or a cyclometalating equivalent (N, C, N), connected back-to-back either directly or via a p-terphenylene or 1,3-phenylene spacer, have been used to construct new diruthenium complexes. These compounds incorporate various terdentate chelates as capping ligands, to allow a double control of the electronic properties of each subcomplex and of the ensemble: via the terminal ligand or through the bridging fragment. Electronic coupling was studied from the intervalence transitions observed in several bimetallic ruthenium complexes of the bis-(cyclometalated) type differing by the substitution of a nitrogen atom by carbon in the terminal terpyridyl unit. The largest metal-metal interaction was found in complexes for which the terminal complexing unit is of the 1,3-di-2-pyridylbenzene type, i.e., with the carbon atom located on the metal-metal C-2 axis of the molecule. Investigations of the mechanism of interaction by extended Huckel calculations showed that the replacement of nitrogen by carbon raises the filled ligand levels, increasing the mixing with ligand orbitals and thus the metal-metal coupling. Finally, the intervalence transition was still observed for a bridging ligand containing three phenylene units as spacers, corresponding to a 24-Angstrom metal-metal distance.