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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

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

Aqueous Titania Interfaces

Water–metal oxide interfaces are central to many phenomena and applications, ranging from material corrosion and dissolution to photoelectrochemistry and bioengineering. In particular, the discovery of photocatalytic water splitting on TiO 2 has motivated intensive studies of water-TiO 2 interfaces for decades. So far, a broad understanding of the interaction of water vapor with several TiO 2 surfaces has been obtained. However, much less is known about liquid water–TiO 2 interfaces, which are more relevant to many practical applications. Probing these complex systems at the molecular level is experimentally challenging and is sometimes possible only through computational studies. This review summarizes recent advances in the atomistic understanding, mostly through computational simulations, of the structure and dynamics of interfacial water on TiO 2 surfaces. The main focus is on the nature, molecular or dissociated, of water in direct contact with low-index defect-free crystalline surfaces. The hydroxyls resulting from water dissociation are essential in the photooxidation of water and critically affect the surface chemistry of TiO 2.

density functional theory↗
DOE OSTI2020

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: Spheres and cylinders

The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Radicals as Exceptional Electron-Withdrawing Groups: Nucleophilic Aromatic Substitution of Halophenols Via Homolysis-Enabled Electronic Activation

While heteroatom-centered radicals are understood to be highly electrophilic, their ability to serve as transient electron-withdrawing groups and facilitate polar reactions at distal sites has not been extensively developed. Here, we report a new strategy for the electronic activation of halophenols, wherein generation of a phenoxyl radical via formal homolysis of the aryl O-H bond enables direct nucleophilic aromatic substitution of the halide with carboxylate nucleophiles under mild conditions. Pulse radiolysis and transient absorption studies reveal that the neutral oxygen radical (O • ) is indeed an extraor-dinarily strong electron-withdrawing group [σ p - (O • ) = 2.79 vs. σ p - (NO 2 ) = 1.27]. Additional mechanistic and computational studies indicate that the key phenoxyl intermediate serves as an open-shell electron-withdrawing group in these reactions, lowering the barrier for nucleophilic substitution by more than 20 kcal/mol relative to the closed-shell phenol form of the substrate. Finally, by using radicals as transient activating groups, this homolysis-enabled electronic activation strategy provides a powerful platform to expand the scope of nucleophile-electrophile couplings and enable previously challenging transformations.

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

How much is surface dopant enough to maximize CO2-to-liquid chemicals conversion at industrially relevant current density?

This invited talk will be presented at the symposium "Waste Feedstock to Fuels and Petrochemicals”, Fuels and Petrochemical Division at the 2024 AIChE Annual Meeting. In this study, we will discuss how minute amount of surface heteroatoms would maximize CO2 reduction to formate/formic acid at industrially relevant current densities and high selectivity which was validated by both experimental and computational studies. The outstanding performance of the best-in-class catalysts in both H-cell and full-cell electrolyzer cell is also demonstrated. Our findings would provide additional design concepts of high performance CO2R electrocatalysts.

Nguyen Phan, Thuy Duong↗
DOE OSTI2022

Intracardiac Electrical Imaging using the 12-lead ECG: A Machine Learning Approach using Synthetic Data

Current state-of-the-art techniques for non-invasive imaging of cardiac electrical phenomena require voltage recordings from dozens of different torso locations and anatomical models built from expensive medical diagnostic imaging procedures. Here this study aimed to assess if recent machine learning advances could alternatively reconstruct electroanatomical maps at clinically relevant resolutions using only the standard 12-lead electrocardiogram (ECG) as input. To that end, a computational study was conducted to generate a dataset of over 16000 detailed cardiac simulations, which was then used to train neural network (NN) architectures designed to exploit both spatial and temporal correlations in the ECG signal. Analysis over a validation set showed average errors in activation map reconstruction below 1.7 msec over 75 intracardiac locations. Furthermore, phenotypical patterns of activation and the morphology of the activation potential were correctly reconstructed. The approach offers opportunities to stratify patients non-invasively, both retrospectively and prospectively, using metrics otherwise only available through invasive clinical procedures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Coupling SCALE with DAKOTA for Axial Burnup Profiles Assessment in Burnup Credit

