Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “formation 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 415 records · Page 23

Efficient Pd on carbon catalyst for ammonium formate dehydrogenation: Effect of surface oxygen functional groups

Formate solution is one of a kind promising liquid organic hydrogen carrier but suffers from kinetic challenges. This study investigates the kinetics of the surfacefunctionalized Pd on carbon catalysts for formate dehydrogenation and the impact of O-functional groups. Here, the fraction of the distinguished O-functional groups was modulated by the different concentrated HNO 3 solutions treatment or by H 2 reduction. This study shows that the O-functional groups play critical roles in dispersing Pd nanoparticles and decreasing the activation energy for dehydrogenation. Density functional theory (DFT) calculation reveals that most O-functional groups enhance formate adsorption on the Pd active site. However, the existence of C=O groups consumes reducing agents and hinders the formation of metallic Pd. The electron transfer from Pd to oxygen functional groups is unfavorable to dehydrogenation. The as-prepared Pd5/re-ACA (reduced activated carbon washed by acid) exhibited significant activity with a higher turnover frequency of 13,511 h -1 than commercial Pd/C.

54 ENVIRONMENTAL SCIENCES↗

Chemical effects in ion mixing of a ternary system (metal-SiO2)

The mixing of Ti, Cr, and Ni thin films with SiO2 by low-temperature (- 196-25 C) irradiation with 290 keV Xe has been investigated. Comparison of the morphology of the intermixed region and the dose dependences of net metal transport into SiO2 reveals that long range motion and phase formation probably occur as separate and sequential processes. Kinetic limitations suppress chemical effects in these systems during the initial transport process. Chemical interactions influence the subsequent phase formation.

Banwell, T.↗

A study of silver acetate under extreme conditions

With the aim of exploring chemical systems that may undergo metallization when irradiated with hard X-rays, we selected silver acetate (AgC 2 H 3 O 2 ) for the subject of this study. X-ray-induced decomposition of silver acetate under ambient and high-pressure conditions was observed in a diamond anvil cell (DAC), leading to the formation of metallic nanograins of silver at ambient pressure and 1.65 GPa. At 4 GPa, no decomposition was observed. The Avrami kinetics equation also provides information about novel structural formation at ambient pressure and 1.65 GPa. By modeling of the XRD data, it was found that the size of the silver nanocrystallites formed at 1.65 GPa pressure steadily increased to ∼5 nm after 600 min of X-ray irradiation as determined by applying the Scherrer equation to the diffraction peak widths. Time-resolved X-ray diffraction (XRD) revealed pressure-dependent kinetics, demonstrating that coupling pressure with irradiation enables controlled photochemical pathways in this model system. Concurrent with previous studies, the application of high pressure (HP) can be considered as a means of controlling X-ray synthetic photochemistry.

36 MATERIALS SCIENCE↗

Gel Precursors as Glass and Ceramic Starting Materials for Space Processing Applications Research

The crystallization kinetics and glass forming ability of reluctant glass formers is investigated. This could ultimately aid the formation of bulk samples of unique glass compositions outside of normal glass forming regions allowing the optimization of certain properties of the glass. One important aspect of processing in space is the containerless undercooling of molten substances. Theoretically, the extent of undercooling can be greatly enhanced by solidifying in the absence of heterogeneous nucleation resulting from contact with crucibles or molds. Techniques were established for the measurement of crystallization kinetics and critical cooling rates. The glass formation ability and crystallization kinetics of Ga2O-43CaO and several Al2O3-CaO compositions were measured. An apparatus was set up to measure the temperature of spherical samples on a thermocouple at large cooling rates. The time and temperature of nucleation is recorded and the probability of nucleation at various cooling rates can be measured.

Downs, R. L.↗

Isotopic Substitution Reveals the Importance of Aluminate Diffusion Dynamics in Gibbsite (Al(OH) 3 ) Crystallization from Alkaline Aqueous Solution

