<p>The aqueous solvation of metal ions is one of the long-standing and complex problems in chemistry, with implications for and applications in a broad range of biochemical and electrochemical systems, where water is the all-pervasive medium.</p><p>This thesis describes computer simulations of Al³⁺(<i>aq</i>), Fe²⁺(<i>aq</i>), Fe³⁺(<i>aq</i>) and Cu²⁺(<i>aq</i>). Various aspects of the solvation of these polyvalent metal ions in water are addressed, at different levels of theory, using Car-Parrinello molecular dynamics, classical molecular dynamics and quantum-mechanical cluster calculations. Polyvalent metal ions are particularly interesting because of their large influence on the solvent structure, dynamics and thermodynamics, as well as on the properties of the individual solvent molecules. Polyvalent metal ions in aqueous solution also constitute a challenging subject for computer simulations since a sophisticated interaction model is needed to incorporate the large many-body effects. </p><p>All the ion-water coordination figures in this thesis are octahedral, except in the Cu²⁺(<i>aq</i>) solution, where the ion is penta-coordinated with four equatorial neighbours in a plane and one axial neighbour located ~0.45 Å further out from the ion. The equatorial ion-water bonds have covalent character, while the axial water molecule is only electrostatically bound. For all the ions, the OD stretching frequencies of the first-shell water molecules are much more downshifted than in liquid water. In the case of Cu²⁺(<i>aq</i>), however, only the OD frequencies of the equatorial water molecules are downshifted with respect to bulk water whereas the OD frequencies of the axial water molecule are slightly upshifted. </p><p>Various limitations of the Car-Parrinello molecular dynamics simulations have been explored and compared, such as finite system-size effects and shortcomings in the electronic structure calculations. The Car-Parrinello simulations are found to give reasonable descriptions of the polyvalent metal ions in aqueous solution.</p>

In this thesis I have performed theoretical studies on the reaction dynamics in few-atom molecules. In particular, I have looked at reaction processes in which highly excited resonant states are involved. When highly excited states are formed, the dynamics becomes complicated and approximations normally used in chemical reaction studies are no longer applicable.To calculate the potential energy curve for some of these states as a function of internuclear distance, a combination of structure calculations and scattering calculations have to be performed, and the reaction dynamics on the potentials has been studied using both time-independent and time-dependent methods.The processes that have been studied and which are discussed in this thesis are ion-pair formation in electron recombination with H3+, dissociative recombination and ion-pair formation of HF+, mutual neutralization in H++F- collisions and dissociative recombination of BeH+. Isotope effects in these reactions have also been investigated. Our calculated cross sections are compared with experimentally measured cross sections for these reactions.

The main focus of this thesis is the investigation of the dissociation dynamics of aryl halides using femtosecond pump-probe spectroscopy. In the monohalogenated aryl halides, iodo-, bromo- and chlorobenzene, the rate of dissociation following excitation at 266 nm in the gas phase increased with increasing mass of the halogen atom. This process was assigned to predissociation of the initially excited singlet (π, π*) state via a repulsive triplet (n, σ*) state due to spin-orbit interaction. In addition to the predissociative mechanism, a direct dissociation channel was observed in iodobenzene. The rate of the predissociation in bromobenzene was found to be faster in the condensed phase than in the gas phase, which can be explained by solvent-induced symmetry perturbations. Ab initio calculations of the potential energy surfaces of the ground state and several low lying excited states in bromobenzene have been performed in order to verify the suggested mechanism. Substituting one of the hydrogen atoms in bromobenzene affected the predissociation rate significantly. In o-, m- and p-dibromobenzene the predissociation rate increased with decreasing distance between the bromine atoms in accordance with an increased spin-orbit interaction introduced by the bromine substituent. The fastest predissociation rate was observed in 1,3,5-tribromobenzene. With chlorine and fluorine substitution, inductive and conjugative effects were found to be of importance. In the o- and m-isomers of the dihalogenated aryl halides, an additional faster dissociation channel was observed. Guided by ab initio calculations of the potential energy surfaces in the dibromobenzene isomers, we ascribed the fast dissociation pathway to predissociation of an initially excited triplet state. Upon methyl group substitution in bromobenzene, the decreased lifetime of the initially excited state was attributed to an incresaed density of coupled states. Another system which has been studied in the condensed phase is diiodomethane. Using Car-Parrinello molecular dynamics simulations we observed a prompt dissociation and subsequent recombination to the isomer, iso-diiodomethane, in acetonitrile solution. Vibrational wavepacket dynamics in the C (1Σ+) state of NaK were studied using a direct ionization probing scheme. A simple analytical expression for the pump-probe signal was developed in order to see what factors that govern direct ionization of the vibrational wavepacket. Our experimental data was consistent with a photoionization transition dipole moment that varies with internuclear distance.

