SEARCH · Engineering Papers

Results for “Transition state”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Development of Physics-Based Transition Models for Unstructured-Mesh CFD Codes Using Deep Learning Models

Predicting transition locations over a vehicle surface is of fundamental importance for many engineering applications. With the transition information, the Reynolds-averaged Navier-Stokes (RANS) computations can turn on the turbulence model at the right locations so that drag, lift and other aerodynamic quantities can be accurately predicted. In contrast to the popularity of RANS-based transition modeling in which transition onset is governed by the turbulence equations, physics-based transition models that account for instability waves within the boundary layer, thus more compliant to flow physics, only gained more attention in recent years. This paper describes the development of a new physics-based transition model based on either the linear stability theory (LST) or parabolized stability equations (PSE). The model is designed to communicate with a structured or unstructured-mesh RANS solver back and forth in order to more accurately compute transition fronts over a three-dimensional body. In the developed model, the Python suite of interface codes in conjunction with the LASTRAC software can be executed autonomously to produce transition onset locations for a given laminar or RANS-computed transitional state. In addition, as a proof of concept, the tool set consists of a deep learning neural network model that has been designed and trained to predict instability wave evolutions inside the boundary layer for various instability wave mechanisms across a selected speed range. A machine-learned intelligent profile interpolation model has also been devised to enable reliable instability-wave spectra predictions with just a few points in the mean flow profiles.

Transition Modeling↗

NASA NTRS1988

Large-angle slewing maneuvers for flexible spacecraft

A new class of closed-form solutions for finite-time linear-quadratic optimal control problems is presented. The solutions involve Potter's solution for the differential matrix Riccati equation, which assumes the form of a steady-state plus transient term. Illustrative examples are presented which show that the new solutions are more computationally efficient than alternative solutions based on the state transition matrix. As an application of the closed-form solutions, the neighboring extremal path problem is presented for a spacecraft retargeting maneuver where a perturbed plant with off-nominal boundary conditions now follows a neighboring optimal trajectory. The perturbation feedback approach is further applied to three-dimensional slewing maneuvers of large flexible spacecraft. For this problem, the nominal solution is the optimal three-dimensional rigid body slew. The perturbation feedback then limits the deviations from this nominal solution due to the flexible body effects. The use of frequency shaping in both the nominal and perturbation feedback formulations reduces the excitation of high-frequency unmodeled modes. A modified Kalman filter is presented for estimating the plant states.

Chun, Hon M.↗

DOE OSTI2024

Proton-Responsive Ligands Promote CO 2 Capture and Accelerate Catalytic CO 2 /HCO 2 – Interconversion

The synthesis and investigation of [Rh(DHMPE) 2 ][BF 4 ] (1) are reported. 1 features proton-responsive 1,2-bis[(dihydroxymethyl)phosphino]ethane (DHMPE) ligands, which readily capture CO 2 from atmospheric sources upon deprotonation. The protonation state of the DHMPE ligand was observed to have a significant impact on the catalytic reactivity of 1 with CO 2 . Deprotonation and CO 2 binding to 1 result in a ∼10-fold rate enhancement in catalytic degenerate CO 2 reduction with formate, monitored by 12 C/ 13 C isotope exchange between H 12 CO 2 – and 13 CO 2 . Studies performed using a similar complex lacking the hydroxyl ligand functionality ([Rh(DEPE) 2 ][BF 4 ] where DEPE = 1,2-bis(diethylphosphino)ethane) do not show the same rate enhancements when base is added. Based upon the cation-dependent activity of the catalyst, Eyring analysis, and cation sequestration experiments, CO 2 binding to 1 is proposed to facilitate preorganization of formate/CO 2 in the transition state via ligand-based encapsulation of Na + or K + cations to lower the activation energy and increase the observed catalytic rate. Incorporation of proton-responsive DHMPE ligands provides a unique approach to accelerate the kinetics of catalytic CO 2 reduction to formate.

Anions↗

