First-principle studies of the formation and diffusion of hydrogen vacancies in magnesium hydride

Ab initio density functional theory (DFT) calculations are performed to study the formation and diffusion of hydrogen vacancies on MgH2(110) surface and in bulk. We find that the formation energies for a single H-vacancy increase slightly from the surface to deep layers. The energies for creating adjacent surface divancacies at two inplane sites and at an inplane and a bridge site are even smaller than that for the formation of a single H-vacancy, a fact that is attributed to the strong vacancy−vacancy interactions. The diffusion of an H-vacancy from an in-plane site to a bridge site on the surface has the smallest activation barrier calculated at 0.15 eV and should be fast at room temperature. The activation barriers computed for H-vacancy diffusion from the surface into sublayers are all less than 0.70 eV, which is much smaller than the activation energy for desorption of hydrogen on the MgH2(110) surface (1.78−2.80 eV/H2). This suggests that surface desorption is more likely than vacancy diffusion to be rate determining, such that finding effective catalyst on the MgH2 surface to facilitate desorption will be very important for improving overall dehydrogenation performance.

First-principle study of adsorption of hydrogen on ti-doped Mg(0001) surface

Ab initio density functional theory (DFT) calculations are performed to study the adsorption of H2 molecules on a Ti-doped Mg(0001) surface. We find that two hydrogen molecules are able to dissociate on top of the Ti atom with very small activation barriers (0.103 and 0.145 eV for the first and second H2 molecules, respectively). Additionally, a molecular adsorption state of H2 above the Ti atom is observed for the first time and is attributed to the polarization of the H2 molecule by the Ti cation. Our results parallel recent findings for H2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.

Atomic hydrogen diffusion in novel magnesium nanostructures : the impact of incorporated subsurface carbon atoms

Ab initio Density Functional Theory (DFT) calculations are performed to study the diffusion of atomic hydrogen on a Mg(0001) surface and their migration into the subsurface layers. A carbon atom located initially on a Mg(0001) surface can migrate into the sub-surface layer and occupy a fcc site, with charge transfer to the C atom from neighboring Mg atoms. The cluster of postively charged Mg atoms surrounding a sub-surface C is then shown to facilitate the dissociative chemisorption of molecular hydrogen on the Mg(0001) surface, and the surface migration and subsequent diffusion into the subsurface of atomic hydrogen. This helps rationalize the experimentally-observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

The role of Ti as a catalyst for the dissociation of hydrogen on a Mg(0001) surface

In this paper, the dissociative chemisorption of hydrogen on both pure and Ti-incorporated Mg(0001) surfaces are studied by ab initio density functional theory (DFT) calculations. The calculated dissociation barrier of hydrogen molecule on a pure Mg(0001) surface (1.05 eV) is in good agreement with comparable theoretical studies. For the Ti-incorporated Mg(0001) surface, the activated barrier decreases to 0.103 eV due to the strong interaction between the molecular orbital of hydrogen and the d metal state of Ti. This could explain the experimentally observed improvement in absorption kinetics of hydrogen when transition metals have been introduced into the magnesium materials.

Phonon anomalies predict superconducting Tc for AlB2-type structures

We show that the well-known Kohn anomaly predicts Tc for ordered AlB2-type structures. We use ab initio Density Functional Theory to calculate phonon dispersions for Mg1-xAlxB2 compositions and identify a phonon anomaly with magnitude that predicts experimental values of Tc for all x. Key features of these anomalies correlate with the electronic structure of Mg1-xAlxB2. This approach predicts Tc for other known AlB2-type structures as well as new compositions. We predict that Mg0.5Ba0.5B2 will show Tc = 63.6 ± 6.6 K. Other forms of the Mg1-xBaxB2 series will also be superconductors when successfully synthesised. Our calculations predict that the end-member composition, BaB2, is likely to show a Tc significantly higher than currently achieved by other diborides although an applied pressure ~16 GPa may be required to stabilise the structure.

