Nonadiabatic forces in ion-solid interactions:The initial stages of radiation damage

The Born-Oppenheimer approximation is the keystone for molecular dynamics simulations of radiation damage processes; however, actual materials response involves nonadiabatic energy exchange between nuclei and electrons. In this work, time dependent density functional theory is used to calculate the electronic excitations produced by energetic protons in Al. We study the influence of these electronic excitations on the interatomic forces and find that they differ substantially from the adiabatic case, revealing a nontrivial connection between electronic and nuclear stopping that is absent in the adiabatic case. These results unveil new effects in the early stages of radiation damage cascades.

Electronic stopping power in gold:The role of d electrons and the H/He anomaly

The electronic stopping power of H and He moving through gold is obtained to high accuracy using time-evolving density-functional theory, thereby bringing usual first principles accuracies into this kind of strongly coupled, continuum nonadiabatic processes in condensed matter. The two key unexplained features of what observed experimentally have been reproduced and understood: (i)The nonlinear behavior of stopping power versus velocity is a gradual crossover as excitations tail into the d-electron spectrum; and (ii)the low-velocity H/He anomaly (the relative stopping powers are contrary to established theory) is explained by the substantial involvement of the d electrons in the screening of the projectile even at the lowest velocities where the energy loss is generated by s-like electron-hole pair formation only. © 2012 American Physical Society.

On the Mechanism of Lewis Acid Catalyzed Glucose Transformations in Ionic Liquids

A complementary computational and experimental study of the reactivity of Lewis acidic CrCl2, CuCl2 and FeCl2 catalysts towards glucose activation in dialkylimidazolium chloride ionic liquids is performed. The selective dehydration of glucose to 5-hydroxymethylfurfural (HMF) proceeds through the intermediate formation of fructose. Although chromium(II) and copper(II) chlorides are able to dehydrate fructose with high HMF selectivity, reasonable HMF yields from glucose are only obtained with CrCl2 as the catalyst. Glucose conversion by CuCl2 is not selective, while FeCl2 catalyst does not activate sugar molecules. These differences in reactivity are rationalized on the basis of in situ X-ray absorption spectroscopy measurements and the results of density functional theory calculations. The reactivity in glucose dehydration and HMF selectivity are determined by the behavior of the ionic liquid-mediated Lewis acid catalysts towards the initial activation of the sugar molecules. The formation of a coordination complex between the Lewis acidic Cr2+ center and glucose directs glucose transformation into fructose. For Cu2+ the direct coordination of sugar to the copper(II) chloride complex is unfavorable. Glucose deprotonation by a mobile Cl- ligand in the CuCl42- complex initiates the nonselective conversion. In the course of the reaction the Cu2+ ions are reduced to Cu+. Both paths are prohibited for the FeCl2 catalyst.

Coordination Properties of Ionic Liquid-Mediated Chromium(II) and Copper(II) Chlorides and Their Complexes with Glucose

The structural and coordination properties of complexes formed upon the interaction of copper(II) and chromium(II) chlorides with diallrylimidazolium chloride (RMlm(+)Cl(-)) ionic liquids and glucose are studied by a combination of density functional theory (DFT) calculations and X-ray absorption spectroscopy (XAS). In the absence of the carbohydrate substrate, isolated mononuclear four-coordinated MeCl42- species (Me = Cu, Cr) dominate in the ionic liquid solution. The organic part of the ionic liquid does not directly interact with the metal centers. The interactions between the RMlm(+) cations and the anionic metal chloride complexes are limited to hydrogen bonding with the basic Cl- ligands and the overall electrostatic stabilization of the anionic metal complexes. Exchange of Cl ligands by a hydroxyl group of glucose is only favorable for CrCl42-. For Cu2+ complexes, the formation of hydrogen bonded complexes between CuCl42- and glucose is preferred. No preference for the coordination of metal chloride species to specific hydroxyl group of the carbohydrate is found. The formation of binuclear metal chloride complexes is also considered. The reactivity and selectivity patterns of the Lewis acid catalyzed reactions of glucose are discussed in the framework of the obtained results.

Selective Hydrogenation of alpha,beta-Unsaturated Aldehydes and Ketones using Novel Manganese Oxide and Platinum Supported on Manganese Oxide Octahedral Molecular Sieves as Catalysts

The selective hydrogenation of ,-unsaturated aldehydes and ketones has been studied using ketoisophorone and cinnamaldehyde as model substrates using manganese oxide octahedral molecular sieve (OMS-2) based catalysts. For the first time, OMS-2 has been shown to be an efficient and selective hydrogenation catalyst. High selectivities for either the CC or CO double bond at approximate to 100% conversion were achieved by using OMS-2 and platinum supported on OMS-2 catalysts. Density functional theory (DFT) calculations showed that the dissociation of H2 on OMS-2 was water assisted and occurred on the surface Mn of OMS-2(001) that had been modified by an adsorbed H2O molecule. The theoretically calculated activation barrier was in good agreement with the experimentally determined value for the hydrogenation reactions, indicating that H2 dissociation on OMS-2 is likely to be the rate-determining step. A significant increase in the rate of reaction was observed in the presence of Pt as a result of the enhancement of H2 dissociative adsorption and subsequent reaction on the Pt or spillover of the hydrogen to the OMS-2 support. The relative adsorption strengths of ketoisophorone and cinnamaldehyde on the OMS-2 support compared with the Pt were found to determine the product selectivity.

