A candidate hydrogen storing material should have high storage capacity and fast dehydrogenation kinetics. On this basis, magnesium hydride (MgH2) is an outstanding compound with 7.66 wt % storage capacity, despite its slow dehydriding kinetics and high desorption temperature. Therefore in this study, bulk and surface alloys of Mg with improved hydrogen desorption characteristics were investigated. In this respect, formation energies of alloyed bulk MgH2 as well as the adsorption energies on alloyed magnesium (Mg) and MgH2 surface structures were calculated by total energy pseudopotential methods. Furthermore, the effect of substitutionally placed dopants on the dissociation of hydrogen molecule (H2) at the surface of Mg was studied via Molecular Dynamics (MD). The results displayed that 31 out of 32 selected dopants contributed to the decrease in formation energy of MgH2 within a range of ~ 37 kJ/mol-H2 where only Sr did not exhibit any such effect. The most favorable elements in this respect came out to be / P, K, Tl, Si, Sn, Ag, Pb, Au, Na, v Mo, Ge and In. Afterwards, a systematical study within adsorption characteristics of hydrogen on alloyed Mg surfaces (via dynamic calculations) as well as calculations regarding adsorption energies of the impurity elements were performed. Accordingly, Mo and Ni yielded lower adsorption energies / -9.2626 and -5.2995 eV for substitutionally alloyed surfaces, respectively. MD simulations presented that Co is found to have a splitting effect on H2 in 50 fs, where the first hydrogen atom is immediately adsorbed on Mg substrate. Finally, charge density distributions were realized to verify the distinguished effects of most 3d and 4d transition metals in terms of their catalyzer effects.

We use the first-principles molecular-dynamics¡@simulation method (MD), which is based on the density functional theory (DFT) with local-density approximation (LDA), to calculate the sticking property of a carbon atom on hydrogen covered C(100) surface. We focused on trajectories and kinetic energy transfer of the gas-phase C atom for four incident angles of =0, £k/8, £k/6 and £k/4. We find that the calculated trajectories and the kinetic energy transfer of the gas-phase atom, Cn, overall are not very sensitive to the change of the incident angle. The insensitivity of the sticking property on the incident angle may be due to a large chemisorption energy, which bends the trajectory of Cn toward the surface, so that Cn is confined to move within a small range.

The (001) surface of silicon has been the topic of our study in this thesis. The clean surface, an-adatom or submonolayer adsorption on the surface, the monolayer adsorption and its stability conditions as well as growth simulation on the surface were investigated using the state of the art techniques. We have used ab initio density functional calculations based on norm-conserving pseudopotentials to investigate the Mg adsorption on the Si(001) surface for 1/4, 1/2 and 1 monolayer (ML) coverages. For both 1/4 and 1/2 ML coverages it has been found that the most favorable site for the Mg adsorption is the cave site between two dimer rows consistent with recent experiments. For the 1 ML coverage (2 Mg atoms per 2X1 unit cell) we have found that the most preferable configuration is when both Mg atoms on 2X1 reconstruction occupy the two shallow sites. We have found that the minimum energy configurations for 1/4 ML coverage is a 2X2 reconstruction while for the 1/2 and 1 ML coverages they are 2X1. Same method was also used to investigate the surface stress and energetics of the clean-, Sb-adsorbed-, and Sb-interdiffused-Si(001) surface. It is found that interdiffusion of Sb into deeper layers of Si(001) leads to a more isotropic surface stress but corresponds to a higher total energy configuration. As a result of competition between stress relief and energy gain, the surface with all the Sb atoms adsorbed on top of Si(001) surface layer is predicted to have a less ordered geometry and roughness in z-direction. We have repeated the similar calculations on the Ge(001) surface for comparison. Finally using empirical molecular dynamics method, we have investigated the crystalline growth of silicon on Si(001) as a function of substrate temperature and incident particle energy. Our results show that the increase of substrate temperature enhances the crystallinity in the film grown on the Si(001) surface, on the other hand, the crystalline growth can be enhanced at low temperature by using higher incidence energy.

