Dinâmica molecular de sistemas iônicos poliatômicos: modelos polarizável e não-polarizável

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito.

Magnetic and orbital blocking in Ni nanocontacts

We address the fundamental question of whether magnetoresistance (MR) of atomic-sized contacts of nickel is very large because of the formation of a domain wall (DW) at the neck. Using ab initio transport calculations we find that, as in the case of nonmagnetic electrodes, transport in Ni nanocontacts depends very much on the orbital nature of the electrons. Our results are in agreement with several experiments in the average value of the conductance. On the other hand, contrary to existing claims, DW scattering does not account for large MR in Ni nanocontacts.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Spin alignment of extra electrons in K-phenanthrene clusters taken from the crystalline tripotassium-intercalated phenanthrene structure

The appearance of ferromagnetic correlations among π electrons of phenanthrene (C14H10) molecules in the herringbone structure is proven for K doped clusters both by ab initio quantum-chemistry calculations and by the direct solution of the many-body Pariser-Parr-Pople Hamiltonian. Magnetic ground states are predicted for one or three additional electrons per phenanthrene molecule. These results are a consequence of the small overlap between the lowest unoccupied molecular orbitals (and lowest unoccupied molecular orbitals + 1) of neutral neighboring phenanthrene molecules, which makes the gain in energy by delocalization similar to the corresponding increase due to the Coulomb interaction.

Understanding the structure of the first atomic contact in gold

We have studied experimentally jump-to-contact (JC) and jump-out-of-contact (JOC) phenomena in gold electrodes. JC can be observed at first contact when two metals approach each other, while JOC occurs in the last contact before breaking. When the indentation depth between the electrodes is limited to a certain value of conductance, a highly reproducible behaviour in the evolution of the conductance can be obtained for hundreds of cycles of formation and rupture. Molecular dynamics simulations of this process show how the two metallic electrodes are shaped into tips of a well-defined crystallographic structure formed through a mechanical annealing mechanism. We report a detailed analysis of the atomic configurations obtained before contact and rupture of these stable structures and obtained their conductance using first-principles quantum transport calculations. These results help us understand the values of conductance obtained experimentally in the JC and JOC phenomena and improve our understanding of atomic-sized contacts and the evolution of their structural characteristics.

Towards Accurate and Efficient Description of Excited States

Thesis (Ph.D.)--University of Washington, 2016-06

Developments in the Theory of X-ray Absorption and Compton Scattering.

Thesis (Ph.D.)--University of Washington, 2016-06

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.

Experimental and theoretical investigation of the structural, chemical, electronic, and high frequency dielectric properties of barium cadmium tantalate-based ceramics

Single-phase Ba(Cd1/3Ta2/3)O-3 powder was produced using conventional solid state reaction methods. Ba(Cd1/3Ta2/3)O-3 ceramics with 2 wt % ZnO as sintering additive sintered at 1550 degreesC exhibited a dielectric constant of similar to32 and loss tangent of 5x10(-5) at 2 GHz. X-ray diffraction and thermogravimetric measurements were used to characterize the structural and thermodynamic properties of the material. Ab initio electronic structure calculations were used to give insight into the unusual properties of Ba(Cd1/3Ta2/3)O-3, as well as a similar and more widely used material Ba(Zn1/3Ta2/3)O-3. While both compounds have a hexagonal Bravais lattice, the P321 space group of Ba(Cd1/3Ta2/3)O-3 is reduced from P (3) under bar m1 of Ba(Zn1/3Ta2/3)O-3 as a result of a distortion of oxygen away from the symmetric position between the Ta and Cd ions. Both of the compounds have a conduction band minimum and valence band maximum composed of mostly weakly itinerant Ta 5d and Zn 3d/Cd 4d levels, respectively. The covalent nature of the directional d-electron bonding in these high-Z oxides plays an important role in producing a more rigid lattice with higher melting points and enhanced phonon energies, and is suggested to play an important role in producing materials with a high dielectric constant and low microwave loss. (C) 2005 American Institute of Physics.