DOE OSTI2022

Deciphering Distinct Overpotential-Dependent Pathways for Electrochemical CO 2 Reduction Catalyzed by an Iron–Terpyridine Complex

[Fe(tpyPY2Me)] 2+ ([Fe] 2+ ) is a homogeneous electrocatalyst for converting CO 2 into CO featuring low overpotentials of <100 mV, near-unity selectivity, and high activity with turnover frequencies faster than 100 000 s -1 . To identify the origins of its exceptional performance and inform future catalyst design, we report a combined computational and experimental study that establishes two distinct mechanistic pathways for electrochemical CO 2 reduction catalyzed by [Fe] 2+ as a function of applied overpotential. Electrochemical data shows the formation of two catalytic regimes at low (η TOF/2 of 160 mV) and high (η TOF/2 of 590 mV) overpotential plateaus. We propose that at low overpotentials [Fe] 2+ undergoes a two-electron reduction, two-proton-transfer mechanism (electrochemical-electrochemical-chemical-chemical, EECC), where turnover occurs through the dicationic iron complex, [Fe] 2+ . Computational analysis supports the importance of the singlet ground-state electronic structure for CO 2 binding and that the rate-limiting step is the second protonation in this low-overpotential regime. When more negative potentials are applied, an additional electron-transfer event occurs through either a stepwise or proton-coupled electron-transfer (PCET) pathway, enabling catalytic turnover from the monocationic iron complex ([Fe] + ) via an electrochemical-chemical-electrochemical-chemical (ECEC) mechanism. Comparison of experimental kinetic data obtained from variable controlled potential electrolysis (CPE) experiments with direct product detection with calculated rates obtained from the energetic span model supports the PCET pathway as the most likely mechanism. Moreover, we build upon this mechanistic understanding to propose the design of an improved ligand framework that is predicted to stabilize the key transition states identified in our study and explore their electronic structures using an energy decomposition analysis. Finally, taken together, this work highlights the value of synergistic computational/experimental approaches to decipher mechanisms of new electrocatalysts and direct the rational design of improved platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Unimolecular and bimolecular formic acid decomposition routes on dispersed Cu nanoparticles

The elementary steps and site requirements in formic acid (HCOOH) dehydrogenation on Cu surfaces remain of keen interest because formate species act as intermediates or spectators in methanol synthesis and water–gas shift reactions. Steady-state and transient kinetic data, isotopic effects, infrared spectra during catalytic and stoichiometric reactions, and theoretical treatments based on density functional theory (DFT) provide evidence for bimolecular reactions, in which saturated bidentate formate (*HCOO*) adlayers, present at 0.25 ML (0.25 *HCOO* per surface Cu atom), react with undissociated species (HCOOH$^\square$) bound at interstices within formate adlayers ($^\square$) to form H-bonded bimolecular HCOOH$^\square$-*HCOO* adducts. The co-existence of vicinal HCOOH$^\square$ and *HCOO* moieties is evident from antisymmetric infrared bands for *HCOO* that become stronger as a result of their H-bonding that perturbs the induced dipole moment of *HCOO* upon vibration, consistent with DFT-derived vibrational frequencies and intensities for such perturbed species. The *HCOO* moiety in this complex undergoes C-H activation via a transition state that is preferentially stabilized through H-bonding with the vicinal HCOOH$^\square$ relative to its *HCOO* precursor. DFT-derived HCOOH dehydrogenation activation barriers and those determined from the evolution of CO 2 from pre-adsorbed *HCOO* species are about 10 kJ mol -1 smaller in the presence of gaseous HCOOH reactants (because of HCOOH$^\square$-*HCOO* interactions) than those for the unimolecular decomposition of bound *HCOO* species. Such bimolecular routes are consistent with measured effects of HCOOH, H 2 , and CO pressures and of H/D isotopic substitution on dehydrogenation turnover rates and represent the predominant channel for the formation of CO 2 and H 2 during catalytic HCOOH dehydrogenation on Cu nanoparticles. A saturated *HCOO* adlayer that retains binding interstices and the presence of HCOOH(g) enable a sequence of elementary steps unavailable for *HCOO* species, thus circumventing unassisted unimolecular routes that exhibit higher activation barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

