As an important chemical raw material, alcohols can be used as fuels, solvents and chemical feedstocks to produce a variety of downstream products. With limited fossil fuel resources, alcohol synthesis from syngas reactions can be a potential alternative to the traditional petroleum based alcohol synthesis. Among many catalysts active for syngas to alcohol processes, alkali promoted Mo2C has shown promising performance. More interestingly, the alkali promoter was found to play an important role in shifting the reaction selectivity from hydrocarbons to alcohols. However, limited understanding of the mechanism of this alkali promoter effect is available due to the complexity of syngas reaction mechanism and low content of alkali added to the catalysts. In this thesis, we performed a comprehensive investigation of the alkali promoter effect with density functional theory (DFT) calculations as our primary tool. We first examine various Mo2C surfaces to determine a representative surface structure active to alkali adsorption. On this particular surface, we develop a syngas reaction network including relevant reaction mechanisms proposed in previous literature. With energetics derived from DFT calculations and a BEP relation, we predict the syngas reaction selectivity and find it to be in excellent agreement with experimental results. The dominant reaction mechanism and selectivity determining steps are determined from sensitivity analysis. We also propose a formation mechanism of alkali promoters on Mo2C catalysts that shows consistency between experimental IR and DFT computed vibrational frequencies. Finally, the effect of alkali promoters on the selectivity determining steps for syngas reactions are investigated from DFT calculations and charge analysis. We are able to rationalize the role of alkali promoters in shifting the reaction selectivity from hydrocarbons to alcohols on Mo2C catalysts.

CO hydrogenation and CO oxidation are two important processes addressing the energy and environmental issues of great interest. Both processes are carried out using metallic catalysts. The objective of this dissertation is to study the catalytic processes that govern these two reactions from a molecular perspective using quantum mechanical calculations. Density Functional Theory (DFT) has proven to be a valuable tool to study adsorption, dissociation, chain growth, reaction pathways etc., on well-defined surfaces. DFT was used to study the CO reduction reactions on promoted cobalt catalyst surfaces and CO oxidation mechanisms on cobalt surfaces. CO hydrogenation via Fischer-Tropsch Synthesis (FTS) is a process used to produce liquid fuels from synthesis gas. The economics of the Fischer-Tropsch process strongly depends on the performance of the catalyst used. The desired properties of a catalyst include selectivity towards middle distillate products such as diesel and jet fuel, higher activity and longer catalyst life. Catalysts are often modified by adding promoters to obtain these desirable properties. Promoters can influence the reaction pathways, reducibility, dispersion, activity and selectivity. In FTS, understanding the effect of promoters in the molecular scale would help in tailoring catalysts with higher activity and desired selectivity. Preventing deactivation of catalyst is important in FTS to increase the catalyst life. Deactivation of Co catalyst can occur by reoxidation, C deposition, sintering, formation of cobalt-support compounds etc. Designing catalyst with resistance to deactivation by the use of promoters is explored in this dissertation. The influence of promoters on the initiation pathways of CO hydrogenation is also explored as a first step towards determining the selectivity of promoted catalyst. The influence of Pt promoter on O removal from the surface of Co catalyst showed that Pt promoter reduced the activation barrier for the removal of O on both flat and stepped Co surfaces. An approximate kinetic model was developed and a volcano plot was established. The turn-over frequency (TOF) calculated based on the activation barriers showed that Pt promoted Co surface had a higher rate than unpromoted Co surface. The effect of Pt and Ru promoters on various pathways of C deposition on Co catalyst was studied to gain a mechanistic understanding. The promoters did not affect the subsurface C formation but they increased the barriers for C-C and C-C-C formation and also decreased the barriers for C-H formation. The promoters also influence the stabilities of C compounds on the Co surface suggesting that Pt and Ru promoters would decrease C deposition on Co catalysts. The effect of Pt promoter on unassisted and H-assisted CO activation pathways on Co catalyst was studied. Pt promoted Co surface followed H-assisted CO activation. Pt promoter decreased the activation barriers for CO activation pathways on Co catalyst thereby increasing the activity of Co catalyst. CO oxidation is a process used to prevent poisoning of fuel cell catalysts and reduce pollution of the atmosphere through exhaust gases containing CO. Expensive catalysts like Pt are widely used for CO oxidation which significantly increases the cost of the process and hence it is necessary to search for alternative lower cost catalysts. Understanding the mechanism of a reaction is the first step towards designing better and efficient catalyst. DFT is helpful in determining the basic mechanism and intermediates of reactions. The mechanism of CO oxidation on CoO catalyst was explored. Four possible mechanisms for CO oxidation on CoO catalyst were studied to determine the most likely mechanism. The mechanism was found to be a two-step process with activation barrier for formation of CO2 larger than the barrier for formation of the intermediate species.