The aqueous solvation of metal ions is one of the long-standing and complex problems in chemistry, with implications for and applications in a broad range of biochemical and electrochemical systems, where water is the all-pervasive medium. This thesis describes computer simulations of Al3+(aq), Fe2+(aq), Fe3+(aq) and Cu2+(aq). Various aspects of the solvation of these polyvalent metal ions in water are addressed, at different levels of theory, using Car-Parrinello molecular dynamics, classical molecular dynamics and quantum-mechanical cluster calculations. Polyvalent metal ions are particularly interesting because of their large influence on the solvent structure, dynamics and thermodynamics, as well as on the properties of the individual solvent molecules. Polyvalent metal ions in aqueous solution also constitute a challenging subject for computer simulations since a sophisticated interaction model is needed to incorporate the large many-body effects. All the ion-water coordination figures in this thesis are octahedral, except in the Cu2+(aq) solution, where the ion is penta-coordinated with four equatorial neighbours in a plane and one axial neighbour located ~0.45 Å further out from the ion. The equatorial ion-water bonds have covalent character, while the axial water molecule is only electrostatically bound. For all the ions, the OD stretching frequencies of the first-shell water molecules are much more downshifted than in liquid water. In the case of Cu2+(aq), however, only the OD frequencies of the equatorial water molecules are downshifted with respect to bulk water whereas the OD frequencies of the axial water molecule are slightly upshifted. Various limitations of the Car-Parrinello molecular dynamics simulations have been explored and compared, such as finite system-size effects and shortcomings in the electronic structure calculations. The Car-Parrinello simulations are found to give reasonable descriptions of the polyvalent metal ions in aqueous solution.

In this thesis I have studied some problems of condensed matter at high pressures and temperatures by means of numerical simulations based on Density Functional Theory (DFT). The stability of MgCO3 and CaCO3 carbonates at the Earth's mantle conditions may play an important role in the global carbon cycle through the subduction of the oceanic crust. By performing ab initio electronic structure calculations, we observed a new high-pressure phase transition within the Pmcn structure of CaCO3. This transformation is characterized by the change of the sp-hybridization state of carbon atom and indicates a change to a new crystal-chemical regime. By performing ab initio Molecular Dynamics simulations we show the new phase to be stable at 250 GPa and 1000K. Thus, the formation of sp3 hybridized bonds in carbonates can explain the stability of MaCO3 and CaCO3 at pressures corresponding to the Earth's lower mantle conditions. We have also calculated phase transition sequence in CaCO3, SrCO3 and BaCO3, and have found that, despite the fact that these carbonates are isostructural and undergo the same type of aragonite to post-aragonite transition, their phase transformation sequences are different at high pressures. The continuous improvement of the high-pressure technique led to the discovery of new composite structures at high pressures and complex phases of many elements in the periodic table have been determined as composite host-guest incommensurate structures. We propose a procedure to accurately describe the structural parameters of an incommensurate phase using ab initio methods by approximating it with a set of analogous commensurate supercells and exploiting the fact that the total energy of the system is a function of structural parameters. By applying this method to the Sc-II phase, we have determined the incommensurate ratio, lattice parameters and Wyckoff positions of Sc-II in excellent agreement with the available experimental data. Moreover, we predict the occurrence of an incommensurate high-pressure phase in Ca from first-principle calculations within this approach. The implementation of DFT in modern electronic structure calculation methods proved to be very successful in predicting the physical properties of a solid at low temperature. One can rigorously describe the thermodynamics of a crystal via the collective excitation of the ionic lattice, and the ab initio calculations give an accurate phonon spectra in the quasi-harmonic approximation. Recently an elegant method to calculate phonon spectra at finite temperature in a self-consistent way by using first principles methods has been developed. Within the framework of self-consistent ab initio lattice dynamics approach (SCAILD) it is possible to reproduce the observed stable phonon spectra of high-temperature bcc phase of Ti, Zr and Hf with a good accuracy. We show that this method gives also a good description of the thermodynamics of hcp and bcc phases of Ti, Zr and Hf at high temperatures, and we provide a procedure for the correct estimation of the hcp to bcc phase transition temperature.

