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.

Effect of reaction mechanism on precursor exposure time in atomic layer deposition of silicon oxide and silicon nitride

Atomic layer deposition (ALD) of highly conformal, silicon-based dielectric thin films has become necessary because of the continuing decrease in feature size in microelectronic devices. The ALD of oxides and nitrides is usually thought to be mechanistically similar, but plasma-enhanced ALD of silicon nitride is found to be problematic, while that of silicon oxide is straightforward. To find why, the ALD of silicon nitride and silicon oxide dielectric films was studied by applying ab initio methods to theoretical models for proposed surface reaction mechanisms. The thermodynamic energies for the elimination of functional groups from different silicon precursors reacting with simple model molecules were calculated using density functional theory (DFT), explaining the lower reactivity of precursors toward the deposition of silicon nitride relative to silicon oxide seen in experiments, but not explaining the trends between precursors. Using more realistic cluster models of amine and hydroxyl covered surfaces, the structures and energies were calculated of reaction pathways for chemisorption of different silicon precursors via functional group elimination, with more success. DFT calculations identified the initial physisorption step as crucial toward deposition and this step was thus used to predict the ALD reactivity of a range of amino-silane precursors, yielding good agreement with experiment. The retention of hydrogen within silicon nitride films but not in silicon oxide observed in FTIR spectra was accounted for by the theoretical calculations and helped verify the application of the model.

Low valence cation doping of bulk Cr2O3: Charge compensation and oxygen vacancy formation

The different oxidation states of chromium allow its bulk oxide form to be reducible, facilitating the oxygen vacancy formation process, which is a key property in applications such as catalysis. Similar to other useful oxides such as TiO2, and CeO2, the effect of substitutional metal dopants in bulk Cr2O3 and its effect on the electronic structure and oxygen vacancy formation are of interest, particularly in enhancing the latter. In this paper, density functional theory (DFT) calculations with a Hubbard + U correction (DFT+U) applied to the Cr 3d and O 2p states, are carried out on pure and metal-doped bulk Cr2O3 to examine the effect of doping on the electronic and geometric structure. The role of dopants in enhancing the reducibility of Cr2O3 is examined to promote oxygen vacancy formation. The dopants are Mg, Cu, Ni, and Zn, which have a formal +2 oxidation state in their bulk oxides. Given this difference in host and, dopant oxidation states, we show that to predict the correct ground state two metal dopants charge compensated with an oxygen vacancy are required. The second oxygen atom removed is termed "the active" oxygen vacancy and it is the energy required to remove this atom that is related to the reduction process. In all cases, we find that substitutional doping improves the oxygen vacancy formation of bulk Cr2O3 by lowering the energy cost.

Ultrafast charge-transfer in organic photovoltaic interfaces: geometrical and functionalization effects

Understanding the microscopic mechanisms of electronic excitation in organic photovoltaic cells is a challenging problem in the design of efficient devices capable of performing sunlight harvesting. Here we develop and apply an ab initio approach based on time-dependent density functional theory and Ehrenfest dynamics to investigate photoinduced charge transfer in small organic molecules. Our calculations include mixed quantum–classical dynamics with ions moving classically and electrons quantum mechanically, where no experimental external parameter other than the material geometry is required. We show that the behavior of photocarriers in zinc phthalocyanine (ZnPc) and C60 systems, an effective prototype system for organic solar cells, is sensitive to the atomic orientation of the donor and the acceptor units as well as the functionalization of covalent molecules at the interface. In particular, configurations with the ZnPc molecules facing on C60 facilitate charge transfer between substrate and molecules that occurs within 200 fs. In contrast, configurations where ZnPc is tilted above C60 present extremely low carrier injection efficiency even at longer times as an effect of the larger interfacial potential level offset and higher energetic barrier between the donor and acceptor molecules. An enhancement of charge injection into C60 at shorter times is observed as binding groups connect ZnPc and C60 in a dyad system. Our results demonstrate a promising way of designing and controlling photoinduced charge transfer on the atomic level in organic devices that would lead to efficient carrier separation and maximize device performance.

Epitaxially Grown Strained Pentacene Thin Film on Graphene Membrane

Organic-graphene system has emerged as a new platform for various applications such as flexible organic photovoltaics and organic light emitting diodes. Due to its important implication in charge transport, the study and reliable control of molecular packing structures at the graphene-molecule interface are of great importance for successful incorporation of graphene in related organic devices. Here, an ideal membrane of suspended graphene as a molecular assembly template is utilized to investigate thin-film epitaxial behaviors. Using transmission electron microscopy, two distinct molecular packing structures of pentacene on graphene are found. One observed packing structure is similar to the well-known bulk-phase, which adapts a face-on molecular orientation on graphene substrate. On the other hand, a rare polymorph of pentacene crystal, which shows significant strain along the c-axis, is identified. In particular, the strained film exhibits a specific molecular orientation and a strong azimuthal correlation with underlying graphene. Through ab initio electronic structure calculations, including van der Waals interactions, the unusual polymorph is attributed to the strong graphene-pentacene interaction. The observed strained organic film growth on graphene demonstrates the possibility to tune molecular packing via graphene-molecule interactions.

Experimental and theoretical investigation of ligand effects on the synthesis of ZnO nanoparticles

ZnO nanoparticles with highly controllable particle sizes(less than 10 nm) were synthesized using organic capping ligands in Zn(Ac)2 ethanolic solution. The molecular structure of the ligands was found to have significant influence on the particle size. The multi-functional molecule tris(hydroxymethyl)-aminomethane (THMA) favoured smaller particle distributions compared with ligands possessing long hydrocarbon chains that are more frequently employed. The adsorption of capping ligands on ZnnOn crystal nuclei (where n = 4 or 18 molecular clusters of(0001) ZnO surfaces) was modelled by ab initio methods at the density functional theory (DFT) level. For the molecules examined, chemisorption proceeded via the formation of Zn...O, Zn...N, or Zn...S chemical bonds between the ligands and active Zn2+ sites on ZnO surfaces. The DFT results indicated that THMA binds more strongly to the ZnO surface than other ligands, suggesting that this molecule is very effective at stabilizing ZnO nanoparticle surfaces. This study, therefore, provides new insight into the correlation between the molecular structure of capping ligands and the morphology of metal oxide nanostructures formed in their presence.