On the Emergence of Ferromagnetism in LaCoO 3 Ultrathin Films

It is well known that the properties of a crystal evolve as it increases in size from a single atomic plane to that of the bulk. Such size-dependent transitions can stem from many different origins and depend on minute changes to crystal bonding and composition. Here, a model example is that of LaCoO 3 , which is non-magnetic in the bulk but can display ferromagnetism at the nanoscale. Here, the evolution of structure-property relationships is studied in the LaCoO 3-δ /SrTiO 3 (001) system as the thickness of LaCoO 3-δ is increased from a single plane to 10 unit cells. In situ synchrotron X-ray studies are performed during and post-deposition to probe changes in the interactions between structure, stoichiometry, and magnetic behavior. Structural quantification indicates that the oxygen octahedral rotation pattern evolves with thickness, due to inherent differences in crystal symmetry between the film and substrate. The change in rotation modifies the required energy barrier for the spin state transition via the Co–O bond length and Co–O–Co bond angle, affecting the appearance of ferromagnetism. Our results highlight the contributions of high spin Co 2+ and/or high spin Co 3+ to respective weak and robust ferromagnetism and the evolution of properties with size in ultrathin LaCoO 3-δ heterostructures.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Spin and charge excitations in the correlated multiband metal Ca 3 Ru 2 O 7

We use Ru $L_3$-edge resonant inelastic x-ray scattering to study the full range of excitations in Ca 3 Ru 2 O 7 from meV-scale magnetic dynamics through to the eV-scale interband transitions. This bilayer 4d-electron correlated metal expresses a rich phase diagram, displaying long-range magnetic order below 56 K followed by a concomitant structural, magnetic, and electronic transition at 48 K. In the low-temperature phase, we observe a magnetic excitation with a bandwidth of ~30 meV and a gap of ~8 meV at the zone center, in excellent agreement with inelastic neutron scattering data. The dispersion can be modeled using a Heisenberg Hamiltonian for a bilayer S = 1 system with single-ion anisotropy terms. At a higher energy loss, dd-type excitations show heavy damping in the presence of itinerant electrons, giving rise to a fluorescencelike signal appearing between the $t_{2g}$ and $e_g$ bands. At the same time, we observe a resonance originating from localized $t_{2g}$ excitations, in analogy to the structurally related Mott insulator Ca 2 RuO 4 . But whereas Ca 2 RuO 4 shows sharp separate spin-orbit excitations and Hund’s-rule driven spin-state transitions, here we identify only a single broad asymmetric feature. These results indicate that local intraionic interactions underlie the correlated physics in Ca 3 Ru 2 O 7 , even as the excitations become strongly mixed in the presence of itinerant electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

An Electron-scattering Time Delay in Black Hole Accretion Disks

Abstract Universal to black hole X-ray binaries, the high-frequency soft lag gets longer during the hard-to-intermediate state transition, evolving from ≲1 to ∼10 ms. The soft lag production mechanism is thermal disk reprocessing of nonthermal coronal irradiation. X-ray reverberation models account for the light-travel time delay external to the disk, but assume instantaneous reprocessing of the irradiation inside the electron-scattering-dominated disk atmosphere. We model this neglected scattering time delay as a random walk within an α -disk atmosphere, with approximate opacities. To explain soft lag trends, we consider a limiting case of the scattering time delay that we dub the thermalization time delay, t th ; this is the time for irradiation to scatter its way down to the effective photosphere, where it gets thermalized, and then scatter its way back out. We demonstrate that t th plausibly evolves from being inconsequential for low mass accretion rates m ̇ characteristic of the hard state, to rivaling or exceeding the light-travel time delay for m ̇ characteristic of the intermediate state. However, our crude model confines t th to a narrow annulus near peak accretion power dissipation, so cannot yet explain in detail the anomalously long-duration soft lags associated with larger disk radii. We call for time-dependent models with accurate opacities to assess the potential relevance of a scattering delay.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Methodology for object-oriented real-time systems analysis and design: Software engineering

