Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗
DOE OSTI2020

BASS: An R Package for Fitting and Performing Sensitivity Analysis of Bayesian Adaptive Spline Surfaces

In this work, we present the R package BASS as a tool for nonparametric regression. The primary focus of the package is fitting fully Bayesian adaptive spline surface (BASS) models and performing global sensitivity analyses of these models. The BASS framework is similar to that of Bayesian multivariate adaptive regression splines (BMARS) from Denison, Mallick, and Smith (1998), but with many added features. The software is built to efficiently handle significant amounts of data with many continuous or categorical predictors and with functional response. Under our Bayesian framework, most priors are automatic but these can be modified by the user to focus on parsimony and the avoidance of overfitting. If directed to do so, the software uses parallel tempering to improve the reversible jump Markov chain Monte Carlo (RJMCMC) methods used to perform inference. We discuss the implementation of these features and present the performance of BASS in a number of analyses of simulated and real data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tensor Decompositions for Count Data that Leverage Stochastic and Deterministic Optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the global maximum likelihood estimator from local minima. Simultaneously, a recent trend in theoretical computer science and numerical linear algebra leverages randomization to solve very large, hard problems. The typical approach is to use randomization for a fast approximation and determinism for refinement to yield effective algorithms with theoretical guarantees. Two popular algorithms for Poisson CPD reflect that emergent dichotomy: CP Alternating Poisson Regression is a deterministic algorithm and Generalized Canonical Polyadic decomposition makes use of stochastic algorithms in several variants. This work extends recent work to develop two new methods that leverage randomized and deterministic algorithms for improved accuracy and performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Aligned carbon nanotubes for 1,000x improvement in computing (CRADA Final Report)

As part of the Cyclotron Road program, the Aligned Carbon team worked to improve and bring to market their technology to manufacture carbon nanotubes (CNTs) for integrated circuits. There are multiple applications which can benefit from the use of CNT based transistors in integrated circuits, including high performance computing, RF data transmission, and biosensing. Each of these applications faces a problem of poor starting CNT material to come to market. Aligned Carbon offers material with the highest ceiling of performance based around a fully aligned, uniform, on-wafer growth of CNTs and our proprietary purification process to remove metallic CNTs leaving only the semi-conducting ones behind. Further, Aligned Carbon has been guided by scalable manufacturing processes that will fit seamlessly into the integrated circuit foundry model of chip fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2025

KIPM Detector Characterization in QUIET

Kinetic Inductance Phonon‑Mediated (KIPM) detectors are superconducting microwave resonators on silicon that sense bursts of energy through tiny shifts in their resonance. In this work, we installed new devices in the QUIET facility's cryogenic test setup and used a network analyzer to drive them across a range of power levels. At each setting, we measured the transmitted signal and fit it to a simple notch‑filter model, allowing us to extract key performance metrics. And despite the increased noise at low power, our model consistently captured the resonance behavior. We observed that the resonant frequency stayed stable across most powers, the device s internal quality factor improved with stronger signals, and its coupling strength peaked mid‑range before leveling off. Overall, these findings confirm the detector s performance, laying groundwork for future work on energy calibration with pulsed optical sources.

Savitala, Akash [Washington U., Seattle; Fermilab]↗
DOE OSTI2026

Dataset_for_Molecular_Motion_Below_the_Glass_Transition_A_Solid-State_NMR_Study_of_Siloxane_Polymer_Dynamics Study

