Magnetic structure of Li2CuO2: From ab initio calculations to macroscopic simulations

The magnetic structure of the edge-sharing cuprate compound Li2CuO2 has been investigated with highly correlated ab initio electronic structure calculations. The first- and second-neighbor in-chain magnetic interactions are calculated to be 142 and -22 K, respectively. The ratio between the two parameters is smaller than suggested previously in the literature. The interchain interactions are antiferromagnetic in nature and of the order of a few K only. Monte Carlo simulations using the ab initio parameters to define the spin model Hamiltonian result in a Nel temperature in good agreement with experiment. Spin population analysis situates the magnetic moment on the copper and oxygen ions between the completely localized picture derived from experiment and the more delocalized picture based on local-density calculations.

Half metallicity in a zigzag double-walled nanotube nanodot : an Ab initio prediction

Ab initio density functional calculations were performed to study finite-length zigzag (7, 0) @ (16, 0) double-walled carbon nanotubes (DWCNTs) with H-termination at the open ends. We find that such a DWCNT nanodot displays a very large magnetic moment at the zigzag edges and the ground state displays symmetric anti-ferromagnetic coupling. When an external electric field is applied along the direction of tube axis, a gap is opened for one spin channel, whereas another spin channel remains metallic, i.e. half metallicity occurs. Our results suggest an important new avenue for the development of CNT-based spintronic materials with enhanced properties.

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

Role of lithium vacancies in accelerating the dehydrogenation kinetics on a LiBH4(010) surface : an Ab initio study

In this work, ab initio density functional calculations were performed to explore the effect of surface lithium vacancies on the initial dehydrogenation kinetics of lithium borohydride. We found that some B−H bonds in neighboring BH4-1 complexes around the vacancy became elongated (weakened). The activation barriers for the recombination of H atoms to form H2 were decreased from 3.64 eV for the stoichiometrically complete LiBH4(010) surface to 1.53 and 0.23 eV in the presence of mono- and di-vacancies, respectively. Our results indicate that the creation of Li vacancies may play a critical role in accelerating the dehydrogenation kinetics of LiBH4.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(0 0 1) and MgH2(1 1 0), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(1 1 0) surface is more stable than MgH2(0 0 1) surface, which is in good agreement with the experimental observation. The H2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved – they are found to be generally high, due to the thermodynamic stability of the MgH2 system, and are larger for the MgH2(0 0 1) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(1 1 0) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates.

Ground and Excited State Intramolecular Proton Transfer in Salicylic Acid: an Ab Initio Electronic Structure Investigation

Energetics of the ground and excited state intramolecular proton transfer in salicylic acid have been studied by ab initio molecular orbital calculations using the 6-31G** basis set at the restricted Hartree-Fock (RHF) and configuration interaction-single excitation (CIS) levels and also using the semiempirical method AM1 at the RHF level as well as with single and pair doubles excitation configuration interaction spanning eight frontier orbitals (PECI = 8). The ab initio potential energy profile for intramolecular proton transfer in the ground state reveals a single minimum corresponding to the primary form, in the first excited singlet state, however, there are two minima corresponding to the primary and tautomeric forms, separated by a barrier of similar to 6 kcal/mol, thus accounting for dual emission in salicylic acid. Electron density changes with electronic excitation and tautomerism indicate no zwitterion formation. Changes in spectral characteristics with change in pH, due to protonation and deprotonation of salicylic acid, are also accounted for, qualitatively. Although the AM1 calculations suggest a substantial barrier for proton transfer in the ground as well as the first excited state of SA, it predicts the transition wavelength in near quantitative accord with the experimental results for salicylic acid and its protonated and deprotonated forms.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Preparation and properties of gallaborane, GaBH6: structure of the gaseous molecule H2Ga(μ-H)2BH2 as determined by vibrational, electron diffraction, and ab initio studies, and structure of the crystalline solid at 110 K as determined by x-ray diffraction``

Gallaborane (GaBH6, 1), synthesized by the metathesis of LiBH4 with [H2GaCl]n at ca. 250 K, has been characterized by chemical analysis and by its IR and 1H and 11B NMR spectra. The IR spectrum of the vapor at low pressure implies the presence of only one species, viz. H2Ga(μ-H)2BH2, with a diborane-like structure conforming to C2v symmetry. The structure of this molecule has been determined by gas-phase electron diffraction (GED) measurements afforced by the results of ab initio molecular orbital calculations. Hence the principal distances (rα in Å) and angles ( α in deg) are as follows: r(Ga•••B), 2.197(3); r(Ga−Ht), 1.555(6); r(Ga−Hb), 1.800(6); r(B−Ht), 1.189(7); r(B−Hb), 1.286(7); Hb−Ga−Hb, 71.6(4); and Hb−B−Hb, 110.0(5) (t = terminal, b = bridging). Aggregation of the molecules occurs in the condensed phases. X-ray crystallographic studies of a single crystal at 110 K reveal a polymeric network with helical chains made up of alternating pseudotetrahedral GaH4 and BH4 units linked through single hydrogen bridges; the average Ga•••B distance is now 2.473(7) Å. The compound decomposes in the condensed phases at temperatures exceeding ca. 240 K with the formation of elemental Ga and H2 and B2H6. The reactions with NH3, Me3N, and Me3P are also described.