Successful application of software engineering methodologies requires an integrated analysis and design life-cycle in which the various phases flow smoothly 'seamlessly' from analysis through design to implementation. Furthermore, different analysis methodologies often lead to different structuring of the system so that the transition from analysis to design may be awkward depending on the design methodology to be used. This is especially important when object-oriented programming is to be used for implementation when the original specification and perhaps high-level design is non-object oriented. Two approaches to real-time systems analysis which can lead to an object-oriented design are contrasted: (1) modeling the system using structured analysis with real-time extensions which emphasizes data and control flows followed by the abstraction of objects where the operations or methods of the objects correspond to processes in the data flow diagrams and then design in terms of these objects; and (2) modeling the system from the beginning as a set of naturally occurring concurrent entities (objects) each having its own time-behavior defined by a set of states and state-transition rules and seamlessly transforming the analysis models into high-level design models. A new concept of a 'real-time systems-analysis object' is introduced and becomes the basic building block of a series of seamlessly-connected models which progress from the object-oriented real-time systems analysis and design system analysis logical models through the physical architectural models and the high-level design stages. The methodology is appropriate to the overall specification including hardware and software modules. In software modules, the systems analysis objects are transformed into software objects.

Schoeffler, James D.↗

DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Real-Time Diagnostics of 2D Crystal Transformations by Pulsed Laser Deposition: Controlled Synthesis of Janus WSSe Monolayers and Alloys

Energetic processing methods such as hyperthermal implantation hold special promise to achieve the precision synthesis of metastable two-dimensional (2D) materials such as Janus monolayers; however, they require precise control. Here, we report a feedback approach to reveal and control the transformation pathways in materials synthesis by pulsed laser deposition (PLD) and apply it to investigate the transformation kinetics of monolayer WS 2 crystals into Janus WSSe and WSe 2 by implantation of Se clusters with different maximum kinetic energies (<42 eV/Se-atom) generated by laser ablation of a Se target. Real-time Raman spectroscopy and photoluminescence are used to assess the structure, composition, and optoelectronic quality of the monolayer crystal as it is implanted with well-controlled fluxes of selenium for different kinetic energies that are regulated with in situ ICCD imaging, ion probe, and spectroscopy diagnostics. First-principles calculations, XPS, and atomic-resolution HAADF STEM imaging are used to understand the intermediate alloy compositions and their vibrational modes to identify transformation pathways. The real-time kinetics measurements reveal highly selective top-layer conversion as WS 2 transforms through WS 2(1–x) Se 2x alloys to WSe 2 and provide the means to adjust processing conditions to achieve fractional and complete Janus WSSe monolayers as metastable transition states. The general approach demonstrates a real-time feedback method to achieve Janus layers or other metastable alloys of the desired composition, and a general means to adjust the structure and quality of materials grown by PLD, addressing priority research directions for precision synthesis with real-time adaptive control.

2D materials↗

DOE OSTI2023

Clean Energy Technology Applications on US Mine Land: Technical Analysis