This dataset contains solid-state 1H and 13C NMR relaxometry data, differential scanning calorimetry (DSC) data, and size exclusion chromatography (SEC/GPC) data supporting the study of sub-glass-transition (sub-Tg) molecular dynamics in a composition- and sequence-controlled series of diphenyl-substituted polysiloxanes (PDMS, 14Ph, 33Ph, 50Ph, 67Ph, and 100Ph; 0–100% diphenylsiloxane content by mole).All solid-state NMR data were acquired on a 200 MHz Bruker Avance III HD spectrometer using a static 7 mm HX probe or a 4 mm HX probe under 4 kHz magic-angle spinning. Raw Bruker TopSpin experiment folders are included for: (1) variable-temperature 1H lineshape measurements used to determine linewidth (FWHM) as a function of temperature across the glass transition; (2) 1H T1 (saturation recovery with solid-echo detection), probing nanosecond-scale dynamics near the 1H Larmor frequency; (3) 1H T1rho (direct spin-lock, 62.5 kHz), probing microsecond-scale segmental dynamics; (4) 13C-detected Lee–Goldburg cross-polarization 1H T1rho (LGCPH T1rho) for 33Ph and 50Ph, resolving aromatic and aliphatic proton environments; and (5) 13C T1 relaxation for 33Ph and 50Ph. Differential scanning calorimetry data (TA Instruments DSC 25, −150 to +120 °C, up to +300 °C for 100Ph, 10 °C/min) are included for all six compositions and support the glass-transition temperatures in Table 1 and Figure 1. Size exclusion chromatography data (Agilent 1200 Series, PL-Gel 300 mixed-C column, THF mobile phase, polystyrene calibration standards) are included for the three synthesized copolymers (33Ph, 50Ph, 67Ph) and support the number-average molecular weights in Table 1. Processed data include per-composition relaxation-time summaries (Excel), curve-fitting and Bloembergen-Purcell-Pound (BPP) model analysis notebooks (Jupyter/Python), and Igor Pro (.pxp) master files used to generate the manuscript's figures.

Bloembergen-Purcell-Pound theory↗
DOE OSTI2025

Predictions for the Detectability of Milky Way Satellite Galaxies and Outer-Halo Star Clusters with the Vera C. Rubin Observatory

We predict the sensitivity of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) to faint, resolved Milky Way satellite galaxies and outer-halo star clusters. We characterize the expected sensitivity using simulated LSST data from the LSST Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) accessed and analyzed with the Rubin Science Platform as part of the Rubin Early Science Program. We simulate resolved stellar populations of Milky Way satellite galaxies and outer-halo star clusters over a wide range of sizes, luminosities, and heliocentric distances, which are broadly consistent with expectations for the Milky Way satellite system. We inject simulated stars into the DC2 catalog with realistic photometric uncertainties and star/galaxy separation derived from the DC2 data itself. We assess the probability that each simulated system would be detected by LSST using a conventional isochrone matched-filter technique. We find that assuming perfect star/galaxy separation enables the detection of resolved stellar systems with $M_V$ = 0 mag and $r_{1/2}$ = 10 pc with >50% efficiency out to a heliocentric distance of ~250 kpc. Similar detection efficiency is possible with a simple star/galaxy separation criterion based on measured quantities, although the false positive rate is higher due to leakage of background galaxies into the stellar sample. When assuming perfect star/galaxy classification and a model for the galaxy-halo connection fit to current data, we predict that 89 +/- 20 Milky Way satellite galaxies will be detectable with a simple matched-filter algorithm applied to the LSST wide-fast-deep data set. Different assumptions about the performance of star/galaxy classification efficiency can decrease this estimate by ~7%-25%, which emphasizes the importance of high-quality star/galaxy separation for studies of the Milky Way satellite population with LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Detection and Bulk Properties of the HR 8799 Planets with High-resolution Spectroscopy

Using the Keck Planet Imager and Characterizer, we obtained high-resolution (R ~ 35,000) K-band spectra of the four planets orbiting HR 8799. We clearly detected H{sub 2}O and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and H{sub 2}O in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward-modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured vsin (i) values of 10.1{sub -2.7}{sup +2.8} km s{sup -1} for HR 8799 d and 15.0{sub -2.6}{sup +2.3} km s{sup -1} for HR 8799 e, and placed an upper limit of <14 km s{sup -1} of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anticorrelated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower-mass planets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

NuSTAR Observations of Abell 2163: Constraints on Non-thermal Emission

