Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “supercell”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-principles study of Al/Al 3 Ni interfaces

Al-Ni alloys have shown promise for high-temperature applications due to the strengthening of Al 3 Ni fibers resistant to coarsening and spheroidization up to 400°C. While the interface between Al and Al 3 Ni phases affects the coarsening rate of Al 3 Ni at elevated temperatures, its characteristics are largely unknown to date. Here, we have constructed various supercells to model this interface and performed a first-principles study based on density functional theory (DFT). We have considered three groups of Al/Al 3 Ni interfaces: experimentally reported orientation relationships from the solidification studies, crystallographically similar Fe-Fe 3 C pearlite interfaces, and the family of low-index (100) termination planes. We have analyzed the correlation between the DFT Al/Al 3 Ni interfacial energies and characteristic features, e.g., excess free volume and the number of broken bonds. We outline the further experimental and computational analysis required to improve the interface modeling of Al/Al 3 Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Point defect energetics in gallium arsenide, a comprehensive density functional theory study

In materials, point defects often control or modify functional properties. To predict the performance of materials intended for application in optoelectronic devices, it is imperative to understand the properties of those point defects. For the first time, all six intrinsic defects of GaAs, a key optoelectronics material, and their charge transition levels are calculated using density functional theory with the HSE06 functional. For comparison, both PBE and r 2 SCAN calculations are also carried out. The HSE06 results are found to be in better agreement with experimental data than previous calculations. In conclusion, the importance of using the exact electron exchange present in hybrid functionals and larger supercells to accurately determine defect levels and ground state defect configurations is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine learning guided prediction of solute segregation at coherent and semi-coherent metal/oxide interfaces

Investigation of semi-coherent metal/oxide interfaces with misfit dislocations using density functional theory (DFT) is computationally intensive to the point of being prohibitive, as it involves several hundreds to many thousands of atoms. In this study, we examined the solute segregation behavior at the Fe/Y 2 O 3 interface—a model interface for cladding applications in nuclear fission reactors—using a combination of DFT calculations and machine learning (ML) approaches. Both coherent and semi-coherent interfaces were considered. ML models were trained on DFT-calculated segregation energies to identify the key chemical, geometric and strain energy related features that govern solute segregation behavior at coherent Fe/Y 2 O 3 interfaces. Furthermore, it was found that ML models when trained on DFT calculated segregation energy of elements at a coherent interface, comprising of about a hundred-atom supercell, can predict the segregation energy of elements at a semi-coherent Fe/Y 2 O 3 interface (with multiple hundreds of atoms) at a fraction of computational cost (1/35th), with an accuracy comparable to DFT calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2025

A microstructural signature of the coesite-quartz transformation: New insights from high-pressure experiments and EBSD

Ultra-high pressure (UHP) metamorphism is difficult to identify in continental crust as few petrological barometers are suitable for dominantly felsic lithologies. In such cases, burial to extreme depths is commonly identified through the preservation of coesite, a high-pressure polymorph of SiO 2 that typically forms at depths exceeding ∼ 100 km (i.e., > 2 GPa pressure). Unfortunately, coesite readily transforms to quartz upon exhumation, meaning that UHP terranes may often be overlooked. While some studies have suggested that quartz may inherit an orientation signature indicative of former coesite, both the specific nature of this signature and the conditions favouring its development remain uncertain. Here, to address this problem, we combine electron backscatter diffraction analysis of natural and experimental samples to explore microstructural evolution across the coesite-quartz phase transformation. We demonstrate that neighbouring domains of quartz commonly feature an 84 ± 4° rotation of [c] axes around the pole of a common {m} plane. This orientation relationship is a product of epitaxy, whereby the {$11\bar{2}2$} Japan twin plane in quartz nucleates on the (010) plane in coesite. In supercell simulations, the nucleation of Japan twins can be explained by the energetically favourable alignment of quartz tetrahedra on parental coesite tetrahedra. Through experiments, we demonstrate that this signature emerges over a broad range of conditions, regardless of the availability of nucleation sites (e.g., grain boundaries) or the density of crystal lattice defects (e.g., dislocations). Overall, our work provides a quantitative and unambiguous tool for identifying UHP terranes from quartz in isolation.

Coesite↗
DOE OSTI2021

Neutron scattering study of tantalum monohydride and monodeuteride