<p>Methods for performing large-scale parallel Molecular Dynamics(MD) simulations are investigated. A perspective on the field of parallel MD simulations is given. Hardware and software aspects are characterized and the interplay between the two is briefly discussed. </p><p>A method for performing <i>ab initio </i>MD is described; the method essentially recomputes the interaction potential at each time-step. It has been tested on a system of liquid water by comparing results with other simulation methods and experimental results. Different strategies for parallelization are explored.</p><p>Furthermore, data-parallel methods for short-range and long-range interactions on massively parallel platforms are described and compared. </p><p>Next, a method for treating electrostatic interactions in MD simulations is developed. It combines the traditional Ewald summation technique with the nonuniform Fast Fourier transform---ENUF for short. The method scales as <i>N log N</i>, where <i>N </i>is the number of charges in the system. ENUF has a behavior very similar to Ewald summation and can be easily and efficiently implemented in existing simulation programs.</p><p>Finally, an outlook is given and some directions for further developments are suggested.</p>

Methods for performing large-scale parallel Molecular Dynamics(MD) simulations are investigated. A perspective on the field of parallel MD simulations is given. Hardware and software aspects are characterized and the interplay between the two is briefly discussed. A method for performing ab initio MD is described; the method essentially recomputes the interaction potential at each time-step. It has been tested on a system of liquid water by comparing results with other simulation methods and experimental results. Different strategies for parallelization are explored. Furthermore, data-parallel methods for short-range and long-range interactions on massively parallel platforms are described and compared. Next, a method for treating electrostatic interactions in MD simulations is developed. It combines the traditional Ewald summation technique with the nonuniform Fast Fourier transform---ENUF for short. The method scales as N log N, where N is the number of charges in the system. ENUF has a behavior very similar to Ewald summation and can be easily and efficiently implemented in existing simulation programs. Finally, an outlook is given and some directions for further developments are suggested.

The catalytic oxidation activity of platinum particles in automobile catalysts is thought to originate from the presence of highly reactive superficial oxide phases which form under oxygen-rich reaction conditions. The thermodynamic stability of platinum oxide surfaces and thin films was studied, as well as their reactivities towards oxidation of carbon compounds by means of first-principles atomistic thermodynamics calculations and molecular dynamics simulations based on density functional theory. On the Pt(111) surface the most stable superficial oxide phase is found to be a thin layer of alpha-PtO2, which appears not to be reactive towards either methane dissociation or carbon monoxide oxidation. A PtO-like structure is most stable on the Pt(100) surface at oxygen coverages of one monolayer, while the formation of a coherent and stress-free Pt3O4 film is favoured at higher coverages. Bulk Pt3O4 is found to be thermodynamically stable in a region around 900 K at atmospheric pressure. The computed net driving force for the dissociation of methane on the Pt3O4(100) surface is much larger than on all other metallic and oxide surfaces investigated. Moreover, the enthalpy barrier for the adsorption of CO molecules on oxygen atoms of this surface is as low as 0.34 eV, and desorption of CO2 is observed to occur without any appreciable energy barrier in molecular dynamics simulations. These results, combined, indicate a high catalytic oxidation activity of Pt3O4 phases that can be relevant in the contexts of Pt-based automobile catalysts and gas sensors.

