Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Organic Photovoltaic Efficiency Predictor: Data-Driven Models for Non-Fullerene Acceptor Organic Solar Cells

In the design of organic solar cells, there has been a need for materials with high power conversion efficiencies. Scharber’s model is commonly used to predict efficiency; however, it exhibits poor performance with new non-fullerene acceptor (NFA) devices, since it was designed for fullerene-based devices. In this work, an empirical model is proposed that can be a more accurate alternative for NFA organic solar cells. Additionally, many screening studies use computationally expensive methods. A model based on using semiempirical simplified time-dependent density functional theory (sTD-DFT) as an alternative method can accelerate the calculations and yield a similar accuracy. The models presented in this paper, termed organic photovoltaic efficiency predictor (OPEP) models, have shown significantly lower errors than previous models, with OPEP/B3LYP yielding errors of 1.53% and OPEP/sTD- DFT of 1.55%. As a result, the proposed computational models can be used for the fast and accurate screening of new high-efficiency NFAs/donor pairs.

14 SOLAR ENERGY↗
DOE OSTI2020

Mechanistic Insights into C–H Borylation of Arenes with Organoiridium Catalysts Embedded in a Microporous Metal–Organic Framework

Organometallic iridium catalysts can be used in conjunction with bispinacolatodiboron (B 2 Pin 2 ) to effect the borylation of a variety of substrates such as arenes, alkanes, heteroarenes, and oxygenates. Recently, efforts have also focused on integrating these catalysts into porous supports, such as metal–organic frameworks (MOFs). While the mechanism of homogeneous borylation systems has been thoroughly investigated experimentally and computationally, analogous studies in MOF-supported iridium catalysts have not been conducted. In this paper, we report the mechanistic investigation of a phenanthroline-iridium catalyst immobilized in the organic linker of Universitetet i Oslo (UiO)-67 (Zr 6 O 4 (OH) 4 (BPDC) 4 (PhenDC) 2 , BPDC = biphenyl-4,4'-dicarboxylate, PhenDC = 1,10-phenanthroline-4,4'-dicarboxylate). By using benzene as a model substrate, variable time normalization analysis (VTNA) of the kinetic data suggested a rate law consistent with zero-order in B 2 Pin 2 , and first-order in arene. A primary kinetic isotope effect (KIE) in the time course of benzene-d 6 borylation also provided complementary information about the role of the arene in the rate-determining step of the reaction. Characterization by techniques such as X-ray absorption spectroscopy (XAS) confirmed the presence of Ir(III), while pair distribution function (PDF) analysis suggested structures containing an Ir–Cl bond, further substantiated by X-ray photoelectron spectroscopy (XPS). Analysis of postcatalysis materials by inductively coupled plasma–optical emission spectroscopy (ICP-OES) revealed low boron accumulation, which may indicate an absence of boron in the resting state of the catalyst. Finally, in comparing borylation of benzene and toluene, a slight selectivity for benzene is observed, which is similar to the analogous homogeneous reaction, indicating the influence of substrate sterics on reactivity.

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

Pathways of Photocatalytic Oxidation of Formic Acid on Dry and Hydrated Anatase TiO 2 Surfaces

The photocatalytic oxidation of formic acid (FA), which is one of the most abundant volatile organic compounds, is a promising air remediation technology inspired by nature. However, the detailed mechanism of this photocatalytic reaction on the surface of TiO 2 , a typical photocatalyst, is not yet well-understood. In this work, we present a computational mechanistic study of the thermal vs photocatalytic oxidation of FA on dry and hydrated anatase TiO 2 (101) surfaces, based on periodic hybrid density functional theory (DFT) calculations, in which the photooxidation is treated as an excited-state process in a constrained triplet spin state. We first compare the adsorption modes of FA on the anatase (101) surface in the ground and excited states, followed by identification of the corresponding reaction intermediates that lead to the formation of CO 2 . We unveil the pivotal role of photogenerated holes localized at surface under-coordinated oxygen sites in mediating the C−H bond cleavage, thereby promoting CO 2 formation through a highly stable intermediate and an exergonic reaction step. Further investigation of the effect of coadsorbed water molecules shows that hydrogen bonding with water stabilizes FA in a monodentate configuration. This is favored over the unreactive bidentate structure that is the most stable under dry conditions, thus providing insight into the experimentally observed increase of the reaction rate in the presence of water.

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

