Engineering topics

Auerbach, Scott M.

Publications and source records attributed to Auerbach, Scott M..

7 records

DOE OSTI2024

Data Science Shows that Entropy Correlates with Accelerated Zeolite Crystallization in Monte Carlo Simulations

We have performed a data science study of Monte Carlo simulation trajectories to understand factors that can accelerate formation of zeolite nanoporous crystals, a process that can take days or even weeks. In previous work, Monte Carlo simulations predicted and experiments confirmed that using a secondary organic structure-directing agent (OSDA) accelerates crystallization of all-silica LTA zeolite, with experiments finding a three-fold speedup [PCCP 24, 142-148 (2022)]. However, it remains unclear what physical factors cause the speed-up. Here, we apply data science to analyze the simulation trajectories to discover what drives accelerated zeolite crystallization in Monte Carlo going from a one-OSDA synthesis (1OSDA) to a two-OSDA version (2OSDA). We encoded simulation snapshots using the Smooth Overlap of Atomic Positions approach, which represents all 2- and 3-body correlations within a given cutoff distance. Principal component analyses failed to discriminate datasets of structures from 1OSDA and 2OSDA simulations, while the Support Vector Machine (SVM) approach succeeded at classifying such structures with an area-under-curve (AUC) score of 0.99 (where AUC = 1 is a perfect classification) with all 3-body correlations, and as high as 0.94 with only 2-body correlations. SVM decision functions reveal relatively broad / narrow histograms for 1OSDA / 2OSDA datasets, suggesting that the two simulations differ strongly in information heterogeneity. Informed by these results, we performed pair (2-body) entropy calculations during crystallization, resulting in entropy differences that semi-quantitatively account for the speedup observed in the previous Monte Carlo simulations. We conclude that altering synthesis conditions in ways that substantially changes the entropy of labile silica networks may accelerate zeolite crystallization, and we discuss possible approaches for achieving such acceleration.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2024

Removing Fluoride from Double Four-Membered Rings Yielding Defect-Free Zeolites under Mild Conditions Using Ozone

We have investigated ozone treatment of as-made LTA zeolites under mild temperature conditions (175 °C) using experiments and periodic DFT as a method of energy savings and engineering defects such as silanol nests in comparison with conventional calcination at 550 °C. We have studied ozone treatment on LTA samples synthesized with 1,2-dimethyl-3-(4-methylbenzyl) imidazolium (denoted as “BULKY”) as the primary organic structure-directing agent (OSDA) and with various amounts of tetramethylammonium (TMA) as a secondary OSDA. Ozone treatment of LTA-BULKY at 175 °C was found to give defect-free, pristine LTA materials as determined by 29 Si NMR, 13 C NMR, Raman spectra, and DFT to assign the spectra. This represents a significant and unexpected finding: that fluoride ions can be completely removed from double four-membered rings (D4Rs) under such mild conditions. Furthermore, ozone treatment of LTA-BULKY-TMA samples removed BULKY but left behind TMA/F, giving a new and more diverse structural landscape of Si environments in LTA. Ab initio MetaDynamics calculations provide pathways with relatively low barriers, explaining how fluoride ions can be removed from D4Rs, leaving behind defect-free LTA materials under mild conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Titrating Controlled Defects into Si-LTA Zeolite Crystals Using Multiple Organic Structure-Directing Agents