Understanding molecular-scale factors governing the precipitation of aluminum hydroxides, such as gibbsite, under alkaline conditions is important for the formation of laterite deposits, as well as aluminum processing. However, mechanisms enabling tetrahedral aluminate ions to assemble into octahedral sites of the gibbsite lattice remain unclear. Formation of oligomeric complexes has been hypothesized as a critical intermediate step. Here, we report a study of gibbsite solubility in highly alkaline solutions using deuterium substitution to probe equilibrium and kinetic factors that could affect oligomeric intermediate formation, including the reactivity and the diffusivity of aluminate ions. When substituting sodium hydroxide with sodium deuteroxide, solution analysis shows a nearly 40% and 50% increase in gibbsite solubility in 2.4 and 3.3 mol·kg –1 total sodium solutions, respectively. Raman spectroscopy indicated that monomeric aluminate ions are the predominant species in solution irrespective of deuteration. However, both 27 Al and 23 Na nuclear magnetic resonance (NMR) spectroscopy revealed significant differences in chemical shifts in deuterated solutions, and analysis of 1 H, 23 Na, and 27 Al diffusion coefficients with pulsed-field gradient, stimulated echo NMR spectroscopy shows a decrease in ion diffusivity with increasing total deuterium in solution, consistent with the notion of increasing strength in the hydrogen/deuterium bonding network. In conclusion, because the relative change in 1 H diffusion coefficients are commensurate with the difference in apparent solubility constants, there is evidence for an oligomeric intermediate whose steady-state concentration is maintained by the collision frequency of aluminate ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of the hydrochlorination of SiCl4

The hyrochlorination of SiC14 and m.g. silicon metal to produce SiHC13, was investigated. Reaction kinetic measurements were carried out to collect additional rate data at 525 C and 550 C. A theoretical study was carried out to provide a kinetic model and a rate equation for the hydrochlorination reaction. Results of this preliminary study show that the rate of formation of SiHC13 follows a pseudo first order kinetics. The rate constants were measured at three temperatures, 550 C, 500 C and 450 C, respectively. The activation energy was determined from the Arrhenius plot to give a value of 13.2 Kcal/mole. The design of a quartz reactor to measure reaction rates and equilibrium conversion of SiHC13 at reaction temperature up to 650 C was completed.

Mui, J. Y. P.↗

Reaction Pathways and Energy Consumption in NH 3 Decomposition for H 2 Production by Low Temperature, Atmospheric Pressure Plasma

Pathways for NH 3 decomposition to N 2 and N 2 H 4 by atmospheric pressure nonthermal plasma are analyzed using a combination of molecular beam mass spectrometry measurements and zero-dimensional kinetic modeling. Experimental measurements show that NH 3 conversion and selectivity towards N 2 formation scale monotonically with the specific energy input into the plasma with ~ 100% selectivity to N 2 formation achieved at specific energy inputs above 0.12 J cm −3 (3.1 eV (molecule NH 3 ) −1 ). The kinetic model recovers these trends, although it underpredicts N 2 selectivity at low specific energy input. These discrepancies can be explained by the underestimation of reaction rate coefficients for reactions that consume N 2 H x species in collisions with H radicals and/or radial nonuniformities in power deposition, gas temperature, and species concentrations that are not represented by the plug flow approximation used in the model. The kinetic model shows that N 2 formation proceeds through N 2 H x decomposition pathways rather than NH x decomposition pathways in low temperature, atmospheric pressure plasma. Higher selectivity toward N 2 production can be achieved by operating at higher NH 3 conversion and with a higher gas temperature. Furthermore, the high energy cost of NH 3 decomposition by atmospheric pressure nonthermal plasma found in this work (25–50 eV (molecule NH 3 converted) −1 ; 17–33 eV (molecule H 2 formed) −1 ) is a result of the energy requirement for electron-impact dissociation of NH 3 and the significant re-formation of NH 3 by three-body recombination reactions between NH 2 and H.

Nonthermal plasma↗

Cation Site Occupancy in Natural Ferroan Double Carbonates via Mössbauer and Fe K-Edge X-ray Absorption Spectroscopy

Double carbonates are minerals with a calcite type structure with alternating cation layers composed of Ca and Mg/Fe coordinated by carbonate groups. While the perfectly ordered AB-stacked crystal is the thermodynamically most stable configuration, natural mineral formation pathways can leave signatures through kinetically trapped disorder such as AB antisite cation substitutions. This study probes the degree of cation ordering in naturally occurring double carbonate samples.in particular, the distribution of Fe A and B crystallographic sites (A site is Ca layer and B site is Mg/Fe layer). Mössbauer, X-ray diffraction, X-ray absorption fine structure (both X-ray absorption near edge structure, XANES, and extended x-ray absorption fine structure, EXAFS), and energy dispersive X-ray spectroscopies were used to collect a comprehensive experimental data set, which we interpret using density functional theory to elucidate the structural effects of cation disorder. Our results show that the A (nominally Ca) site can host a relatively high Fe fraction. We discuss the implications in terms of mineral formation.