Optimal Design of Sustainable Ammonia-Based Food–Energy–Water Systems with Nitrogen Management

As the basis for virtually any form of nitrogen fertilizers, ammonia plays a vital role in agriculture; in addition, there has been an increased interest in its use as a carbon-free energy carrier. However, ammonia is also associated with two major environmental concerns: CO 2 emissions from the conventional production process and nitrogen pollution from the excessive use of ammonia-based fertilizers. To mitigate these environmental impacts, we develop an optimization framework for the design of a sustainable ammonia-based agricultural system that synergistically integrates the production of ammonia from renewable resources and effective measures for nitrogen management. The proposed model captures the effect of intermittency by incorporating both design and detailed operational decisions. Here, by applying a multiscale time representation that reduces the problem size and a tailored surrogate model that accurately approximates model nonlinearity, we are able to achieve optimal solutions within reasonable computation times. A computational case study is conducted using real-world data from a local farm in Morris, Minnesota, and the results indicate the trade-off between cost and nitrogen loss. Importantly, we show that practicing effective nitrogen management can significantly reduce the nitrogen loss with only a small increase in net present cost.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Two-State Reactivity in Iron-Catalyzed Alkene Isomerization Confers σ-Base Resistance

A low-coordinate, high spin (S = 3/2) organometallic iron(I) complex is a catalyst for the isomerization of alkenes. A combination of experimental and computational mechanistic studies supports a mechanism in which alkene isomerization occurs by the allyl mechanism. Importantly, while substrate binding occurs on the S = 3/2 surface, oxidative addition to an η 1 -allyl intermediate only occurs on the S = 1/2 surface. In conclusion, since this spin state change is only possible when the alkene substrate is bound, the catalyst has high immunity to typical σ-base poisons due to the antibonding interactions of the high spin state.

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

Alkane Dehydrogenation and H/D Exchange by a Cationic Pincer-Ir(III) Hydride: Cooperative C–H Addition and β -H Elimination Modes Induce Anomalous Selectivity

We report that the cationic iridium complex ( iPr PCP)IrH + catalyzes the transfer-dehydrogenation of alkanes to give alkenes and hydrogen isotope exchange (HIE) of alkanes and arenes. Contrary to established selectivity trends found for C–H activation by transition metal complexes, strained cycloalkanes, including cyclopentane, cycloheptane, and cyclooctane, undergo C–H addition much more readily than n-alkanes, which in turn are much more reactive than cyclohexane. Aromatic C–H bonds also undergo H/D exchange much less rapidly than those of the strained cycloalkanes, but much more favorably than cyclohexane. The order of reactivity toward dehydrogenation correlates qualitatively with the reaction thermodynamics, but the magnitude is much greater than can be explained by thermodynamics. Accordingly, the cycloalkenes corresponding to the strained cycloalkanes undergo hydrogenation much more readily than cyclohexene, despite the less favorable thermodynamics of such hydrogenations. Here, computational (DFT) studies allow rationalization of the origin of reactivity and the unusual selectivity. Specifically, the initial C–H addition is strongly assisted by β-agostic interactions, which are particularly favorable for the strained cycloalkanes. Subsequent to α-C–H addition, the H atom of the β-agostic C–H bond is transferred directly to the hydride ligand of ( iPr PCP)IrH + to give a dihydrogen ligand.

02 PETROLEUM↗
DOE OSTI2025

P–P Coupling with and without Terminal Metal–Phosphorus Intermediates