Origin of the Increase of Activity and Selectivity of Nickel Doped by Au, Ag, and Cu for Acetylene Hydrogenation

Activity and selectivity are both important issues in heterogeneous catalysis and recent experimental results have shown that Ni catalysts doped by gold exhibit high activity for the hydrogenation of acetylene with good selectivity of ethylene formation. To unravel the underlying mechanism for this observation, the general trend of activity and selectivity of Ni surfaces doped by Au, Ag, and Cu has been investigated using density functional theory calculations. Complete energy profiles from C2H2 to C2H4 on Ni(111), Au/Ni(111), Ag/Ni(111) and Cu/Ni(111) are obtained and their turnover frequencies (TOFs) are computed. The results show that acetylene adsorption on Ni catalyst is strong which leads to the low activity while the doping of Au, Ag, and Cu on the Ni catalyst weakens the acetylene adsorption, giving rise to the increase of activity. The selectivity of ethylene formation is also quantified by using the energy difference between the hydrogenation barriers and the absolute value of the adsorption energies of ethylene. It is found that the selectivity of ethylene formation increases by doping Au and Ag, while those of Cu/Ni and Ni are similar.

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.

yambo: An ab initio tool for excited state calculations

yambo is an ab initio code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the GW approximation for the self-energy. Optical properties are evaluated either by solving the Bethe-Salpeter equation or by using the adiabatic local density approximation. yambo is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. yambo relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible 110 procedures and is interfaced to several publicly available density functional ground-state codes.

Implementation and testing of Lanczos-based algorithms for Random-Phase Approximation eigenproblems

The treatment of the Random-Phase Approximation Hamiltonians, encountered in different frameworks, like time-dependent density functional theory or Bethe-Salpeter equation, is complicated by their non-Hermicity. Compared to their Hermitian Hamiltonian counterparts, computational methods for the treatment of non-Hermitian Hamiltonians are often less efficient and less stable, sometimes leading to the breakdown of the method. Recently [Gruning et al. Nano Lett. 8 (2009) 28201, we have identified that such Hamiltonians are usually pseudo-Hermitian. Exploiting this property, we have implemented an algorithm of the Lanczos type for Random-Phase Approximation Hamiltonians that benefits from the same stability and computational load as its Hermitian counterpart, and applied it to the study of the optical response of carbon nanotubes. We present here the related theoretical grounds and technical details, and study the performance of the algorithm for the calculation of the optical absorption of a molecule within the Bethe-Salpeter equation framework. (C) 2011 Elsevier B.V. All rights reserved.

Quasiparticle calculations of the electronic properties of ZrO2 and HfO2 polymorphs and their interface with Si

Quasiparticle calculations are performed to investigate the electronic band structures of various polymorphs of Hf and Zr oxides. The corrections with respect to density-functional-theory results are found to depend only weakly on the crystal structure. Based on these bulk calculations as well as those for bulk Si, the effect of quasiparticle corrections is also investigated for the band offsets at the interface between these oxides and Si assuming that the lineup of the potential at the interface is reproduced correctly within density-functional theory. On the one hand, the valence-band offsets are practically unchanged with a correction of a few tenths of electron volts. On the other hand, conduction-band offsets are raised by 1.3-1.5 eV. When applied to existing calculations for the offsets at the density-functional-theory level, our quasiparticle corrections provide results in good agreement with the experiment.

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

The failure of generalized gradient approximations (GGAs) and meta-GGAs for the two-center three-electron bonds in He-2(+), (H2O)(2)(+), and (NH3)(2)(+)

The radical cations He-2(+) (H2O)(2)(+), and (NH3)(2)(+) with two-center three-electron A-A bonds are investigated at the configuration interaction (CI), accurate Kohn-Sham (KS), generalized gradient approximation (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A(2)(+) systems remain a difficult case for density functional theory (DFT). All methods tested consistently overestimate the stability of A(2)(+): the corresponding D-e errors decrease for more diffuse valence densities in the series He-2(+) > (H2O)(2)(+) > (NH3)(2)(+). Upon comparison to the energy terms of the accurate Kohn-Sham solutions, the approximate exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A(2)(+). These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe nondynamical correlation. However, in the hemibonded A(2)(+) systems the typical molecular (left-right, nondynamical) correlation of the two-electron pair bond is absent. The nondynamical correlation built into the exchange functionals is then spurious and yields too low energies.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 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. 

Crystal and magnetic structure of Mn3IrSi

A new ternary Ir-Mn-Si phase with stoichiometry Mn₃IrSi has been synthesized and found to crystallize in the cubic AlAu₄-type structure, space group P2₁3 with Z=4, which is an ordered form of the beta-Mn structure. The unit cell dimension was determined by x-ray powder diffraction to a=6.4973(3) Angstrom. In addition to the crystal structure, we have determined the magnetic structure and properties using superconducting quantum interference device magnetometry and Rietveld refinements of neutron powder diffraction data. A complex noncollinear magnetic structure is found, with magnetic moments of 2.97(4)u(B) at 10 K only on the Mn atoms. The crystal structure consists of a triangular network built up by Mn atoms, on which the moments are rotated 120degrees around the triangle axes. The magnetic unit cell is the same as the crystallographic and carries no net magnetic moment. The Neel temperature was determined to be 210 K. A first-principles study, based on density functional theory in a general noncollinear formulation, reproduces the experimental results with good agreement. The observed magnetic structure is argued to be the result of frustration of antiferromagnetic couplings by the triangular geometry.