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

Homoleptic An(IV) Dimethyl Sulfoxide Complexes of Th–Pu Enable Insight into Solution Speciation and Redox

A family of 8- and 9- coordinate homoleptic dimethyl sulfoxide An(IV) complexes is used to evaluate the solvation properties of tetravalent actinide ions in dimethyl sulfoxide (DMSO) vs aqueous media. These compounds are prepared from aqueous mixtures of DMSO and yield crystalline solids, whose solid-state electronic absorption spectra are compared to solution phase spectra, providing insight into the solution speciation. The thermodynamics of the equilibria between the 8- and 9-coordinate compounds were determined experimentally and computationally, showing a trend that correlates with ionic radius. Cyclic voltammetry data for the +IV/+III couples of Np and Pu in the DMSO electrolyte indicate that the +IV ions are significantly stabilized compared to aqueous media due to the more donating nature of DMSO compared to water. As a result, computational studies of these systems indicate that further reactivity may take place upon reduction to the +III oxidation state.

Actinides↗
DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

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

Polarized, V-Shaped, and Conjoined Biscoumarins: From Lack of Dipole Moment Alignment to High Brightness

Eleven conjoined coumarins possessing a chromeno[3,4-c]chromene-6,7-dione skeleton have been synthesized via the reaction of electron-rich phenols with esters of coumarin-3-carboxylic acids, catalyzed by either Lewis acids or 4-dimethylaminopyridine. Furthermore, Michael-type addition to angular benzo[f]coumarins is possible, leading to conjugated helical systems. Arrangement of the electron-donating amino groups at diverse positions on this heterocyclic skeleton makes it possible to obtain π-expanded coumarins with emission either sensitive to, or entirely independent of, solvent polarity with large Stokes shifts. Computational studies have provided a rationale for moderate solvatochromic effects unveiling the lack of collinearity of the dipole moments in the ground and excited states. Depending on the functional groups present, the obtained dyes are highly polarized with dipole moments of ~14 D in the ground state and ~20–25 D in the excited state. Strong emission in nonpolar solvents, in spite of the inclusion of a NO 2 group, is rationalized by the fact that the intramolecular charge transfer introduced into these molecules is strong enough to suppress intersystem crossing yet weak enough to prevent the formation of dark twisted intramolecular charge transfer states. Photochemical transformation of the dye possessing a chromeno[3,4-c]pyridine-4,5-dione scaffold led to the formation of a spirocyclic benzo[g]coumarin.

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

Mapping the Electronic Structure of the Uranium(VI) Dinitride Molecule, UN 2

A combined anion photoelectron spectroscopic and relativistic electronic-structure computational study of the electronic structure of the UN 2 molecule is reported. Since the photoelectron spectrum of uranium dinitride’s negative ion, UN 2 - , (UN 2 plus an additional electron) directly reflects the electronic structure of neutral UN 2 , we have measured and relied upon the photoelectron spectrum of the UN 2 – anion as a means of mapping the electronic structure of neutral UN 2 . In addition to the electron affinity of the UN 2 ground state (the energy released when attaching an electron to UN 2 molecule), energy levels of the UN 2 excited states (electronic states with higher energies than the ground state) were well characterized by the close interplay between experiment and high-level theory. We found that both electron attachment and electronic excitation significantly bend the UN 2 molecule and elongate its U≡N bond. Implications for the activation of UN 2 are discussed.

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

Charge Tethering Drives Intermediate-Range Order and Slow Dynamics in Zwitterionic Liquids

We have synthesized and characterized a family of zwitterionic liquids (ZwLs) based on the poly(ethylene oxide) imidazolium cation and the alkyl sulfonate anion. To better rationalize ZwLs′ structural behavior and networkforming properties, we have computationally studied them in contrast to isostructural but chemically untethered ionic liquid (IL) analogues. Perhaps surprisingly, it is the molecular liquids and not the ILs that show the most selfassembly and intermediate-range order associated with scattering prepeaks; the ZwLs display also the slowest dynamics. This is consistent with prior work showing the importance of networks in considering the viscoelastic relaxation of liquids. We find that, because of their significantly large molecular dipoles, these zwitterions are versatile network-building blocks that result in highly viscous liquids. From a structural perspective, networks in these ZwLs sit somewhere in between what would be true chemical connectivity and the Coulombic networks observed in ionic liquids.

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

First-Principles Grand-Canonical Simulations of Water Adsorption in Proton-Exchanged Zeolites

