Filtering by
- All Subjects: Density Functional Theory
- All Subjects: Molecular Dynamics
- Creators: Muhich, Christopher
- Creators: Singh, Arunima
Anatase exists often in its reduced form (TiO2-x), enabling it to perform redox reactions through the absorption and release of oxygen into/from the crystal lattice. These processes result in structural changes, induced by defects in the material, which can theoretically be observed using advanced characterization methods. In situ electron microscopy is one of such methods, and can provide a window into these structural changes. However, in order to interpret the structural evolution caused by defects in materials, it is often necessary and pertinent to use atomistic simulations to compare the experimental images with models.
In this thesis project, we modeled the defect structures in anatase, around oxygen vacancies and at surfaces, using molecular dynamics, benchmarked with density functional theory. Using a “reactive” forcefield designed for the simulation of interactions between anatase and water that can model and treat bonding through the use of bond orders, different vacancy structures were analyzed and simulated. To compare these theoretical, generated models with experimental data, the “multislice approach” to TEM image simulation was used. We investigated a series of different vacancy configurations and surfaces and generated fingerprints for comparison with TEM experiments. This comparison demonstrated a proof of concept for a technique suggesting the possibility for the identification of oxygen vacancy structures directly from TEM images. This research aims to improve our atomic-level understanding of oxide materials, by providing a methodology for the analysis of vacancy formation from very subtle phenomena in TEM images.
Project 1 Abstract: Ethylene Oxychlorination
The current two-step process for the industrial process of vinyl chloride production involves CuCl2 catalyzed ethylene oxychlorination to ethylene dichloride followed by thermal cracking of the latter to vinyl chloride. To date, no industrial application of a one-step process is available. To close this gap, this work evaluates a wide range of self-prepared supported CeO2 and EuOCl catalysts for one-step production of vinyl chloride from ethylene in a fixed-bed reactor at 623 773 K and 1 bar using feed ratios of C2H4:HCl:O2:Ar:He = 3:3 6:1.5 6:3:82 89.5. Among all studied systems, CeO2/ZrO2 and CeO2/Zeolite MS show the highest activity but suffer from severe combustion of ethylene, forming COx, while 20 wt.% EuOCl/γ-Al2O3 leads to the best vinyl chloride selectivity of 87% at 15.6% C2H4 conversion with complete suppression of CO2 formation and only 4% selectivity to CO conversion for over 100 h on stream. Characterization by XRD and EDX mapping reveals that much of the Eu is present in non-active phases such as Al2Eu or EuAl4, indicating that alternative synthesis methods could be employed to better utilize the metal. A linear relationship between conversion and metal loading is found for this catalyst, indicating that always part of the used Eu is available as EuOCl, while the rest forms inactive europium aluminate species. Zeolite-supported EuOCl slightly outperforms EuOCl/γ Al2O3 in terms of total yield, but is prone to significant coking and is unstable. Even though a lot of Eu seems locked in inactive species on EuOCl/γ Al2O3, these results indicate possible savings of nearly 16,000 USD per kg of catalyst compared to a bulk EuOCl catalyst. These very promising findings constitute a crucial step for process intensification of polyvinyl chloride production and exploring the potential of supported EuOCl catalysts in industrially-relevant reactions.
Project 2 Abstract: Alkyne Semihydrogenation
Despite strongly suffering from poor noble metal utilization and a highly toxic selectivity modifier (Pb), the archetypal catalyst applied for the three-phase alkyne semihydrogenation, the Pb-doped Pd/CaCO3 (Lindlar catalyst), is still being utilized at industrial level. Inspired by the very recent strategies involving the modification of Pd with p-block elements (i.e., S), this work extrapolates the concept by preparing crystalline metal phosphides with controlled stoichiometry. To develop an affordable and environmentally-friendly alternative to traditional hydrogenation catalysts, nickel, a metal belonging to the same group as Pd and capable of splitting molecular hydrogen has been selected. Herein, a simple two-step synthesis procedure involving nontoxic precursors was used to synthesize bulk nickel phosphides with different stoichiometries (Ni2P, Ni5P4, and Ni12P5) by controlling the P:Ni ratios. To uncover structural and surface features, this catalyst family is characterized with an array of methods including X-ray diffraction (XRD), 31P magic-angle nuclear magnetic resonance (MAS-NMR) spectroscopy, and X-ray photoelectron spectroscopy (XPS). Bulk-sensitive techniques prove the successful preparation of pure phases while XPS analysis unravels the facile passivation occurring at the NixPy surface that persists even after reductive treatment. To assess the characteristic surface fingerprints of these materials, Ar sputtering was carried out at different penetration depths, reveling the presence of Ni+ and P-species. Continuous-flow three-phase hydrogenations of short-chain acetylenic compounds display that the oxidized layer covering the surface is reduced under reaction conditions, as evidenced by the induction period before reaching the steady state performance. To assess the impact of the phosphidation treatment on catalytic performance, the catalysts were benchmarked against a commercial Ni/SiO2-Al2O3 sample. While Ni/SiO2-Al2O3 presents very low selectivity to the alkene (the selectivity is about 10% at full conversion) attributed to the well-known tendency of naked nickel nanoparticles to form hydrides, the performance of nickel phosphides is highly selective and independent of P:Ni ratio. In line with previous findings on PdxS, kinetic tests indicate the occurrence of a dual-site mechanism where the alkyne and hydrogen do not compete for the same site.
This work is the subject of a publication of which I am a co-author, as cited below.
D. Albani; K. Karajovic; B. Tata; Q. Li; S. Mitchell; N. López; J. Pérez-Ramírez. Ensemble Design in Nickel Phosphide Catalysts for Alkyne Semi-Hydrogenation. ChemCatChem 2019. doi.org/10.1002/cctc.201801430
Graph neural networks (GNN) offer a potential method of bypassing the Kohn-Sham equations in density functional theory (DFT) calculations by learning both the Hohenberg-Kohn (HK) mapping of electron density to energy, allowing for calculations of much larger atomic systems and time scales and enabling large-scale MD simulations with DFT-level accuracy. In this work, we investigate the feasibility of GNNs to learn the HK map from the external potential approximated as Gaussians to the electron density 𝑛(𝑟), and the mapping from 𝑛(𝑟) to the energy density 𝑒(𝑟) using Pytorch Geometric. We develop a graph representation for densities on radial grid points and determine that a k-nearest neighbor algorithm for determining node connections is an effective approach compared to a distance cutoff model, having an average graph size of 6.31 MB and 32.0 MB for datasets with 𝑘 = 10 and 𝑘 = 50 respectively. Furthermore, we develop two GNNs in Pytorch Geometric, and demonstrate a decrease in training losses for a 𝑛(𝑟) to 𝑒(𝑟) of 8.52 · 10^14 and 3.10 · 10^14 for 𝑘 = 10 and 𝑘 = 20 datasets respectively, suggesting the model could be further trained and optimized to learn the electron density to energy functional.
Using DFT calculations and GAMESS computational software, porphine and its derivatives were analyzed for unique sites to accept the adsorbates As(III), As(V) and P(V) in order to compare resulting adsorption energies and determine if any of these molecules prefer arsenic oxyanions over phosphate. Pure porphine preferred As(III) over P(V) with a resulting adsorption energy of -0.7974 eV. Of the functionalized porphyrins tested, carboxyl porphyrin preferred As(V) over P(V) with a total adsorption energy of -0.7345 eV. Ethyl, methyl, chlorine and amino porphyrin all preferred As(III), with energies of -0.7934, -0.8239, -0.7602, and -0.8508 eV, respectively. Of the metalated porphyrins tested, copper and vanadium porphyrin preferred As(V) over P(V) with adsorption energies of -0.7645 and -2.0915 eV. Chromium, iron and magnesium porphyrin all preferred As(III) over P(V) with energies of -0.5993, -1.4539, and - 1.0790 eV, respectively.
The homogenous mixture of 1-butyl-3-methylimidazolium iodide and water is examined extensively through MD as well as physical characterization of properties. Molecular ordering within the liquid mixture is related to macroscopic properties. These mixtures are then used as the basis of an electrolyte with unusual characteristics, specifically a wide liquid temperature range with an extremely low lower bound combined with relatively low viscosity allowing excellent performance in the MET sensor. Electrolyte performance is further improved by the addition of fullerene nanoparticles, which dramatically increase device sensitivity. The reasons behind this effect are explored by testing the effect of graphene surface size and through MD simulations of fullerene and a silica nanoparticle (for contrast) in [BMIM][I]/water mixtures.
Immiscible ionic liquid/water systems are explored through MD studies of particles at IL/water interfaces. By increasing the concentration of hydrophobic nanoparticles at the IL/water interface, one study discovers the formation of a commingled IL/water/particle pseudo-phase, and relates this discovery to previously-observed unique behaviors of these interfaces, particularly spontaneous particle transport across the interface. The other study demonstrates that IL hydrophobicity can influence the deformation of thermo-responsive soft particles at the liquid/liquid interface.
A proof of principle demonstration of a resonator that can switch from a high-Q “on state” to a low-Q “off state” at reduced temperatures is demonstrated in (Al1-xFex)2O3 and La(Al1-xFex)O3. The Fe3+ ions are in a high spin state (S=5/2) and undergo electron paramagnetic resonance absorption transitions that increase the microwave loss of the system. Transitions occur between mJ states with a corresponding change in the angular momentum, J, by ±ħ (i.e., ΔmJ=±1) at small magnetic fields. The paramagnetic ions also have an influence on the dielectric and magnetic properties, which I explore in these systems along with another low loss complex perovskite material, Ca[(Al1-xFex)1/2Nb1/2]O3. I describe what constitutes an optimal microwave loss switchable material induced from EPR transitions and the mechanisms associated with the key properties.
As a first step to modeling the properties of high-performance microwave host lattices and ultimately their performance at microwave frequencies, a first-principles approach is used to determine the structural phase stability of various complex perovskites with a range of tolerance factors at 0 K and finite temperatures. By understanding the correct structural phases of these complex perovskites, the temperature coefficient of resonant frequency can be better predicted.
A strong understanding of these parameters is expected to open the possibility to produce new types of high-performance switchable filters, time domain MIMO’s, multiplexers, and demultiplexers.
Overall, this work not only provides insights into the structure-property relationship of 2D pentagonal materials and opens up a new route of studying 2D materials by combining geometry and computational materials science, but also shows the potential applications of 2D pentagonal materials in electronic and magnetic devices.