Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data 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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A machine learning pipeline for membrane segmentation of cryo-electron tomograms

We describe how to use several machine learning techniques organized in a learning pipeline to segment and identify cell membrane structures from cryo electron tomograms. These tomograms are difficult to analyze with traditional segmentation tools. The learning pipeline in our approach starts from supervised learning via a special convolutional neural network trained with simulated data. It continues with semi-supervised reinforcement learning and/or a region merging technique that tries to piece together disconnected components belonging to the same membrane structure. A parametric or non-parametric fitting procedure is then used to enhance the segmentation results and quantify uncertainties in the fitting. Domain knowledge is used in generating the training data for the neural network and in guiding the fitting procedure through the use of appropriately chosen priors and constraints. We demonstrate that the approach proposed here works well for extracting membrane surfaces in two real tomogram datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Visible-light Phase Curves from the Second Year of the TESS Primary Mission

We carried out a systematic study of full-orbit phase curves for known transiting systems in the northern ecliptic sky that were observed during Year 2 of the TESS primary mission. We applied the same methodology for target selection, data processing, and light-curve fitting as we did in our Year 1 study. Out of the 15 transiting systems selected for analysis, seven—HAT-P-7, KELT-1, KELT-9, KELT-16, KELT-20, Kepler-13A, and WASP-12—show statistically significant secondary eclipses and day–night atmospheric brightness modulations. Small eastward dayside hot-spot offsets were measured for KELT-9b and WASP-12b. KELT-1, Kepler-13A, and WASP-12 show additional phase-curve variability attributed to the tidal distortion of the host star; the amplitudes of these signals are consistent with theoretical predictions. We combined occultation measurements from TESS and Spitzer to compute dayside brightness temperatures, TESS-band geometric albedos, Bond albedos, and phase integrals for several systems. The new albedo values solidify the previously reported trend between dayside temperature and geometric albedo for planets with 1500 K < T {sub day} < 3000 K. For Kepler-13Ab, we carried out an atmospheric retrieval of the full secondary eclipse spectrum, which revealed a noninverted temperature–pressure profile, significant H{sub 2}O and K absorption in the near-infrared, evidence for strong optical atmospheric opacity due to sodium, and a confirmation of the high geometric albedo inferred from our simpler analysis. We explore the implications of the phase integrals (ratios of Bond to geometric albedos) for understanding exoplanet clouds. We also report updated transit ephemerides for all of the systems studied in this work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

We present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. We inverted waveforms from 192 Global Centroid Moment Tensor earthquakes (MW 5.5-7.0) recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region. Spatial coverage of the available data is highly uneven due to earthquakes clustered along plate boundaries and sparse coverage of open seismic networks in the region. We considered three possible starting models: the SPiRaL global model (Simmons et al., 2021); MEC-1 (Kaviani et al., 2020); and CSEM2.0 (Noe et al., 2023). Because the SPiRaL model provides good fits to the observed waveforms measured by the time-bandwidth product of selected windows in several period bands, provides all the necessary parameters and covers the entire domain we used it for the starting model with the period band 50-100 seconds. Inversion iterations proceeded using time-frequency phase misfits in six stages and 54 total iterations reducing the minimum period to 30 seconds. Our final model, MESWA, provides improved waveform fits compared to the starting model for both the data used in the inversion and an independent validation data set of 66 events. Two metrics of waveform fit (the time-frequency phase misfit used in the optimization and normalized L2 misfit) were both reduced by nearly 60% for both data sets and MESWA provides significantly larger misfit reductions relative to the SPiRaL model than the MEC-1 or CSEM models. We also find that MESWA provides a larger time-bandwidth product of selected windows indicating that more information content of the observed waveforms is explained by MESWA than the other models. Our new model reveals tectonic features imaged by other studies and methods but in a new holistic model of shear and compressional wavespeeds (v S and v P , respectively) with anisotropy covering the crust and uppermost mantle of a larger domain. MESWA has smaller scale-length features and tends to sharpen some features relative to the SPiRaL starting model. Examples include: low crustal v S in the TurkishIranian Plateau, Zagros Mountains, Afghan Central Blocks and Sulaiman Fold Belt; low mantle vSfollowing divergent (Gulf of Aden, Red Sea) and transform (Dead Sea Fault) margins of the Arabian Plate; low and high v S in the mantle beneath the Arabian Shield and Platform, respectively. Low vS is imaged below Cenozoic volcanic centers of the Arabian Peninsula, the so-called Mecca-Madina-Nafud (MMN) Line. Positive anisotropy (v SH > v SV ) is inferred for asthenospheric depths across the region except where up/downwelling may influence fabric alignment (e.g. Afar, Red Sea, Arabian Shield). Elevated vS tracks Makran subduction under southeast Iran. MESWA resembles the SPiRaL model in its long-wavelength structure, but enhances shorter wavelengths features on the order of 200 km and smaller. The resulting model could be used for as a starting model for further improvements, say using waveforms from in-country seismic networks that are not openly available or smaller-scale studies targeting shorter period waveforms. The model also could be used for source characterization and moment tensor inversion to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2024