The topic of this Thesis is magnetization dynamics on atomic length scales. A computational scheme, Atomistic Spin Dynamics, based on density functional theory, the adiabatic approximation and the atomic moment approximation is presented. Simulations are performed for chemically disordered systems, antiferromagnets and ferrimagnets and also systems with reduced dimensionality The autocorrelation function of the archetypical spin glass alloy CuMn is sampled in simulations following a quenching protocol. The aging regime can be clearly identified and the dependence of the relaxation on the damping parameter is investigated. The time evolution of pair correlation and autocorrelation functions has been studied in simulations of the dilute magnetic semiconductor GaMnAs. The dynamics reveal a substantial short ranged magnetic order even at temperatures at or above the ordering temperature. The dynamics for different concentrations of As antisites are discussed. Antiferromagnets offer opportunities for ultrafast switching, this is studied in simulations of an artificial antiferromagnet. For the right conditions, the cooperative effect of applied field torque and and the torque from the other sublattice enables very fast switching. The dynamics of bcc Fe precessing in a strong uniaxial anisotropy are investigated. It is demonstrated that the magnetization can shrink substantially due to a spin wave instability. The dynamics of a two-component model ferrimagnet at finite temperature are investigated. At temperatures where the magnetic system is close to the magnetic and angular momentum compensations points of the ferrimagnet, the relaxation in a uniaxial easy exis anisotropy resembles results in recent experiments on ferrimagnetic resonance. The different cases of uniaxial or colossal magnetic anisotropy in nanowires at different temperatures are compared. The magnon softening in a ferromagnetic monolayer is investigated, giving results that compare well with recent experiments. The effect of lattice relaxation can be treated in first principles calculations. Subsequent simulations captures the softening of magnons caused by reduced dimensionality and temperature. / Felaktigt tryckt som Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 706

The ionomer-membrane interface in a membrane electrode assembly connects the catalyst and membrane and allows hydrated protons to move between the catalyst and membrane. The continuous operation of the polymer membrane electrolyte fuel cell at high temperature and/or in frequent freeze/thaw cycles leads to membrane degradation and delamination of the interface, which lower the proton conductivity. In this dissertation, we modeled the chemical degradation and proton conductivity of perfluorosulfonic acid (PFSA) ionomers by ab initio calculations and macroscopic modeling. All ab initio calculations were performed using Gaussian 03 suites of program by employing B3LYP/6-311++G** method/basis set. The macroscopic modeling involves nonequilibrium thermodynamics. The results show that PFSA membranes can degrade both via side-chain and backbone in the presence of hydroxyl radical. The energetics of homolytic bond cleavage show that the C–S bond in the side-chain is the weakest link and breaks exothermally in the presence of hydroxyl radical. The C–S bond in the membrane fragment radical can break at low activation energy. The side-chain degradation also leads to the split of the backbone into two parts. The backbone degradation starts with the reaction of –COOH impurities in the backbone with the hydroxyl radical, which has the lowest activation energy, and follows an “unzipping mechanism”. The reactions in this mechanism are exothermic. The channels in the interface were modeled as cylindrical pores and the anionic charges were fixed on the pore wall. The analytical expression of proton conductivity was derived from the evolution equations for mass and momentum of hydronium ions by using an order of magnitude analysis. The results show that the conductivity increases with increasing water content and pore radius. The conductivity usually increases on decreasing the separation distance between sulfonates on the length and decreases with decreasing sulfonates separation distance on the circumference. The conductivity of the two pores, one of the interface and the other of the membrane, is closer to the conductivity of the pore with the lowest conductivity and its magnitude depends on the relative radius and length of the pores.