Controlling defects in zeolites is crucial for tuning their adsorption and catalytic properties. In this work, we have performed an integrated zeolite synthesis, spectroscopy, and density functional theory study to test the limit of F– as a charge-balancing agent that mitigates defects in siliceous zeolites. We focused on the synthesis of siliceous zeolite LTA at 150 °C with 1,2-dimethyl-3-(4-methylbenzyl) imidazolium as the primary organic structure-directing agent (OSDA) and tetramethyl ammonium (TMA) as the secondary OSDA. By varying the amount of TMA in the synthesis gel, positive charges were titrated into the resulting as-made Si-LTA. Surprisingly, we found that greater TMA concentration does not induce more F– to enter into the zeolite. 29Si solid-state MAS NMR, Raman spectroscopy, and density functional theory suggest that this system has surpassed its capacity for F– to balance OSDA charge, and additional positive charge is balanced by Si–O– framework defects. The number of defects in the as-made Si-LTA can be precisely titrated by the amount of TMA in the zeolite structures. For the Si-LTA synthesized without TMA, framework defects formed in the early stage of crystal growth were found to heal during later crystallization, leading to defect-free Si-LTA. However, for the Si-LTA synthesized with TMA, the defects formed in early stages do not heal. A DFT thermodynamic analysis explains that crowding of Si-LTA pores by TMA impedes defect healing; this prediction is corroborated by synthesis experiments at an elevated temperature (170 °C). These results indicate that F– can have a limited capacity to balance OSDA charge in zeolite synthesis, opening up a third route to zeolite synthesis intermediate between the fluoride and hydroxide routes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Monte carlo simulations and experiments of all-silica zeolite LTA assembly combining structure directing agents that match cage sizes

In this work, we investigated the influence of organic structure-directing agents (OSDAs) on the formation rates of all-silica zeolite LTA using both simulations and experiments, to shed light on the crystallization process. We compared syntheses using one OSDA with a diameter close to the size of the large cavity in LTA, and two OSDAs of diameters matching the sizes of both the small and large LTA cavities. Reaction-ensemble Monte Carlo (RxMC) simulations predict a speed up of LTA formation using two OSDAs matching the LTA pore sizes; this qualitative result is confirmed by experimental studies of crystallization kinetics, which find a speedup in all-silica LTA crystallization of a factor of 3. Analyses of simulated rings and their Si–O–Si angular energies during RxMC crystallizations show that all ring sizes in the faster crystallization exhibit lower angular energies, on average, than in the slower crystallization, explaining the origin of the speedup through packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Shape-Selective Synthesis of Alkylcyclopentenyl Cations in Zeolites and Spectroscopic Distinction of Constitutional Isomers

Alkylcyclopentenyl cations belong to the long-lived intermediates that make up the “hydrocarbon pool” during the catalytic conversion of methanol on zeolites, and recent works show that such cations contribute to olefin and aromatics formation. From liquid phase chemistry, two types of alkylcyclopentenyl cations are known and distinguished by the substituent at the central carbon (C-2) of the allylic system: a more stable type with a methyl group at the C-2 and a less stable type with hydrogen at the C-2. The following three linked objectives are pursued: (i) IR spectroscopic distinction of the different substitution patterns at the allylic system of alkylcyclopentenyl cations, (ii) the role of the zeolite framework in determining the substitution pattern, and (iii) identification of alkylcyclopentenyl cations from precursors relevant to methanol-to-olefins conversion. UV–vis and IR spectroscopy are applied in situ to characterize alkylcyclopentenyl cations produced by adsorption of a pentaalkylcyclopentadiene or by adsorption and thermally induced cyclization of 2,6-dimethyl-2,4,6-octatriene. Prior knowledge of electronic spectra is combined with DFT-computed and experimental IR spectra to establish the frequency of the C–H vibration of the C-2 hydrogen-substituted type, and a characteristic red shift of the asymmetric allylic stretching vibration (Δν ≈ -20 to -30 cm –1 ) after replacing hydrogen by methyl at the central carbon of the allylic system. Although DFT demonstrates that both types of ions fit into medium- and large-pore zeolites and that the C-2 methyl-substituted type is thermodynamically favored even in the pores of the considered zeolites, formation of alkylcyclopentenyl cations is found by UV–vis and IR spectra to be shape-selective. The bulkier C-2 methyl-substituted type is detected in large-pore zeolites (MOR, BEA) and in the intersections of medium-pore zeolites (MFI), whereas in channels of medium size (TON), the less bulky C-2 hydrogen-substituted type is exclusively formed. Finally, the type of ion formed and its overall size are dictated by the zeolite framework and are independent of the precursor; the same type of alkylcyclopentenyl cation as found through cyclization of dimethyloctatriene could be generated from ethene.

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↗

DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