This paper presents a computational study that demonstrates the application of the SCALE code system in conjunction with the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) for the analysis of key factors influencing the evaluation of burnup credit (BUC) in pressurized water reactors (PWRs). The primary objective of this analysis is to characterize the model by utilizing parameterization, uncertainty quantification, and optimization studies. Using this approach, we can comprehensively assess the system and conduct informed predictive studies. This study highlights the effectiveness of the SCALE code system integrated within the DAKOTA framework in terms of efficiency and capability. With the coupling of the burnup code ORIGAMI with the CSAS or TSUNAMI-3D sequence embedded in a DAKOTA analysis, we can characterize the factors that influence the k eff of PWR 17x17 spent nuclear fuel (SNF) in the GBC-32 computational benchmark cask for the assessment of BUC in criticality safety analysis. The coupling methodology used in this study is not exclusive to BUC analysis. However, the choice to apply this methodology to the BUC problem is particularly significant because of the diverse range of aspects it encompasses in nuclear criticality safety analyses. This problem presents a unique opportunity to explore and address multiple facets of such analyses related to BUC and illustrates the capability of the SCALE code system with DAKOTA. This analysis makes use of historical reference data for the axial burnup profile, where the entire space within the bounds is considered. Both SCALE and DAKOTA are currently integrated in the Nuclear Energy Advanced Modeling Simulation (NEAMS) Workbench code system, which has a user-friendly graphical interface that simplifies the setup of simulations and configuration of input parameters as well as the visualization of simulation results.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

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

Reductive Elimination From Tetra‐Alkyl Cuprates [Me n Cu(CF 3 ) 4− n ] − ( n =0–4): Beyond Simple Oxidation States

Abstract In recent years, the electronic structures of organocuprates in general and the complex [Cu(CF 3 ) 4 ] − in particular have attracted significant interest. A possible key indicator in this context is the reactivity of these species. Nonetheless, this aspect has received only limited attention. Here, we systematically study the series of tetra‐alkyl cuprates [Me n Cu(CF 3 ) 4− n ] − and their unimolecular reactivity in the gas phase, which includes concerted formal reductive eliminations as well as radical losses. Through computational studies, we characterize the electronic structures of the complexes and show how these are connected to their reactivity. We find that all [Me n Cu(CF 3 ) 4− n ] − ions feature inverted ligand fields and that the distinct reactivity patterns of the individual complexes arise from the interplay of different effects.

Zimmer, Bastian↗
DOE OSTI2025

Electrocatalytic Ammonia Oxidation by Pyridyl-Substituted Ferrocenes

Ammonia (NH 3 ) is a promising carbon-free fuel when prepared from sustainable resources. First-row transition metal electrocatalysts for ammonia oxidation are an enabling technology for sustainable energy production. We describe electrocatalytic ammonia oxidation using robust molecular complexes based on Earth-abundant iron. Electrochemical studies of ferrocenes with covalently attached pyridine arms reveal facile ammonia oxidation in DMSO (2.4 M NH 3 ) with modest overpotentials (η = 770–820 mV) and turnover frequencies (125–560 h –1 ). Experimental and computational studies indicate that the pendant pyridyl base serves as an H-bond acceptor with an N–H bond of ammonia that transfers a proton to the pyridine following oxidation by the attached ferrocenium moiety in a proton-coupled electron transfer (PCET) step. This generates an amidyl (•NH 2 ) radical stabilized via H-bonding to a pendant pyridinium moiety that rapidly dimerizes to hydrazine (H 2 N–NH 2 ), which is easily oxidized to nitrogen (N 2 ) at the glassy carbon working electrode. This report identifies a general strategy to oxidize ammonia via H-bonding to a base (B:), thereby activating [B···H-NH 2 ] toward PCET by a proximal oxidant to form [BH···NH 2 ] +/• radical cations, which are susceptible to dimerization to form easily oxidized hydrazine.

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

The electron affinity of the uranium atom

The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U – . Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin–orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th – , the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both scenarios, the ground state of the anion corresponds to electron attachment to the 6d orbital.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Pressure induced modification of electronic and magnetic properties of MnCrNbAl and MnCrTaAl