Terminal metal–phosphorus (M–P) complexes are of significant contemporary interest as potential platforms for P-atom transfer (PAT) chemistry. Decarbonylation of metal–phosphaethynolate (M–PCO) complexes has emerged as a general synthetic approach to terminal M–P complexes. M–P complexes that are stabilized by strong M–P multiple bonds are kinetically persistent and isolable. In the absence of strong M–P stabilization, the formation of diphosphorus-bridged complexes (i.e., M–P–P–M species) is often interpreted as evidence for the intermediacy of reactive, unobserved M–P species. Here, we demonstrate that while diphosphorus complexes can arise from reactive M–P species, P–P coupling can also proceed directly from M–PCO species without the intermediacy of M–P complexes. Photochemical decarbonylations of a pincer-supported Ni (II)–PCO complex at 77 K afford a spectroscopically observed terminal Ni–P complex, which is best described as a triplet, Ni(II)-metallophosphinidene with two unpaired electrons localized on the atomic phosphorus ligand. Thermal annealing of this transient Ni–P complex results in rapid dimerization to afford the corresponding P 2 2– -bridged dinickel complex. Unexpectedly, the same P 2 2– -bridged dinickel complex can also be accessed via a thermally promoted process in the absence of light. The analysis of reaction kinetics, isotope-labeling studies, and computational results indicate that the thermal P–P coupling process proceeds via a noncanonical mechanism that avoids terminal M–P intermediates. Together, these results represent the first observation of P–P coupling from characterized terminal M–P species and demonstrate that terminal M–P intermediates are not required to obtain P–P coupling products. These observations provide critical mechanistic understanding of the activation modes relevant to P-atom transfer.

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

Development, Application, and Mechanistic Interrogation of a Dual Ni Catalysis Approach to Photoredox-Based C(sp 3 )–C(sp 3 ) Cross-Coupling

The installation of alkyl substituents such as methyl groups is a crucial tactic for the synthesis and diversification of medicinal compounds with improved pharmacological profiles. However, strategies that leverage methyl radical for C(sp 3 ) methylation remain underdeveloped due to the challenge of obtaining cross-selectivity between fleeting aliphatic radicals and alkyl electrophiles. Here, we report the development, application, and interrogation of a conceptually novel mechanistic framework for C(sp 3 )–Me bond formation using two distinct Ni catalysts capable of cross-coupling sterically and electronically diverse alkyl halides (chlorides and bromides) with methyl radical generated photocatalytically from benzaldehyde dimethyl acetal. Furthermore, by modifying the alkyl substituents on the acetal coupling partner, we demonstrate cross-couplings beyond methylation to access an array of 1°–1° and 1°–2° alkyl–alkyl bonds. Experimental and computational mechanistic studies provide support for cooperativity between an in situ-generated (bpy)Ni I (X) catalyst that facilitates XAT and inner-sphere C–C bond formation and a (Tp*)Ni II (acac) cocatalyst that captures methyl radical and engages in concurrent outer-sphere S H 2 coupling.

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

Photoinduced chemomimetic biocatalysis for enantioselective intermolecular radical conjugate addition

Exploiting nature’s catalysts for non-natural transformations that are inaccessible to chemocatalysis is highly desirable but challenging. On the one hand, the widespread nicotinamide-dependent oxidoreductases have not been utilized for single-electron-transfer-induced bimolecular cross-couplings; on the other, the addition of catalytic asymmetric radical con- jugate to terminal alkenes remains a challenge owing to strong racemic background reaction and unselective termination of prochiral radical species. Here we report a chemomimetic biocatalysitic approach for construction of alpha-carbonyl stereocentres via an unnatural intermolecular conjugate addition of N-(acyloxy)phthalimides-derived radicals with acceptor-substituted terminal alkenes, by combination of visible-light excitation and nicotinamide-dependent ketoreductases (KREDs). Based on protein crystal structure, we engineered KREDs via a semi-rational mutagenesis strategy to improve reaction outcomes with a small and high-quality variants library. Here, mechanistic investigations combining wet experiments, crystallographic studies and computational simulations demonstrate that the repurposed biocatalyst can suppress racemic background reaction and unselected side reactions, yielding enantioselectivity that is challenging to achieve by chemocatalysis.

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

Impact of processing humidity on ionomer film structure and performance in hydroxide exchange membrane electrolyzers

Hydroxide exchange membrane electrolyzers (HEMELs) enable hydrogen production using low-cost, earth-abundant materials. Improving electrode fabrication is integral to enhancing device performance, and ionomer-responsible for transporting hydroxide and mechanically supporting the catalyst-is a major component. Here, we use experiments and computation to study the effects of relative humidity (RH) during the drying process of poly(aryl piperidinium) ionomer films on HEMEL electrodes. Broadly, the drying environments determine the physical structure and electrochemical traits of the ionomer network. High RH drying yields a highly porous network with excessive water uptake, structural defects, washout, and 64% reduction in hydroxide conductivity. Extremely low RH drying produces an overly compact pore network that hinders hydroxide mobility. In contrast, moderately low RH drying (9% RH) creates an ionomer film with well-balanced traits: excellent mechanical stability and connectivity needed for catalyst retention and hydroxide transport, which improves HEMEL performance by 40% at 1.8 V compared to suboptimal RHs. This research advances HEMEL manufacturing by providing a simple, scalable, and low-cost approach to optimize electrode ionomer films.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pinning dislocations in colloidal crystals with active particles that seek stacking faults

