Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Accounting”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Controlling Heterogeneous Catalysis with Organic Monolayers on Metal Oxides

A key theme of heterogeneous catalysis research is achieving control of the environment surrounding the active site to precisely steer the reactivity toward desired reaction products. One method toward this goal has been the use of organic ligands or self-assembled monolayers (SAMs) on metal nanoparticles. Metal-bound SAMs are typically employed to improve catalyst selectivity but often decrease the reaction rate as a result of site blocking from the ligands. Recently, the use of metal oxide-bound organic modifiers such as organophosphonic acid (PA) SAMs has shown promise as an additional method for tuning reactions on metal oxide surfaces as well as modifying oxide-supported metal catalysts. In this Account, we summarize recent approaches to enhance catalyst performance with oxide-bound monolayers. These approaches include (1) modification of metal oxide catalysts to tune surface reactions, (2) formation of SAMs on the oxide component of supported metal catalysts to modify sites at the metal–support interface, and (3) enhancement of catalyst performance (e.g., stability) through modification of sites remote from the active sites. Further, both the headgroups and organic tail groups of PA SAMs or other ligands can influence reactions on metal oxide surfaces. Binding of the headgroup can selectively poison certain active sites, altering the selectivity in a manner analogous to metal-bound ligands (at the expense of active site quantity). Moreover, tail groups can be functionalized to interact favorably with reactants and intermediates, for instance through dipole–dipole interactions. On supported metal catalysts like Pt/Al 2 O 3 , PA SAMs can selectively form on the oxide support. This selective deposition allows for modification of the metal–support interface with minimal blockage of metal sites. PA headgroups were shown to provide tunable acid sites at the interface, dramatically improving hydrodeoxygenation rates of various alcohols. Additionally, organic tail functionality was used to activate or stabilize specific reactants at the interface, such as with the use of amine-functionalized PAs to stabilize chemisorption of CO 2 during the reverse water gas shift reaction. PAs have also been found to affect the electronic properties of bulk metal sites through long-range electron withdrawal via the oxide, providing an additional avenue to tune catalytic behavior. Finally, organic modifiers were shown to enhance catalytic performance without directly modifying the active site. For instance, in biphasic liquid environments the modification of catalyst particles with hydrophobic or hydrophilic SAMs shifts the selectivity of multipath reactions on the basis of the hydrophobicities of different intermediates and products. As another “long-range” effect, the deposition of ligands on oxide supports improved catalyst stability through both improved resistance to sintering and suppression of active site poisoning. The recent contributions discussed in this Account demonstrate the versatility and significant potential for the approach of modifying catalysts with oxide-bound organic monolayers.

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

Chemistry, Structure, and Function of Lone Pairs in Extended Solids

The lone pair has been a known feature of the electronic structure of molecules for over 100 years. Beginning with the pioneering work of Lewis and others that was later developed into useful guidelines for predicting molecular structure, lone pairs and their steric consequences are now taught at the very earliest stages of a chemistry education. In the crystalline solid state, lone pairs have perhaps had a less visible yet equally consequential role, with a significant impact on a range of properties and functionalities. Important properties associated with s 2 electron-derived lone pairs include their role in creating conditions favorable for ion transport, in the formation and correlation of local dipoles and the resulting polar behavior leading to ferroics and multiferroics, in increasing the refractive index of glass, in reducing the thermal conductivity of thermoelectric materials, and in breaking local symmetry permitting second-harmonic light generation.. In recent years, the role of the lone pair in developing the electronic structure of some topological quantum materials has also been recognized. While structural distortions due to lone pairs have traditionally been characterized through their crystallography, recent advances in scattering and spectroscopy have revealed the presence of local lone pair-driven distortions that do not correlate over long length scales. The role of these crystallographically “hidden” lone pairs, their detection, and their impact on properties have become a growing body of work in the literature. Hidden lone pairs are an effective argument for considering a role for lone pairs that goes beyond their being objects that occupy space in the coordination polyhedra of cations. This Account introduces the chemistry of lone pairs in extended crystalline solids, including a discussion of when they are stereochemically active, how they manifest in the structure, and how their chemistry can be tuned by the chemical environment around them. Eventually, all of these factors work in unison to help develop and tune properties of interest. Certain specific examples of structure–property relationships in materials that are driven by lone pair behavior are described here, including the potential impact of lone pairs on the optical and electronic properties of hybrid halide perovskite compounds that are relevant to their photovoltaic applications. Here, we highlight the role of lone pairs in the dielectric behavior of geometrically frustrated pyrochlores, the temperature-dependent optoelectronic behavior of halide perovskites, the polar phase transitions in lead-free ferroelectric perovskites, and the compositional insulator-to-metal transition in ruthenium pyrochlores. The theme underpinning this Account is that the lone pair can be considered to be a powerful design element for a broad range of material function.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Emerging Electrochemical Techniques for Probing Site Behavior in Single-Atom Electrocatalysts

