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.

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.

Synergetic effects of the Cu/Pt{110} surface alloy: enhanced reactivity of water and carbon monoxide

We have used synchrotron-based high-resolution X-ray photoelectron spectroscopy in combination with ab initio density functional theory calculations to investigate the characteristics of water and CO adsorption on the bimetallic Cu/Pt{110}-(2 x 1) surface at a Cu coverage near 0.5 ML. Cu fills the troughs of the reconstructed clean surface forming nanowires, which are stable up to 830 K. Their presence dramatically influences the adsorption of water and CO. Water adsorption changes from intact to partially dissociated while the desorption temperature of CO on this surface increases by up to 27 K with respect to the clean Pt{110} surface. Ab initio calculations and experimental valence band spectra reveal that the Cu 3d-band is narrowed and shifted upward with respect to bulk Cu surfaces. This and electron donation to surface Pt atoms cause the increase in the bond strength between CO and the Pt surface atoms. The pathway for water dissociation occurs via Cu surface atoms. The heat of adsorption of water bonding to Cu surface atoms was calculated to be 0.82 eV, which is significantly higher than on the clean Pt{110} surface; the activation energy for partial dissociation is 0.53 eV (not corrected for zero point energy).

Inelastic neutron scattering in spectroscopic studies of hydrogen on carbon-supported catalysts-experimental spectra and computed spectra of model systems

Inelastic neutron scattering spectroscopy has been used to observe and characterise hydrogen on the carbon component of a Pt/C catalyst. INS provides the complete vibration spectrum of coronene, regarded as a molecular model of a graphite layer. The vibrational modes are assigned with the aid of ab initio density functional theory calculations and the INS spectra by the a-CLIMAX program. A spectrum for which the H modes of coronene have been computationally suppressed, a carbon-only coronene spectrum, is a better representation of the spectrum of a graphite layer than is coronene itself. Dihydrogen dosing of a Pt/C catalyst caused amplification of the surface modes of carbon, an effect described as H riding on carbon. From the enhancement of the low energy carbon modes (100-600 cm(-1)) it is concluded that spillover hydrogen becomes attached to dangling bonds at the edges of graphitic regions of the carbon support. (C) 2003 Elsevier Science B.V. All rights reserved.

Ab initio and tight-binding calculations of noncollinear magnetism in manganese clusters

Combining ab initio and tight-binding calculations, we have studied the noncollinear magnetism in manganese clusters. The oscillations in the per-atom moments observed experimentally are reproduced theoretically. The tendency of antiferromagnetic coupling between near neighbors leads to noncollinear coupling between atoms within the clusters. For clusters containing 12, 13, 15, 19, and 23 atoms, the geometrical structures were optimized from ab initio calculations with collinear coupled spin moments among different atomic sites. For larger clusters such as Mn-36 and Mn-55, the geometries are taken as portions of an fcc structure. Although the local atomic moments have high values close to 4 mu(B), the net moments lie in the range of 0.4-1.2 mu(B)/atom. Taking the noncollinear coupling into account brings the calculated magnetic moments much closer to the experimental results.

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.

Relationship between phonons and thermal expansion in Zn(CN)2 and Ni(CN)2 from inelastic neutron scattering and ab initio calculations

Zn(CN)2 and Ni(CN)2 are known for exhibiting anomalous thermal expansion over a wide temperature range. The volume thermal expansion coefficient for the cubic, three dimensionally connected material, Zn(CN)2, is negative (alpha(V) = −51  10(-6) K-1) while for Ni(CN)2, a tetragonal material, the thermal expansion coefficient is negative in the two dimensionally connected sheets (alpha(a) = −7  10(-6) K-1), but the overall thermal expansion coefficient is positive (alpha(V) = 48  10(-6) K-1). We have measured the temperature dependence of phonon spectra in these compounds and analyzed them using ab initio calculations. The spectra of the two compounds show large differences that cannot be explained by simple mass renormalization of the modes involving Zn (65.38 amu) and Ni (58.69 amu) atoms. This reflects the fact that the structure and bonding are quite different in the two compounds. The calculated pressure dependence of the phonon modes and of the thermal expansion coefficient, alpha(V), are used to understand the anomalous behavior in these compounds. Our ab initio calculations indicate that phonon modes of energy approx. 2 meV are major contributors to negative thermal expansion (NTE) in both the compounds. The low-energy modes of approx.8 and 13 meV in Zn(CN)2 also contribute significantly to the NTE in Zn(CN)2 and Ni(CN)2, respectively. The measured temperature dependence of the phonon spectra has been used to estimate the total anharmonicity of both compounds. For Zn(CN)2, the temperature-dependent measurements (total anharmonicity), along with our previously reported pressure dependence of the phonon spectra (quasiharmonic), is used to separate the explicit temperature effect at constant volume (intrinsic anharmonicity).

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

Strongly nucleophilic RhI centre in square-planar complexes with terdentate (κ3) 2,2′:6′,2′′-terpyridine ligands: crystallographic, electrochemical and density functional theoretical studies