Boglaienko, Daria (ORCID:0000000160943877)↗

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗

Soot formation and precursor chemistry in Counterflow flames of aviation fuel surrogates

To meet market demands, the aviation sector is interested in utilizing drop-in Synthetic Aviation Turbine Fuels (SATF), either as neat fuels or in blends with conventional Jet A. SATF currently approved in standard specifications may have lower aromatic content with significant fractions of normal, branched, and cyclo-alkanes. Fundamental studies on soot formation from aviation fuels (Jet A, SATF) and their surrogate components are essential to understand how fuel composition influences soot and soot precursor formation. Here, this study reports new measurements of polycyclic aromatic hydrocarbons (PAH) and soot in counterflow diffusion flames (CDFs) of aviation fuel surrogates. Both intrusive and non-intrusive diagnostics are employed to determine the profiles of temperature, gas phase species, PAHs (up to C16), and soot volume fraction (SVF) in CDFs of iso-octane and surrogate mixtures. These measurements shed light on the transition of soot precursors to primary soot particles. In addition to serving as a common surrogate component in Jet A surrogate mixtures, iso-octane is a template species for larger, less volatile branched alkanes found in SATF mixtures. The newly developed Lawrence Livermore National Laboratory (LLNL) PAH and soot model successfully captures temperature, precursor species, and SVF profiles for the mixtures and conditions discussed in this work. Finally, a high-fidelity surrogate for Jet A is proposed that matches targeted physical and chemical properties well, while leveraging the wide range of candidate fuel molecules available in the LLNL detailed chemical model. The proposed surrogate formulation is validated against newly acquired measurements of the surrogate and literature measurements of Jet A. These new experiments and simulations provide critical insights into the PAH and soot formation from aviation fuels. Reaction pathways which require further investigation are highlighted, such that future work may bridge the remaining quantitative gaps in predicting soot formation from aviation fuel surrogates and surrogate components.

Aviation fuels↗

Continuum kinetic investigation of the impact of bias potentials in the current saturation regime on sheath formation

In this work, we examine sheath formation in the presence of bias potentials in the current saturation regime for pulsed power fusion experiments. It is important to understand how the particle and heat fluxes at the wall may impact the wall material and affect electrode degradation. Here, simulations are performed using the 1X-1V Boltzmann–Poisson system for a proton–electron plasma in the presence of bias potentials ranging from 0 to 10 kV. The results indicate that the sheath near the anode remains generally the same as that of a classical sheath without the presence of a bias potential. However, the sheath near the cathode becomes more prominent with a larger potential drop, a significant decrease in the electron density, and larger sheath lengths. The spatially constant current density increases to a saturation value with the increasing bias potential. For high bias potentials, the material choice needs to consider that the anode faces significantly larger particle and heat fluxes compared to the cathode. In general, the results trend with theory with differences attributed to the simplified assumptions in the theory and the kinetic effects considered in the simulations. Due to the significant computational cost of a well resolved 1X-2V simulation, only one such simulation is performed for the 5 kV case showing higher current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rational Regulation of Layer‐by‐Layer Processed Active Layer via Trimer‐Induced Pre‐Swelling Strategy for Efficient and Robust Thick‐Film Organic Solar Cells

Thick-film (>300 nm) organic solar cells (OSCs) have garnered intensifying attention due to their compatibility with commercial roll-to-roll printing technology for the large-scale continuous fabrication process. However, due to the uncontrollable donor/acceptor (D/A) arrangement in thick-film condition, the restricted exciton splitting and severe carrier traps significantly impede the photovoltaic performance and operability. For this work, combined with layer-by-layer deposition technology, a twisted 3D star-shaped trimer (BTT-Out) is synthesized to develop a trimer-induced pre-swelling (TIP) strategy, where the BTT-Out is incorporated into the buried D18 donor layer to enable the fabrication of thick-film OSCs. The integrated approach characterizations reveal that the exceptional configuration and spontaneous self-organization behavior of BTT-Out trimer could pre-swell the D18 network to facilitate the acceptor's infiltration and accelerate the formation of D/A interfaces. This enhancement triggers the elevated polarons formation with amplified hole-transfer kinetics, which is essential for the augmented exciton splitting efficiency. Furthermore, the regulated swelling process can initiate the favorable self-assembly of L8-BO acceptors, which would ameliorate carrier transport channels and mitigate carrier traps. As a result, the TIP-modified thin-film OSC devices achieve the champion performance of 20.3% (thin-film) and 18.8% (thick-film) with upgraded stability, among one of the highest performances reported of thick-film OSCs.

