SEARCH · Engineering Papers
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.
A simplified theory for heavy ion scattering shows good agreement with heavy ion absorption experiments. Theoretical implications on the complete coupled channel reaction equations are discussed.
Quantum computing leverages the principles of quantum mechanics in novel ways to tackle complex chemistry problems that cannot be accurately addressed using traditional quantum chemistry methods. However, the high computational cost and available number of physical qubits with high fidelity limit its application to small chemical systems. This work employed a quantum-classical framework which features a quantum active space-embedding approach to perform simulations of chemical reactions that require up to 14 qubits. This framework was applied to prototypical example metal hydrogenation reactions: the coupling between hydrogen and Li 2 , Li 3 , and Li 4 clusters. Particular attention was paid to the computation of barriers and reaction energies. The predicted reaction profiles compare well with advanced classical quantum chemistry methods, demonstrating the potential of the quantum embedding algorithm to map out reaction profiles of realistic gas-phase chemical reactions to ascertain qualitative energetic trends. Additionally, the predicted potential energy curves provide a benchmark to compare against both current and future quantum embedding approaches.
The past couple of decades have seen tremendous advances in nuclear structure and reaction theory. Innovative theory frameworks for describing the nuclear many-body system, increasingly powerful computers, and opportunities for confronting theory predictions with data on unstable nuclei, have been driving the field. An important goal is to move from phenomenological ingredients in reaction calculations to predictive theories based on microscopic frameworks. We discuss ongoing efforts aimed at integrating microscopic descriptions of nuclear structure into reaction predictions for medium-mass and heavy nuclei. This contribution highlights areas where Eric Bauge, a champion for building bridges, has made important contributions by encouraging and enabling collaborations between communities with complementary expertise.
Hyoscyamine 6β-hydroxylase (H6H) is an Fe(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase that catalyzes the last two steps in the biosynthesis of scopolamine, a prolifically administered anti-nausea drug. After its namesake first reaction, H6H couples the newly installed C6-bonded oxygen to C7 to form the epoxide of scopolamine. Oxoiron(IV) (ferryl) intermediates initiate both reactions by cleaving C–H bonds, but it remains unclear how the enzyme switches target site and promotes (C6)O–C7 coupling in preference to C7 hydroxylation in the second step. In one possible epoxidation mechanism, the C6 oxygen would – analogously to mechanisms proposed for the Fe/2OG halogenases and, in the preceding paper, N-acetylnorloline synthase (LolO) – coordinate as alkoxide to the C7–H-cleaving ferryl intermediate to enable alkoxyl coupling to the ensuing C7 radical. Here we provide structural and kinetic evidence that H6H instead exploits the distinct spatial dependencies of competitive C–H-cleavage (C6 vs C7) and C–O-coupling (oxygen rebound vs cyclization) steps to promote the two-step sequence without substrate coordination or repositioning for the epoxidation step. Structural comparisons of ferryl-mimicking vanadyl complexes of wild-type H6H and a variant that preferentially hydroxylates C7 of 6-hydroxyhyoscyamine suggest that only a modest (~ 10°) shift in the Fe–O–H(C7) approach angle is sufficient to change the outcome. Finally, the observation that, in wild-type H6H, 2 H 2 O solvent also increases the C7-hydroxylation:epoxidation ratio by ~ 8-fold implies that the latter outcome requires cleavage of the alcohol O-H bond, which, unlike in the LolO oxacyclization, is not accomplished in advance of C–H cleavage.
The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.
Electrochemical conversion of CO2 to liquid products is limited by parasitic reactions that consume products and destabilize operation. Here, we show that salt-free formic acid synthesis in a forward-bias bipolar-membrane CO2 electrolyzer is governed by a coupled reaction-transport instability driven by parasitic anode formic acid oxidation. Operando mass spectrometry reveals that this process proceeds through CO-mediated poisoning, identifying catalyst tolerance to CO as a key descriptor of anode selectivity. Among the catalysts examined, PtRu/C preferentially promotes hydrogen oxidation over formic acid oxidation. Coupling this catalyst with transport-interface engineering improves product removal, enhances Faradaic and energy efficiency, and enables stable operation at 200 mA/cm2 for 190 h with a voltage decay of 0.64 mV/h. Techno-economic analysis indicates an 18% reduction in minimum selling price and highlights product concentration as the dominant cost lever. These results establish a general framework for suppressing product oxidation in liquid-product electrolyzers.
The discovery of chirality-induced spin selectivity (CISS) revolutionized our understanding of the capabilities of chiral molecules, revealing that chiral molecules can function as spin filters, aligning the spin orientation of electrons when they transmit through them. Recently, CISS has been exploited to direct energy conversion, especially the oxygen evolution reaction (OER). However, despite the remarkable progress that has been achieved, the effect of CISS in influencing the intermediate species formation and changing the rate-determining step (RDS) is still vague. To understand those key reaction mechanism steps, electrocatalysts with distinct magnetic characteristics, ferromagnetic CoFe2O4 and paramagnetic Co3O4, were synthesized. The results show that spin-polarized charge carriers retain their spin alignment when coupled with ferromagnetic CoFe2O4, akin to the behavior observed under a magnetic field. The Tafel analysis and kinetic isotope studies (kinetic isotope effect) suggest that in the absence of chiral molecules, the initial electron transfer step, the formation of O* species, governs the rate-determining step (RDS). However, introducing chiral molecules shifts the RDS to a combination of the first and second electron transfers, leading to the formation of OOH*. This conclusion was further supported by in situ infrared spectroscopy, which shows that l-methionine-modified CoFe2O4 (l-CoFe2O4) promotes the formation of OOH*, a key intermediate for O2 generation. This study highlights the critical role of CISS in affecting the OER mechanism and intermediate species formation.
Materials with hierarchical architectures that combine soft and hard material domains with coalesced interfaces possess superior properties compared with their homogeneous counterparts. These architectures in synthetic materials have been achieved through deterministic manufacturing strategies such as 3D printing, which require an a priori design and active intervention throughout the process to achieve architectures spanning multiple length scales. Here we harness frontal polymerization spin mode dynamics to autonomously fabricate patterned crystalline domains in poly(cyclooctadiene) with multiscale organization. This rapid, dissipative processing method leads to the formation of amorphous and semi-crystalline domains emerging from the internal interfaces generated between the solid polymer and the propagating cure front. The size, spacing and arrangement of the domains are controlled by the interplay between the reaction kinetics, thermochemistry and boundary conditions. Small perturbations in the fabrication conditions reproducibly lead to remarkable changes in the patterned microstructure and the resulting strength, elastic modulus and toughness of the polymer. Furthermore, this ability to control mechanical properties and performance solely through the initial conditions and the mode of front propagation represents a marked advancement in the design and manufacturing of advanced multiscale materials. Drawing inspiration from biological systems in which structural complexity develops through dissipative reaction–diffusion processes, this study explores a transformative synthetic manufacturing strategy aimed at harnessing the principles underpinning morphogenic growth, unlocking new avenues for advanced materials design and fabrication. Synthetic coupled reaction-transport processes offer a versatile yet relatively underexplored method to manipulate the spatial attributes of synthetic materials10. Here we introduce an innovative manufacturing approach based on frontal ring-opening metathesis polymerization (FROMP) that draws parallels with morphogenic growth and development, enabling the formation of patterned microstructures within polymeric materials.
The thermodynamics of organic chemistry under mild aqueous conditions was examined in order to begin to understand its influence on the structure and operation of metabolism and its antecedents. Free energies were estimated for four types reactions of biochemical importance carbon-carbon bond cleavage and synthesis, hydrogen transfer between carbon groups, dehydration of alcohol groups, and aldo-keto isomerization. The energies were calculated for mainly aliphatic groups composed of carbon, hydrogen, and oxygen. The energy values showed that (1) when carbon-carbon bond cleavage involves two different types of functional groups, transfer of the shared electron-pair to the more reduced carbon group is energetically favored over transfer to the more oxidized carbon group, and (2) the energy of carbon-carbon bond transformation is strongly dependent on the type of functional group that donates the shared electron-pair during cleavage, and the group that accepts the shared electron-pair during synthesis, and (3) the energetics of C-C bond transformation is determined primarily by the half-reaction energies of the couples: carbonyl/carboxylic acid, carboxylic acid/carbon dioxide, alcohol/carbonyl, and hydrocarbon/alcohol. The energy of hydrogen-transfer between carbon groups was found to depend on the functional group class of both the hydrogen-donor and hydrogen-acceptor. From these and other observations we concluded that the chemistry of the origin of metabolism (and to a lesser degree modem metabolism) is strongly constrained by the (1) limited disproportionation energy of organic substrates that can be dissipated in a few irreversible reactions, (2) the energy-dominance of few half-reaction couples in carbon-carbon bond transformation that establishes whether a chemical reaction is energetically irreversible, reversible or unfeasible, and (3) the dependence of the transformation-energy on the oxidation state of carbon groups (functional group type) which is contingent on prior reactions in the synthetic pathway.
The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.
In this paper, we present an optimal open-loop slew trajectory algorithm developed at GSFC for the so-called "Yardstick design" of the James Webb Space Telescope (JWST). JWST is an orbiting infrared observatory featuring a lightweight, segmented primary mirror approximately 6 meters in diameter and a sunshield approximately the size of a tennis court. This large, flexible structure will have significant number of lightly damped, dominant flexible modes. With very stringent requirements on pointing accuracy and image quality, it is important that slewing be done within the required time constraint and with minimal induced vibration in order to maximize observing efficiency. With reaction wheels as control actuators, initial wheel speeds as well as individual wheel torque and momentum limits become dominant constraints in slew performance. These constraints must be taken into account when performing slews to ensure that unexpected reaction wheel saturation does not occur, since such saturation leads to control failure in accurately tracking commanded motion and produces high frequency torque components capable of exciting structural modes. A minimum-time constraint is also included and coupled with reaction wheel limit constraints in the optimization to minimize both the effect of the control torque on the flexible body motion and the maneuver time. The optimization is on slew command parameters, such as maximum slew velocity and acceleration, for a given redundant reaction wheel configuration and is based on the dynamic interaction between the spacecraft and reaction wheel motion. Analytical development of the slew algorithm to generate desired slew position, rate, and acceleration profiles to command a feedback/feed forward control system is described. High-fidelity simulation and experimental results are presented to show that the developed slew law achieves the objectives.
Heterobimetallic complexes are studied for their ability to mimic biological systems as well as active sites in heterogeneous catalysts. While specific interest in early/late heterobimetallic systems has fluctuated, they serve as important models to fundamentally understand metal–metal bonding. Specifically, the polarized metal–metal multiple bonds formed in highly reduced early/late heterobimetallic complexes exemplify how each metal modulates the electronic environment and reactivity of the complex as a whole. In this Perspective, we chronicle the development of phosphinoamide-supported group IV/cobalt heterobimetallic complexes. This combination of metals allows access to a low valent Co −I center, which performs a rich variety of bond activation reactions when coupled with the pendent Lewis acidic metal center. Conversely, the low valent late transition metal is also observed to act as an electron reservoir, allowing for redox processes to occur at the d 0 group IV metal site. Most of the bond activation reactions carried out by phosphinoamide-bridged M/Co −I (M = Ti, Zr, Hf) complexes are facilitated by cleavage of metal–metal multiple bonds, which serve as readily accessible electron reservoirs. Comparative studies in which both the number of buttressing ligands as well as the identity of the early metal were varied to give a library of heterobimetallic complexes are summarized, providing a thorough understanding of the reactivity of M/Co −I heterobimetallic systems
A series of detailed studies comparing various vibration dissociation coupling models, reaction systems and rates, and radiative heating models has been conducted for the nonequilibrium stagnation region of an AFE/AOTV vehicle. Atomic and molecular nonequilibrium radiation correction factors have been developed and applied to various absorption coefficient step models, and a modified vibration dissociation coupling model has been shown to yield good vibration/electronic temperature and concentration profiles. While results indicate sensitivity to the choice of vibration dissociation coupling model and to the nitrogen electron impact ionization rate, by proper combinations accurate flowfield and radiative heating results can be obtained. These results indicate that nonequilibrium effects significantly affect the flowfield and the radiative heat transfer. However, additional work is needed in ionization chemistry and absorption coefficient modeling.
High-loading lithium-sulfur (Li-S) cells operated with lean electrolyte are limited by polysulfide crossover to Li metal and by transport-limited liquid-solid conversion that forms passivating Li 2 S films. Here, we show that a cathode-facing separator coating of carboxylated multiwalled carbon nanotubes acts as a cathode-confined polysulfide reservoir with intermediate binding. Early in discharge it captures newly generated polysulfides at the separator interface, suppressing shuttle reactions. As polysulfides are consumed, the reservoir buffers concentration gradients and feeds reactants back to the cathode, shifting Li 2 S deposition from burst-like film growth to progressive, three-dimensional, porous formation. Synchrotron XRD and S K-edge XANES, together with Scharifker–Hills nucleation analysis and depth-of-discharge EIS/DRT, substantiate this coupled transport–reaction control. With 4.3 mg S cm -2 and E/S = 5, cells reach 4.2 mAh cm -2 and retain 90% capacity over 100 cycles at 20 °C.
The high-pressure carbon monoxide (HiPco) technique for producing single-wall carbon nanotubes (SWNTs) is analyzed with the use of a chemical reaction model coupled with flow properties calculated along streamlines, calculated by the FLUENT code for pure carbon monoxide. Cold iron pentacarbonyl, diluted in CO at about 30 atmospheres, is injected into a conical mixing zone, where hot CO is also introduced via three jets at 30 degrees with respect to the axis. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Then iron nucleates and forms clusters that catalyze the formation of SWNTs by a disproportionation reaction (Boudouard) of CO on Fe-containing clusters. Alternative nucleation rates are estimated from the theory of hard sphere collision dynamics with an activation energy barrier. The rate coefficient for carbon nanotube growth is estimated from activation energies in the literature. The calculated growth was found be about an order of magnitude greater than measured, regardless of the nucleation rate. A study of cluster formation in an incubation zone prior to injection into the reactor shows that direct dimer formation from Fe atoms is not as important as formation via an exchange reaction of Fe with CO in FeCO.
The high-pressure carbon monoxide (HiPco) technique for producing single wall carbon nanotubes (SWNT) is analyzed using a chemical reaction model coupled with properties calculated along streamlines. Streamline properties for mixing jets are calculated by the FLUENT code using the k-e turbulent model for pure carbon monixide. The HiPco process introduces cold iron pentacarbonyl diluted in CO, or alternatively nitrogen, at high pressure, ca. 30 atmospheres into a conical mixing zone. Hot CO is also introduced via three jets at angles with respect to the axis of the reactor. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Cluster reaction rates are from Krestinin, et aI., based on shock tube measurements. Another model is from classical cluster theory given by Girshick's team. The calculations are performed on streamlines that assume that a cold mixture of Fe(CO)5 in CO is introduced along the reactor axis. Then iron forms clusters that catalyze the formation of SWNTs from the Boudouard reaction on Fe-containing clusters by reaction with CO. To simulate the chemical process along streamlines that were calculated by the fluid dynamics code FLUENT, a time history of temperature and dilution are determined along streamlines. Alternative catalyst injection schemes are also evaluated.
Solid-state synthesis involves a web of coupled chemical reactions and physical changes that unfold across multiple scales. Efforts to fine-tune its parameters have historically followed heuristic, trial-driven workflows that demand significant time and resources. In this study, we aimed to open this black box by employing multiscale in situ synchrotron imaging and diffraction. Using LiNi 0.5 Mn 0.3 Co 0.2 O 2 battery positive electrode material as a model system and Ba-based sintering aids, we reveal dopant segregation, intergranular mass transport, and porosity evolution as key drivers of single-crystalline particle formation. Notably, we uncovered a dynamic competition between particle-level grain coalescence and atomic-scale cation disordering, both of which are thermally activated yet have opposing impacts on battery performance. These findings highlight the coupled, multiscale nature of structure development and offer a mechanistic basis for optimizing the solid-state synthesis process. This framework provides a path toward more controlled, efficient, and scalable production of high-performance battery positive electrode materials.
The $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction rate is crucial in determining the carbon-to-oxygen abundance ratio in stellar nucleosynthesis. Measuring this reaction’s cross section at stellar energies is challenging due to its extremely small value, approximately 10 -17 barn at E c.m. = 300 keV. To address this, R-matrix calculations are employed to extrapolate data to lower energies, requiring a comprehensive understanding of each contribution to the cross section. The dominant contributions to the cross section at stellar energies arise from electric dipole (E1) and electric quadrupole (E2) transitions to the ground state of 16 O, along with a significant cascade contribution. Traditionally, these contributions have been separated using the γ-ray angular distribution. In this work, we propose a novel technique using the energy distribution of the 16 O recoils at the focal plane. This method involves a neural network trained on detailed Monte Carlo simulations of the energy distribution of recoils transported through the recoil mass separator ERNA. This approach enables the simultaneous determination of all three contributions with errors around 10% in the energy range E c.m. = 1.0–2.2 MeV. In conclusion, by employing this new technique, we aim to significantly improve the accuracy of determining the cross section of the $^{12}$C($\alpha $,$\gamma $)$^{16}$O reaction at astrophysical energies.