Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.

Crystal and molecular structures of three 2-sulphur-substituted cyclohexanones studied by X-ray crystallography and by ab initio molecular orbital calculations

The conformational features of three 2-sulphur-substituted cyclohexanone derivatives, which differ in the number of sulphur-bound oxygen atoms, i.e. zero (I), one (II) and two (III), were investigated by single crystal X-ray crystallography and geometry optimized structures determined using Hartree-Fock method. In each of (I)-(III) an intramolecular S center dot center dot center dot O(carbonyl) interaction is found with the magnitude correlated with the oxidation state of the sulphur atom, i.e. 2.838(3) angstrom in (I) to 2.924(2) angstrom in (II) to 3.0973(18) angstrom in (III). There is an inverse relationship between the strength of this interaction and the magnitude of the carbonyl bond. The supramolecular aggregation patterns are primarily determined by C-H center dot center dot center dot O contacts and are similarly influenced by the number of oxygen atoms in the molecular structures. Thus, a supramolecular chain is found in the crystal structure of (I). With an additional oxygen atom available to participate in C-H center dot center dot center dot O interactions, as in (II), a two-dimensional array is found. Finally, a three-dimensional network is found for (III). Despite there being differences in conformations between the experimental structures and those calculated in the gas-phase, the S center dot center dot center dot O interactions persist. The presence of intermolecular C-H center dot center dot center dot O interactions involving the cyclohexanone-carbonyl group in the solid-state, disrupts the stabilising intramolecular C-H center dot center dot center dot O interaction in the energetically-favoured conformation. (I): C(12)H(13)NO(3)S, triclinic space group P (1) over bar with a = 5.392(3) angstrom b = 10.731(6) angstrom, c = 11.075(6) angstrom, alpha = 113.424(4)degrees, beta = 94.167(9)degrees, gamma = 98.444(6)degrees, V = 575.5(6) angstrom(3), Z = 2, R(1) = 0.052; (II): C(12)H(13)NO(4)S, monoclinic P2(1)/n, a = 7.3506(15) angstrom, b = 6.7814(14) angstrom, c = 23.479(5) angstrom, beta = 92.94(3)degrees, V = 1168.8(4) angstrom(3), Z = 4, R(1) = 0.046; (III): C(12)H(13)NO(5)S, monoclinic P2(1)/c, a = 5.5491(11) angstrom, b = 24.146(3) angstrom, c = 11.124(3) angstrom, beta = 114.590(10)degrees, V = 1355.3(5) angstrom(3), Z = 4, R(1) = 0.051.

Ab Initio correlated all electron dirac-fock calculations for Eka-Francium Fluoride (E119F)

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

Ab Initio correlated all electron Dirac-Fock calculations for eka-francium fluoride (E119F)

Resultados obtidos por cálculos Dirac-Fock correlacionados de 4 componentes para o fluoreto do elemento E119 (Eka-Frâncio) com base estável e precisa, livre de prolapso variacional, são reportados neste trabalho. No nível CCSD(T), a distância de equilíbrio R_e, frequência harmônica ω_e e energia de dissociação D_e são 2,432 Å, 354,97 cm_-1 e 116,92 kcal mol_-1, respectivamente. Também são reportados base livre de prolapso variacional de 4 componentes para o elemento 119, uma curva analítica de energia potencial precisa e o espectro vibracional a partir dos dados obtidos no nível CCSD(T). Nossos resultados sugerem que a molécula E119F deva ser menos iônica que seus fluoretos alcalinos homólogos mais leves, em contraste com o senso químico comum baseado nas propriedades periódicas - era de se esperar nesta molécula a ligação química mais iônica possível. Também encontramos que a correção do tipo modelo de carga para negligenciar as integrais do tipo SS resulta em erros insignificantes e acelera os cálculos cerca de 3 vezes no nível CCSD(T) e cerca de 4 vezes no nível DFT/B3LYP.

Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

Hydrogen spillover mechanism on a pd-doped mg surface as revealed by ab initio density functional calculation

The hydrogenation kinetics of Mg is slow, impeding its application for mobile hydrogen storage. We demonstrate by ab initio density functional theory (DFT) calculations that the reaction path can be greatly modified by adding transition metal catalysts. Contrasting with Ti doping, a Pd dopant will result in a very small activation barrier for both dissociation of molecular hydrogen and diffusion of atomic H on the Mg surface. This new computational finding supports for the first time by ab initio simulationthe proposed hydrogen spillover mechanism for rationalizing experimentally observed fast hydrogenation kinetics for Pd-capped Mg materials.

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.

The location of Ti atom in sodium alanate : an ab initio spin-polarised study

In this work, ab initio spin-polarised Density Functional Theory (DFT) calculations are performed to study the interaction of a Ti atom with a NaAlH4(001) surface. We confirm that an interstitially located Ti atom in the NaAlH4 subsurface is the most energetically favoured configuration as recently reported (Chem. Comm. (17) 2006, 1822). On the NaAlH4(001) surface, the Ti atom is most stable when adsorbed between two sodium atoms with an AlH4 unit beneath. A Ti atom on top of an Al atom is also found to be an important structure at low temperatures. The diffusion of Ti from the Al-top site to the Na-bridging site has a low activation barrier of 0.20 eV and may be activated at the experimental temperatures (∼323 K). The diffusion of a Ti atom into the energetically favoured subsurface interstitial site occurs via the Na-bridging surface site and is essentially barrierless.

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.

Vacancy mediated desorption of hydrogen from a sodium alanate surface : An ab initio spin-polarized study

Ab initio spin-polarized density functional theory calculations are performed to explore the effect of single Na vacancy on NaAlH4(001) surface on the initial dehydrogenation kinetics. The authors found that two Al–H bond lengths become elongated and weakened due to the presence of a Na vacancy on the NaAlH4(001) surface. Spontaneous recombination from the surface to form molecular hydrogen is observed in the spin-polarized ab initio molecular dynamics simulation. The authors’ results indicate that surface Na vacancies play a critical role in accelerating the dehydrogenation kinetics in sodium alanate. The understanding gained here will aid in the rational design and development of complex hydride materials for hydrogen storage

Structural and electronic properties of diazonium functionalized (4, 4) single walled carbon nanotube : an ab initio study

In this work, ab initio density functional theory (DFT) calculations are performed to study the structural and electronic properties of diazonium reagent functionalized (4, 4) single-walled carbon nanotube (SWCNT). We find the aryl group covalently bonds with SWCNT and prefers to be perpendicular to the side wall of nanotube. It has a rotational barrier of 0.35 eV around the formed aryl-tube bond axis and should be thermodynamically stable at room temperature. Additionally, new peaks appeared around the Fermi energy in the density of state (DOS) due to the weak band dispersion. Increasing of the coverage of the functional group will result in significant upshift of the Fermi level.