LDA or GGA? A combined experimental inelastic neutron scattering and ab initio lattice dynamics study of alkali metal hydrides

In a previous work, we carried out inelastic neutron scattering (INS) spectroscopy experiments and preliminary first principles calculations on alkali metal hydrides. The complete series of alkali metal hydrides, LiH, NaH, KH, RbH and CsH was measured in the high-resolution TOSCA INS spectrometer at ISIS. Here, we present the results of ab initio electronic structure calculations of the properties of the alkali metal hydrides using both the local density approximation (LDA) and the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) parameterization. Properties calculated were lattice parameters, bulk moduli, dielectric constants, effective charges, electronic densities and inelastic neutron scattering (INS) spectra. We took advantage of the currently available computer power to use full lattice dynamics theory to calculate thermodynamic properties for these materials. For the alkali metal hydrides (LiH, NaH, KH, RbH and CsH) using lattice dynamics, we found that the INS spectra calculated using LDA agreed better with the experimental data than the spectra calculated using GGA. Both zero-point effects and thermal contributions to free energies had an important effect on INS and several thermodynamic properties.

Hydrogen adsorption on boron doped graphene: an ab initio study

(i) The electronic and structural properties of boron doped graphene sheets, and (ii) the chemisorption processes of hydrogen adatoms on the boron doped graphene sheets have been examined by ab initio total energy calculations. In (i) we find that the structural deformations are very localized around the boron substitutional sites, and in accordance with previous studies (Endo et al 2001 J. Appl. Phys. 90 5670) there is an increase of the electronic density of states near the Fermi level. Our simulated scanning tunneling microscope (STM) images, for occupied states, indicate the formation of bright (triangular) spots lying on the substitutional boron (center) and nearest-neighbor carbon (edge) sites. Those STM images are attributed to the increase of the density of states within an energy interval of 0.5 eV below the Fermi level. For a boron concentration of similar to 2.4%, we find that two boron atoms lying on the opposite sites of the same hexagonal ring (B1-B2 configuration) represents the energetically most stable configuration, which is in contrast with previous theoretical findings. Having determined the energetically most stable configuration for substitutional boron atoms on graphene sheets, we next considered the hydrogen adsorption process as a function of the boron concentration, (ii). Our calculated binding energies indicate that the C-H bonds are strengthened near boron substitutional sites. Indeed, the binding energy of hydrogen adatoms forming a dimer-like structure on the boron doped B1-B2 graphene sheet is higher than the binding energy of an isolated H(2) molecule. Since the formation of the H dimer-like structure may represent the initial stage of the hydrogen clustering process on graphene sheets, we can infer that the formation of H clusters is quite likely not only on clean graphene sheets, which is in consonance with previous studies (Hornekaer et al 2006 Phys. Rev. Lett. 97 186102), but also on B1-B2 boron doped graphene sheets. However, for a low concentration of boron atoms, the formation of H dimer structures is not expected to occur near a single substitutional boron site. That is, the formation (or not) of H clusters on graphene sheets can be tuned by the concentration of substitutional boron atoms.

Estudo da adsorção de hidrogênio e sulfeto na superfície de paládio: aspectos experimentais (eletroquímica) e teóricos (ab initio e Teoria do Funcional da Densidade)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

Ab initio studies of hydrogen desorption from low index magnesium hydride surface

The low index Magnesium hydride surfaces, MgH2(001) and MgH2(110), have been studied by ab intio Density Functional Theory (DFT) calculations. It was found that the MgH2(110) surface is more stable than MgH2(001) surface, which is in good agreement with the experimental observation. The H-2 desorption barriers vary depending on the crystalline surfaces that are exposed and also the specific H atom sites involved-they are found to be generally high, due to the thermodynamic stability of the MgH2, system, and are larger for the MgH2(001) surface. The pathway for recombinative desorption of one in-plane and one bridging H atom from the MgH2(110) surface was found to be the lowest energy barrier amongst those computed (172 KJ/mol) and is in good agreement with the experimental estimates. (c) 2006 Elsevier B.V. All rights reserved.

Thiothionyl bromide (S=SBr2): The elusive isomer of dibromodisulfane (BrSSBr)

CCSD(T) with a series of correlation consistent basis up to quadruple-zeta is used to investigate the structures, vibrational spectra, relative stability, heats of formation, and barrier to isomerization of S=SBr2 and BrSSBr. It represents the most accurate and detailed characterization of these molecules to date. We show that the frequency mode at 302 cm(-1), detected in various studies and assigned to impurities by some authors, and to the anti-symmetric SBr stretch in BrSSBr by others, thus in fact corresponds to the anti-symmetric SBr stretch in the elusive S=SBr2 species; it thus corroborates and complements an earlier partial IR spectra study attributable to S=SBr2. Including corrections for relativistic and core-valence correlation effects, we also predict 26.33 (12.74) kcal/mol for Delta H-f (298.15 K) of S=SBr2 (BrSSBr). For the S=SBr2 -> BrSSBr reaction, our best estimates for the Gibbs free energy and the enthalpy of the reaction at 298.15 K are -13.71 and -13.44 kcal/mol, respectively. For a value of Delta G(#) equal to 23.52 kcal/mol, we estimate a TST rate constant, at 298.15 K, of 3.57 x 10(-5) s(-1). (c) 2007 Elsevier B.V. All rights reserved.

Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.