Single-atom catalysts (SACs) have aroused tremendous interest over the past decade, particularly in the community of energy and environment-related electrocatalysis. A rapidly growing number of recent publications have recognized it as a promising candidate with maximum atomic utilization, distinct activity, and selectivity in comparison to bulk catalysts and nanocatalysts. However, the complexity of localized coordination environments and the dispersion of isolated sites lead to significant difficulties when it comes to gaining insight into the intrinsic behavior of electrocatalytic reactions. Furthermore, the low metal loadings of most SACs make conventional ensemble measurements less likely to be accurate on the subnanoscale. Thus, it remains challenging to probe the activity and properties of individual atomic sites by available commercial instruments and analytical methods. In spite of this, continuing efforts have lately focused on the development of advanced measurement methodologies, which are very useful to the fundamental understanding of SACs. There have recently been a number of in situ/operando techniques applied to SACs, such as electron microscopy, spectroscopy, and other analysis methods, which support relevant functions to identify the active sites and reaction intermediates and to investigate the dynamic behavior of localized structures of the catalytic sites. This Account aims to present recent electrochemical probing techniques which can be used to identify single-atomic catalytic sites within solid supports. First, we describe the basic principles of molecular probe methods for the study and analysis of electrocatalytic site behavior. In particular, the in situ probing technique enabled by surface interrogation scanning electrochemical microscopy (SISECM) can measure the active site density and kinetic rate with high resolution. An alternative electrochemical probing technique is further demonstrated on the basis of single-entity electrochemistry, which allows the unique electrochemical imaging of the size and catalytic rate of single atoms, molecules, and clusters. Next, the merits and limitations of different electrochemical techniques are then discussed, along with perspectives for future prospects. Apart from this, we further showcase the powerful capability of emerging electrochemical probing techniques for determining significant effects and properties of SACs for various electrocatalytic reactions, including oxygen reduction and evolution, hydrogen evolution, and nitrate reduction. Overall, electrochemical techniques with atomic resolution have greatly increased opportunities for observing, measuring, and understanding the surface and interface chemistry during energy conversion. In the future, it is anticipated that the development of electrochemical probing techniques will be advanced with innovative perspectives on the behavior and features of SACs. We hope that this Account can contribute in several ways to promoting the fundamental knowledge and technical progress of emerging electrochemical measurements for studying SACs.

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

The Influence of Nonlocal Metal–Dielectric Environments on Physical and Chemical Processes