Powder samples of TaH 0.89 and TaD0.96 are synthesized under a hydrogen (deuterium) pressure of 2.8 GPa and a temperature of 250 °C, then quenched to the liquid nitrogen temperature, recovered to ambient pressure and studied by neutron diffraction (ND) and inelastic neutron scattering (INS). The ND study shows that both hydrogen and deuterium atoms occupy tetrahedral interstitial sites in a distorted body centered cubic (bcc) crystal structure of metal atoms, while the ordering scenarios in TaH 0.89 and TaD 0.96 are different. Hydrogen and deuterium atoms are ordered in a layered fashion, forming long period superstructures with space groups P$\bar{4}$ and P222, respectively, so that the unit cells of the crystal structures of TaH 0.89 and TaD 0.96 are $\sqrt{2}$ × $\sqrt{2}$ × 7 and $\sqrt{2}$ × $\sqrt{2}$ × 8 supercells of the initial cubic unit cell. The INS study demonstrates a pronounced “soft” (trumpet-like) anharmonicity of the potential well for H and D atoms.

08 HYDROGEN↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Optical vibrational spectra and proposed crystal structure of ε-UO 3

ε-UO 3 is an exotic polymorph in the uranium trioxide system with an undetermined crystal structure and limited optical vibrational spectroscopic data. To improve understanding of this compound, we synthesize and investigate the crystal structure and optical vibrational spectra of ε-UO 3 . Infrared spectra collected for ε-UO 3 are in good agreement with previously published results, and our studies extend the available data into the low-energy (600–100 cm –1 ) regime. For the first time, Raman spectra are presented for ε-UO 3 using both 785 and 532 nm excitation wavelengths. Previous reports suggest an impurity phase may be present in ε-UO 3 produced by calcination of U 3 O 8 ; however, spectral center-of-mass calculations, principal component analyses, and Raman spectroscopic mapping employed to investigate this possibility indicate that the product of U 3 O 8 calcined in O 3 (g) in this work is likely phase-pure. A possible novel structure solution for ε-UO 3 is determined via Rietveld refinement of powder X-ray diffraction data and is triclinic, P-1, with a = 4.01 Å, b = 3.85 Å, c = 4.18 Å, and α = 98.26°, β = 90.41°, γ = 120.46° (R wp = 8.30%). The asymmetric unit of ε-UO 3 consists of U(VI) in hexagonal bipyramidal coordination with displaced equatorial oxygen. Further analysis reveals that the structure of ε-UO3 is best described by a 2 × 1 × 2 supercell structure in P-1 with a = 8.03 Å, b = 3.86 Å, c = 8.37 Å with α = 98.26°, β = 90.41°, and γ = 120.46°, although a higher-symmetry structure is possible. Optical vibrational spectroscopic and structural measurements of ε-UO 3 presented here furthers our understanding of this complex uranium oxide and clarifies the origin of reported structural similarity to U 3 O 8 .

-UO3↗
DOE OSTI2022

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the broadband vibration isolation performance of nonlocal total-internal-reflection metasurfaces

The concept of a nonlocal elastic metasurface has been recently proposed and experimentally demonstrated in Zhu et al. (2020). When implemented in the form of a total-internal-reflection (TIR) interface, the metasurface can act as an elastic wave barrier that is impenetrable to deep subwavelength waves over an exceptionally wide frequency band. The underlying physical mechanism capable of delivering this broadband subwavelength performance relies on an intentionally nonlocal design that leverages long-range connections between the units forming the fundamental supercell. This paper explores the design and application of a nonlocal TIR metasurface to achieve broadband passive vibration isolation in a structural assembly made of multiple dissimilar elastic waveguides. The specific structural system comprises shell, plate, and beam waveguides, and can be seen as a prototypical structure emulating mechanical assemblies of practical interest for many engineering applications. The study also reports the results of an experimental investigation that confirms the significant vibration isolation capabilities afforded by the embedded nonlocal TIR metasurface. We report these results are particularly remarkable because they show that the performance of the nonlocal metasurface is preserved when applied to a complex structural assembly and under non-ideal incidence conditions of the incoming wave, hence significantly extending the validity of the results presented in Zhu et al. (2020). Results also confirm that, under proper conditions, the original concept of a planar metasurface can be morphed into a curved interface while still preserving full wave control capabilities.

42 ENGINEERING↗
DOE OSTI2023

Accurate prediction of short-range order and its effect on thermodynamic, structural, and electronic properties of disordered alloys: exemplified in archetypical Cu 3 Au

Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solid breeder material crystal structure evolution due to Li burn up–Loss of crystal stability