The synthesis and characterization of Pt core/ Cu shell (Pt@Cu) dendrimer-encapsulated nanoparticles (DENs) having full and partial Cu shells deposited via electrochemical underpotential deposition (UPD) is described. Pt DENs containing averages of 55, 147, and 225 Pt atoms immobilized on glassy carbon electrodes served as the substrate for the UPD of a Cu monolayer. This results in formation of Pt@Cu DENs. Evidence for this conclusion is based on results from the analysis of cyclic voltammograms (CVs) for the UPD and stripping of Cu on Pt DENs, and from experiments showing that the Pt core DENs catalyze the hydrogen evolution reaction before Cu UPD, but that after Cu UPD this reaction is inhibited. Results obtained by in-situ electrochemical X-ray absorption spectroscopy (XAS) confirm the core@shell structure. Calculations from density functional theory (DFT) show that the first portion of the Cu shell deposits onto the (100) facets, while Cu deposits lastly onto the (111) facets. The DFT-calculated energies for Cu deposition on the individual facets are in good agreement with the peaks observed in the CVs of the Cu UPD on the Pt DENs. Finally, structural analysis of Pt DENs having just partial Cu shells by in-situ XAS is consistent with the DFT-calculated model, confirming that the Cu partial shell selectively decorates the (100) facets. These results are of considerable significance because site-selective Cu deposition has not previously been shown on nanoparticles as small as DENs. In summary, the application of UPD as a synthetic route and characterization tool for core@shell DENs having well defined structures is established. A study of the degradation mechanism and degradation products of Pd DENs is provided as well. These DENs consisted of an average of 147 atoms per dendrimer. Elemental analysis and UV-vis spectroscopy indicate that there is substantial oxidation of the Pd DENs in air-saturated solutions, less oxidation in N₂-saturated solution, and no detectable oxidation when the DENs are in contact with H₂. Additionally, the stability improves when the DEN solutions are purified by dialysis to remove Pd²⁺-complexing ligands such as chloride. For the air- and N₂-saturated solutions, most of the oxidized Pd recomplexes to the interiors of the dendrimers, and a lesser percentage escapes into the surrounding solution. The propensity of Pd DENs to oxidize so easily is a likely consequence of their small size and high surface energy. Calculations from density functional theory (DFT) show that the first portion of the Cu shell deposits onto the (100) facets, while Cu deposits lastly onto the (111) facets. The DFT-calculated energies for Cu deposition on the individual facets are in good agreement with the peaks observed in the CVs of the Cu UPD on the Pt DENs. Finally, structural analysis of Pt DENs having just partial Cu shells by in-situ XAS is consistent with the DFT-calculated model, confirming that the Cu partial shell selectively decorates the (100) facets. These results are of considerable significance because site-selective Cu deposition has not previously been shown on nanoparticles as small as DENs. In summary, the application of UPD as a synthetic route and characterization tool for core@shell DENs having well defined structures is established. A study of the degradation mechanism and degradation products of Pd DENs is provided as well. These DENs consisted of an average of 147 atoms per dendrimer. Elemental analysis and UV-vis spectroscopy indicate that there is substantial oxidation of the Pd DENs in air-saturated solutions, less oxidation in N2-saturated solution, and no detectable oxidation when the DENs are in contact with H2. Additionally, the stability improves when the DEN solutions are purified by dialysis to remove Pd2+-complexing ligands such as chloride. For the air- and N2-saturated solutions, most of the oxidized Pd recomplexes to the interiors of the dendrimers, and a lesser percentage escapes into the surrounding solution. The propensity of Pd DENs to oxidize so easily is a likely consequence of their small size and high surface energy. / text