Water appears by design or as impurities in many important reactive systems. For those catalyzed by porous solid acids, such as widely used zeolites, experimentally quantifying the amount or elucidating the structure of adsorbed water clusters at reaction conditions is challenging, while computational studies (e.g., first-principles molecular dynamics simulations) often need to assume the loading to examine solvation effects. Furthermore, we perform first-principles grand-canonical simulations to predict water adsorption to H-ZSM-5 zeolites under specified experimental conditions. Presampling with inexpensive force fields and a pool-based parallelization algorithm are used to improve simulation efficiency, while molecular dynamics is used to sample configurations involving hydronium species. We observe that H + exchange dramatically increases the hydrophilicity of zeolite MFI and an appreciable amount of water is present at very low relative humidities. At all conditions examined, the zeolitic protons are found to dissociate readily from the surface basic sites and become mobile by participating in the hydrogen-bonded chains of adsorbed water molecules.

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

Mechanism of Stoichiometrically Governed Titanium Oxide Brownian Tree Formation on Stepped Au(111)

Previously observed formation of substoichiometric titanium oxide dendritic structures across terraces of Au(111) is computationally studied and shown to follow the classical fractal formation mechanism of diffusion-limited aggregation (DLA). Globally optimized gas-phase oxide cluster structures are sampled in a variety of landing formations on gold surfaces and shown to favor isomers driving polymerization to Brownian tree fractal structures. Mobility of Ti 3 O 5 monomers is shown to be extremely high, with diffusion barriers of 0.21 eV or less. Through bonding stabilization, polymerization of these monomers is energetically favorable and irreversible on the 111 terrace but geometrically impossible to propagate along the step edge. Simulated scanning tunneling microscopy (STM) images show strong similarity to experiment. By contrast, observation of Ti 3 O 6 aggregating as wires along step edges is explained by the affinity of oxygen to step edges and statistical arguments for aggregation entropy at the step, in addition to low barriers for monomer diffusion and polymerization.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electronic Structure of Excitons in Hematite Fe 2 O 3

We present a computational study of the electronic structure of excitons in hematite Fe 2 O 3 using density functional theory (DFT) and time-dependent DFT theory. Excitons resulting from solar light absorption are precursors to photogenerated charge carriers, which are the active species in solar-to-fuel redox conversions. Upon vertical excitation, the lowest energy exciton is found to have a multideterminantal wave function, describing an electron excited from O 2p states to Fe 3d states of the nearest basal plane, possibly supporting a recent XUV-derived exciton structure with a radius of a single Fe–O bond. Subsequently, the exciton self-traps into stable electron–hole pair structures whose wave functions are well-approximated by single excitations conveniently labeled “HOMO–LUMO”. Self-trapping yields electron–hole pair structures whereby the electron state (electron polaron-like) is separated from the hole state (hole polaron-like) by 3, 5, 7, 9, ... basal planes in structures referred to as Exc-3, Exc-5, Exc-7, Exc-9, ... The natural transition orbitals (electron–hole states) exhibit a strong “localized” character. The hole state is best described as a (FeO 6 ) + octahedral moiety carrying ~0.65 h + charge, and overall ~70% of the hole is assigned to O 2p atomic states, mostly on the four (4) equatorial O atoms of the moiety. The electron state is best described as a (FeFe) – moiety carrying ~0.80 e – charge in Fe 3d atomic states. The lattice distortions around the hole site exhibit Fe–O bond shortening due to the removal of electron density from an orbital state with O–Fe–O antibonding character. The electron site shows Fe–Fe distance shortening due to excess electron density added to an orbital state with in-phase Fe–Fe interactions. Furthermore, these excitonic structures with increasing electron–hole separation can be viewed as the onset of exciton separation into isolated electrons and hole polarons.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Selective Adsorption of Thiol-Containing Molecules on Copper Sulfide Surfaces via Molecule–Surface Disulfide Bridges

Recent results in the fields of nanoenhanced agriculture and expanding interest in prebiotic chemistry have placed increased emphasis on understanding the chemically selective interaction of small molecules with the surfaces of metal sulfides. Here, we present an integrated experimental and computational study of the interaction of thiol-containing molecules with copper sulfide (covellite) surfaces in aqueous media. In situ Fourier-transform infrared (FTIR) measurements and ex situ X-ray photoelectron spectroscopy (XPS) measurements show that molecules bearing free thiol groups, including glutathione and cysteine, bind strongly to CuS (covellite) nanoparticles and to CuS (001) single crystals, while control studies show that similar molecules lacking the free thiol group exhibit much less binding. Additional experiments show that these thiol-containing molecules interact transiently with CuO nanoparticle surfaces but are readily removed by rinsing. The FTIR and XPS experiments demonstrate that adsorption of molecular thiols to CuS surfaces occurs in a chemically selective manner. Further experimental studies and density functional calculations show that the preferred mode of binding is through the surface S atoms, forming a Solid–S–S–Molecule disulfide linkage. While the role of disulfide linkages in controlling structure and function of proteins and other biomolecules is widely known, the formation of surface disulfide linkages as a motif for covalent molecular binding at surfaces has not been established previously.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Role of Surface Hydroxyls in Dehydration and Dehydrogenation of Formic Acid on Fe 3 O 4 (001)