36 MATERIALS SCIENCE↗

Self-organization of photoionized plasmas via kinetic instabilities

Abstract Self-organization in an unmagnetized collisionless plasma (in this paper) refers to formation of transient coherent structures such as collective oscillations (electrostatic waves) or magnetic fields resulting from so-called kinetic effects in the plasma. This topical review provides a comprehensive analysis of the self-organization of strong-field photoionized, non-equilibrium plasmas through kinetic instabilities. The authors propose and demonstrate a novel experimental platform that enables the formation of dense plasmas with known highly anisotropic and non-thermal electron velocity distribution functions on a timescale on the order of an inverse electron plasma frequency. We then show that such plasmas are highly susceptible to a hierarchy of kinetic instabilities, including two-stream, current filamentation and Weibel, that convert a fraction of the electron kinetic energy into electric and/or magnetic energy stored in self-organized structures. The electrostatic waves so produced are measured using a collective light (Thomson) scattering technique with femtosecond resolution as the kinetic instabilities aided by collisions eventually thermalize the plasma electrons. In addition, we describe a novel experimental technique that has made it possible to map the temporal evolution of the wavenumber spectrum of the thermal Weibel instability with picosecond resolution, which leads to the formation of quasi-static coherent magnetic fields with different topologies in photoionized plasmas. Finally, the paper summarizes the important results and discusses future directions on this topic.

Zhang, Chaojie (ORCID:0000000180353014)↗

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

Hydrogenation and C-S bond activation pathways in thiophene and tetrahydrothiophene reactions on sulfur-passivated surfaces of Ru, Pt, and Re nanoparticles

