Improved description of chemical barriers with generalized gradient approximations (GGAs) and meta-GGAs

The performance of exchange and correlation (xc) functionals of the generalized gradient approximation (GGA) type and of the meta-GGA type in the calculation of chemical reactions is related to topological features of the electron density which, in turn, are connected to the orbital structure of chemical bonds within the Kohn-Sham (KS) theory. Seventeen GGA and meta-GGA xc functionals are assessed for 15 hydrogen abstraction reactions and 3 symmetrical S(N)2 reactions. Systems that are problematic for standard GGAs characteristically have enhanced values of the dimensionless gradient argument s(sigma)(2) with local maxima in the bonding region. The origin of this topological feature is the occupation of valence KS orbitals with an antibonding or essentially nonbonding character. The local enhancement of s(sigma)(2) yields too negative exchange-correlation energies with standard GGAs for the transition state of the S(N)2 reaction, which leads to the reduced calculated reaction barriers. The unwarranted localization of the effective xc hole of the standard GGAs, i.e., the nondynamical correlation that is built into them but is spurious in this case, wields its effect by their s(sigma)(2) dependence. Barriers are improved for xc functionals with the exchange functional OPTX as x component, which has a modified dependence on s(sigma)(2). Standard GGAs also underestimate the barriers for the hydrogen abstraction reactions. In this case the barriers are improved by correlation functionals, such as the Laplacian-dependent (LAP3) functional, which has a modified dependence on the Coulomb correlation of the opposite- and like-spin electrons. The best overall performance is established for the combination OLAP3 of OPTX and LAP3.

STO and GTO field-induced polarization functions for H to Kr

Field-induced polarization (FIP) functions were proposed over two decades ago to improve the accuracy of calculated response properties, and the FIP functions in GTO form for H and C to F were tested on small molecules, with encouraging results. The concept of FIP,is now extended to all atoms up to Kr. New simplifying approximations for the description of asymptotic highest occupied atomic orbitals. (HOAOs) are introduced in this study. They provide the basis for STO and GTO exponents of a complete set of FIP functions from H to Kr, which are both listed for the convenience of the users. Tests on the polarizabilities of a series of atoms and molecules demonstrate that addition of the FIP basis functions to a series' of standard basis sets drastically improves the performance of all these basis sets compared to converged results. Moreover, the byproduct of this study (approximate asymptotic HOAOs) provides information for the construction of accurate basis sets for long-range ground state properties. (C) 2003 Wiley Periodicals, Inc.

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.

Magnetic properties of copper hexadecaphthalocyanine (F16CuPc) thin films and powders

The structural and magnetic properties of F16CuPc thin films and powder, including x-ray diffraction (XRD), superconducting quantum interference device (SQUID) magnetometry, and theoretical modelling of exchange interactions are reported. Analysis of XRD from films, with thickness ranging between 100 and 160 nm, deposited onto Kapton and a perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride (PTCDA) interlayer shows that the stacking angle (defined in the text) of the film is independent of the thickness, but that the texture is modified by both film thickness and substrate chemistry. The SQUID measurements suggest that all samples are paramagnetic, a result that is confirmed by our theoretical modelling including density functional theory calculations of one-dimensional molecular chains and Green's function perturbation theory calculations for a molecular dimer. By investigating theoretically a range of different geometries, we predict that the maximum possible exchange interaction between F16CuPc molecules is twice as large as that in unfluorinated copper-phthalocyanine (CuPc). This difference arises from the smaller intermolecular spacing in F16CuPc. Our density functional theory calculation for isolated F16CuPc molecule also shows that the energy levels of Kohn-Sham orbitals are rigidly shifted similar to 1 eV lower in F16CuPc compared to CuPc without a significant modification of the intramolecular spin physics, and that therefore the two molecules provide a suitable platform for independently varying magnetism and charge transport. 

Log gf values for astrophysically important transitions Fe II

Aims. In a recent measurement, Meléndez & Barbuy (2009, A&A, 497, 611) report accurate log gf values for 142 important astrophysical lines with wavelengths in the range 4000 Å to 8000 Å. Their results include both solar and laboratory measurements. In this paper, we describe a theoretical study of these lines. Methods. The CIV3 structure codes, combined with our "fine-tuning" extrapolation process, are used to undertake a large-scale CI calculation involving the lowest 262 fine-structure levels belonging to the 3d4s, 3d, 3d4s, 3d4p, and 3d4s4p configurations. Results. We find that many of the 142 transitions are very weak intercombination lines. Other transitions are weak because the dominant configurations in the two levels differ by two orbitals. Conclusions. The comparison between our log gf values and the experimental values generally shows good agreement for most of these transitions, with our theoretical values agreeing slightly more closely with the solar than with the laboratory measurements. A detailed analysis of the small number of transitions for which the agreement between theory and experiment is not as good shows that such disagreements largely arise from severe cancellation due to CI mixing.

SOLVING THE INHOMOGENEOUS SCHRODINGER-EQUATION ON A BIDIRECTIONAL PIPELINE OF TRANSPUTERS

In this paper we describe the design of a parallel solution of the inhomogeneous Schrodinger equation, which arises in the construction of continuum orbitals in the R-matrix theory of atomic continuum processes. A prototype system is described which has been programmed in occam2 and implemented on a bi-directional pipeline of transputers. Some timing results for the prototype system are presented, and the development of a full production system is discussed.

Oscillator strengths and transition probabilities for the W xlv ion

