Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plasma-based global pathway analysis to understand the chemical kinetics of plasma-assisted combustion and fuel reforming

The Global Pathway Analysis (GPA) algorithm helps analyze the chemical kinetics of complex combustion systems by identifying important global reaction pathways connecting a source species to a sink species through various important intermediate species (i.e., hub species). Here, the present work aims to extend GPA algorithm to plasma-assisted combustion and fuel reforming systems to identify the dominant global pathways in such systems at various conditions. In addition, the present study extends the ability of GPA algorithm to identify reaction cycles involving the excitation of high-concentration species (e.g., O 2 , N 2 , and fuel) to their vibrational and electronic states and the subsequent de excitation to their ground state, based on their significance on the reactivity of plasma-assisted systems in terms of gas heating and radical production. Provisions are made in the GPA algorithm to evaluate the reactivity of identified re action pathways and cycles based on the element-flux transfer (i.e., dominance), heat release, and radical production rate. The newly developed Plasma-based Global Pathway Analysis (PGPA) algorithm is then used to analyze the plasma assisted combustion of ammonia and reforming of methane. The PGPA analyses elucidated the significance of vibrational-translational cycles on the reactivity of NH 3 /air mixtures. Further, analyses on the production of NO ascribed the early reforming of NH 3 to N 2 and H 2 in impeding the production of NO during plasma-assisted NH 3 ignition. Lastly, the enhanced reforming of CH 4 /N 2 mixtures using plasma has been attributed to electron impact dissociation of CH 4 when compared to thermal reforming. In contrast, conventional path-Flux analysis (PFA) was found to require significant manual effort and pre-analysis intuitions from expert knowledge, making it arduous to provide valuable in sights into plasma chemistry. The user-friendly and automated nature of PGPA thus provides a valuable tool for assessing the kinetics of plasma-assisted systems helpful in analyzing and, further, a foundation in reducing plasma-assisted chemistry, without the needs of expert knowledge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

A Python Multiscale Thermochemistry Toolbox (pMuTT) for thermochemical and kinetic parameter estimation

Estimating the thermochemical properties of systems is important in many fields such as material science and catalysis. The Python multiscale thermochemistry toolbox (pMuTT) is a Python software library developed to streamline the conversion of ab-initio data to thermochemical properties using statistical mechanics, to perform thermodynamic analysis, and to create input files for kinetic modeling software. Its open-source implementation in Python leverages existing scientific codes, encourages users to write scripts for their needs, and allows the code to be expanded easily. The core classes developed include a statistical mechanical model in which energy modes can be included or excluded to suit the application, empirical models for rapid thermodynamic property estimation, and a reaction model to calculate kinetic parameters or changes in thermodynamic properties. In addition, pMuTT supports other features, such as Brønsted–Evans–Polanyi (BEP) relationships, coverage effects, and ab-initio phase diagrams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗
DOE OSTI2020

Kinetic study of hydrogen transport in graphite under molten fluoride salt environment

Here, in this work, a kinetic model, which describes the reactions occurring during the hydrogen charging process on a graphite surface, was optimized on Electrochemical Impedance Spectroscopic (EIS) to study the entry of hydrogen into graphite in high temperature molten FLiNaK salt. The results show that the rate constants of k 1 increases with decreasing applied potential (becoming more negative), and k 2 and k 3 are independent of the applied potential. At the same charging potential, all the three rate constants are independent of the water concentration in the FLiNaK salt. The surface coverage of absorbed hydrogen increases with decreasing the charging potential and increasing the moisture content in the salt. The adsorption efficiency of hydrogen increases with increasingly charging potential and decreases with increasing moisture content of the FLiNaK salt. These results demonstrate that EIS is a promising technique for measuring the kinetic parameters for hydrogen transport in graphite under molten fluoride salt environment.

EIS↗
DOE OSTI2025

Geologic hydrogen as an emerging fuel: experimental insights, thermodynamics, kinetics, and reactive transport modeling

Geologic hydrogen (GeoH 2 ) is emerging as a viable clean energy source. It is largely produced through serpentinization, a geological process in which ultramafic rocks react with water under suitable temperature and pressure. Here, this review synthesizes the current understanding of H 2 generation by serpentinization, with an emphasis on reaction mechanisms, kinetics, and thermodynamics, as well as on modeling flow and transport of reacting fluids in geological formations. We describe the role of mineral assemblages, such as olivine and pyroxene, fluid-rock interactions, and catalytic surfaces, in influencing GeoH 2 yield and reaction rates. By integrating models of reaction kinetics, subsurface reactive flow and transport, and the serpentinization process, and by accounting for the thermodynamic state of the system, this review aims to guide future GeoH 2 research and to evaluate the potential of natural hydrogen as a sustainable clean energy source.