Methane and carbon dioxide adsorption on edge-functionalized graphene: a comparative DFT study

With a view towards optimizing gas storage and separation in crystalline and disordered nanoporous carbon-based materials, we use ab initio density functional theory calculations to explore the effect of chemical functionalization on gas binding to exposed edges within model carbon nanostructures. We test the geometry, energetics, and charge distribution of in-plane and out-of-plane binding of CO2 and CH4 to model zigzag graphene nanoribbons edge-functionalized with COOH, OH, NH2, H2PO3, NO2, and CH3. Although different choices for the exchange-correlation functional lead to a spread of values for the binding energy, trends across the functional groups are largely preserved for each choice, as are the final orientations of the adsorbed gas molecules. We find binding of CO2 to exceed that of CH4 by roughly a factor of two. However, the two gases follow very similar trends with changes in the attached functional group, despite different molecular symmetries. Our results indicate that the presence of NH2, H2PO3, NO2, and COOH functional groups can significantly enhance gas binding, making the edges potentially viable binding sites in materials with high concentrations of edge carbons. To first order, in-plane binding strength correlates with the larger permanent and induced dipole moments on these groups. Implications for tailoring carbon structures for increased gas uptake and improved CO2/CH4 selectivity are discussed. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4736568]

Spin-glass state in nanoparticulate (La0.7Sr0.3MnO3)(1-x) (BaTiO3)(x) solid solutions: Experimental and density-functional studies

We report the transition from robust ferromagnetism to a spin- glass state in nanoparticulate La0.7Sr0.3MnO3 through solid solution with BaTiO3. The field- and temperature-dependent magnetization and the frequency-dependent ac magnetic susceptibility measurements strongly indicate the existence of a spin- glass state in the system, which is further confirmed from memory effect measurements. The breaking of long-range ordering into short-range magnetic domains is further investigated using density-functional calculations. We show that Ti ions remain magnetically inactive due to insufficient electron leakage from La0.7Sr0.3MnO3 to the otherwise unoccupied Ti-d states. This results in the absence of a Mn-Ti-Mn spin exchange interaction and hence the breaking of the long-range ordering. Total-energy calculations suggest that the segregation of nonmagnetic Ti ions leads to the formation of short-range ferromagnetic Mn domains.

Carbon monoxide adsorption on molybdenum phosphides: Fourier transform infrared spectroscopic and density functional theory studies

Molybdenum phosphide (MoP) and supported molybdenum phosphide (MoP/gamma-Al2O3) have been prepared by the temperature-programmed reduction method. The surface sites of the MoP/gamma-Al2O3 catalyst were characterized by carbon monoxide (CO) adsorption with in situ Fourier transform infrared (FT-IR) spectroscopy. A characteristic IR band at 2037 cm(-1) was observed on the MoP/gamma-Al2O3 that was reduced at 973 K. This band is attributed to linearly adsorbed CO on Mo atoms of the MoP surface and is similar to IR bands at 2040-2060 cm(-1), which correspond to CO that has been adsorbed on some noble metals, such as platinum, palladium, and rhodium. Density functional calculations of the structure of molybdenum phosphides, as well as CO chemisorption on the MoP(001) surface, have also been studied on periodic surface models, using the generalized gradient approximation (GGA) for the exchange-correlation functional. The results show that the chemisorption of CO on MoP occurred mainly on top of molybdenum, because the bonding of CO requires a localized mininum potential energy. The adsorption energy obtained is DeltaH(ads) approximate to -2.18 eV, and the vibrational frequency of CO is 2047 cm-1, which is in good agreement with the IR result of CO chernisorption on MoP/gamma-Al2O3.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

Density functional/Monte Carlo study of ring-opening polymerization

Density functional calculations of the structure, potential energy surface and reactivity for organic systems closely related to bisphenol-A-polycarbonate (BPA-PC) provide the basis for a model describing the ring-opening polymerization of its cyclic oligomers by nucleophilic molecules. Monte Carlo simulations using this model show a strong tendency to polymerize that is increased by increasing density and temperature, and is greater in 3D than in 2D. Entropy in the distribution of inter-particle bonds is the driving force for chain formation. (C) 2002 Elsevier Science B.V. All rights reserved.