In this paper we present oscillator strengths and transition probabilities for W xlv transitions between levels arising from configurations 3d¹⁰4s²,4p²,4d², 3d¹⁰4k4l (k = s,p,d,f and l = p,d,f), 3d⁹4s²4l (l = p,d,f) and 3d⁹4s4p². The model used to calculate these contained all configurations which can be constructed from the available orbitals (up to n = 4), with either a 3d¹⁰ or 3d⁹ core. The calculations were performed with the configuration interaction CIV3 program with the inclusion of relativistic effects achieved through the use of the Breit-Pauli approximation. We compare our ab initio energy levels, oscillator strengths and transition rates with other experimental and theoretical values available in the literature. There is generally good agreement when only levels with 3d¹⁰ cores are considered. The literature is sparse for levels in which the 3d-subshell is opened: for the majority of the fine-structure lines considered, there is either no comparison data available or substantial differences are found. This paper also investigates how the inclusion of relativistic effects can result in a significant redistribution of the oscillator strength from the LS calculations.

Density functional study of the C atom adsorption on the α-Fe 2O3 (001) surface

The adsorption of C atoms on the α-Fe₂O₃ (001) surface was studied based on density function theory (DFT), in which the exchange-correlation potential was chosen as the PBE (Perdew, Burke and Ernzerhof) generalized gradient approximation (GGA) with a plane wave basis set. Upon the optimization on different adsorption sites with coverage of 1/20 and 1/5 ML, it was found that the adsorption of C atoms on the α-Fe ₂O₃ (001) surface was chemical adsorption. The coverage can affect the adsorption behavior greatly. Under low coverage, the most stable adsorption geometry lied on the bridged site with the adsorption energy of about 3.22 eV; however, under high coverage, it located at the top site with the energy change of 8.79 eV. Strong chemical reaction has occurred between the C and O atoms at this site. The density of states and population analysis showed that the s, p orbitals of C and p orbital of O give the most contribution to the adsorption bonding. During the adsorption process, O atom shares the electrons with C, and C can only affect the outermost and subsurface layers of α-Fe₂O₃; the third layer can not be affected obviously. Copyright © 2008 Chinese Journal of Structural Chemistry.

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.

Site symmetry dependence of repulsive interactions between chemisorbed oxygen atoms on Pt{100}-(1x1)

Ab initio total energy calculations using density functional theory with the generalized gradient approximation have been performed for the chemisorption of oxygen atoms on a Pt{100}-(1 x 1) slab. Binding energies for the adsorption of oxygen on different high-symmetry sites are presented. The bridge site is the most stable at a coverage of 0.5 ML, followed by the fourfold hollow site. The atop site is the least stable. This finding is rationalized by analyzing the ''local structures'' formed upon oxygen chemisorption. The binding energies and heats of adsorption at different oxygen coverages show that pairwise repulsive interactions are considerably stronger between oxygen atoms occupying fourfold sites than those occupying bridge sites. Analysis of the partial charge densities associated with Bloch states demonstrates that the O-Pt bond is considerably more localized at the bridge site. These effects cause a sharp drop in the heats of adsorption for oxygen on hollow sites when the coverage is increased from 0.25 to 0.5 ML. Mixing between oxygen p orbitals and Pt d orbitals can be observed over the whole metal d-band energy range.

Gamma-ray Spectra and Enhancement Factors for Positron Annihilation with Core Electrons

Many-body theory is developed to calculate the γ spectra for positron annihilation in noble-gas atoms. Inclusion of electron-positron correlation effects and core annihilation gives spectra in excellent agreement with experiment [K. Iwata et al., Phys. Rev. Lett. 79, 39 (1997)]. The calculated correlation enhancement factors γnl for individual electron orbitals nl are found to scale with the ionization energy I_nl (in eV), as γnl=1+ √A/Inl+(B/I_nl)β, where A≈40  eV, B≈24  eV, and β≈2.3.

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.

Electron-impact ionization of diatomic molecules using a configuration-average distorted-wave method

Electron-impact ionization cross sections for diatomic molecules are calculated in a configuration-average distorted-wave method. Core bound orbitals for the molecular ion are calculated using a single-configuration self-consistent-field method based on a linear combination of Slater-type orbitals. The core bound orbitals are then transformed onto a two-dimensional (r,θ) numerical lattice from which a Hartree potential with local exchange is constructed. The single-particle Schrödinger equation is then solved for the valence bound orbital and continuum distorted-wave orbitals with S-matrix boundary conditions. Total cross section results for H2 and N2 are compared with those from semiempirical calculations and experimental measurements.

Synthesis, characterization, electrochemical behavior and in vitro protein tyrosine kinase inhibitory activity of the cymene-halogenobenzohydroxamato [Ru(eta(6)-cymene)(bha)Cl] complexes

The ruthenium(II)-cymene complexes [Ru(eta(6)-cymene)(bha)Cl] with substituted halogenobenzohydroxamato (bha) ligands (substituents = 4-F, 4-Cl, 4-Br, 2,4-F-2, 3,4-F-2, 2,5-F-2, 2,6-F-2) have been synthesized and characterized by elemental analysis, IR, H-1 NMR, C-13 NMR, cyclic voltammetry and controlled-potential electrolysis, and density functional theory (DFT) studies. The compositions of their frontier molecular orbitals (MOs) were established by DFT calculations, and the oxidation and reduction potentials are shown to follow the orders of the estimated vertical ionization potential and electron affinity, respectively. The electrochemical E-L Lever parameter is estimated for the first time for the various bha ligands, which can thus be ordered according to their electron-donor character. All complexes exhibit very strong protein tyrosine kinase (PTK) inhibitory activity, even much higher than that of genistein, the clinically used PTK inhibitory drug. The complex containing the 2,4-difluorobenzohydroxamato ligand is the most active one, and the dependences of the PTK activity of the complexes and of their redox potentials on the ring substituents are discussed. (C) 2012 Elsevier B.V. All rights reserved.