OpenPATH - Leveraging Technology to Measure Travel Behavior

Shifting transportation to more sustainable modes is a key piece of the decarbonization puzzle. However, mobility behavior and travel patterns are difficult to influence because they are difficult to measure. OpenPATH provides a tool to capture longitudinal behaviors through a smartphone application. Agencies interested in gathering data about a population's travel behavior can set up a deployment of the app customized to the needs of their community. Partners can choose between simple mode and purpose labels or surveys for each trip to balance the level of user engagement with the associated burden. The labels, trip surveys, and an initial demographic survey can all be tailored to the specific context of the deployment. The OpenPATH tool is unique in its open-source nature, ability to gather detailed longitudinal travel data, and design allowing direct engagement with travelers. A valuable technological advancement, this tool enables partners to measure the way changes in the transportation landscape impact their community. The suite of tools includes both public and administrator dashboards. The public dashboard supports continuous data analysis through charts presenting trip information updated daily. The administrator dashboard displays geospatial data and supports data export. Example applications have included e-bike programs; gathering valuable metrics on increased access to opportunities and reduction in VMT, and studies aimed at understanding existing mobility behavior to see where advancements such as electric vehicles could fit into these habits. OpenPATH collects travel data in association with an initial demographic survey, enabling detailed insight into the behavior patterns or impact of a certain program on different populations.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Quantitative XPS of plutonium: Evaluation of the Pu4f peak shape, relative sensitivity factors and estimated detection limits

We report that high-resolution X-ray photoelectron spectra have been acquired from sputter cleaned and in situ oxidized samples of α-plutonium and a face-centred cubic δ-plutonium–gallium alloy. The differences in the Pu4f peak shape in these two materials have been investigated, and the poorly screened satellite peaks have been quantified. Curve fitting parameters for the Pu4p 3/2 , Pu4d 5/2 , Pu4f and Pu6p 3/2 photoelectron peaks are reported, and relative sensitivity factors have been determined. The Pu4f curve fit model has been applied to data acquired using different spectrometers and alloys. Examples of quantification of the plutonium spectra are provided and minimum detection limits are calculated for common impurities in plutonium metal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Search for vector-like leptons coupling to first- and second-generation Standard Model leptons in pp collisions at s = 13 TeV with the ATLAS detector

A search for pair production of vector-like leptons coupling to first- and second-generation Standard Model leptons is presented. The search is based on a dataset of proton-proton collisions at s$$ \sqrt{s} $$ = 13 TeV recorded with the ATLAS detector during Run 2 of the Large Hadron Collider, corresponding to an integrated luminosity of 140 fb−1. Events are categorised depending on the flavour and multiplicity of leptons (electrons or muons), as well as on the scores of a deep neural network targeting particular signal topologies according to the decay modes of the vector-like leptons. In each of the signal regions, the scalar sum of the transverse momentum of the leptons and the missing transverse momentum is analysed. The main background processes are estimated using dedicated control regions in a simultaneous fit with the signal regions to data. No significant excess above the Standard Model background expectation is observed and limits are set at 95% confidence level on the production cross-sections of vector-like electrons and muons as a function of the vector-like lepton mass, separately for SU(2) doublet and singlet scenarios. The resulting mass lower limits are 1220 GeV (1270 GeV) and 320 GeV (400 GeV) for vector-like electrons (muons) in the doublet and singlet scenarios, respectively.

Aad, G↗
DOE OSTI2023

Isospin-breaking effects in the three-pion contribution to hadronic vacuum polarization

Isospin-breaking (IB) effects are required for an evaluation of hadronic vacuum polarization at subpercent precision. While the dominant contributions arise from the e + e - → π + π - channel, also IB in the subleading channels can become relevant for a detailed understanding, e.g., of the comparison to lattice QCD. Here, we provide such an analysis for e + e - → 3π by extending our dispersive description of the process, including estimates of final-state radiation (FSR) and ρ–ω mixing. In particular, we develop a formalism to capture the leading infrared-enhanced effects in terms of a correction factor η3π that generalizes the analog treatment of virtual and final-state photons in the 2π case. The global fit to the e+e- → 3π data base, subject to constraints from analyticity, unitarity, and the chiral anomaly, gives for the total 3π contribution to the anomalous magnetic moment of the muon, of which ${a}^{FSR}_\mu [3\pi]$ = 0.51(1) x 10 -10 and ${a}^{\rho - \omega}_\mu [3\pi]$ = -2.68(70) x 10 -10 can be ascribed to IB. We argue that the resulting cancellation with ρ–ω mixing in e + e - → 2π can be understood from a narrow-resonance picture, and provide updated values for the vacuum-polarization-subtracted vector-meson parameters M ω = 782.70(3) MeV, M Φ = 1019.21(2) MeV, Γ ω = 8.71(3) MeV, and Γ Φ = 4.27(1) MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