Mg-based nanocomposites with high capacity and fast kinetics for hydrogen storage

Magnesium and its alloys have shown a great potential in effective hydrogen storage due to their advantages of high volumetric/ gravimetric hydrogen storage capacity and low cost. However, the use of these materials in fuel cells for automotive applications at the present time is limited by high hydrogenation temperature and sluggish sorption kinetics. This paper presents the recent results of design and development of magnesium-based nanocomposites demonstrating the catalytic effects of carbon nanotubes and transition metals on hydrogen adsorption in these materials. The results are promising for the application of magnesium materials for hydrogen storage, with significantly reduced absorption temperatures and enhanced ab/desorption kinetics. High level Density Functional Theory calculations support the analysis of the hydrogenation mechanisms by revealing the detailed atomic and molecular interactions that underpin the catalytic roles of incorporated carbon and titanium, providing clear guidance for further design and development of such materials with better hydrogen storage properties.

Oxoketene-oxoketene, imidoylketene-imidoylketene and oxoketenimine-imidoylketene rearrangements. 1,3-shifts of phenyl groups

Dibenzoylketene 5 undergoes degenerate 1,3-shifts of the phenyl group between acyl and ketene carbon atoms, thus interconverting it with 6 and 7. This 1,3-shift takes place in the gas phase under flash vacuum thermolysis (FVT) conditions, but not in solution at 110-145 degrees C. Imidoyl(benzoyl)ketene 13 undergoes degenerate 1,3-shift of the phenyl group on FVT, thus interconverting it with 14, but the ketenimine isomer 15 is not formed, and none of these shifts take place in the solid state at 250 degrees C. Imidoyl(p-toluoyl)ketene 21 undergoes a 1,3-p-tolyl shift, interconverting it with ketene 22 but not with ketenimine 23. The imidoyl(p-toluoyl)ketene rotamer 25 cyclizes to 4-toluoyloxyquinoline 28 and 4-quinolone 29. The cyclization of imidoyl(benzoyl)ketene 13 to 4-benzoyloxyquinoline 18, and of 25 to 28 involves 1,3-C-to-O shifts of benzoyl (toluoyl) groups. Calculations of the transition states for the transformations at the B3LYP/6-31G** level of theory are in agreement with the observed reaction preferences.

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 H-2 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 H-2 molecules, respectively). Additionally, a molecular adsorption state of H-2 above the Ti atom is observed for the first time and is attributed to the polarization of the H-2 molecule by the Ti cation. Our results parallel recent findings for H-2 adsorption on Ti-doped carbon nanotubes or fullerenes. They provide new insight into the preliminary stages of hydrogen adsorption onto Ti-incorporated Mg surfaces.

Vibrational dynamics and structure of natural feldspars