Solid breeder materials for fusion power plants are subjected to extreme conditions of radiation damage, high temperature and designed component burn-up. Here, the loss of lattice stability of Li 2 TiO 3 due to Li burn-up was investigated using first principles-based approaches to aid evaluation of upper limit of component lifetime. The lattice stability is analyzed using calculated phonon dispersion in the Li 2 TiO 3 supercell, while Li burn-up is modeled by the introduction of Li-vacancies. It has been determined that the studied ceramic can be structurally stable up to ~40% of Li atoms burn-up. At 50% burn-up, the negative frequency branch in the phonon dispersion spectrum appears indicating the loss of lattice dynamical stability. Moreover, the structure obtained by introduction of unstable frozen phonons with minima energy amplitude also results in a phonon dispersion with negative branches in the structure obtained. Therefore, the system completely loses stability when approximately half of Li atoms are consumed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anionic Oxygen Redox in the High-Lithium Material Li 8 SnO 6

Lithiated transition metal oxide cathode materials that combine oxygen-anion redox [O 2 – /(O 2 ) n– ] with cation redox offer substantial capacities and higher voltage than cathodes without oxygen redox, but oxygen-anion redox is usually accompanied by the formation of peroxo or superoxo species that may cause oxygen release. To clarify the relationship between the oxygen release and peroxide and superoxide dimers without the complication of transition-metal redox, we performed density functional theory calculations to study the layered Li 8 SnO 6 cathode because it has only oxygen-anion redox during delithiation. Features with O–O distances of ~1.5 Å (corresponding to the formation of a peroxide) were observed at 75% lithium concentration, and further delithiation to 62.5% lithium concentration shortened the O–O bond distance to ~1.3 Å (corresponding to the formation of superoxide). The solely anionic oxygen redox stabilizes the delithiated structure of Li 8–x SnO 6 , and it does not cause structural disorder or release oxygen gas. Ab initio molecular dynamics calculations in the supercell at 300 K indicate that oxygen dimerization may occur at 75% lithium concentration. On the (0001) surface with a high lithium concentration (98%), oxygen dimerization is unfavorable. The calculations by HSE06 show that oxygen redox in Li 8 SnO 6 offers an initial voltage (vs Li/Li + ) greater than 4.2 V. The GGA+U calculations indicate that the Li + -polaron in Li 95 Sn 12 O 72 migrates with a migration barrier of only 0.43 eV. Here, this work also quantifies a solely anionic oxygen redox mechanism in Li 8 SnO 6 , and it provides a deeper understanding of anionic oxygen redox in Li-excess oxide cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

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

Identifying Order and Disorder in Double Four-Membered Rings via Raman Spectroscopy during Crystallization of LTA Zeolite

Fluoride (F – ) has been essential for the synthesis of low-defect siliceous zeolites. It has been hypothesized that F – balances the positive charges from organic structure-directing agents (OSDAs) and stabilizes key building units during zeolite crystallization such as the double four-membered ring (D4R). However, due to the lack of characterization techniques for investigating medium-range structures, including rings and cages formed during zeolite crystallization, the roles of F – in stabilizing building units and maintaining local charge balance during zeolite assembly are not yet fully understood. Here, the crystallization of siliceous Linde type A (LTA) zeolite in the presence of F – was investigated using Raman spectroscopy and periodic density functional theory (DFT) calculations. We have discovered that the F – -filled double four-membered ring (F – /D4R) and the empty D4R exhibit rather distinct Raman features. Both the F – /D4R and empty D4R are formed in the LTA zeolite synthesized in the presence of F – using 1,2-dimethyl-3-(4-methylbenzyl) imidazolium as the OSDA. Observed Raman bands of the F – /D4R and empty D4R, along with predictive DFT calculations on LTA supercells, reveal an ordered distribution of these two D4R units in the final as-made LTA zeolite. The discovery of these distinct Raman signatures of F – /D4R and empty D4R units opens an interesting new window for studying defects in the D4R during zeolite formation. In particular, we have observed variation in Raman intensities of F – /D4R and empty D4R bands during LTA crystallization; periodic DFT calculations indicate that the observed Raman behavior is consistent with empty D4R units containing one or two Si vacancies surrounded by Q 3 Si—consistent also with solid-state nuclear magnetic resonance measurement. These defects appear to heal during further crystallization, leading to the formation of defect-free LTA zeolite crystals. Altogether, our results provide deeper understanding on the roles of F – in charge balancing and stabilizing intact D4R units during zeolite formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