Au cours de cette thèse, nous avons étudiés théoriquement les propriétés structurelles et électroniques des super-réseaux d'oxyde de manganèse en structure perovskiteö au moyen de calculs ab initio.Les oxydes de manganèse au lanthane, donnés avec la formule générique La(1-x)A(x)MnO(3) (LAMO) (A un élément divalent), constituent une classe importante d'oxydes de manganèse en raison de leurs diverses propriétés, telles que l'effet de magnétorésistance colossale, leur riche diagramme de phase en fonction du dopage, de la température ou de champs externes, et leur grande température Curie. Ces propriétés peuvent être exploitées dans de nombreuses applications technologiques potentielles telles que les valves de spin ou les injecteurs de spin. Le contrôle des propriétés de ces matériaux peut se faire par dépôt sous forme de films minces ou comme blocs de construction dans des super-réseaux. Lorsque x = 1/3, le La(1-x)A(x)MnO(3) massif est ferromagnétique et métallique grâce au mécanisme de double échange dans la electrons 3d de Mn. Lorsque Mn est dans un état de valence mixte, les orbitales eg (dx2-y2 and dz2) sont partiellement occupées et peuvent se délocaliser sur les atomes de Mn voisins, seulement si ceux-ci sont alignés ferromagnétiquement. Dans des films très minces, puisque la direction perpendiculaire au substrat, c, n'a que quelques cellules unitaires d'épaisseur, seules les interactions dans le plan (ab) sont importantes pour les propriétés thermodynamiques. En agissant sur la géométrie de la couche LAMO, on peut ainsi maximiser l'occupation de l'orbitales dx2-y2 et augmenter l'échange magnétique et la température Curie associée.Notre but était donc de concevoir de nouveaux matériaux avec un ordre orbital 3d spécifique afin d'assurer les propriétés magnétiques souhaitées.Dans cette thèse, nous avons travaillé sur deux types de super-réseaux. Le premier était constitué de deux oxydes de manganèse antiferromagnétiques, non dopés, LaMnO(3) et SrMnO(3), c'est-à-dire des super-réseaux [LaMnO(3)]n/[SrMnO(3)]m. Nous avons étudié l'état fondamental magnétique pour différentes valeurs n et m afin d'expliquer les résultats expérimentaux surprenants. Le deuxième type de super-réseaux que nous avons étudiés est composé de couches métalliques LAMO en alternance avec des couches isolantes. En effet, les super-réseaux avec des interfaces métal-isolant ont un grand potentiel dans les applications de valves de spin. Ainsi, nous avons d'abord considéré des super-réseaux entre composés ferromagnétiques-métalliques et ferroélectriques-isolantes [LAMO](3)/[BTiO(3)](3) (A = Sr ou Ba, B = Ba ou Pb). Dans ces super-réseaux, les propriétés magnétiques sont malheuresement réduites en raison de la délocalisation d'électrons dz2 à l'interface entre Mn et Ti. Dans de tels super-réseaux , nous avons clarifié le rôle de la polarisation des couches ferroélectriques et le rôle des mouvements antiferrodistortifs dans les couches de manganite. Enfin, de manierè à 'éviter que la délocalisation ait lieu à l'interface, nous avons conçu un autre super-réseau avec interface métal-isolant dans lequel nous avons remplacé l'isolant (BTiO(3)) par un oxyde simple (BO): [LAMO]n/[BO]p superlattices (A = Sr ou Ba, B = Ba, Sr ou Mg et n = 3 ou 6, p = 6 ou 2). Dans ces nouveaux super-réseaux, nous avons réussi à promouvoir les occupations des orbitales dx2-y2 dans les interfaces assurant un fort moment magnetique à l'interface et a priori une fort temperature de Curie. Nous avons également montré une faible corrélation entre la conductivité électrique et l'ordre orbitaire. / This thesis theoretically studies structural and electronic properties of perovskite manganese oxide superlattices by means of ab-initio calculations.Lanthanum manganese oxides, given with the generic formula La1−xAxMnO3 (LAMO) (A a divalent element), are an important class of perovskite manganese oxides due to their various exotic properties, such as giant and colossal magnetoresistance effect, rich phase diagrams with respect to doping, temperature or external fields, and intinsic large Curie temperature. These properties can be exploited in many potential technological applications such as spin valves or spin injectors. Controlling the properties of these materials can be done through deposition as thin films or as building blocks in superlattices. When x = 1/3, bulk La1−xAxMnO3 is ferromagnetic and metallic due to the double-exchange mechanism in the Mn 3d shell. When Mn is in a mixed valence state, the eg orbitals (dx2-y2 and dz2) are partially occupied, and can delocalize on neighboring Mn atoms only if the latter are ferromagnetically aligned. In very thin films, since the direction perpendicular to the substrate, c, is only a few unit cell thick, only in-plane (ab) interactions are important for the thermodynamic properties. By acting on the LAMO layer geometry, one can thus maximize the dx2-y2 occupancy and increase the magnetic exchange and related Curie temperature.Our aim was thus to design new materials with desired 3d orbital order so that to ensure desired magnetic properties.In this thesis, we worked on two types of superlattices. The first one was made of two undoped antiferromagnetic manganese oxides LaMnO3 and SrMnO3, i.e. [LaMnO3]n/[SrMnO3]m superlattices. We investigated the magnetic ground state for different n and m values in order to explain suprising experimental results. The second type of superlattices is composed of metallic LAMO layers with alternated insulating layers. Indeed, the superlattices with metal-insulator interfaces have a great potential in spin valves applications. Thus, we first considered the ferromagnetic-metallic/ferroelectric-insulating [LAMO]3/[BTiO3]3 superlattices (A= Sr or Ba, B=Ba or Pb) where the magnetic properties are reduced due to delocalization of dz2 electrons at the interface from Mn to Ti. In such superlattices, we clarified the role of the polarization of the ferroelectric layers and the role of the antiferrodistortive motions in the manganite layers. In order to prevent the delocalization taking place at the interface, we designed another superlattice with metallic-insulator interface where we replaced the insulator (BTiO3) with a simple oxide (BO); that is [LAMO]n/[BO]p superlattices (A= Sr or Ba, B=Ba, Sr or Mg and n=3 or 6, p=6 or 2). Within this new superlattice, we successully managed to promote dx2-y2 orbital occupancies at the interfaces which ensures a large magnetic moment at the interfaces and an expected large Curie temperature. We also showed the weak correlation between electrical conductivity and orbital ordering.