The present study gives a contribution to the knowledge on the Na-feldspar and plagioclases, extending the database of the Raman spectra of plagioclases with different chemical compositions and structural orders. This information may be used for the future planetary explorations by “rovers”, for the investigation of ceramics nanocrystal materials and for the mineralogical phase identification in sediments. Na-feldspar and plagioclase solid solution have been investigated by Raman spectroscopy in order to determine the relationships between the vibrational changes and the plagioclase crystal chemistry and structure. We focused on the Raman micro-spectroscopy technique, being a non-destructive method, suited for contactless analysis with high spatial resolution. Chemical and structural analyses have been performed on natural samples to test the usefulness of Raman spectroscopy as a tool in the study of the pressure-induced structural deformations, the disordering processes due to change in the Al-Si distribution in the tetrahedral sites and, finally, in the determination of the anorthitic content (Anx) in plagioclase minerals. All the predicted 39 Ag Raman active modes have been identified and assigned to specific patterns of atomic vibrational motion. A detailed comparison between experimental and computed Raman spectra has been performed and previous assignments have been revised, solving some discrepancies reported in recent literature. The ab initio calculation at the hybrid HF/DFT level with the WC1LYP Hamiltonian has proven to give excellent agreement between calculated and experimentally measured Raman wavenumbers and intensities in triclinic minerals. A short digression on the 36 infrared active modes of Na-feldspar has been done too. The identification of all 39 computed Raman modes in the experimentally measured spectra of the fully ordered Na-feldspar, known as low albite, along with the detailed description of each vibrational mode, has been essential to extend the comparative analysis to the high pressure and high temperature structural forms of albite, which reflect the physical–chemical conditions of the hosting rocks. The understanding of feldspar structure response to pressure and temperature is crucial in order to constrain crustal behaviour. The compressional behaviour of the Na-feldspar has been investigated for the first time by Raman spectroscopy. The absence of phase transitions and the occurrence of two secondary compression mechanisms acting at different pressures have been confirmed. Moreover, Raman data suggest that the internal structural changes are confined to a small pressure interval, localized around 6 GPa, not spread out from 4 to 8 GPa as suggested by previous X-rays studies on elasticity. The dominant compression mechanisms act via tetrahedral tilting, while the T-O bond lengths remain nearly constant at moderate compressional regimes. At the spectroscopic level, this leads to the strong pressure dependencies of T-O-T bending modes, as found for the four modes at 478, 508, 578 and 815 cm-1. The Al-Si distribution in the tetrahedral sites affects also the Raman spectrum of Na-feldspar. In particular, peak broadening is more sensitive than peak position to changes in the degree of order. Raman spectroscopy is found to be a good probe for local ordering, in particular being sensitive to the first annealing steps, when the macroscopic order parameter is still high. Even though Raman data are scattered and there are outliers in the estimated values of the degree of order, the average peak linewidths of the Na-feldspar characteristic doublet band, labelled here as υa and υb, as a function of the order parameter Qod show interesting trends: both peak linewidths linearly increase until saturation. From Qod values lower than 0.6, peak broadening is no more affected by the Al-Si distribution. Moreover, the disordering process is found to be heterogeneous. SC-XRD and Raman data have suggested an inter-crystalline inhomogeneity of the samples, i.e., the presence of regions with different defect density on the micrometric scale. Finally, the influence of Ca-Na substitution in the plagioclase Raman spectra has been investigated. Raman spectra have been collected on a series of well characterized natural, low structural plagioclases. The variations of the Raman modes as a function of the chemical composition and the structural order have been determined. The number of the observed Raman bands at each composition gives information about the unit-cell symmetry: moving away from the C1 structures, the number of the Raman bands enhances, as the number of formula units in the unit cell increases. The modification from an “albite-like” Raman spectrum to a more “anorthite-like” spectrum occurs from sample An78 onwards, which coincides with the appearance of c reflections in the diffraction patterns of the samples. The evolution of the Raman bands υa and υb displays two changes in slope at ~An45 and ~An75: the first one occurs between e2 and e1 plagioclases, the latter separates e1 and I1 plagioclases with only b reflections in their diffraction patterns from I1 and P1 samples having b and c reflections too. The first variation represents exactly the e2→e1 phase transitions, whereas the second one corresponds in good approximation to the C1→I1 transition, which has been determined at ~An70 by previous works. The I1→P1 phase transition in the anorthite-rich side of the solid solution is not highlighted in the collected Raman spectra. Variations in peak broadening provide insights into the behaviour of the order parameter on a local scale, suggesting an increase in the structural disorder within the solid solution, as the structures have to incorporate more Al atoms to balance the change from monovalent to divalent cations. All the information acquired on these natural plagioclases has been used to produce a protocol able to give a preliminary estimation of the chemical composition of an unknown plagioclase from its Raman spectrum. Two calibration curves, one for albite-rich plagioclases and the other one for the anorthite-rich plagioclases, have been proposed by relating the peak linewidth of the most intense Raman band υa and the An content. It has been pointed out that the dependence of the composition from the linewidth can be obtained only for low structural plagioclases with a degree of order not far away from the references. The proposed tool has been tested on three mineralogical samples, two of meteoric origin and one of volcanic origin. Chemical compositions by Raman spectroscopy compare well, within an error of about 10%, with those obtained by elemental techniques. Further analyses on plagioclases with unknown composition will be necessary to validate the suggested method and introduce it as routine tool for the determination of the chemical composition from Raman data in planetary missions.