This dissertation includes work done on graphene and related materials, examining their electronic and topological properties using first-principles methods. Ab-initio computational methods, like density functional theory (DFT), have become increasingly popular in condensed matter and material science. Motivated by the search for novel materials that would help us devise fast, low-power, post-CMOS transistors, we explore the properties of some of these promising materials. We begin by studying graphene and its interaction with dielectric oxides. Graphene has recently inspired a flurry of research activity due to its interesting electronic and mechanical properties. For the device community, graphene's high charge carrier mobility and continuous gap tunability can have immense use in novel transistors. In Chapter 3 we examine the properties of graphene placed on two oxides, namely quartz and alumina. We find that oxygen-terminated quartz is a useful oxide for the purpose of graphene based FETs. Inspired by a recent surge of interest in topological insulators, we then explore the topological properties of two-dimensional materials. We conduct a theoretical study to examine the relationship between crystal space group symmetry and the electric polarization of a two-dimensional crystal. We show that the presence of symmetry restricts the polarization values to a small number of distinct groups. There groups in turn are topologically inequivalent, making polarization a topological index. We also conduct density functional theory calculations to obtain actual polarization values of materials belonging to C3 symmetry and show that our results are consistent with our theoretical analysis. Finally we prove that any transformation from one class of polarization to another is a topological phase transition. In Chapter 5 we use density functional theory to examine the electronic properties of graphene intercalation compounds. Bilayer pseudospin field effect transistor (BiSFET) has been proposed as an interesting low-power, efficient post-CMOS switch. In order to implement this device we need bilayer graphene with reduced interlayer interaction. One way of achieving that is by inserting foreign molecules between the layers, a process which is called intercalation. In this chapter we examine the electronic properties of bilayer graphene intercalated with iodine monochloride and iodine monobromide molecules. We find that intercalation of graphene indeed makes it promising for the implementation of BiSFET, by reducing interlayer interaction. As an interesting side problem, we also use hybrid, more extensive approaches in DFT, to examine the electronic and optical properties of dilute nitrides. Dilute nitrides are highly promising and interesting materials for the purposes of optoelectronic applications. Together, we hope this work helps in elucidating the electronic properties of promising material systems as well as act as a guide for experimentalists. / text

Atomically two-dimensional (2-D) graphene, as well as the hexagonal boron nitride dielectric have been and are continuing to be widely investigated for the next generation nanoelectronic devices. More recently, other 2-D materials and electronic systems including the surface states of topological insulators (TIs) and monolayers of transition metal dichalcogenides (TMDs) have also attracted considerable interest. In this work I have focused on these latter two material systems on possible device applications. TIs are characterized by an insulating bulk band gap and metallic Dirac surface states which are spin-polarized. Here, the electronic structures of bulk and thin film TIs are studied using ab-initio density functional theory (DFT). Band inversion, an essential characteristic of TIs, is shown in the bulk band structures. Properties of TI surface bands in thin film such as the critical film thickness to induce a gap, the thickness dependent gap size, and the localization length of surface states are reported. Effects of crystalline dielectric materials on TI surface states are also addressed by ab-initio calculations. I discuss the sensitivity of Dirac point degeneracy and linear band dispersion of TI with respect to different dielectric surface terminations as well as different relative atom positions of the dielectric and TI. Additionally, this work presents research on exciton condensation in TI using a tight-binding model combined with self-consistent non-local Hartree-Fock mean-field theory. Possibility of exciton condensation in the TI Bi₂Se₃ thin film is assessed. Non-equilibrium Green's function (NEGF) simulations with the atomistic tight-binding (TB) Hamiltonian are carried out to explore the performance of metal-oxide-semiconductor field-effect-transistor (MOSFET) and tunnel field-effect-transistor (TFET) based on the Bi₂Se₃ TI thin film. How the high dielectric constant of Bi₂Se₃ affects the performance of MOSFET and TFET is presented. Bulk TMDs such as MoS₂, WS₂ and others are the van der Waals-bonded layered material, much like graphite, except monolayer (and Bulk) TMDs have a large band gap in-contrast to graphene (and graphite). Here, the performance of nanoscale monolayer MoS₂ n-channel MOSFETs are examined through NEGF simulations using an atomistic TB Hamiltonian. N- and p-channel MOSFETs of various monolayer TMDs are also compared by the same approach. I correlate the performance differences with the band structure differences. Finally, ab-initio calculations of adatom doping effects on the monolayer MoS₂ is shown. I discuss the most stable atomic configurations, the bonding type and the amount of charge transfer from adatom to the monolayer MoS₂. / text