Spin-gapless semiconductor (SGS) is a new class of material that has been studied recently for potential applications in spintronics. This material behaves as an insulator for one spin channel, and as a gapless semiconductor for the opposite spin. In this work, we present results of a computational study of two quaternary Heusler alloys, MnCrNbAl and MnCrTaAl that have been recently reported to exhibit spin-gapless semiconducting electronic structure. In particular, using density functional calculations we analyze the effect of external pressure on electronic and magnetic properties of these compounds. It is shown that while these two alloys exhibit nearly SGS behavior at optimal lattice constants and at negative pressure (expansion), they are half-metals at equilibrium, and magnetic semiconductors at larger lattice constant. At the same time, reduction of the unit cell volume has a detrimental effect on electronic properties of these materials, by modifying the exchange splitting of their electronic structure and ultimately destroying their half-metallic/semiconducting behavior. Thus, our results indicate that both MnCrNbAl and MnCrTaAl may be attractive for practical device applications in spin-based electronics, but a potential compression of the unit cell volume (e.g. in thin-film applications) should be avoided.

Materials Science↗
DOE OSTI2024

Structural insights into the unexpected agonism of tetracyclic antidepressants at serotonin receptors 5-HT 1e R and 5-HT 1F R

Serotonin [5-hydroxytryptamine (5-HT)] acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT 1e R have confirmed roles in native tissue and are validated drug targets. Despite 5-HT 1e R’s therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT 1e R’s pharmacology in relation to the highly homologous 5-HT 1F R, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT 1e R/5-HT 1F R agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed antimigraine properties. Using cryo-EM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT 1e R distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT 1e R and 5-HT 1F R contribute to the agonist activity of these antidepressants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A new approach for simulating inhomogeneous chemical kinetics

Abstract In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling $$\alpha$$ α - and $$\gamma$$ γ -radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method’s efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach’s capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method’s suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Study of self-assembly behavior and ionic conductivity of conjugated liquid crystals with T-shaped facial-polyphilic structure

The unique self-assembly of liquid crystals (LCs), combined with their potential application as organic semiconductors, has become a focus of recent research. Here, a joint experimental and computational study of the self-assembly and ionic conduction was carried out on a series of T-shaped conjugated LCs consisting of three incompatible components. By extending the EOn side-chain length, several experimental evaluations confirmed a decrease of the order-disorder transition temperature, while coarse-grained simulations revealed a structural evolution from a smectic phase to a columnar phase. Ionic conductivity of these molecules was achieved by adding Li salt, leading to a maximum conductivity of 1.1 × 10 −3 siemens per centimeter observed at 120°C. All-atom simulations were performed to examine the Li-ion solvation environment and to evaluate the intrachain and interchain Li-ion hopping mechanisms. The molecule with a long EOn side chain was found to generate a densely distributed network of Li-ion solvation sites, which can facilitate effective interchain hopping to promote ion transport.