Since the first non-thermal reports of inverse Compton (IC) emission from the intracluster medium (ICM) of galaxy clusters at hard X-ray energies, we have yet to unambiguously confirm IC emission in observations with newer facilities. RXTE detected IC emission in one of the hottest known clusters, Abell 2163 (A2163), a massive merging cluster with a giant radio halo—the presumed source of relativistic electrons IC scattering CMB photons to X-ray energies. The cluster's redshift (z ~ 0.2) allows its thermal and non-thermal radio emission to fit the field of view of the Nuclear Spectroscopic Telescope Array (NuSTAR), permitting a deep observation capable of confirming or ruling out the RXTE report. The IC flux provides constraints on the average magnetic field strength in a cluster. To determine the global diffuse IC emission in A2163, we fit its global NuSTAR spectrum with four models: single- (1T) and two-temperature (2T), 1T+power-law component (T+IC), and multi-temperature+power law (9T+IC). Each represents different characterizations of the thermal ICM emission, with power-law components added to represent IC emission. We find that the 3–30 keV spectrum can be described by purely thermal emission, with a global average temperature of kT = (11.8 ± 0.2) keV. The IC flux is constrained to <4.0 × 10 –12 erg s –1 cm –2 using the 1T+IC model and <1.6 × 10 –12 erg s –1 cm –2 with the more physical 9T+IC model, both to 90% confidence levels. Here, combining these limits with 1.4 GHz diffuse radio data from the VLA, we find the average magnetic field strength to be >0.22 μG and >0.35 μG, respectively, providing the strongest constraints on these values in A2163 to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Vertical Structure of Local Disk Galaxies Revealed by DESI Imaging Data

The vertical structure of galactic disks is an important probe of the disk assembly history. Here we investigate the vertical structure of a sample of 79 local disk galaxies within 50 Mpc using data from the Dark Energy Spectroscopic Instrument Legacy Imaging Surveys. Vertical luminosity profiles as a function of radius in the g, r, and z bands are extracted and fitted using a single-component sech2 model to determine scale height and its radial variation (i.e., flaring). Our measurements indicate local galactic disks are overall thin with negligible flaring. The median scale heights at 1 R e are 0.21, 0.22, and 0.22 kpc in the g, r, and z bands, respectively, while the median radial gradients of scale height are −0.006, 0.003, and 0.001 in these bands. These values are consistent with that of the geometric thin disk of the Milky Way represented by metal-rich, low-[α/Fe] populations, confirming the weak flaring of the geometric thin disk. A clear positive correlation of scale height with stellar mass is observed down to the low-mass end of 10 7 M ⊙ . These results provide a homogeneous benchmark for the vertical structure of nearby disk galaxies and establish the cosmological representativeness of the Milky Way thin disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Matrix Multiply Performance of GPUs on Exascale-class HPE/Cray Systems

The computation of dense matrix-matrix products (GEMMs) is central to many modeling and simulation workloads as well as AI/ML deep learning campaigns. In fact, millions of dollars are spent annually on computing GEMMs, and large model training demands are increasing exponentially. Specialized processors such as GPUs are designed to perform well for these operations. However, the performance of GEMMs on GPUs can exhibit complex behaviors depending on many factors, making it challenging to optimize the performance of GEMMs on these processors. In this study we undertake an examination of GEMM performance on several leading GPU models taken from product lines of GPUs to be deployed in forthcoming exascale computing systems. We show results to illustrate the many factors that can affect performance of GEMMs on GPUs. We then present data collected from a large number of test runs for an example GEMM operation to show the dependence behaviors of GEMM rate on matrix dimensions. Finally, we show results from machine learning-based performance models using novel feature engineering methods to fit the measured performance, providing a potential basis for GEMM performance tuning and autotuning methods for GPUs. Recommendations are also given for how to achieve high GEMM performance on modern GPUs.

Melesse Vergara, Veronica↗
DOE OSTI2026

In-Band Scattering and Absorption of Infrared Blocking Foam Filters for Millimeter-wave Cameras