Moradi, Rasoul [Univ. of Southern California, Los ↗
DOE OSTI2014

Kinetics and Mechanisms of Cadmium Carbonate Heteroepitaxial Growth at the Calcite (101¯4) Surface

Elucidating the kinetics and mechanisms of heteroepitaxial nucleation and growth at mineral-water interfaces is essential to understanding surface reactivity in geochemical systems. In the present work, the formation of heteroepitaxial cadmium carbonate coatings at calcite-water interfaces was investigated by exposing calcite (10-14) surfaces to Cd-bearing aqueous solutions. In situ atomic force microscopy (AFM) was employed as the primary technique. The AFM results indicate that the heteroepitaxial growth of cadmium carbonate proceeds via three different mechanisms depending on the initial supersaturation of the aqueous solution: advancement of existing steps, nucleation and growth of three-dimensional (3D) islands, and nucleation and spread of two-dimensional (2D) nuclei. The 3D islands and 2D nuclei exhibit different morphologies and growth kinetics. The effects of supersaturation on heteroepitaxial growth mechanisms can be interpreted in terms of the free energy barrier for nucleation. At low initial supersaturation, where 3D nucleation dominates, it is hypothesized, from the growth rate and morphology of the 3D islands observed with AFM, that the crystallization of the overgrowth follows a non-classical pathway involving the formation of a surface precursor that is not fully crystalline, whereas high supersaturation favors the formation of crystalline 2D nuclei whose morphology is based on the atomic structure of the calcite substrate. Cross-sectional transmission electron microscopy (TEM) images reveal that the atomic structure of the interface between the cadmium carbonate coating and calcite shows perfect, dislocation-free epitaxy.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2017

Neon diffusion kinetics and implications for cosmogenic neon paleothermometry in feldspars

Observations of cosmogenic neon concentrations in feldspars can potentially be used to constrain the surface exposure duration or surface temperature history of geologic samples. The applicability of cosmogenic neon to either application depends on the temperature-dependent diffusivity of neon isotopes. Here in this work, we investigate the kinetics of neon diffusion in feldspars of different compositions and geologic origins through stepwise degassing experiments on single, proton-irradiated crystals. To understand the potential causes of complex diffusion behavior that is sometimes manifest as nonlinearity in Arrhenius plots, we compare our results to argon stepwise degassing experiments previously conducted on the same feldspars. Many of the feldspars we studied exhibit linear Arrhenius behavior for neon whereas argon degassing from the same feldspars did not. This suggests that nonlinear behavior in argon experiments is an artifact of structural changes during laboratory heating. However, other feldspars that we examined exhibit nonlinear Arrhenius behavior for neon diffusion at temperatures far below any known structural changes, which suggests that some preexisting material property is responsible for the complex behavior. In general, neon diffusion kinetics vary widely across the different feldspars studied, with estimated activation energies (E a ) ranging from 83.3 to 110.7 kJ/mol and apparent pre-exponential factors (D 0 ) spanning three orders of magnitude from 2.4 ×10 -3 to 8.9 × 10 -1 cm 2 s -1 . Finally, as a consequence of this variability, the ability to reconstruct temperatures or exposure durations from cosmogenic neon abundances will depend on both the specific feldspar and the surface temperature conditions at the geologic site of interest.

58 GEOSCIENCES↗
DOE OSTI2022

Kinetic and equilibrium reactions on natural and laboratory generation of thermogenic gases from Type II marine shale

The phenomenon that laboratory pyrolysis experiments produce much wetter gases than those in natural reservoirs is a long-recognized and debated problem in the investigation of natural gases in sedimentary basins. In this study, we explore the discrepancy by pyrolyzing a type II kerogen from the Woodford Shale in Oklahoma, compared with the previous results on the produced natural gases from the Arkoma Basin generated from the same source rock (Liu et al., 2019) with the discussion of gas and isotopic compositions at bulk and position-specific (PS) levels. An improved GC-pyrolysis-GC IRMS method is applied for the determination of PS δ 13 C of propane produced in the pyrolysis of the Woodford Shale at Easy %R o from 0.76 to 3.27. Kinetic and thermodynamic considerations of the chemical and isotopic compositions of the natural and laboratory pyrolysis gases suggest that the generation of light hydrocarbons involves uni-directional cracking reactions, exchange reactions with water, and likely reversible reactions among light hydrocarbons and other H-containing volatiles. After the gas generation in the unconventional Woodford Shale reservoirs, the C 1 -C 4 gases might have approached close to chemical equilibrium of C 1 -C 3 and isotope equilibrium of C 2 -C 1 and C 3 -C 1 pairs at their peak temperatures. The capping H for the generation of C 1 -C 4 in the Woodford Shale gases appears to have experienced at least partial exchange with the water, while that in the pyrolysis gases is only originated from organic-bound compounds with large kinetic isotope effects (KIE). Our findings indicate that elevated compound-specific and PS δ 13 C values of propane in the wet-gas cracking stage are significantly influenced by the breakdown of the thermally stable compounds (e.g., remaining kerogen, residues). A first synthesis of PS δ 13 C and δ 2 H isotopic compositions of propane from this study and the literature data suggests relatively similar isotopic structures of propane precursors in kerogens. Finally, this study demonstrates that PS isotope analysis of propane can contribute to identifying various geological (e.g., maturation, wet-gas cracking, H exchange, diffusion) and biodegradation processes.

58 GEOSCIENCES↗
DOE OSTI2020

Interpenetration and kinetic effects in converging, high-energy plasma jets

In this work, we report on numerical simulations of laser-driven convergent plasma fusion targets. These “inverted corona” fusion targets are useful for the study of counter-streaming and converging rarefied plasma flows, and previous experiments have demonstrated their potential as neutron sources. The scheme consists of a fuel layer lined along the interior surface of a hollow plastic shell that is laser-ablated and expands inward towards the target center. The plasma streams generated in these targets are initially nearly collisionless as they converge, leading to wide interaction length scales and long interaction time scales as the jets interpenetrate. Such kinetic effects impact mixing of constituent ions - a phenomenon not properly captured by single-fluid hydrodynamic simulations. Here we conduct numerical simulations using two different methods: (1) single-fluid simulations in HYDRA, and (2) kinetic-ion, fluid-electron hybrid particle-in-cell (PIC) simulations in the code Chicago. It is shown that the initially nearly collisionless plasma fronts interpenetrate deeply and lead to broader interaction regions in space and time resulting in significant beam-beam fusion. The two approaches make different, testable predictions for the effect of fuel-layer thickness on neutron yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Thermally-activated non-Schmid glide of screw dislocations in W using atomistically-informed kinetic Monte Carlo simulations

Thermally-activated 1/2<111> screw dislocation motion is the controlling plastic mechanism at low temperatures in body-centered cubic (bcc) crystals. Dislocation motion proceeds by nucleation and propagation of atomic-sized kink pairs in close-packed planes. The atomistic character of kink pairs can be studied using techniques such as molecular dynamics (MD). However, MD’s natural inability to properly sample thermally-activated processes as well as to capture {110} screw dislocation glide calls for the development of other methods capable of overcoming these limitations. Here we develop a kinetic Monte Carlo (kMC) approach to study single screw dislocation dynamics from room temperature to 0.5T m and at stresses 0 < σ < 0.9σ P , where Tm and σ P are the melting point and the Peierls stress. The method is entirely parameterized with atomistic simulations using an embedded atom potential for tungsten. To increase the physical fidelity of our simulations, we calculate the deviations from Schmid’s law prescribed by the interatomic potential used and we study single dislocation kinetics using both projections. We calculate dislocation velocities as a function of stress, temperature, and dislocation line length. We find that considering non-Schmid effects has a strong influence on both the magnitude of the velocities and the trajectories followed by the dislocation. We finish by condensing all the calculated data into effective stress and temperature dependent mobilities to be used in more homogenized numerical methods.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gap junctions fine-tune ganglion cell signals to equalize response kinetics within a given electrically coupled array

Retinal ganglion cells (RGCs) summate inputs and forward a spike train code to the brain in the form of either maintained spiking (sustained) or a quickly decaying brief spike burst (transient). We report diverse response transience values across the RGC population and, contrary to the conventional transient/sustained scheme, responses with intermediary characteristics are the most abundant. Pharmacological tests showed that besides GABAergic inhibition, gap junction (GJ)–mediated excitation also plays a pivotal role in shaping response transience and thus visual coding. More precisely GJs connecting RGCs to nearby amacrine and RGCs play a defining role in the process. These GJs equalize kinetic features, including the response transience of transient OFF alpha (tOFFα) RGCs across a coupled array. We propose that GJs in other coupled neuron ensembles in the brain are also critical in the harmonization of response kinetics to enhance the population code and suit a corresponding task.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Carbidic Mo is the sole kinetically-relevant active site for catalytic methane dehydroaromatization on Mo/H-ZSM-5

Systematic variation of Mo/Al = {0.10, 0.25, 0.35} in Mo/H-ZSM-5 catalysts effects ~20× change of Mo-to-H + ratio during methane dehydroaromatization (DHA) reactions at steady state and on deactivating catalysts. Mo and Brønsted-acid site count, respectively enumerated by ethane hydrogenolysis probe reaction and dimethyl ether chemical titration, evolve disparately during DHA catalysis, indicating the two sites reduce in number by distinct mechanisms. Here, the deactivation of Mo/H-ZSM-5 catalytic activity is uniquely attributed to the loss of active, or accessible, Mo sites, evinced by (i) the non-selective nature of deactivation (i.e. that deactivation occurs without modification of residual active sites), (ii) invariance in product distribution from Mo/H + = 0.10–2.1, and (iii) linear correspondence between Mo-catalyzed ethane hydrogenolysis rate and active site count during DHA on deactivating samples. DHA forward rate – the intrinsic kinetic descriptor of catalyst activity – normalizes by Mo content at steady state and on deactivating catalysts with Mo/Al = 0.10–0.35 for Mo-based contact times in the range τ Mo = 0.79–44 mol Mo s mol –1 C , unequivocally establishing Mo aggregates as the sole kinetically-relevant active site on Mo/H-ZSM-5. Brønsted-acid site content over a ~40× range in τ H+ = 1.4–54 mol H+ s mol –1 C has no discernable correlation with or effect on DHA rate, methane conversion, or product distribution – demonstrating H + does not catalyze any rate- or selectivity-determining steps in the benzene formation pathway and suggesting that the benefits of zeolitic acid sites are only to disperse Mo during catalyst synthesis and, conceivably, to catalyze equilibrated reaction steps.

02 PETROLEUM↗
DOE OSTI2021

H 2 O-assisted O 2 reduction by H 2 on Pt and PtAu bimetallic nanoparticles: Influences of composition and reactant coverages on kinetic regimes, rates, and selectivities

Hydrogen peroxide (H 2 O 2 ) can replace hazardous oxidants in industrial processes but is currently too expensive for many such applications. While direct synthesis of H 2 O 2 (H 2 + O 2 → H 2 O 2 ) may reduce costs in comparison to incumbent technology, current catalysts lack the requisite stability and selectivity. Here, we examine the direct synthesis of H 2 O 2 on bimetallic Pt 1 Au x (0 ≤ x ≤ 230) and Pt catalysts at steady-state in pure water and relate kinetic parameters for H 2 O 2 and H 2 O formation to possible active site structures informed by complementary characterization methods. X-ray photoelectron spectra show significant Pt surface enrichment compared to the bulk composition. Analysis of infrared spectra of mixed monolayers of 12 CO* and 13 CO* indicate that Pt and Au form substitutional surface alloys. The Pt 1 Au x nanoparticles with the greatest mole fractions of Au predominantly expose Pt monomers (i.e., isolated Pt atoms), yet Pt atoms exposed upon all these nanoparticles possess electronic structures distinct from bulk Pt. Despite these differences, rate measurements are consistent with product formation through proton-electron transfer pathways for all Pt 1 Au x catalysts. In situ XAS indicate that Pt remains metallic during H 2 O 2 synthesis. Under the most oxidizing conditions, selectivities toward H 2 O 2 increase strongly with the Au to Pt ratio from 2% for monometallic Pt to 85% for Pt 1 Au 170 . However, selectivities are similar among all catalysts within reducing conditions. Comparisons of apparent activation enthalpies for the formation of H 2 O 2 and H 2 O across these catalysts and the range of conditions suggest that Pt monomers within Au provide the greatest selectivities for H 2 O 2 formation, because these active sites present high barriers for O-O bond rupture. Further, selectivities decrease with increasing ratios of H 2 to O 2 pressures, because Pt atoms aggregate and form oligomers that readily dissociate dioxygen intermediates. The combined use of spectroscopy, kinetics, and concepts employed in reaching these conclusions take inspiration from the legacy of Prof. Michel Boudart, and specifically his elegant methods for interrogating bimetallic catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Adsorption kinetics, thermodynamics, and isotherm studies for functionalized lanthanide-chelating resins

Conventional ion exchange resins are widely utilized to remove metals from aqueous solutions, but their limited selectivity precludes dilute ion extraction. This research investigated the adsorption performance of ligand-functionalized resins towards rare earth elements (REE). Functionalized resin particles were synthesized by grafting different ligands (diethylenetriaminepentaacetic dianhydride (DTPADA), phosphonoacetic acid (PAA), or N,N-bis(phosphonomethyl)glycine (BPG)) onto pre-aminated polymeric adsorbents (diameter-0.6 mm). Lanthanide uptake trends were evaluated for the functionalized resins using batch adsorption experiments with a mixture of three REEs (Nd, Gd, and Ho at 0.1–1000 mg/L each). Resin physical-chemical properties were determined by measuring their surface area, ligand concentrations, and acidity constants. The aminated supports contained 4.0 mmol/g primary amines, and ligand densities for the functionalized resins were 0.33 mmol/g (PAA), 0.22 mmol/g (BPG), and 0.42 mmol/g (DTPADA). Kinetic studies revealed that the functionalized resins followed pseudo-second order binding kinetics with rates limited by intraparticle diffusion. Capacity estimates for total REE adsorption based on Langmuir q Max were 0.12 mg/g (amine; ≈ 0.77 µmol/g), 5.0 mg/g (PAA; ≈ 32.16 µmol/g), 3.0 mg/g (BPG; ≈ 19.30 µmol/g), and 2.9 mg/g (DTPADA; ≈ 18.65 µmol/g). Lastly, attaching ligands to the aminated resins greatly improved their REE binding strength and adsorption efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

A multidimensional unified gas-kinetic scheme for radiative transfer equations on unstructured mesh

In order to extend the unified gas kinetic scheme (UGKS) to solve radiative transfer equations in a complex geometry, a multidimensional asymptotic preserving implicit method on unstructured mesh is constructed in this paper. With an implicit formulation, the CFL condition for the determination of the time step in UGKS can be much relaxed, and a large time step is used in simulations. Differently from previous direction-by-direction UGKS on orthogonal structured mesh, on unstructured mesh the interface flux transport takes into account multi-dimensional effect, where gradients of radiation intensity and material temperature in both normal and tangential directions of a cell interface are included in the flux evaluation. The multiple scale nature makes the UGKS be able to capture the solutions in both optically thin and thick regions seamlessly. In the optically thick region the condition of cell size being less than photon's mean free path is fully removed, and the UGKS recovers a solver for diffusion equation in such a limit on unstructured mesh. For a distorted quadrilateral mesh, the UGKS goes to a nine-point scheme for the diffusion equation, and it naturally reduces to the standard five-point scheme for a orthogonal quadrilateral mesh. Numerical computations covering a wide range of transport regimes on unstructured and distorted quadrilateral meshes will be presented to validate the current approach. - Highlights: • A unified gas-kinetic scheme (UGKS) on unstructured mesh has been developed in all flow regimes. • The scheme solves diffusion equation in optical thick region with large cell size. • The implicit UGKS removes the CFD condition. The efficiency of the scheme has been much improved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The role of amorphous ZIF in ZIF-8 crystallization kinetics and morphology

Understanding the composition and structure of amorphous precursor phases is fundamental for elucidating two-step crystallization mechanisms and designing shape- and size-controlled nanomaterials. However, that understanding is largely lacking for metal–organic framework compounds despite their growing significance as functional materials. Here, in this study, we report the crystallization of zeolite imidazolate frameworks-8 (ZIF-8, Zn(C 4 H 5 N 2 ) 2 ) via an amorphous ZIF (am-ZIF) solid precursor phase with a rough stoichiometric composition of Zn(C 4 H 5 N 2 ) 1.78 (C 4 H 6 N 2 ) 0.17 (CH 3 COO) 0.22 . The formation of am-ZIF is attributed to the incomplete deprotonation of 2-Methylimidazole (HmIm) and the involvement of the hydrogen bond between CH 3 COO– and –HN, which can further transform into the dense Dia(Zn) structure with a diamondoid crystal topology in pure water. Taking am-ZIF as a precursor, the tunable dissolution and recrystallization kinetics of am-ZIF into ZIF-8, due to the addition of EtOH and CTAB, allows the selective fabrication of dodecahedral, cubic, and hollow ZIF-8. Overall, an in-depth understanding of the differences in composition and structure of am-ZIF from ZIF-8 and the resulting crystallization kinetics suggests a novel approach to designing metal–organic frameworks with controlled crystal morphology.

36 MATERIALS SCIENCE↗