Plasmonic nanolayers and laminar metallic/dielectric multilayers were originally developed for optical cloaking applications and lensing applications that could potentially image objects whose size was below the diffraction limit. These assemblies were initially formed from gold or silver nanorods grown within an alumina mesh. However, more recently, assemblies with similar properties have also been prepared by sequential thin-layer deposition of alternating layers of gold and magnesium fluoride (MgF 2 ). These metal/dielectric composite materials enable control of the dielectric constant in the directions perpendicular to the layers and balance the real and imaginary dielectric constants of the assembly such that the speed and the amplitude of the waves traveling through the assembly are not attenuated. In this Account, we will also focus on a few of the applications ranging from surface wetting to fluorescence quenching to enhancement of photochemical reactions. First, we will share an introduction to processes used to create these materials, which are combinations of low refractive index metals and transparent higher index materials arranged in a scalable repeating fashion. Two fabrication methods were employed: an electrochemical deposition of Ag nanorods into an anodized alumina matrix which produced materials with an anisotropic negative refractive index material within the plane of the film and lamellar metal/dielectric layers in which the negative index perpendicular to the growth direction. These alternating layers of plasmonic metals and dielectric materials were ultimately chosen to prepare films for further testing, because of their relative ease of fabrication. We will continue with a discussion of a few of the applications of both of these nonlocal dielectric composite materials including more specialized plasmonic, composite, and hyperbolic metamaterials including fluorescence quenching, photochemical reactions, and surface wetting. In each of these applications, the unique response caused by the enhancement of the electric field and the interface between hyperbolic materials and plasmonic materials as they interact photophysically with their near neighbors is presented. In each of the applications, the enhanced electric field extends from the composite substrate layer to interact with its near neighbors and beyond. Furthermore, the presence of this extended interaction can be observed in the form of decreased emission lifetime, enhancement of photochemical reaction rates, and changes in the surface energies measured by contact angle goniometry. In this Account, all of these situations will be addressed. Finally, we will conclude with a summary and vision for the future as well as a discussion of the unique challenges and opportunities available as research active faculty at an HBCU.

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

Nontraditional Approaches To Enable High-Energy and Long-Life Lithium–Sulfur Batteries

In this study, lithium-sulfur (Li-S) batteries are promising for automotive applications due to their high theoretical energy density (2600 Wh/kg). In addition, the natural abundance of sulfur could mitigate the global raw material supply chain challenge of commercial lithium-ion batteries that use critical elements, such as nickel and cobalt. However, due to persistent polysulfide shuttling and uncontrolled lithium dendrite growth, Li-S batteries using nonencapsulated sulfur cathodes and conventional ether-based electrolytes suffer from rapid cell degradation upon cycling. Despite significant improvements in recent decades, there is still a big gap between lab research and commercialization of the technology. To date, the reported cell energy densities and cycling life of practical Li-S pouch cells remain largely unsatisfactory. Traditional approaches to improving Li-S performance are primarily focused on confining polysulfides using electronically conductive hosts. However, these micro- and mesoporous hosts suffer from limited pore volume to accommodate high sulfur loading and the associated volume change during cycling. Moreover, they fail to balance adsorption-conversion of polysulfides during charge-discharge, leading to the formation of massive dead sulfur. Such hosts are themselves electrochemically inactive, which decreases the practical energy density. In contrast, a series of nontraditional approaches, paired with advances in multiscale mechanistic understanding, have recently demonstrated exciting performance outcomes not only in conventional coin cells but also in practical pouch cells. In this Account, we first introduce our novel cathode design strategies to overcome polysulfide shuttling and sluggish redox kinetics in thick S cathodes via selenium-sulfur chemistry and cathode host engineering. Next, we gain a mechanistic understanding of Li-S batteries in various types of electrolytes via a series of spectroscopic, nuclear magnetic resonance, and electrochemical methods. Meanwhile, a novel cathode solid electrolyte interphase encapsulation strategy via nonviscous highly fluorinated ether-based electrolyte is introduced. The established selection rule by investigating how solvating power retards the shuttle effect and induces robust cathode/solid-electrolyte interphase formation is also included. We then discuss how the synergistic interactions between rational cathode structures and electrolytes can be exploited to tailor the reaction pathways and kinetics of S cathodes under high mass loading and lean electrolyte conditions. In addition, a novel interlayer design to simultaneously overcome degradation processes (polysulfide shuttling and lithium dendrite formation) and accelerate redox reaction kinetics is presented. Finally, this Account concludes with an overview of the challenges and strategies to develop Li-S pouch cells with high practical energy density, long cycle life, and fast-charging capability.