There is growing interest in functional, adaptive devices built from colloidal subunits of micron size or smaller. A colloidal material with dynamic mechanical properties could facilitate such microrobotic machines. Here we study via computer simulation how active interstitial particles in small quantities can be used to modify the bulk mechanical properties of a colloidal crystal. Passive interstitial particles are known to pin dislocations in metals, thereby increasing resistance to plastic deformation. We extend this tactic by employing anisotropic active interstitials that travel super-diffusively and bind strongly to stacking faults associated with partial dislocations. We find that: (1) interstitials that are effective at reducing plasticity compromise between strong binding to stacking faults and high mobility in the crystal bulk. (2) Reorientation of active interstitials in the crystal depends upon rotational transitions between high-symmetry crystal directions. (3) The addition of certain active interstitial shapes at concentrations as low as 60 per million host particles (0.006%) can create a shear threshold for dislocation migration. This work demonstrates how active materials in a dense matrix can locally sense their environment and lead to bulk property changes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulation of the photodetachment spectra of the nitrate anion (NO 3 ₋ ) in the B ~ 2 E' energy range and non-adiabatic electronic population dynamics of NO 3

The photodetachment spectrum of the nitrate anion (NO 3 ₋ ) in the energy range of the NO 3 second excited state is simulated from first principles using quantum wave packet dynamics. The prediction at 10 K and 435 K relies on the use of an accurate full-dimensional fully coupled five state diabatic potential model utilizing an artificial neural network. The ability of this model to reproduce experimental spectra was demonstrated recently for the lower energy range. Analysis of the spectra indicates a weaker Jahn–Teller coupling compared to the first excited state. The detailed non-adiabatic dynamics is studied by computing the population dynamics. An ultra-fast non-statistical radiationless decay is found only among the Jahn–Teller components, which is followed by a slow statistical non-radiative decay among the different state manifolds. The latter is reproduced perfectly by a simple first order kinetics model. The dynamics in the second excited state is not affected by the presence of a conical intersection with the first excited state manifold.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High throughput exploration of the oxidation landscape in high entropy alloys

We report high entropy alloys (HEAs) have gained interest for structural applications in extreme environments. With a potentially vast chemical and phase space, there are significant opportunities to discover superior performing alloys. Crucial for most high-temperature applications is understanding and mitigating the oxidation behavior of these chemically complex alloys. Most experimental and computational HEA studies have focused on a limited set of compositions and only a fraction these compositions have been characterized for oxidation. We present a high-throughput framework that utilizes density-functional theory (DFT) in concert with a combined machine-learning model and grand-canonical linear programming for assessing phase stability, phase-fraction, chemical activity and high-temperature survivability of arbitrary HEAs. This framework considers temperature dependent contributions to the Gibbs energy of the competing phases arising from short-range order and vibrational entropy. We demonstrate the effectiveness of the framework by assessing the thermodynamic stability, oxidation behavior, chemical activity, and phase decomposition of body-centered cubic Mo-W-Ta-Ti-Zr refractory HEAs. A total of 51 compositions were analyzed and ranked in order of their survivability based on the Pareto-front analysis. Oxidation was performed at 1373 K on fours samples in air showing the difference in oxidation behavior determined experimentally through scale thickness and their mass changes. The insights on oxidation behavior presented in this work will enable the fast assessment of technologically useful HEAs needed for future structural application in extreme conditions

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental equation of state of 11 lanthanide nitrides (NdN to LuN) and pressure induced phase transitions in NdN, SmN, EuN, and GdN