As the United States transitions toward a clean energy economy, an opportunity exists for redeveloping the more than 17,000 mine land sites located across the nation with clean energy technologies, which have a combined potential for generating more than 85 GW of clean electricity. This report provides an overview of the potential of demonstrating and deploying clean energy projects on current and former mine land. Clean energy project refers to a project that demonstrates one or more of the following technologies: solar; microgrids; geothermal; direct air capture; fossil-fueled electricity generation with carbon capture, utilization, and sequestration; energy storage, including pumped storage hydropower and compressed air energy storage; and advanced nuclear technologies. The report discusses the following technologies and their potential for creating jobs and generating tax revenue that would result in direct and indirect benefits to the local economy: Solar photovoltaics (PV) is being developed on current and former mine land in various parts of the world, including the United States. This approach is attractive because it requires limited infrastructure investment and would utilize the bare surfaces of mines and tailing ponds. Solar resource availability may be greater in the southern regions, including the Interior and Appalachian Basins and the southwestern United States. However, since some mine land sites include areas of significant change in elevation, the deployment of PV on mine land may require sophisticated planning to account for shading and irradiance, or may require regrading of the areas. PV does not create significant environmental risks and generally does not face public resistance; Geothermal systems are often spatially and genetically associated with ore deposits, and in some cases, they have been discovered while in search for epithermal mineral resources. Numerous diverse geothermal applications have been employed at mine land around the world, including power generation, mineral extraction from geothermal brines, process heating, direct use for other mining operations, and direct use for non-mining operations and subsurface energy storage, including geothermal heat pumps. Case studies highlighting these applications provide key lessons relating to identifying drivers and barriers to geothermal resource deployment and can be used to create screening tools for identifying the types and locations of mine land most amenable to utilizing geothermal resources; Carbon capture, utilization, and sequestration technologies include direct air capture (DAC) and enhanced weathering. DAC technologies include air contactors, regeneration systems, and CO 2 compression systems. Captured CO 2 can be converted to valuable feedstocks or possibly injected into abandoned subsurface mines where it would be absorbed by alkaline rock waste and mine tailings or by the porous minerals along the walls of the mine. DAC systems can be coupled with energy sources such as wind, solar, grid, or geothermal. Many DAC systems require a source of water or steam; however, some are expected to be net producers of water. Local impacts of DAC systems are expected to be low, and are related to land footprint, material disposal, and upstream impacts of energy and material production; Compressed air energy storage is an established energy storage technology in salt caverns. It has the potential for implementation in underground mines by pressurizing and storing a large amount of air using electrical compressors when excess electricity is available. When a need for discharge emerges, the air is used to spin turbines and produce the necessary volume of electricity. Abandoned or unused mine openings, including shafts, adits, access tunnels, and mined workings of any orientation, offer potential for vast amounts of compressed air energy storage if the site characteristics meet operational requirements; Pumped hydropower storage can be implemented in surface and subsurface mines. In surface mine applications, both reservoirs may be located in a mine pit or artificial reservoirs made of excavated materials. In subsurface mines, the lower reservoir may be implemented by waterproofing and flooding mine shafts and tunnels. The water is then pumped from the lower reservoir to the upper reservoir during periods of low load and high production, and it is discharged through the turbines during periods of peak demand. The potential environmental damages associated with acidity of mine water or the presence of toxic chemicals incentivizes the development of closed-loop technologies, in which water circulates inside the pumped hydropower facility without being discharged into the external water basins; Advanced nuclear energy technologies include small modular reactors, which can be deployed locally to produce electricity and heat. Such units require seismic stability and a supply of cooling water, but population constraints may exist in some areas. Therefore, remote mine land could represent an optimal location for siting advanced nuclear energy technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Mechanism for Acetone and Crotonaldehyde Production during Steam Reforming of Ethanol over La 0.7 Sr 0.3 MnO 3–x Perovskite: Evidence for a Shared C4 Aldol Addition Intermediate