Nous avons étudié des magmas ferrifères silicatés magnésiens à la pression du manteau terrestre en utilisant la dynamique moléculaire (First Principles Molecular Dynamics). Les résultats de l’équation d’état que nous avons obtenus à partir de nos simulations ont été utilisés pour créer un modèle chimique et minéralogique pour les zones de très basse vitesse sismique (ULVZ, anomalies régionales dans le manteau proche de la limite noyau-manteau). De plus, nous avons étudié le comportement du Ni, du Co et du Fe dans ces magmas et établi la dépendance du spin en fonction de la concentration, de la pression, de la température et du degré de polymérisation du magma silicaté. Nous avons montré qu’une baisse du spin moyen peut être corrélée au changement de pente (kink) observé précédemment pour les coefficients de partage du Ni et du Co. Nous avons analysé la structure du magma pour toutes les compositions étudiées en fonction de la pression. Nos résultats donnent un nouvel aperçu de la coordination des éléments majeurs et traces dans les magmas silicatés de différents degrés de polymérisation. Nous interprétons l’anomalie de coordination Ni-O en fonction de la pression comme un changement d’état de spin. L’effet de la polymérisation du magma silicaté sur les coefficients de partage du Co, du Ni et du W entre le métal et le magma silicaté a été étudié par expériences multi-enclumes en conditions isobares et isothermes. Nous avons réalisé des simulations FPMD de magmas à des degrés de polymérisation similaires aux expériences afin d’expliquer le caractère de plus en plus lithophile du W lorsque le degré de polymérisation du magma silicaté diminue. Nous proposons une explication structurale pour expliquer l’affinité décroissante apparente du W dans les magmas silicatés dépolymérisés. / We explore Fe-bearing Mg-silicate melts through the pressure regime of the Earth’s mantle using First Principles Molecular Dynamics (FPMD). The equation of state results we obtained from our simulations are used to create a chemical and mineralogical model for Ultra-Low Velocity Zones (anomalous region on the mantle side of the core-mantle boundary). Furthermore we study the behaviour of Ni, Co, and Fe in these melts, and asses their spin-crossover dependencies on their concentration, pressure, temperature, and the degree of polymerization of the silicate melts. We show that a decrease in the average spin can be correlated with the previously observed kink in the partitioning coefficient of Ni and Co. We investigate the melt structure of all the compositions studied as a function of pressure. Our results provide new insight into the coordination of major and trace elements in silicate melts with different degrees of polymerization. We interpret the anomalous Ni-O coordination trend with pressure as the result of the spin state change. The effect of silicate melt polymerization on the partitioning of Co, Ni, and W between a metal and silicate melt, is investigated at isobaric and isothermic conditions using multi-anvil experiments. We have performed FPMD simulations of melts with similar degrees of polymerization as the experiments in order to explain the increasing lithophile character of W with the decrease in polymerization of the silicate melt. We propose a structural explanation for tungsten’s apparent increased affinity for depolymerized silicate melts.

The catalytic oxidation activity of platinum particles in automobile catalysts is thought to originate from the presence of highly reactive superficial oxide phases which form under oxygen-rich reaction conditions. The thermodynamic stability of platinum oxide surfaces and thin films was studied, as well as their reactivities towards oxidation of carbon compounds by means of first-principles atomistic thermodynamics calculations and molecular dynamics simulations based on density functional theory. On the Pt(111) surface the most stable superficial oxide phase is found to be a thin layer of alpha-PtO2, which appears not to be reactive towards either methane dissociation or carbon monoxide oxidation. A PtO-like structure is most stable on the Pt(100) surface at oxygen coverages of one monolayer, while the formation of a coherent and stress-free Pt3O4 film is favoured at higher coverages. Bulk Pt3O4 is found to be thermodynamically stable in a region around 900 K at atmospheric pressure. The computed net driving force for the dissociation of methane on the Pt3O4(100) surface is much larger than on all other metallic and oxide surfaces investigated. Moreover, the enthalpy barrier for the adsorption of CO molecules on oxygen atoms of this surface is as low as 0.34 eV, and desorption of CO2 is observed to occur without any appreciable energy barrier in molecular dynamics simulations. These results, combined, indicate a high catalytic oxidation activity of Pt3O4 phases that can be relevant in the contexts of Pt-based automobile catalysts and gas sensors.