Liu, Ziwei [Cornell Univ., Ithaca, NY (United Stat↗
DOE OSTI2026

Concerted Proton–Electron Transfer Minimizes Substituent Effects on Adsorbed Phthalocyanine Electrocatalysis

Molecularly modified electrodes (MMEs) are potent electrocatalysts, but few principles exist for their rational design. Electrocatalysis by soluble molecules depends strongly on substituents that tune the catalyst redox potential (E 1/2 ), but it is unclear if this parameter similarly impacts MME catalysis. Herein, we employ the hydrogen evolution reaction (HER) as a test case for comparing carbon-adsorbed cobalt phthalocyanine (CoPc/C) and cobalt hexadecafluoro-phthalocyanine (CoFPc/C). By correlating HER activity and voltammetric data to total Co surface concentration across a wide range of catalyst loadings, we find that only 5–25% of adsorbed Co sites contribute to the Co(II/I) redox wave and that this subpopulation poorly correlates with catalytic activity. Instead, in the low-loading limit, catalytic activity correlates linearly with the majority Co(II/I)-silent Co population, revealing per-site turnover frequency (TOF) values for HER. Despite a 230 mV difference in Co(II/I) redox potentials, CoPc/C and CoFPc/C display TOF values differing by less than a factor of 3 when compared over a wide potential range. Mechanistic studies point to an inner-sphere concerted proton–electron transfer step as rate-determining, suggesting that the Co–H bond dissociation free energy (BDFE) rather than the Co(II/I) E 1/2 is thermodynamically relevant. Computational studies indicate that the fluoro-substituents lead to compensatory changes in Co(II/I) E 1/2 and Co(I) basicity, leaving the Co–H BDFE largely unchanged between CoPc and CoFPc and thereby manifesting in similar catalytic rates. Furthermore, these results highlight the limited effect of E 1/2 -tuning on MME catalytic activity and motivate the development of methods to directly alter active site–substrate BDFE.

Catalysts↗
DOE OSTI2026

Kinetic and modeling studies of the mechanism of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2

Since its discovery over 15 years ago, the reversible dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 has remained one of the more intriguing hydrogen-cycling systems. While the mechanism of this reaction has been the subject of a good deal of speculation and computational studies, prior to this work it had not been probed through kinetic studies. Previous reports of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 have not included kinetic studies. The present studies have shown that the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 is suppressed by hydrogen pressure indicating that the rate-limiting step in this process involves hydrogen elimination. Computational modeling of kinetic data obtained from monitoring the hydrogen elimination from Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 under static vacuum over a range of temperatures supports that the dehydrogenation occurs through a reversible three-step process in which the elimination of hydrogen from the [B 3 H 10 ] − intermediate is rate limiting. A mechanism involving the low energy transfer of neighboring BH3 groups is proposed to account for the formation of [B 3 H 8 ] − at relatively low temperatures.

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

New High-Entropy Perovskite Oxides with Increased Reducibility and Stability for Thermochemical Hydrogen Generation

This project aims to design, synthesize, and test a transformative class of High-Entropy Perovskite Oxides (HEPOs) as redox oxides to enable thermochemical hydrogen generation with improved stability, kinetics, and efficiency. These developed HEPOs are expected to demonstrate improved kinetics with oxygen surface exchange coefficient ( k > 7.5×10 -4 cm/s) in Budget Period (BP) 1, retain its structural stability in a broad range of oxygen non-stoichiometry (Δδ > 0.15) at a low operating reduction temperature of T red < 1400°C in BP 2, and deliver a H 2 yield of over 400 µmol per gram of oxide and high stability with less than 20% degradation after at least 50 cycles in BP 3. This project is feasible due to the unique thermodynamic properties (simultaneously increased reducibility and phase stability) and kinetic characters (stability against particle coarsening and potentially enhanced oxygen transport and surface reaction kinetics) of such HEPOs, and it is enabled by a unique active learning computational design approach. Computational studies have been conducted to investigate the oxygen vacancy formation in complex perovskite systems. * Accurate prediction of V O .. concentration with disordered A-site cations in Fe-based high-entropy perovskite oxides * Combined MC/DFT computation elucidates the mechanism of Co preference on the redox due to the strain introduced by local distortion. In this project, we explored a large number (~150) of perovskite compositions, which are listed in Tables 2 – 4). * All perovskite specimens have been synthesized through a high-throughput high-energy ball milling process, followed by the conventional sintering process. * XRD, SEM/EDS and TGA were performed to confirm the crystal structure, phase homogeneity and oxygen non-stoichiometry for compositionally complex perovskite oxides (CCPOs). * 110 compositions show single-phase from XRD * Unusual aliovalent doping effects in medium-entropy perovskite compositions. * V-shape relation between Δδ vs. x (La 1-x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3 * The highest reported hydrogen production for the CCPOs made in this project ( T re = 1350 ºC 30 min, T Ox = 1100 ºC 30 min) * B-site mixing (La 0.8 Sr 0.2 )(Mn 0.2 Fe 0.2 Co 0.4 Al 0.2 )O 3 : 89.97 ± 2.73 mmol H2 /mol oxide (395 ± 10 μmol/g oxide ) (i) No phase transformation during reactions when Co molar ratio is less than 61% (ii) Balance between intrinsic kinetics (oxygen surface exchange) and thermodynamics (oxygen non-stoichiometry) (iii) Preference of Co identified by in-situ XPS * A-site mixing (La 1/6 Pr 1/6 Nd 1/6 Gd 1/6 Ba 1/6 Sr1/6)MnO 3 : 98.48 mmol H2 /mol oxid e (~415 μmol/g oxide )

08 HYDROGEN↗