Here a mechanistic study was conducted on the catalytic conversion of ethanol over La 0.7 Sr 0.3 MnO 3–x perovskite catalysts in the presence and absence of water. The study sought insights into the path of C–C coupling toward acetone and crotonaldehyde and also into clarifying whether the lack of previous reports of C–C coupling over La 0.7 Sr 0.3 MnO 3–x (100) could be due to a “pressure gap”. Several types of experiments were performed at 400–800 K: flow experiments with a torr range reactant gas flown over La 0.7 Sr 0.3 MnO 3–x powders; ultra-high vacuum experiments with continuous gas exposures to a La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample; and torr range continuous gas exposures to a La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample. When ethanol and water were flown over La 0.7 Sr 0.3 MnO 3–x powders at 400–800 K, the products detected were ethene, acetaldehyde, acetone, crotonaldehyde, CO, CO 2 , and H 2 . Acetone was catalytically produced over both the La 0.7 Sr 0.3 MnO 3–x powder and the La 0.7 Sr 0.3 MnO 3–x (100) single-crystal sample at temperatures of 700–800 K when reaction conditions were on the order of 1 Torr of reactant gas and with an excess of water relative to ethanol (1 ethanol/9 water). Isotopic labeling with deuterium was used to gain insights into the C–C coupling reaction mechanism and paths in species with three and four carbons (C 3 and C 4 species). Additionally, steady-state isotopic transient kinetic analysis (SSITKA) experiments + simulations using carbon labeling of the ethanol feed were performed. Three mechanistic paths were considered for the C–C coupling step: the first two paths, A and B, involve coupling between two intermediates which are both in oxygen vacancies; and the third path, C, involves coupling between one intermediate in an oxygen vacancy and one intermediate outside of an oxygen vacancy. The results suggest that the dominant path to the C 3 product, acetone, depends on the conditions. The less active path (attributed to path A or B) occurs at 600–700 K and involves coupling between two irreversibly bound species. The more active path (attributed to path C) requires an excess of water, becomes dominant at 600–800 K, and involves coupling between one irreversibly bound species and one reversibly bound species. Based on these various observations from experiments and simulations, an elementary step is proposed for acetone formation involving a previously unreported C 4 transition state that is formed after aldol addition. Density functional theory calculations were performed based on this hypothesis, and it confirmed that this specific and previously unreported aldol addition path to acetone does exist and that this path consistent with the experimental data. In this path, C–C formation occurs to create a C 4 intermediate that is bound to an oxygen vacancy, then a hydrogen transfer with C–C bond breaking occurs that results in the production of the acetone molecule. The proposed mechanism is also consistent with the experimental observation that acetone formation has a greater than first-order dependence on the water vapor pressure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

A self-adjusting platinum surface for acetone hydrogenation

We show that platinum displays a self-adjusting surface that is active for the hydrogenation of acetone over a wide range of reaction conditions. Reaction kinetics measurements under steady-state and transient conditions at temperatures near 350 K, electronic structure calculations employing density-functional theory, and microkinetic modeling were employed to study this behavior over supported platinum catalysts. The importance of surface coverage effects was highlighted by evaluating the transient response of isopropanol formation following either removal of the reactant ketone from the feed, or its substitution with a similarly structured species. The extent to which adsorbed intermediates that lead to the formation of isopropanol were removed from the catalytic surface was observed to be higher following ketone substitution in comparison to its removal, indicating that surface species leading to isopropanol become more strongly adsorbed on the surface as the coverage decreases during the desorption experiment. This phenomenon occurs as a result of adsorbate–adsorbate repulsive interactions on the catalyst surface which adjust with respect to the reaction conditions. Reaction kinetics parameters obtained experimentally were in agreement with those predicted by microkinetic modeling when the binding energies, activation energies, and entropies of adsorbed species and transition states were expressed as a function of surface coverage of the most abundant surface intermediate (MASI, C 3 H 6 OH*). It is important that these effects of surface coverage be incorporated dynamically in the microkinetic model (e.g., using the Bragg–Williams approximation) to describe the experimental data over a wide range of acetone partial pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

NASA NTRS1990

Velocity-Resolved Far-Infrared Spectra of [Fe 2]: Evidence for Mixing and Clumping in SN 1987A