L'avenir de la spintronique repose sur le développement de matériaux ayant des propriétés magnétiques remarquables. L'objectif de cette thèse est de comprendre la physique des deux matériaux fonctionnels proposés pour des applications spintroniques qui utilisent des simulations de la densité fonctionnelle.Nous nous sommes intéressés dans une première partie au ferrite de gallium pour lequel il a été montré que les propriétés dépendaient de la concentration de fer.Les spectres optiques ont été calculés et comparés aux spectres expérimentaux suggérant des niveaux élevés de désordre. Dans la deuxième partie, nous avons montré une polarisation de spin à l’interface hybride formée entre la phthalocyanine de manganèse et la surface de cobalt,en accord avec les expériences de photoémission.La formation de la spinterface a été expliquée par différents mécanismes d'hybridation dans chaque canal de spin.Cette polarisation de spin est coordonnée avec des moments magnétiques induits sur les sites moléculaires. / The future of the spintronics technology requires developing functional materials with remarkable magnetic properties. The aim of this thesis is to understand the physics of functional materials proposed for spintronic applications using ab-initio density functional simulations. We investigated the properties of two different functional materials. We first studied the magnetoelectric gallium ferrite GFO. The dependence of the different properties on the iron concentration has been demonstrated and discussed. The optical spectra were calculated and compared to the experimental once suggesting high levels of iron disorder. In the second part, we demonstrated a highly spin polarized hybrid interface formed between manganese phthalocyanine and cobalt surface in agreement with photoemission experiments. The formation of this spinterface was described by different hybridization mechanisms in each spin channel. This high spin polarization is coordinated with induced magnetic moments on the molecular sites.