Models of gas-liquid solubilities

A recent method for phase equilibria, the AGAPE method, has been used to predict activity coefficients and excess Gibbs energy for binary mixtures with good accuracy. The theory, based on a generalised London potential (GLP), accounts for intermolecular attractive forces. Unlike existing prediction methods, for example UNIFAC, the AGAPE method uses only information derived from accessible experimental data and molecular information for pure components. Presently, the AGAPE method has some limitations, namely that the mixtures must consist of small, non-polar compounds with no hydrogen bonding, at low moderate pressures and at conditions below the critical conditions of the components. Distinction between vapour-liquid equilibria and gas-liquid solubility is rather arbitrary and it seems reasonable to extend these ideas to solubility. The AGAPE model uses a molecular lattice-based mixing rule. By judicious use of computer programs a methodology was created to examine a body of experimental gas-liquid solubility data for gases such as carbon dioxide, propane, n-butane or sulphur hexafluoride which all have critical temperatures a little above 298 K dissolved in benzene, cyclo-hexane and methanol. Within this methodology the value of the GLP as an ab initio combining rule for such solutes in very dilute solutions in a variety of liquids has been tested. Using the GLP as a mixing rule involves the computation of rotationally averaged interactions between the constituent atoms, and new calculations have had to be made to discover the magnitude of the unlike pair interactions. These numbers have been seen as significant in their own right in the context of the behaviour of infinitely-dilute solutions. A method for extending this treatment to "permanent" gases has also been developed. The findings from the GLP method and from the more general AGAPE approach have been examined in the context of other models for gas-liquid solubility, both "classical" and contemporary, in particular those derived from equations-of-state methods and from reference solvent methods.

Formation of Polycyclic Aromatic Hydrocarbons and Nitrogen Containing Polycyclic Aromatic Compounds in Titan's Atmosphere, the Interstellar Medium and Combustion

Several different mechanisms leading to the formation of (substituted) naphthalene and azanaphthalenes were examined using theoretical quantum chemical calculations. As a result, a series of novel synthetic routes to Polycyclic Aromatic Hydrocarbons (PAHs) and Nitrogen Containing Polycyclic Aromatic Compounds (N-PACs) have been proposed. On Earth, these aromatic compounds originate from incomplete combustion and are released into our environment, where they are known to be major pollutants, often with carcinogenic properties. In the atmosphere of a Saturn's moon Titan, these PAH and N-PACs are believed to play a critical role in organic haze formation, as well as acting as chemical precursors to biologically relevant molecules. The theoretical calculations were performed by employing the ab initio G3(MP2,CC)/B3LYP/6-311G** method to effectively probe the Potential Energy Surfaces (PES) relevant to the PAH and N-PAC formation. Following the construction of the PES, Rice-Ramsperger-Kassel-Markus (RRKM) theory was used to evaluate all unimolecular rate constants as a function of collision energy under single-collision conditions. Branching ratios were then evaluated by solving phenomenological rate expressions for the various product concentrations. The most viable pathways to PAH and N-PAC formation were found to be those where the initial attack by the ethynyl (C2H) or cyano (CN) radical toward a unsaturated hydrocarbon molecule led to the formation of an intermediate which could not effectively lose a hydrogen atom. It is not until ring cyclization has occurred, that hydrogen elimination leads to a closed shell product. By quenching the possibility of the initial hydrogen atom elimination, one of the most competitive processes preventing the PAH or N-PAC formation was avoided, and the PAH or N-PAC formation was allowed to proceed. It is concluded that these considerations should be taken into account when attempting to explore any other potential routes towards aromatic compounds in cold environments, such as on Titan or in the interstellar medium.