Understanding the role of surface structure and hydroxylation in catalytic reactions on metal oxide surfaces is important for developing a mechanistic insight into the complex interface processes. Here, we investigate the reactivity of formic acid on reconstructed Fe 3 O 4 (001) using a combination of X-ray photoelectron spectroscopy, infrared reflection absorption spectroscopy, temperature-programmed reaction spectroscopy, low energy electron diffraction, and electronic structure calculations. We find that formic acid initially dissociates at low temperatures (< 80 K) into bidentate formate and a hydroxyl up to an initial dosed coverage of two HCOOH per Fe 3 O 4 (001) unit cell. At higher temperatures (> 450 K), formate largely decomposes along the dehydration pathway, producing CO and H 2 O, with dehydrogenation to CO 2 being a minority side reaction. As a first step, water formation leads to surface oxygen extraction via the Mars-van Krevelen mechanism. Computational studies reveal formate embedded in oxygen vacancies as a key intermediate in the CO formation mechanism. CO formation proceeds via two reaction pathways with desorption that peaks at 530 K on the hydroxyl-rich surface and 560 K on the hydroxyl-deficient surface. Atomic hydrogen coadsorption experiments and ab initio calculations reveal that the presence of surface hydroxyls reduces the CO formation barrier. Furthermore, these results highlight the complex interactions between substrate and intermediate species occurring during reactions on metal oxide surfaces.

Adsorption↗
DOE OSTI2025

Investigation of Electrocatalytic CO 2 Reduction on MXene Materials via First-Principles Simulations

Computational studies of CO 2 reduction to yield various products were carried out on the basal plane and edges of different MXene materials. The impact of vacancies upon Mo 2 TiC 2 T x , W 2 TiC 2 T x , and Ti 3 C 2 T x (T x = O and OH) was also examined for both edge and basal sites. Initial calibrations were carried out to generate surfaces with optimal oxide/hydroxide ratios and proper termination sites upon which various vacancy sites were explored to ensure an accurate model of the surfaces’ resting states. From this work, Mo 2 TiC 2 O was determined to exhibit the lowest theoretical overpotential for methane as determined by volcano plot analyses. At a large enough vacancy concentration, the CO 2 reduction reaction (CO 2 RR) is predicted to outcompete the hydrogen evolution reaction (HER) as the predominant reaction on the surface. In conclusion, when examining the edge of Mo 2 TiC 2 O, stronger CO 2 binding was exhibited to such an extent that the reaction was predicted to terminate after the generation of formate on the edge.