Through an extensive data analysis of powder X-ray diffraction data obtained at pressures up to at least 78 GPa, we report the experimental equations of state for all lanthanide nitrides between NdN and LuN, excluding the radioactive Pm. By fitting the obtained unit cell volumes as a function of pressure with the third order Birch–Murnaghan equation of state, we find that the bulk modulus increases with an increasing lanthanide number from K 0 = 146(12) GPa for NdN to 182(7) GPa in EuN. Hereafter, the bulk modulus reaches a plateau for the rest of the series except for TmN which has a lower bulk modulus. We find that the first derivative of the bulk modulus is around 4 for all compounds except TbN, which displays a significantly different compression behavior. In addition, we find a B1 to B10 pressure-induced phase transition in NdN, SmN, EuN, and GdN at increasingly higher pressures. In fact, we observe that the onset pressure of the phase transition increases linearly with Ln atomic number. From TbN and onwards, we do not observe any sign of a B1 to B10 transition indicating that the transition pressure exceeds the studied pressure range. Therefore, we believe that, for the heavier lanthanides, the linear relationship between the onset pressure and the lanthanide number does not hold and even higher pressures are needed to observe the transition. This coherent study of the series of lanthanide nitrides offers a unique opportunity for benchmark studies of computational methods applied to compounds with 4f electrons.

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

Modeling of carbon and tungsten transient dust influx in tokamak edge plasma

The paper presents computer simulation studies of burst injection of carbon and tungsten dust particles in DIII-D-like edge plasmas. The injection causes a large transient influx of the low- and high-Z impurities associated with the dust ablation in the plasmas. The dust transport and the effects of the ablated impurities on the edge plasma dynamics in a modern mid-size tokamak geometry are investigated for low- and high-power plasma discharge conditions. The core plasma contamination with dust-ablated impurities and the factors affecting it are evaluated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2026

Prospects for detecting UF6 hydrolysis intermediates by mass spectrometry, resonance Raman, and NQR spectroscopy

Simulations are performed to consider several spectrometric and spectroscopic candidates for elucidating the mechanism of the hydrolysis of uranium hexafluoride (UF6). This study is among the first to benchmark the def-mTZVP basis sets for actinide-containing molecules, and it is shown to be a suitable basis set for surveying geometrical structures and vibrational spectra when used in conjunction with density functional theory. An experiment is proposed coupling mass spectrometric ion selection with vibrational spectroscopy, and supporting infrared and Raman spectral simulations demonstrate that ionization blue-shifts bands and can change their qualitative features. Ultraviolet resonance Raman is shown to have good prospects for discriminating U–O–U bridged intermediates, evidence for which was observed recently [L. E. McNamara et al., J. Phys. Chem. A 130, 775–786 (2026)]. As a more speculative approach, we also consider the prospect of addressing 233U or 235U nuclei by quadrupole resonance (NQR) spectroscopy for assigning early stage intermediate complexes. In doing so, we provide a first order-of-magnitude estimate for the collision-induced NQR signal for the UF6 dimer, which is of interest for the interpretation of a fast decoherence time observed in liquid-phase nuclear magnetic resonance. Having provided critical insights into the formation and stability of UF6 hydrolysis intermediates in previous studies, this computational spectroscopy survey is expected to help guide and expedite future laboratory campaigns.

Lutz, Jesse J. [Center for Computing Research, San↗
DOE OSTI2020

Crystallographic evidence of Watson–Crick connectivity in the base pair of anionic adenine with thymine

Utilizing an ionic liquid strategy, we report crystal structures of salts of free anionic nucleobases and base pairs previously studied only computationally and in the gas phase. Reaction of tetrabutylammonium ([N 4444 ] + ) or tetrabutylphosphonium ([P 4444 ] + ) hydroxide with adenine (HAd) and thymine (HThy) led to hydrated salts of deprotonated adenine, [N 4444 ][Ad]·2H 2 O, and thymine, [P 4444 ][Thy]·2H 2 O, as well as the double salt cocrystal, [P 4444 ] 2 [Ad][Thy]·3H 2 O·2HThy. The cocrystal includes the anionic [Ad-(HThy)] base pair which is a stable formation in the solid state that has previously not even been suggested. It exhibits Watson–Crick connectivity as found in DNA but which is unusual for the free neutral base pairs. The stability of the observed anionic bases and their supramolecular formations and hydrates has also been examined by electronic structure calculations, contributing to more insight into how base pairs can bind when a proton is removed and highlighting mechanisms of stabilization or chemical transformation in the DNA chains.

59 BASIC BIOLOGICAL SCIENCES↗