25 ENERGY STORAGE↗
DOE OSTI2025

Population Balance Models for Catalytic Depolymerization: From Elementary Steps to Multiphase Reactors

Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.

Manis, Lela K. [University of Illinois at Urbana-C↗
DOE OSTI2020

Electrochemical Oxidation of Organic Molecules at Lower Overpotential: Accessing Broader Functional Group Compatibility with Electron-Proton Transfer Mediators

CONSPECTUS: Electrochemical organic oxidation reactions are highly appealing because protons are often effective terminal electron acceptors, thereby avoiding undesirable stoichiometric oxidants. These reactions are often plagued by high overpotentials, however, that greatly limit their utility. Single-electron transfer (SET) from organic molecules generates high energy radical-cations. Formation of such intermediates often requires electrode potentials far above the thermodynamic potentials of the reaction and frequently causes decomposition and/or side reactions of ancillary functional groups. In this Account, we show how electrocatalytic electron-proton transfer mediators (EPTMs) address this challenge. EPTMs bypass the formation of radical-cation intermediates by supporting mechanisms that operate at electrode potentials much lower (=1 V) than analogous direct electrolysis reactions. The stable aminoxyl radical TEMPO (2,2,6,6-tetramethylpiperidine N-oxyl) is an effective mediator for electrochemical alcohol oxidation, and we have employed such processes for applications ranging from pharmaceutical synthesis to biomass conversion. A complementary electrochemical alcohol oxidation method employs a cooperative Cu/TEMPO mediator system that operates at 0.5 V lower electrode potential than the TEMPO-only mediated process. This difference, which arises from a different catalytic mechanism, rationalizes the broad functional group tolerance of Cu/TEMPO-based aerobic alcohol oxidation catalysts. Aminoxyl mediators address long-standing challenges in the "Shono oxidation", an important method for a-C–H oxidation of tertiary amides and carbamates. Shono oxidations are initiated by a high-potential SET step that limits their utility. Aminoxyl-mediated Shono-type oxidations have been developed that operate at much lower potentials and tolerate diverse functional groups. Analogous reactivity underlies a-C–H cyanation of secondary cyclic amines, a new method that enables efficient diversification of piperidine-based pharmaceutical building blocks and preparation of non-natural amino acids. Electrochemical oxidations of benzylic C–H bonds are commonly initiated by SET to generate radical cations, but such methods are again plagued by large overpotentials. Mediated electrolysis methods that promote hydrogen-atom-transfer (HAT) from benzylic C–H bonds to Fe-oxo species and phthalimide N-oxyl (PINO) support C–H oxygenation, iodination, and oxidative-coupling reactions. A complementary method merges photochemistry with electrochemistry to achieve amidation of C(sp3)–H bonds. This unique process operates at much lower overpotentials compatible with diverse functional groups. These results have broad implications for organic electrochemistry, highlighting the importance of "overpotential" considerations and the prospects for expanding synthetic utility by using mediators to bypass high-energy outer-sphere electron-transfer mechanisms. Principles demonstrated here for oxidation are equally relevant to electrochemical reductions.

Wang, Fei↗
DOE OSTI2026

Quantifying Temperature Dependence of Pu(IV) Absorbance Spectra for Advanced Online Monitoring of Nuclear Processes

This article presents a systematic study of Pu(IV) absorbance spectral features as a function of temperature to develop an understanding of this parameter’s effect on chemometric models that can be used as online monitoring tools to support nuclear processing. The descriptive and predictive models that provide real-time feedback of these processes are usually constructed with data collected in conditions typical of a laboratory environment, which can differ drastically from a processing environment. To assess the impact of temperature on Pu(IV) absorbance spectra, 11 samples of Pu(IV) were synthesized with varying HNO 3 concentrations ranging from 0.6 to 9.5 M and heated between 15 and 45 °C. Ultraviolet (UV)–visible (vis)–near-infrared (NIR) absorption spectra collected at different HNO 3 concentrations and temperatures revealed that features associated with Pu(IV) are sensitive to temperature at all HNO 3 concentrations and that changes in features depend on HNO 3 concentration. The contributions of temperature and HNO 3 concentration to variation in Pu(IV) spectral features were evaluated using the principal component analysis of spectra that were baseline-corrected with an asymmetric least-squares method. Furthermore, predictive modeling for HNO 3 concentration with partial least-squares regression of UV–vis–NIR spectra highlighted the importance of accounting for temperature in the calibration set to optimize model performance. This methodology constitutes a new, systematic approach to account for the effect of temperature on the absorption spectra of metal ions and is useful for process monitoring applications in many industries.

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

Modeling Short-Range Order in Disordered Rocksalt Cathodes by Pair Distribution Function Analysis

Pair distribution function (PDF) analysis is a powerful technique for the characterization of short-range order (SRO) in disordered materials. Accurate interpretation of experimental PDF data is critically reliant on the development of structural models that can account for local variations in site occupancies and bond lengths. To this end, we outline an approach to model SRO using first-principles calculations based on the cluster-expansion formalism. These methods are validated on neutron scattering data from two disordered rocksalt oxyfluorides, Li 1.3 Mn 0.4 Ti 0.3 O 1.7 F 0.3 and Li 1.3 Mn 0.4 Nb 0.2 Ti 0.1 O 1.7 F 0.3 . For each composition, we demonstrate that an average structure without any SRO fails to reproduce several key features in the experimental PDF. To pinpoint the origin of the suspected SRO in these materials, configurational and displacive effects were separately investigated using two disparate models. Special quasi-random structures were relaxed using density functional theory to account for local changes in bond lengths while maintaining a near-random ionic configuration. This leads to slightly improved accuracy but still misrepresents asymmetry in the first few peaks of the PDF. Monte Carlo simulations were performed to model configurational SRO on a fixed lattice, which by itself is shown to have a minimal influence on the PDF. Instead, we find that it is the bond length relaxations within environments created by SRO which controls the details of the PDF, thereby highlighting the subtle but important coupling between configurational and displacive SRO in disordered materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Methane Exhaust Measurements at Gathering Compressor Stations in the United States

Unburned methane entrained in exhaust from natural gas-fired compressor engines (“combustion slip”) can account for a substantial portion of station-level methane emissions. A novel in-stack, tracer gas method was coupled with Fourier transform infrared (FTIR) species measurements to quantify combustion slip from natural gas compressor engines at 67 gathering and boosting stations owned or managed by nine “study partner” operators in 11 U.S. states. The mean methane emission rate from 63 four-stroke, lean-burn (4SLB) compressor engines was 5.62 kg/h (95% CI = 5.15–6.17 kg/h) and ranged from 0.3 to 12.6 kg/h. The mean methane emission rate from 39 fourstroke, rich-burn (4SRB) compressor engines was 0.40 kg/h (95% CI = 0.37– 0.42 kg/h) and ranged from 0.01 to 4.5 kg/h. Study results for 4SLB engines were lower than both the U.S. EPA compilation of air pollutant emission factors (AP-42) and Inventory of U.S. Greenhouse Gas Emissions and Sinks (GHGI) by 8 and 9%, respectively. Study results for 4SRB engines were 43% of the AP-42 emission factor and 8% of the GHGI emission factor, the latter of which does not distinguish between engine types. Total annual combustion slip from the U.S. natural gas gathering and boosting sector was modeled using measured emission rates and compressor unit counts from the U.S. EPA Greenhouse Gas Reporting Program. Modeled results [328 Gg/y (95% CI = 235–436 Gg/y) of unburned methane] would account for 24% (95% CI = 17–31%) of the 1391 Gg of methane emissions for “Gathering and Boosting Stations”, or 6% of the net emissions for “Natural Gas Systems” (5598 Gg) as reported in the 2020 U.S. EPA GHGI. In conclusion, gathering and boosting combustion slip emissions reported in the 2020 GHGI (374 Gg) fall within the uncertainty of this model.

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

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗
DOE OSTI2021

Modeled Emission of Hydroxyl and Ozone Reactivity from Evaporation of Fragrance Mixtures

Volatile chemical products (VCPs) account for increasing fractions of organic carbon emitted to the atmosphere, particularly in urban areas. Fragrances are potentially reactive components that are added to many VCPs. To better constrain these emissions, 11 commercially available liquid fragrance mixtures were characterized for their composition and their evaporation modeled. Here, emissions of mass, hydroxyl reactivity, and ozone reactivity were estimated by modeling under four different scenarios. Fragrance compounds were generally less than one-half the mass of fragrance mixtures, with the balance comprised of solvents and plasticizers and unresolved mass thought to be dominated by plasticizers. The results showed that terpenes and terpenoids account for nearly all of the emitted mass and reactivity while only comprising ~10% w/w on average of the liquid fragrance mixtures. Most of the reactivity is emitted within hours, with ozone reactivity evolving more rapidly than OH reactivity and comprised almost entirely of terpenes. Limonene, a common fragrance constituent, dominates the reactivity of emitted carbon. Generally, 20–40% of the potential hydroxyl reactivity contained in the fragrance mixture does not evaporate on time scales sufficient to have an impact on local or regional air quality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Photolysis Controls Atmospheric Budgets of Biogenic Secondary Organic Aerosol

Secondary organic aerosol (SOA) accounts for a large fraction of tropospheric particulate matter. While SOA production rates and mechanisms have been extensively investigated, loss pathways remain uncertain. Most large-scale chemistry and transport models account for mechanical deposition of SOA, but not chemical losses, such as photolysis. There is also a paucity of laboratory measurements of SOA photolysis, which limits how well photolytic losses can be modeled. Here we show, through a combined experimental and modeling approach, that photolytic loss of SOA mass significantly alters SOA budget predictions. Using environmental chamber experiments, we find that SOA produced from several biogenic volatile organic compounds (VOCs) undergoes photolysis-induced mass loss at rates between 0 and 2.2% ± 0.4% of nitrogen dioxide (NO 2 ) photolysis, equivalent to average atmospheric lifetimes as short as 10 hours. We incorporate our photolysis rates into a regional chemical transport model to test the sensitivity of predicted SOA mass concentrations to photolytic losses. The addition of photolysis at rates consistent with our experimental conditions causes a ~50% reduction in biogenic SOA loadings over the Amazon, indicating that photolysis exerts a substantial control over the atmospheric SOA lifetime, with a likely dependence upon SOA molecular composition and thus production mechanisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

What Factors Drive the Changes in Water Withdrawals in the U.S. Agriculture and Food Manufacturing Industries between 1995 and 2010?

Climate change and increasing world population will directly impact the global food supply chain linkages. In the United States, agricultural production requires less irrigated water than before but it still accounts for a third of total water withdrawals. To better understand the evolution of its water use, we perform a structural decomposition analysis of water withdrawals across eight different crops and six livestock categories and differentiate the trends over 1995–2005 vs 2005–2010 to account for the role of the economic crisis in the second period. Based on USGS data, the results show that both periods experienced an overall decline in water withdrawals in the production of all crops except oilseeds. This trend is driven by a decrease in water intensity, reflecting greater efficiency of irrigation systems, and by reduced local per capita income in the second period. However, increased foreign demand for water-intensive sectors like oilseeds from NAFTA and Asian partners mitigated the decline. Results indicate also a decreasing water use in livestock production partially due to a shift from red to white meat consumption in the country. Arguably, recent tariff wars and border closures have greatly reduced the virtual water embodied in American exports.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

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

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. Here, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