Transition metal alloys are an important class of materials in heterogeneous catalysis due in no small part to the often greatly enhanced activity and selectivity they exhibit compared to their monometallic constituents. A host of experimental and theoretical studies have demonstrated that, in many cases, these synergistic effects can be attributed to atomic-scale features of the catalyst surface. Realizing the goal of designing -- rather than serendipitously discovering -- new alloy catalysts thus depends on our ability to predict their atomic configuration under technologically relevant conditions. This dissertation presents original research into the development and use of computational tools to accomplish this objective. These tools are all based on a similar strategy: For each of the alloy systems examined, cluster expansion (CE) Hamiltonians were constructed from the results of density functional theory (DFT) calculations, and then used in Metropolis Monte Carlo (MC) simulations to predict properties of interest. Following a detailed description of the DFT+CE+MC simulation scheme, results for the AuPd/Pd(111) and AuPt/Pt(111) surface alloys are presented. These two systems exhibit considerably different trends in their atomic arrangement, which are explicable in terms of their interatomic interactions. In AuPd, a preference for heteronuclear, Au-Pd interactions results in the preferential formation of Pd monomers and other small ensembles, while in AuPt, a preference for homonuclear interactions results in the opposite. AuPd/Pd(100) and AuPt/Pt(100) were similarly examined, revealing not only the effects of the same heteronuclear/homonuclear preferences in this facet, but also a propensity for the formation of second nearest-neighbor pairs of Pd monomers, in close agreement with experiment. Subsequent simulations of the AuPd/Pd(100) surface suggest the application of biaxial compressive strain as a means increasing the population of this catalytically important ensemble of atoms. A method to incorporate the effects of subsurface atomic configuration is also presented, using AuPd as an example. This method represents several improvements over others previously reported in the literature, especially in terms of its simplicity. Finally, we introduce the dimensionless scaled pair interaction, whereby the finite-temperature atomic configuration of any bimetallic surface alloy may be predicted from a small number of relatively inexpensive calculations. / text

Complex oxides exhibiting a wide variety of novel functional properties such as ferromagnetism and ferroelectricity have been extensively studied during the past decades. Recent advances in the field of oxide heteroepitaxy have made it possible to create and control hybrid oxide heterostructures with abrupt epitaxial interfaces. The oxide heteroepitaxy with the capability of controlling interface composition, strain, length scales, etc. has opened the totally new and exciting scientific avenue and has offered potential device applications to be explored. Epitaxial integration of functional oxides on semiconductor such as Si (001) and Ge(001) is of great interest, as it potentially leads to further technological development of these interesting oxide systems. In this dissertation, using density functional theory we explore physics and chemistry of novel oxide heterostructures and issues related to the integration of functional oxides on semiconductors. Oxide materials that are studied in this dissertation include polar LaAlO₃, high-k dielectric SrTiO₃, photocatalytic anatase TiO₂ and CoO, and strongly correlated magnetic oxide LaCoO₃. / text