Expanded closed-cell polymer foams are widely used as thermal infrared (IR) blocking filters in millimeter-wave cameras, particularly for Cosmic Microwave Background observations. Precise knowledge of their millimeter-wave properties is essential for optimizing sensitivity. We present broadband (150 GHz - 2 THz) transmittance spectroscopy of Styroace-II and several Zotefoam filters, fitting their spectra with a radiative transfer model incorporating dielectric absorption and Rayleigh, Mie, and higher-order scattering. For a typical 5~cm thick filter stack at 280~GHz, Styroace-II exhibits ${\sim}10\%$ scattering with absorption estimated as ${\lesssim}5\%$ by effective-medium theory, while Zotefoam HD30 offers superior performance at ${\sim}3\%$ scattering and absorption likewise bounded to ${{\lesssim}0.3\%}$. Each model component is constrained at the ${\sim}0.1\%$ transmittance level for millimeter wavelengths. We observe batch-to-batch scattering variability of up to 2 percentage points in foams with multiple tested batches. Less commonly used Zotefoam formulations (LD15 and LD24) can further reduce in-band scattering to ${<}1\%$ while maintaining negligible in-band absorption and likely comparable IR blocking due to shared polyethylene absorption features and similar cell sizes. Based on this work, a filter constructed from the best measured LD24 batch has replaced the Styroace-II filter in a Simons Observatory 220/280 GHz Small Aperture Telescope.