Density functional study of reactions of phenoxides with polycarbonate

Density functional calculations with simulated annealing have been used to study the reactions of chains of bisphenol A polycarbonate (BPA-PC) with sodium phenoxide (NaOPh), diphenyl carbonate (DPC), and tetraphenylphosphonium phenoxide (PPh4OPh). These calculations extend our work on the reactions of LiOPh, NaOPh, and phenol with the cyclic tetramer of BPA-PC. We study, in particular, chain growth catalyzed by NaOPh and PPh4OH. The energy barriers for reactions with PPh4OPh are somewhat larger than those involving LiOPh and NaOPh, but they are significantly lower than those involving phenol (HOPh), due in part to the collective rearrangement of phenyl groups in the reacting molecules. We discuss in the Appendix the bonds between alkali metal atoms (Na in the present calculations) and other atoms (here oxygen) that are analogous to the more familiar "hydrogen bonds".

Examining the redox and formate mechanisms for water-gas shift reaction on Au/CeO2 using density functional theory

We perform DFT calculations to investigate the redox and formate mechanisms of water-gas-shift (WGS) reaction on Au/CeO2 catalysts. In the redox mechanism, we analyze all the key elementary steps and find that the OH cleavage is the key step. Three possible pathways of OH cleavage are calculated: (1) OHad '' + *'--> H-ad' + O-ad"; (2) H-ad' + OHad '' --> H-2(g) + O-ad '' + *'; and (3) OHad" + OHad '' --> 2O(ad '') + H-2(g) (*': the free adsorption sites on the oxides; ad': adsorption on the metal; ad": adsorption on the oxide, respectively). In the formate mechanism, we identify all the possible pathways for the formation and decomposition of surface formates in the WGS reaction. It is found that there is a shortcoming in the redox and formate mechanisms which is related to surface oxygen reproduction. Four possible pathways for producing surface oxygen are studied, and all the barriers of the four pathways are more than 1 eV. Our results suggest that the processes to reproduce surface oxygen in the reaction circle are not kinetically easy. (C) 2008 Elsevier B.V. All rights reserved.

Methyl chemisorption on Ni(111) and C-H-M multicentre bonding:a density functional theory study

Density functional theory has been used to study the adsorption of CH3 on Ni(111). CH3 is found to adsorb strongly at all four high symmetry sites of the Ni(111) surface. Calculated adsorption energies of CH3 on the different sites are in the following order: hcp approximate to fcc>bridge>top. The bonding and structures of CH3 on the different sites are analysed in detail. An important factor, namely three-centre bonding between carbon, hydrogen and nickel which contributes to the 'soft' C-H vibrational frequency of CH3 on Ni(111), and may determine the preferred chemisorption site, is stressed. (C) 1999 Elsevier Science B.V. All rights reserved.

A density functional theory study on the interaction between chemisorbed CO and S on Rh(111)

Density functional theory calculations are carried out for Rh(111)-p(2 x 2)-CO, Rh(111)-p(2 x 2)-S, Rh(111)-p(2 x 2)-(S + CO), Rh(111)-p(3 x 3)-CO, Rh(111)-p(3 x 3)-S and Rh(111)-p(3 x 3)-(S + CO), aiming to shed some light on the S poisoning effect. Geometrical structures of these systems are optimized and chemisorption energies are determined. The presence of S does not significantly influence the geometrical structure and chemisorption energy of CO and vice versa, which strongly suggests that the interaction between CO and S on the Rh(111) surface is mainly short-range in nature. The long range electronic effect for the dramatic attenuation of the CO methanation activity by sulfur is likely to be incorrect. It is suggested that an ensemble effect may be dominant in the catalytic deactivation. (C) 1999 Elsevier Science B.V. All rights reserved.

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.