We present approx. 400 km/s resolution profiles of the 17.94 and 25.99 micron [Fe II] transitions from SN 1987A at t approx. 400 days after core collapse. These observations used the facility cooled grating spectrometer aboard NASA's Kuiper Airborne Observatory. The two profiles are similar and have FWHM line widths of approx. 2700 km/s. The higher signal-to-noise 18 micron profile is somewhat asymmetric, falling off more steeply on the redshifted side than on the blue. Gaussian fits to the profiles yield an average centroid velocity of 280 +/- 140 km/s relative to the Large Magellanic Cloud. The wings of the profiles extend to velocities is approx. greater than 3000 km/s. This shows that a significant fraction of the iron has been mixed outward into the hydrogen-rich envelope, which has a minimum expansion velocity of 2100-2400 km/s. Both profiles also contain an unresolved 3-5 sigma emission feature on the redshifted wing at nu(LSR) approx. + 3900 km/s. We interpret this feature as emission from a high-velocity clump of material containing approx. 3% of the total iron mass. The total line flux of the 26 micron ground-state transition yields an optically thin, singly ionized iron mass of 0.026 solar mass, relatively independent of the assumed temperature. This is significantly less than the 0.06 Me of Fe+ determined from the decline of the optical light curve and the ionization of measured nickel lines, implying that the iron transitions still have appreciable optical depth. However, because of the small change in the 26 micron line flux from our measurement at 250 days, and the similarity of our profiles to the 1.26 micron [Fe II] profile, most of the emission is believed to originate from optically thin material with a temperature of 4406 +/- 400 K. A comparison of the data with spherically symmetric models indicates a power-law density exponent of -3.2 +/- 1.1 and a minimum expansion velocity of 650 +/- 650 km/s for this optically thin component. The [Fe II] line fluxes and profiles also imply that the remainder of the material has high optical depth and is distributed in clumps throughout the ejecta, rather than being concentrated at low velocities in the center of a smooth density distribution.

Haas, Michael R.↗

DOE OSTI2024

Linear Free Energy Relationships Associated with Hydride Transfer From [(6,6'-R 2 -bpy)Re(CO) 3 H]: A Cautionary Tale in Identifying Hydrogen Bonding Effects in the Secondary Coordination Sphere

Here, six rhenium hydride complexes, [(6,6'-R 2 -bpy)Re(CO) 3 H] (bpy = 2,2'-bipyridine, R = OEt, OMe, NHMe, Me, F, Br), were synthesized. These complexes insert CO 2 to form rhenium formate complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 {OC(O)H}]. All the rhenium formate species were characterized using X-ray crystallography, which revealed that the bpy ligand is not coplanar with the metal coordination plane containing the two nitrogen donors of the bpy ligand but tilted. A solid-state structure of [(6,6'-Me 2 -bpy)Re(CO) 3 H] determined using MicroED also featured a tilted bpy ligand. The kinetics of CO 2 insertion into complexes of the type [(6,6'-R 2 -bpy)Re(CO) 3 H] were measured experimentally and the thermodynamic hydricities of [(6,6'-R 2 -bpy)Re(CO) 3 H] species were determined using theoretical calculations. A Brønsted plot constructed using the experimentally determined rate constants for CO 2 insertion and the calculated thermodynamic hydricities for [(6,6'-R 2 -bpy)Re(CO) 3 H] revealed a linear free energy relationship (LFER) between thermodynamic and kinetic hydricity. This LFER is different to the previously determined relationship for CO 2 insertion into complexes of the type [(4,4'-R 2 -bpy)Re(CO) 3 H]. At a given thermodynamic hydricity, CO 2 insertion is faster for complexes containing a 6,6'-substituted bpy ligand. This is likely in part due to the tilting observed for systems with 6,6'-substituted bpy ligands. Notably, the 6,6'-(NHMe) 2 -bpy ligand could in principle stabilize the transition state for CO 2 insertion via hydrogen bonding. This work shows that if only the rate of CO 2 insertion into [(6,6'-(NHMe) 2 -bpy)Re(CO) 3 H] is compared to [(4,4'-R 2 -bpy)Re(CO) 3 H] systems, the increase in rate could be easily attributed to hydrogen bonding, but in fact all 6,6'-substituted systems lead to faster than expected rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