Thomas, Alex [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2024

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗
DOE OSTI2023

The SOUX AGN sample: optical/UV/X-ray SEDs and the nature of the disc

We use the SOUX sample of ~700 active galactic nucleus (AGN) to form average optical-ultraviolet (UV)-X-rays spectral energy distributions (SEDs) on a two-dimensional (2D) grid of M BH and L 2500 . We compare these with the predictions of a new AGN SED model, qsosed, which includes prescriptions for both hot and warm Comptonization regions as well as an outer standard disc. This predicts the overall SED fairly well for 7.5 < log(M BH /M ⊙ ) < 9.0 over a wide range in L/L Edd , but at higher masses the outer disc spectra in the model are far too cool to match the data. We create optical-UV composites from the entire Sloan Digital Sky Survey sample and use these to show that the mismatch is due to there being no significant change in spectral shape of the optical-UV continuum across several decades of M BH at constant luminosity. We show for the first time that this cannot be matched by standard disc models with high black hole spin. These apparently fit, but are not self-consistent as they do not include the General Relativistic effects for the emission to reach the observer. At high spin, increased gravitational redshift compensates for almost all of the higher temperature emission from the smaller inner disc radii. The data do not match the predictions made by any current accretion flow model. Either the disc is completely covered by a warm Comptonization layer whose properties change systematically with L/L Edd , or the accretion flow structure is fundamentally different to that of the standard disc models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Full-waveform tomography reveals iron spin crossover in Earth’s lower mantle

Abstract Three-dimensional models of Earth’s seismic structure can be used to identify temperature-dependent phenomena, including mineralogical phase and spin transformations, that are obscured in 1-D spherical averages. Full-waveform tomography maps seismic wave-speeds inside the Earth in three dimensions, at a higher resolution than classical methods. By providing absolute wave speeds (rather than perturbations) and simultaneously constraining bulk and shear wave speeds over the same frequency range, it becomes feasible to distinguish variations in temperature from changes in composition or spin state. We present a quantitative joint interpretation of bulk and shear wave speeds in the lower mantle, using a recently published full-waveform tomography model. At all depths the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in ferropericlase provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by an enrichment in silica towards the core-mantle boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean.

58 GEOSCIENCES↗
DOE OSTI2023

The importance of flavor in SMEFT Electroweak Precision Fits

Effective field theory tools are essential for exploring non-Standard Model physics at the LHC in the absence of the discovery of new light particles. Predictions for observables are typically made at the lowest order in the QCD and electroweak expansions in the Standard Model effective field theory (SMEFT) and often ignore the effects of flavor. Here, we present results for electroweak precision observables (EWPOs) at the next-to-leading order QCD and electroweak expansions (NLO) of the SMEFT with an arbitrary flavor structure for the fermion operators. Numerical NLO SMEFT fits to EWPOs have a strong dependence on the assumed flavor structures and we demonstrate this using various popular assumptions for flavor symmetries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Effect of Sample Mass, Confinement, and Preheating Time on the Thermal Response of LLM‐105: Experiments and Kinetic Analysis

Various small-scale experiments were performed to provide data for developing a model to predict the thermal response of LLM-105 over a wide range of conditions. The thermal decomposition of LLM-105 was studied as a function of sample mass, confinement of volatile products, and preheating time in both isothermal and ramped heating experiments. The thermal decomposition of LLM-105 is a two-step process, as shown by the two exothermic peaks in the heat flow profiles, which were fitted to two nth-order autocatalytic reaction models with a similar activation energy of ∼289 kJ/mol. The magnitude and shape of these peaks varied with sample mass and confinement. Increasing sample mass enhanced the second exotherm with respect to the first one, while increasing the level of confinement promoted a transition from a sublimation-dominated regime towards thermal decomposition. The effect of LLM-105 particle size on the rate of weight loss was evident for open-pan experiments, where bigger particles sublimed at lower temperatures than smaller particles. Thermal response and solid residue composition of LLM-105 samples were analyzed following preheating for different durations. Longer preheating times caused a shift of the second exotherm to lower temperatures and a decrease in the reaction enthalpy, confirming that LLM-105 decay is a consecutive reaction mechanism, probably autocatalytic. In conclusion, the kinetic model derived from ramped experiments was validated against the measured LLM-105 fraction remaining and the enthalpy remaining of the solid residue as a function of preheating times and showed good agreement.

Chemistry - Chemical explosives↗
DOE OSTI2020

Oxide-catalyzed self- and cross-condensation of cycloketones. Kinetically relevant steps that determine product distribution

Metal oxides with acid-base properties can be used to catalyze aldol condensation and potentially obtain valuable chemicals from biomass-derived oxygenates. Understanding the reaction mechanisms could enable the development of effective catalysts with controllable selectivity. In this work, we investigate the reaction parameters that affect the kinetics and product distribution of self- and cross-aldol condensation of cyclopentanone (CPO) and cyclohexanone (CHO) over CeO 2 and ZrO 2 catalysts by combining detailed analysis of experiments and DFT calculations. It is well-known that either the deprotonation of an α-C–H that generates an enolate or the C–C coupling that produces the dimers can be rate-limiting step. Indeed, the detailed kinetics analysis shows that the observed data can be fitted by unimolecular or bimolecular Langmuir-Hinshelwood models. Specifically, on the ZrO 2 catalyst, the unimolecular α-C–H activation is the rate limiting step for the CPO self-condensation, while the bimolecular C–C coupling limits the CHO self-condensation. By contrast, on CeO 2 , the bimolecular C–C coupling step seems to be rate-limiting for both CPO and CHO self-condensation reactions. The adsorption parameters derived from kinetics fitting indicate that when the two reactants are co-fed, the surface is preferentially covered by CHO. So, the condensation rate of [CPO]-activated products ([CPO]CPO + [CPO]CHO) decreases in the presence of CHO, due to the competitive adsorption. However, despite a higher surface coverage, CHO is not an effective electrophile because it imposes stronger steric constraints to the formation of the C–C bond than CPO. Therefore, for a given enolate, the product distribution shows that [CPO]CPO > [CPO]CHO and, similarly [CHO]CPO > [CHO]CHO. Likewise, the yields of [CPO]- and [CHO]-activated products in mixed CPO/CHO feeds are similar because the higher surface coverage and stability of the [CHO] enolate is counterbalanced by the inhibition of C–C coupling by steric constraints when CHO is involved. In summary, the overall rank of products observed experimentally can be explained in terms of a combination of the following effects: relative surface coverage of the reactants, steric constraints for the C–C coupling, stability of the enolate, and effectiveness of the electrophile to accept an electron.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