NO x adsorption with CO and C 2 H 4 on Pd/SSZ-13: Experiments and modeling

A transient monolith model containing microkinetic schemes for NO uptake and release over Pd/SSZ-13 without and with reductants CO or C2H4 is presented. The scheme involves three cationic Pd species (Z – [PdOH] + , Z – Pd 2+ Z – , Z – Pd + ) as the active sites during uptake of NO, CO and C 2 H 4 , and their desorption and conversion at higher temperature. Kinetic parameters are estimated through a combination of density functional theory (DFT) estimates and a fit of uptake, desorption and conversion data. The tuned model is validated at different uptake temperatures, ramp rates, and flowrates. A “degree of uptake control” parameter is defined that helps to identify the step(s) that are uptake controlling. The model helps to interpret the data features and is used to identify operating conditions to meet application-relevant performance metrics, including NO trapping efficiency and NO release temperature. Furthermore, the model demonstrates higher NO uptake on Pd/SSZ-13 compared to Pd/ZSM-5.

42 ENGINEERING↗
DOE OSTI2022

Recurrent neural networks for short-term and long-term prediction of geothermal reservoirs

Accurate prediction of geothermal reservoir responses to alternative energy production scenarios is critical for optimizing the development of the underlying resources. While the conventional physics-based models offer a comprehensive prediction tool, data-driven models provide an efficient alternative to build fit-for-purpose predictive models by extracting and using the statistical patterns in the collected data to make predictions. The recurrent neural network (RNN) is a data-driven model that is commonly applied to predict time series sequences. This paper presents a variant of RNN that also utilizes the efficiency of convolutional neural networks (CNN) for the prediction of energy production from geothermal reservoirs. Specifically, a CNN–RNN architecture is developed that takes historical well controls as input (features) and their corresponding production response data as output (labels) to learn an input-output mapping that can predict the future well production responses/performance for any given future well control inputs. The model is paired with a labeling scheme to handle real field disturbances that create data gaps. In addition to the model structure, we introduce a thorough workflow for applying the model, which includes data pre-processing, feature selection, as well as different training strategies for short-term and long-term prediction. Finally, the performance and accuracy of the model are evaluated by applying it to multiple datasets, including a field reservoir model.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Characterization of the indoor near-field aerosol transmission in a model commercial office building

To evaluate the exposure potential of infectious aerosols containing SARS-CoV-2 in an office building setting, synthetic test aerosols were used to experimental study airborne particle transmission in a multizone small office test building at the Oak Ridge National Laboratory. Nine measurement points in a single zone using active aerosol impactors report that the coefficient of variation of the time-averaged concentration is <10% in two campaigns and < 15% in one campaign, so a nearly well-mixed condition was noted. To understand the effect of HVAC system operation on the dynamic concentration of aerosols in office spaces, an aerosol transport model that includes factors such as outside air (OA) ratio, filtration, return air fraction, transport loss in air ducts, and particle deposition has been developed. The results of model fitting demonstrate strong agreement with experimental data. Overall, our investigation finds the return air fraction effects outweigh other mechanisms for the aerosol recirculation in this study, and the impact of air change rate (ACR) is more important than the small particle deposition for aerosol removal. Because ACR dominates the aerosol transport, the full model can be simplified to just one factor, the ACR, while maintaining an acceptable representation of the experimental data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Fouling modeling and prediction approach for heat exchangers using deep learning

In this article, we develop a generalized and scalable statistical model for accurate prediction of fouling resistance using commonly measured parameters of industrial heat exchangers. This prediction model is based on deep learning where a scalable algorithmic architecture learns non-linear functional relationships between a set of target and predictor variables from large number of training samples. Here, the efficacy of this modeling approach is demonstrated for predicting fouling in an analytically modeled cross-flow heat exchanger, designed for waste heat recovery from flue-gas using room temperature water. The performance results of the trained models demonstrate that the mean absolute prediction errors are under 10 –4 KW –1 for flue-gas side, water side and overall fouling resistances. The coefficients of determination (R 2 ), which characterize the goodness of fit between the predictions and observed data, are over 99%. Even under varying levels of measurement noise in the inputs, we demonstrate that predictions over an ensemble of multiple neural networks achieves better accuracy and robustness to noise. We find that the proposed deep-learning fouling prediction framework learns to follow heat exchanger flow and heat transfer physics, which we confirm using locally interpretable model agnostic explanations around randomly selected operating points. Overall, we provide a robust algorithmic framework for fouling prediction that can be generalized and scaled to various types of industrial heat exchangers.