Thiophene-H2 reactions proceed via sulfur removal and hydrogenation routes on dispersed metal nanoparticles that become decorated by refractory S-adlayers during catalysis. The identity and kinetic relevance of the required elementary steps are described here based on rates measured at S-chemical potentials set by H2S/H2 ratios similar to those prevalent during practical catalysis on Re, ReSx, Ru, and Pt catalysts. Free energies for S adatom formation (from H2S decomposition and H2 evolution) are strongly exothermic (< -50 kJ mol-1 on Pt(111) and < -150 kJ mol-1 on Re and Ru(0001)), but strong repulsions between S adatoms cause adsorption free energies to increase significantly with coverage on all three surfaces, preventing complete monolayer formation. These adlayers, composed of unreactive S-atoms (S') that cover 1/3–2/3 ML leave residual interstitial spaces (*) that bind S-atoms (S*), intermediates, and transition states reversibly, as required for catalytic turnovers. The number and binding properties of these interstices depend on the identity and chemical state of the nanoparticle bulk phase, which influences S'-binding and coverages and cause large differences in direct desulfurization and hydrogenation turnover rates (per exposed metal atom) on dispersed Re, ReSx, Ru, and Pt. The identity and kinetic relevance of elementary steps for desulfurization (to C4¬ hydrocarbons) and hydrogenation (to tetrahydrothiophene; THT) are similar among these catalysts; they involve the kinetically-relevant formation of a thiophene-derived intermediate (monohydrothiophene on Re and ReSx; dihydrothiophene on Ru and Pt) that either cleaves its C-S bond or “over-hydrogenates” to THT in one surface sojourn. THT then undergoes C-S bond cleavage in secondary reactions that correct such over-hydrogenation to form the more unsaturated species that cleave C-S bonds. THT/C4 product ratios are insensitive to H2S/H2 ratios and thiophene pressure, even though active interstitial spaces are covered by kinetically-detectable coverages of S* and thiophene; therefore, primary and secondary reactions must involve the same active surfaces. The observed increase in THT/C4 ratios with H2 pressure shows that THT formation transition states involve a larger number of H-atoms than for C-S cleavage. The requirement for bound species with intermediate unsaturation (between THT and thiophene) for C-S bond cleavage is reminiscent of the H-shuttling required in C-C and C-O hydrogenolysis, reactions that involve the partial dehydrogenation of alkanes and alkanols, respectively, to weaken such bonds and to increase the formation entropy of the relevant transition states via the evolution of H2(g). These mechanistic details challenge prevalent paradigms about different site requirements for hydrogenation and desulfurization pathways and about how metal-sulfur bond energies act as descriptors of reactivity; in fact, such binding energies merely act to define the refractory S-adlayers that enable the formation of weakly-binding interstices that reversibly bind intermediates and transition states, thus allowing catalytic turnovers.

Yik, Edwin↗

Nickel/gallium modified HZSM-5 for ethane aromatization: Influence of metal function on reactivity and stability

Bimetallic catalysts often outperform their monometallic counterparts due to the synergistic effect. Here we reported that the Ni/Ga co-functionalized HZSM-5 catalyst, specifically with Ni 1 /Ga 1 stoichiometric ratio, shows significantly improved selectivity (above 80 %) and stability than the Ni/HZSM-5 catalyst and higher activity (with identical selectivity) than the Ga/HZSM-5 catalyst in ethane aromatization. According to the transient experiment, the Ga and Ni 1 Ga 1 /HZSM-5 catalysts show an induction period of a few minutes. However, for Ni/HZSM-5 and Ni 1 Ga 1/3 /HZSM-5, the induction period increased to >10 min and >30 min, respectively, where up to 100 % of ethane was converted into methane and coke through cracking/hydrogenolysis. Furthermore, based on the generalized power-law equation model, a modified model for deactivation kinetic was proposed. Consequently, the influence of the metal function on the kinetic parameters was discussed. We anticipate the formation of Ni 3 Ga intermetallic alloy and its synergy with the exchanged Ga δ+ to be responsible for optimal BTX formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Clumped 13 CH 2 D and 12 CHD 2 compositions of methyl groups from wood and synthetic monomers: Methods, experimental and theoretical calibrations, and initial results

Methyl groups are found in numerous biogenic and synthetic materials including geologically preserved materials such as wood. The carbon and hydrogen isotope compositions of methyl groups are used as tracers in biogeochemical cycles, as paleothermometers, and to determine the hydrogen isotopic composition of ancient rain. Here we present analyses of resolved 13 C–D ( 13 CH 2 D) and D–D ( 12 CHD 2 ) clumped isotope compositions of methyl groups as new variables for the study of methyl groups in the present and past. We first present chemical methods to extract, purify, and derivatize methyl groups from methoxyl (R–O–CH 3 ) groups as CH 3 F and CH 3 Cl, and high-resolution mass spectrometric techniques to determine the clumped isotope compositions of these species. We achieve precisions for 13 C–D clumping of ±0.25‰ and D–D clumping of ±2.5‰. We anchor our clumped isotopic measurements to a thermodynamic reference frame by first calculating the theoretical temperature dependences of 13 C–D and D–D clumping in CH 3 Cl, then placing our measurements onto this reference frame through experimental internal isotopic equilibration of CH 3 Cl at 200 °C. Finally, we provide and analyze an initial dataset of clumped 13 C–D and D–D compositions of methyl groups from various commercial/synthetic monomers and environmental woods. We observe ranges in clumped isotope compositions of ~11‰ in 13 C–D and ~48‰ in D–D, and systematic differences within these ranges between methyl groups from commercial monomers and wood. Specifically, commercial clumped 13 C–D compositions are between 0 and 3‰, which correspond to apparent equilibrium temperatures between 170 °C and the infinite temperature limit. In contrast, the clumped 13 C–D compositions of wood methoxyl groups are distinctively high (9.50–11.25‰) and 3–6‰ higher than would be expected if formed in internal isotopic equilibrium at Earth-surface temperatures. Commercial/synthetic methyl and wood methoxyl clumped D–D compositions are also distinct: -5 to +13‰ in commercial monomers vs. -35 to -8‰ in wood—such negative values cannot result from formation in isotopic equilibrium and require kinetic processes to have occurred. Overall, these results indicate that wood methoxyl groups are formed out of isotopic equilibrium and that clumped isotope compositions of methyl groups may be useful tracers of methyl group sources and sinks in the environment. For instance, isotopic clumping in methyl groups may be useful for understanding controls on isotopic clumping in methane produced by methylotrophic methanogens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