Lund, Colton [Argonne National Laboratory (ANL), A↗
DOE OSTI2025

Comparing Classical and Machine Learning Force Fields for Modeling Deformation of Metal–Organic Frameworks Relevant for Direct Air Capture

Deformation of metal–organic frameworks (MOFs) induced by adsorbate molecules can affect adsorption properties such as capacity and selectivity, but most computational studies of MOFs assume framework rigidity to simplify calculations. Although flexible force fields (FFs) for MOFs have been parametrized for specific materials, the generality of FFs for reliably modeling adsorbate-induced deformation to accuracy nearing that of density functional theory (DFT) has not been established. This work confirms using DFT calculations that adsorbate-induced deformation can affect CO 2 and H 2 O adsorption energies in a considerable fraction of MOFs promising for direct air capture (DAC). We then benchmark the efficacy of several general-purpose FFs in describing adsorbate-induced deformation for DAC against DFT. Our results show that current classical FFs are insufficient for describing MOF deformation, especially in cases of interest for DAC where strong interactions exist between adsorbed molecules and MOF frameworks. Some emerging machine learning force fields (MLFFs) we tested, particularly CHGNet, MACE-MP-0, and Equiformer V2, appear to be more promising than the classical FF for emulating the deformation behavior described by DFT. The best performing FF (CHGNet), however, fails to achieve the accuracy required for practical predictions with a mean absolute adsorption energy error of 0.124 eV.

adsorption↗
DOE OSTI2025

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗
DOE OSTI2025

Low Frequency Vibrational Modes of Two-Dimensional Lead-Free Metal Halide Double Perovskites

2D layered double perovskites of (S-MPA) 4 AgBiI 8 (MPA-AgBiI 8 ) and (S-MPA) 4 CuBiI 8 (MPA-CuBiI 8 ) (S-MPA, S-β-methylphenethylammonium) were synthesized with a hydrothermal method. The crystal structure of MPA-AgBiI 8 was determined using single-crystal X-ray diffraction (scXRD). Powder XRD (pXRD) data suggest that the crystal structure of MPA-CuBiI 8 is more complex than that of MPA-AgBiI 8 . UV-Vis electronic absorption spectra of these perovskites reveal a bandgap of 2.03 eV for both. Short exciton lifetime from time-resolved photoluminescence (TRPL) results and low PL intensity of the Cu-based perovskite suggest a high density of trap states within the bandgap. Low frequency Raman spectra of both materials show distinct peaks and a slightly higher frequency for the Cu-based perovskite than the Ag-based perovskite. Density functional theory (DFT) calculations were conducted to simulate the low frequency Raman spectra and help explain the different phonon modes, which are collective vibrations of metal halide bond bending and stretching within the Ag- and Cu-centered octahedra coupled with MPA libration and twisting within the inorganic layer. The DFT theory also quantified octahedral distortions in the two perovskites. Furthermore, this combined experimental and computational study provides new insights into the low frequency vibrations of 2D perovskites.

Halogens↗
DOE OSTI2020

Atomic-Scale Structure and Catalysis on Positively Charged Bimetallic Sites for Generation of H 2

In this paper, we report that a cationic bimetallic site consisting of one Pd and three Zn atoms (Pd 1 Zn 3 ) supported on ZnO (Pd 1 Zn 3 /ZnO) exhibits an extraordinarily high catalytic activity for the generation of H 2 through methanol partial oxidation (MPO) that is 2–3 orders of magnitude higher than that of a metallic Pd–Zn site on Pd–Zn nanoalloy (Pd–Zn/ZnO). Computational studies uncovered that the positively charged Pd atom of the subnanometer Pd 1 Zn 3 bimetallic site largely decreases the activation barrier for dehydrogenation of methanol as compared to a metallic Pd atom of Pd–Zn alloy, thus switching the rate-determining step of MPO from methanol dehydrogenation over a Pd–Zn alloy with high barrier to the O 2 dissociation step on a cationic Pd 1 Zn 3 site with a low barrier, which is supported by our kinetics studies. The significantly higher catalytic activity and selectivity for H 2 production over a cationic bimetallic site suggest a new approach to design bimetallic catalysts.

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

Thermal Stability of Au Rhombic Dodecahedral Nanocrystals Can Be Greatly Enhanced by Coating Their Surface with an Ultrathin Shell of Pt

Rhombic dodecahedral nanocrystals have been considered particularly difficult to synthesize because they are enclosed by {110}, a low-index facet with the greatest surface energy. Recently, we demonstrated the use of seed-mediated growth for the facile and robust synthesis of Au rhombic dodecahedral nanocrystals (AuRD). While the unique shape and surface structure of AuRD are desirable for potential applications in plasmonics and catalysis, respectively, their high surface energy makes them highly susceptible to thermal degradation. Here we demonstrate that it is feasible to greatly improve the thermal stability with some sacrifice to the plasmonic properties of the original AuRD by coating their surface with an ultrathin shell made of Pt. Our in situ electron microscopy analysis indicates that the ultrathin Pt coating can increase the thermal stability from 60 up to 450 °C, a trend that is also supported by the results from a computational study.

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

Defect Engineering of WO 3 by Rapid Flame Reduction for Efficient Photoelectrochemical Conversion of Methane into Liquid Oxygenates

Photoelectrochemical (PEC) conversion is a promising way to use methane (CH 4 ) as a chemical building block without harsh conditions. However, the PEC conversion of CH 4 to value-added chemicals remains challenging due to the thermodynamically favorable overoxidation of CH 4 . Here, we report WO 3 nanotube (NT) photoelectrocatalysts for PEC CH 4 conversion with high liquid product selectivity through defect engineering. By tuning the flame reduction treatment, we carefully controlled the oxygen vacancies of WO 3 NTs. The optimally reduced WO 3 NTs suppressed overoxidation of CH4 showing a high total C1 liquid selectivity of 69.4% and a production rate of 0.174 μmol cm –2 h –1 . Scanning electrochemical microscopy revealed that oxygen vacancies can restrain the production of hydroxyl radicals, which, in excess, could further oxidize C1 intermediates to CO 2 . Additionally, band diagram analysis and computational studies elucidated that oxygen vacancies thermodynamically suppress overoxidation. In conclusion, this work introduces a strategy for understanding and controlling the selectivity of photoelectrocatalysts for direct conversion of CH 4 to liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