Desde sua descoberta por S. Iijima em 1991, os nanotubos de carbono têm sido considerados um dos materiais nanoestruturados mais promissores para o desenvolvimento de novos dispositivos eletrônicos em escala nanoscópica. Devido _a sua alta razão entre a área superficial e o volume, esse material se destaca para aplicações como sensores de gás. No presente trabalho, estudamos através de simulações computacionais, a possibilidade de nanotubos de carbono com defeitos de nitrogênio (os chamados nanotubos CNx), poderem ser usados como sensores de moléculas gasosas. Na primeira parte do trabalho foram realizados cálculos de estrutura eletrônica baseados na Teoria do Funcional da Densidade (DFT) para diferentes sistemas formados pelo nanotubo e pela molécula. Através de cálculos de energia de ligação, foi possível identificar quais gases poderiam ou não serem adsorvidos à superfície do nanotubo. Dentre as moléculas investigadas, o monóxido de carbono e a amônia mostraram ser as mais facilmente adsorvidas ao nanotubo. Na segunda parte, foram realizados cálculos das propriedades de transporte utilizando o formalismo das funções de Green fora do equilíbrio (NEGF) recursivo. Foi possível concluir que os nanotubos estudados poderiam ser usados para detectar o monóxido de carbono e a amônia. Todavia, em relação à seletividade, os resultados indicaram que não parece possível distinguir essas duas moléculas, caso o sistema fosse inserido em um ambiente contendo uma mistura desses gases. Ainda, foram feitas simulações de nanotubos contendo defeitos aleatoriamente distribuídos, de forma a levar em conta os fatores de desordem característicos de sistemas mais realistas. / Since their discovery by S. Iijima in 1991, carbon nanotubes have been considered as one of the most promising nanostructured materials for the development of new nanoscopic electronic devices. Due to its high surface area to volume ratio, this material stands out as a candidate for possible gas sensoring applications. In this thesis, we have studied, by means of computational simulations, the possibility of using carbon nanotubes containing nitrogen defects (the so-called CNx nanotubes) as gas sensors. In the first part, we have performed electronic structure calculations based on Density Functional Theory (DFT) of several systems to address the possible binding of different molecules to the nanotube surface. Our results indicate that, among the molecules which were investigated, carbon monoxide and ammonia adsorb more easily to the nanotube surface. In the second part of this thesis, we have performed calculations of the transport properties by means of non-equilibrium Green\'s function formalism (NEGF). The results have shown that the nitrogen-defect carbon nanotubes could be used to detect, mainly carbon monoxide and ammonia molecules. On the other hand, when dealing with the selectivity of this system, it seems to be not possible to distinguish these gases, in the case of inserting the system in a environment containing a mixture of these molecules. Finally, we have simulated carbon nanotubes with defects randomly distributed along its length, in order to take into account disordering factors usually found in more realistic nanosensors.