Rh-I-terpyridine complexes have been unambiguously formed for the first time. The 2,21:6',2"-terpyridine (tpy), 4'-chloro-2,2':6',2"-terpyridine (4'-Cl-tpy) and 4'-(tert-butyldimethylsilyl-ortho-carboranyl)-2,2':6',2"-terpyridine (carboranyl-tpy) ligands were used for successful syntheses and characterisation of the corresponding Rh-I complexes with halide coligands, [Rh(X)(4'-Y-terpyridine)] (X = Cl, Y = H, Cl, carboranyl; X = Br, Y = H). All four neutral Rh-tpy complexes are square planar, with Rh-X bonds in the plane of the 4'-Y-terpyridine ligands. Full characterisation of these dark blue, highly air-sensitive compounds was hampered by their poor solubility in various organic solvents. This is mainly due to the formation of pi-stacked aggregates, as evidenced by the crystal structure of [Rh(Cl)(tpy)]; in addition, [Rh(Cl)(carboranyl-tpy)] merely forms discrete dimers. The (bonding) properties of the novel Rh-I-terpyridine complexes have been studied with single-crystal X-ray diffraction, (time-dependent) density functional theoretical (DFT) calculations, far-infrared spectroscopy, electronic absorption spectroscopy and cyclic voltammetry. From DFT calculations, the HOMO of the studied Rh-I-terpyridine complexes involves predominantly the metal centre, while the LUMO resides on the terpyridine ligand. Absorption bands of the studied complexes in the visible region (400-900 nm) can be assigned to MLCT and MLCT/XLCT transitions. The relatively low oxidation potentials of [Rh(X)(tpy)] (X = Cl, Br) point to a high electron density on the metal centre. This makes the Rh-I-terpyridine complexes strongly nucleophilic and (potentially) highly reactive towards various (small) substrate molecules containing carbon-halide bonds.

Synthesis, chemical, electrochemical characterization of oxomolybdenum(V) complexes of 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole ligand: crystal structure of ligand and oxomolybdenum(VI) complex and density functional theory (DFT) calculation

This work reports the ligational behavior of the neutral bidentate chelating molecule 2-(3,5-dimethyl pyrazol-1-yl) benzothiazole towards the oxomolybdenum(V) center. Both mononuclear complexes of the type (MoOX3L)-O-V and binuclear complexes of the formula (Mo2O4X2L2)-O-V (where X = Cl, Br) are isolated in the solid state. The complexes are characterized by elemental analyses, various spectroscopic techniques (UV-Vis IR), magnetic susceptibility measurement at room temperature, and cyclic voltammetry for their redox behavior at a platinum electrode in CH3CN. The mononuclear complexes (MoOX3L)-O-V are found to be paramagnetic while the binuclear complexes Mo2O4X2L2 are diamagnetic. Crystal and molecular structure of the ligand and the dioxomolybdenum complex (MoO2Br2L)-O-VI (obtained from the complex MoOBr3L during crystallization) have been solved by single crystal X-ray diffraction technique. Relevant DFT calculations of the ligand and the complex (MoO2Br2L)-O-VI are also carried out.

The interplay of secondary Hg⋯S, Hg⋯N and Hg⋯π bonding interactions in supramolecular structures of phenylmercury(ii) dithiocarbamates

Three new phenylmercury(II) and one mercury(II) dithiocarbamate complexes viz. PhHg S2CN(PyCH2) Bz (1), PhHg S2CN(PyCH2)CH3 (2), PhHg S2CN(Bz)CH3 (3), and [Hg (NCS2(PyCH2)Bz)(2)] (4) (Py = pyridine; Bz = benzyl) have been synthesized and characterized by elemental analyses, IR, electronic absorption, H-1 and C-13 NMR spectroscopy. The crystal structures of 1, 2 and 3 showed a linear S-Hg-C core at the centre of the molecule, in which the metal atom is bound to the sulfur atom of the dithiocarbamate ligand and a carbon atom of the aromatic ring. In contrast the crystal structure of 4 showed a linear S-Hg-S core at the Hg(II) centre of the molecule. Weak intermolecular Hg center dot center dot center dot N (Py) interactions link molecules into a linear chain in the case of 1, whereas chains of dimers are formed in 2 through intermolecular Hg center dot center dot center dot N (Py) and Hg center dot center dot center dot S interactions. 3 forms a conventional face-to-edge dimeric structure through intermolecular Hg center dot center dot center dot S secondary bonding and 4 forms a linear chain of dimers through face-to-face Hg center dot center dot center dot S secondary bonding. In order to elucidate the nature of these secondary bonding interactions and the electronic absorption spectra of the complexes, ab initio quantum chemical calculations at the MP2 level and density functional theory calculations were carried out for 1-3. Complexes 1 and 2 exhibited photoluminescent properties in the solid state as well as in the solution phase. Studies indicate that Hg center dot center dot center dot S interactions decrease and Hg center dot center dot center dot N interactions increase the chances of photoluminescence in the solid phase

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.