In the first part of this thesis, we present results from extensive molecular dynamics simulations of the collapse transitions of hydrophobic polymers in explicit water. The focus is to understand the roles that curvature and interactions associated with the polymer-water “interface” have on collapse thermodynamics. We show that model hydrophobic polymers can have parabolic, protein-like, temperature-dependent free energies of unfolding. Analysis of the water structure shows that the polymer-water interface can be characterized as soft and weakly dewetted. We also show that an appropriately defined surface tension for the polymer-water interface is independent of the attractive polymer-water interactions. This helped us to develop a perturbation model for predicting the effect of attractions on polymer collapse thermodynamics. In the second part, we explore connections between structure, thermodynamics, and dynamics of inhomogeneous fluids. First, we use molecular dynamics simulations and classical density functional theory (DFT) to study the hard-sphere fluid at approximately 103 equilibrium state points, spanning different confining geometries and particle-boundary interactions. We provide strong empirical evidence that both excess entropy and a new generalized measure of available volume for inhomogeneous fluids correlate excellently with self-diffusivity, approximately independent of the degree of confinement. Next, we study via simulations how tuning particle-wall interactions to flatten or enhance the particle layering of a model confined fluid impacts its self-diffusivity, viscosity, and entropy. Interestingly, interactions that eliminate particle layering can significantly reduce confined fluid mobility, whereas those that enhance layering can have the opposite effect. Excess entropy helps to understand and predict these trends. Finally, we explore the relationships between the effective interparticle interactions, static structure, and tracer diffusivity of a solute in a mixture. We show that knowledge of these relationships can allow one to “tune” the effective interparticle interactions of the solute in a way that increases its tracer diffusivity. One interesting consequence is that the mobility of a hard-sphere solute can be increased by adding a soft-repulsion to its interaction, effectively making it bigger. / text

The chromogenic properties of tungsten trioxide (WO3) have been studied by photoluminescence spectroscopy at 4.2 K in order to characterize the electronic structure of this material and see how this relates to optical responses during chromogenic coloration. Transition processes between electron energy states are often the cause of optical phenomena and it is important to identify such processes in order to understand the chromogenic coloration of tungsten oxide films. Much research work has been devoted to characterize the physical and chemical mechanisms that are responsible for this coloration and this is of fundamental importance to understand the chromogenic behavior. The latest research shows that oxygen vacancies could play an important role in certain coloration processes, but it is still a matter of debate whether these are important for the overall response. This work aims to identify specific transitions that are related to oxygen vacancies by measuring photoluminescence from films with controlled vacancy content. The main goal of the project was to set up an experiment that could measure photoluminescence at liquid helium temperature. This was done by installing and integrating the components included in this experimental set-up. The films had been prepared prior to this work and were deposited on a nanocrystalline CaF2 substrate, which is a material that has a very large band gap and was therefore expected to fully transparent in the UV range. However it was found that the substrate inelastically scattered the UV excitation light, which produced strong signals that overshadowed the photoluminescence and prevented an effective characterization of the electronic structure in the films. Instead, suggestions were given on how to minimize uncertainty factors and overcome the difficulties met in this work. It was also found that the films attain a lasting blue coloration by exposure to UV light in vacuum, and that this might be due to oxygen being desorbed from the film during experiments in vacuum.

In this work, we have implemented an orbital-free density functional theory (OF-DFT) solver using the finite element method. In OF-DFT, the total ground state energy is minimized directly with respect to the electron density, rather than via orbitals like in the standard Kohn-Sham approach. For this to be possible, one needs an approximation of a universal density functional of the non-interacting kinetic energy. Presently available approximations allow for computation with very low computational expense, but which gives inaccurate energies. A stable OF-DFT code can be used as a testbed for new kinetic energy functionals and provide the necessary tool for investigating the accuracy of OF-DFT calculations for complex systems. We have implemented Thomas-Fermi theory with and without nuclear cusp condition, as well as additional exchange terms of Dirac and Amaldi. The program uses an extended version of the steepest descent in order to find the minimizing density in the variational principle. Our results include convergence tests for the hydrogen atom, weak bonding in the H2 molecule, and accurate results for the lightest noble gases (He, Ne, Ar). For heavier atoms (Kr, Xe, Rn), the results are less accurate. In addition, we consider hydrogen in the simple cubic structure without the cusp condition, which is a first attempt to use the code for periodic systems. Lastly, we discuss some possible improvements for the iterative process towards the minimizing density, as well as other possible directions for future development.