42 ENGINEERING↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

The maximum carboxylation rate of Rubisco affects CO2 refixation in temperate broadleaved forest trees

Mesophyll resistance to CO 2 diffusion (rm) and the maximum carboxylation rate of Rubisco (Vcmax) affect photosynthetic rates, and can potentially also influence the percentage of respiratory and photorespiratory CO 2 being refixated (Pr) by mesophyll cells. Here we investigated how various leaf anatomical traits (e.g. leaf mass per area [LMA] and leaf dry matter content [LDMC]) influenced rm in leaves of mature forest trees. We further explored how rm and Vcmax in turn affected Pr, and if these traits varied among species and leaves along a light gradient. Photosynthetic CO 2 response of leaves grown in high-, medium-, and low-light environments was measured, from Pinus sylvestris [Scots pine], Picea abies [Norway spruce], Quercus robur [English oak], and Betula pendula [Silver birch] in southern Sweden. A modified version of the Farquhar-von Caemmerer-Berry model was fitted to the leaf gas exchange data to estimate Vcmax, rm and Pr. We found that of all leaf traits measured, only LMA for Q. robur was significantly higher in leaves from high-light environments. When comparing species, both rm and LMA were significantly higher in the conifers, and rm had a negative correlation with Vcmax. We found that Pr was similar between different species and functional groups, with an average of 73.2% (and SD of ±10.4) across all species. There was a strong, positive correlation between Pr and Vcmax in broadleaves, and we hypothesise that this effect might derive from a higher CO 2 drawdown near Rubisco in leaves with high Vcmax.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Behavior of astatine and bismuth in non-conventional solvents: Extraction into imidazolium-based ionic liquid and methyl anthranilate with active pharmaceuticals binary mixtures from nitric acid media

Astatine is one of the least chemically studied elements and its behavior in the presence of non- conventional solvents has not been investigated before. This work considers both hydrophobic ionic liquids and binary mixtures as alternatives to conventional solvents. The study is based on the extraction of astatine and bismuth (target material required to produce astatine) into imidazolium-based ionic liquid, and binary mixtures formed by active pharmaceuticals (ibuprofen, lidocaine) and a food grade ingredient (methyl anthranilate). It is shown that both the ionic liquid and binary mixtures can successfully extract At from nitric acid media, but extraction of Bi into the ionic liquid is very inefficient, resulting in a good separation factor for these two elements in the entire studied acidity range. Extraction of At into binary mixtures is very efficient, having distribution ratio values as high as 1000, while the behavior of Bi under these conditions depends on the composition of the mixture. Furthermore, a mathematical model has been developed to fit both At and Bi experimental data and applied to determine corresponding thermodynamic extraction constants.

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

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

High-Pressure Investigation of 2,4,6-Trinitro-3-bromoanisole (TNBA): Structural Determination and Piezochromism

Understanding phase transitions in energetic materials is crucial for developing predictive models of detonation. 2,4,6-Trinitro-3-bromoanisole (TNBA), an energetic material, was studied in its single-crystal form up to pressures of 45 GPa in a diamond anvil cell. The material was characterized by using X-ray, Raman, and optical transmission measurements. From single-crystal X-ray diffraction, the ambient structure of TNBA was determined which crystallizes in the P2 1 /c space group having four molecular units per unit cell. The X-ray data up to 9.2 GPa were fitted to a third-order Birch–Murnaghan equation of state by using the parameters K o = 13.2(2.4) GPa and K p = 5.1(1.4). Between 6.8 and 7.3 GPa, a phase transition was inferred in TNBA from concurrent fading of X-ray diffraction, disappearance of Raman peaks, increase in sample fluorescence, and discontinuous color change. The new phase was consistent with an amorphous state of at least partially intact molecules judging from the presence of higher-order Raman modes and irreversibility of the Raman spectra upon release. Piezochromism was observed with the translucent yellow TNBA gradually darkening and becoming opaque black at ~25 GPa. This correlated to the absorption edge gradually shifting to the red in the visible spectrum. Signs of two possible additional structural transitions were detected in the 32.4–41.0 GPa range as suggested by a jump in the absorption edge, the irreversible changes in the absorption spectrum upon release to ambient pressure, and by the lack of Raman modes in recovered samples. Here, the crystal and electronic structures of TNBA were also investigated up to 10 GPa by using DFT calculations and crystal structure prediction (CSP) simulations. In agreement with the experimentally observed transition at 7 GPa, the simulations at 10 GPa found a bevy of polymorphs lower in enthalpy and higher in density than P2 1 /c. The lowest calculated enthalpy structure was determined to be P2 1 2 1 2 1 , being in a different space group than the ambient experimental result.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