Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A simple method for floating graphene oxide films facilitates nanoscale investigations of ion and water adsorption

Graphene oxide (GO) is a promising material for separations. Nanoscale GO thin films at the air/water interface are excellent experimental models to understand molecular-scale interactions of ions and water with GO. However, the characteristics of GO, such as functional groups and flake size, also affect the thin film properties making it difficult to make systematic studies with GO thin films. This paper reports a simple, reliable, and quick method of preparing ultra-thin GO films, irrespective of their origin, and demonstrates the new opportunities possible with the utilization of this method. The total amount of GO used to form the thin film is significantly less compared to previous examples in the literature, minimizing the dissolved GO in the subphase. X-ray reflectivity (XR) studies show that the majority of the GO film has 1.5 nm thickness over a macroscopic area (~100 cm 2 ) with very small roughness. Sum frequency generation (SFG) spectroscopy measurements show that H 2 O and D 2 O interact differently with GO films, a property that was not observed before. SFG data show that functional groups vary significantly between different commercially available GO samples. The differences are also characterized with XR at high resolution. X-ray fluorescence near total reflection (XFNTR) measurements show that these differences strongly affect ion adsorption and interfacial water behavior near GO, which are vital properties in separation applications. The results pave the way for future studies to elucidate the complex separation mechanisms with GO.

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

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

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

Probing the active sites of oxide encapsulated electrocatalysts with controllable oxygen evolution selectivity

Electrocatalysts encapsulated by nanoscopic overlayers can control the rate of redox reactions at the outer surface of the overlayer or at the buried interface between the overlayer and the active catalyst, leading to complex behavior in the presence of two competing electrochemical reactions. This study investigated oxide encapsulated electrocatalysts (OECs) comprised of iridium (Ir) thin films coated with an ultrathin (2–10 nm thick) silicon oxide (SiO x ) or titanium oxide (TiO x ) overlayer. The performance of SiO x |Ir and TiO x |Ir thin film electrodes towards the oxygen evolution reaction (OER) and Fe(II)/Fe(III) redox reactions were evaluated. An improvement in selectivity towards the OER was observed for all OECs. Overlayer properties, namely ionic and electronic conductivity, were assessed using a combination of electroanalytical methods and molecular dynamics simulations. SiO x and TiO x overlayers were found to be permeable to H 2 O and O 2 such that the OER can occur at the MO x |Ir (M = Ti, Si) buried interface, which was further supported with molecular dynamics simulations of model SiO 2 coatings. In contrast, Fe(II)/Fe(III) redox reactions occur to the same degree with TiO x overlayers having thicknesses less than 4 nm as bare electrocatalyst, while SiO x overlayers inhibit redox reactions at all thicknesses. This observation is attributed to differences in electronic transport between the buried interface and outer overlayer surface, as measured with through-plane conductivity measurements of wetted overlayer materials. These findings reveal the influence of oxide overlayer properties on the activity and selectivity of OECs and suggest opportunities to tune these properties for a wide range of electrochemical reactions.

08 HYDROGEN↗
DOE OSTI2025

Structural incorporation into goethite fractionates rare earth elements

Up to 20% of rare earth elements (REEs) in ion adsorption deposits (IADs) are associated with iron oxide minerals, primarily goethite. Often termed “non-extractable”, goethite-hosted REEs are thought to be structurally incorporated into the mineral lattice. The large mismatch in size and charge density between REEs (ionic radius, r = 0.86–1.03 Å) and Fe 3+ (r = 0.65 Å), however, makes direct substitution energetically unfavorable. To determine REE compatibility with and incorporation into goethite on the atomic level, we used X-ray pair distribution function analysis (PDF) and LIII-edge extended X-ray absorption fine structure (EXAFS) spectroscopy. The compatibility of REEs with goethite and the precursor ferrihydrite (FH) is Lu ≥ Yb ≫ Dy > Nd for both phases. Nd and Dy primarily form secondary amorphous phases, with <30% Nd and Dy incorporated into goethite. In the FH precursor at pH 6.8, Yb and Lu assumed a local REE-OOH like structure with next nearest neighbor Fe. The Yb, Lu, and Nd-FH samples were also matured at ambient conditions for 100 days; despite the presence of only ∼5% goethite, Lu and Yb were 42% and 100% in goethite-like structural environments, respectively, whereas the PDF and EXAFS of Nd showed little evidence of any incorporation. Using ab initio molecular dynamics (AIMD) to model the EXAFS, we determined the presence of protonated Fe vacancies, edge-sharing with structural Lu and Yb, likely helped accommodate these REEs in the goethite structure. Incorporation into Fe oxyhydroxides thus potentially fractionates the REEs during weathering associated with formation of lateritic and ion adsorption deposits.

Mergelsberg, Sebastian T. [Pacific Northwest Natio↗
DOE OSTI2020

Precipitation of calcium–alumino–silicate–hydrate gels: The role of the internal stress

Concrete gains its strength from the precipitation of a calcium–alumino–silicate–hydrate (C–A–S–H) colloidal gel, which acts as its binding phase. However, despite concrete’s ubiquity in the building environment, the atomic-scale mechanism of C–A–S–H precipitation is still unclear. Here, we use reactive molecular dynamics simulations to model the early-age precipitation of a C–A–S–H gel. We find that, upon gelation, silicate and aluminate precursors condensate and polymerize to form an aluminosilicate gel network. Notably, we demonstrate that the gelation reaction is driven by the existence of a mismatch of atomic-level internal stress between Si and Al polytopes, which are initially experiencing some local tension and compression, respectively. The polymerization of Si and Al polytopes enables the release of these competitive stresses.

Zhao, Cheng↗
DOE OSTI2022

On the meaning of Berry force for unrestricted systems treated with mean-field electronic structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin–orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction—even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson–Fischer point (which can lead to numerical instabilities). Here, we do not address magnetic fields here.

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

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

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

Continuum shock mixture models for Ni+Al multilayers: Inert mesoscale simulations

Mesoscale modeling of shock waves in Ni+Al multilayers poses significant challenges that are due, in part, to shock-induced chemical reactions. Current modeling approaches utilize reactive molecular dynamics (MD), but they are limited to resolving domains of only a few hundred nanometers. In contrast, actual multilayer superlattices can be tens of micrometers thick, and they exhibit non-ideal (i.e., wavy) interfaces. The second part of our research builds upon previous work developing physically based, thermodynamically complete equations of state for various Ni and Al intermetallic compositions. Here, we introduce a novel workflow for high-fidelity mesoscale simulations of Ni+Al multilayers using a continuum hydrocode. By increasing the simulation domain size beyond MD limitations (e.g., 2 × 6 μm 2 ) and incorporating explicit interfacial roughness, we investigate the shock response of Ni+Al multilayers at previously unexplored scales. Our experimental design encompasses nine multilayer geometries with varying roughness amplitudes and tilt angles (θ = 15°, 30°, and 45°), alongside 19 flyer impact velocities ranging from 0.3 to 3.0 km/s, resulting in a total of 171 high-fidelity simulations. The bulk shock state from inert 2D mesoscale simulations aligns with the law of mixtures, while temperature and pressure fluctuations strongly correlate with multilayer geometry types. A new metric dubbed the “hot spot probability integral” shows a greater dependence on a tilt angle than interfacial roughness.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2026

Epitaxial registry and transport properties in V4O7 thin films

We report the growth of V4O7 thin films deposited simultaneously on amorphous SiO2 and crystalline c-cut Al2O3 substrates. x-ray diffraction shows that films grown on SiO2 are polycrystalline with no preferred orientation, while films grown on Al2O3 exhibit a well-defined out-of-plane orientation and an in-plane registry consistent with epitaxial growth. Transmission electron microscopy confirms the crystallographic relationship between the V4O7 film and the Al2O3 substrate. Atomic force microscopy indicates substantially lower surface roughness for films grown on Al2O3 (~6 nm) compared to those grown on SiO2 (~22 nm), and Raman spectroscopy confirms stabilization of the same V4O7 phase on both substrates. Electrical transport measurements reveal a metal–insulator transition near 244 K for both substrates, with thermal hysteresis not exceeding ~1 K. Although the transition temperature remains essentially unchanged, films grown on Al2O3 exhibit higher conductivity over the entire temperature range, exceeding that of films grown on SiO2 by approximately two orders of magnitude at 100 K and by a factor of five at 300 K. These results indicate that the conductivity differences are consistent with variations in microstructural connectivity associated with crystallographic order.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Crystallographic registry in epitaxial V3O5 thin films

Crystallographic orientation plays a central role in determining transport properties in correlated oxides with low-hysteresis metal–insulator transitions. Here, V3O5 thin films were grown on Al2O3 (0001) and their crystallographic registry was evaluated using x-ray diffraction and cross-sectional transmission electron microscopy. θ–2θ measurements show a single out-of-plane orientation defined by the (202) reflection. Azimuthal φ-scans of the V3O5 (310) reflection exhibit twelve discrete maxima over 360°, indicating a finite set of in-plane orientations imposed by the substrate symmetry. High-resolution TEM and fast Fourier transform analysis confirm crystallographic coherence at the film–substrate interface. The electrical conductivity exhibits a metal–insulator transition at TMIT ≈ 423 K with no measurable thermal hysteresis and a total change of ~1.6 orders of magnitude between 300 and 480 K. These results establish crystallographic registry in epitaxial V3O5 thin films and provide a structurally well-defined system for future studies of orientation-dependent transport.\r\nT

36 MATERIALS SCIENCE↗
DOE OSTI2023

Small molecules disaggregate alpha-synuclein and prevent seeding from patient brain-derived fibrils

The amyloid aggregation of alpha-synuclein within the brain is associated with the pathogenesis of Parkinson’s disease (PD) and other related synucleinopathies, including multiple system atrophy (MSA). Alpha-synuclein aggregates are a major therapeutic target for treatment of these diseases. We identify two small molecules capable of disassembling preformed alpha-synuclein fibrils. The compounds, termed CNS-11 and CNS-11g, disaggregate recombinant alpha-synuclein fibrils in vitro, prevent the intracellular seeded aggregation of alpha-synuclein fibrils, and mitigate alpha-synuclein fibril cytotoxicity in neuronal cells. Furthermore, we demonstrate that both compounds disassemble fibrils extracted from MSA patient brains and prevent their intracellular seeding. They also reduce in vivo alpha-synuclein aggregates in C. elegans . Both compounds also penetrate brain tissue in mice. A molecular dynamics–based computational model suggests the compounds may exert their disaggregating effects on the N terminus of the fibril core. These compounds appear to be promising therapeutic leads for targeting alpha-synuclein for the treatment of synucleinopathies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

CSGL: chemical synthesis graph learning for molecule representation

Abstract Motivation Molecule representation learning (MRL) translates molecules into a real vector space, serving as input to downstream tasks in biology, chemistry, and computer science. This article introduces a chemical synthesis graph learning (CSGL) framework, which enhances MRL by considering both the atomic structures of molecules and their roles in chemical reactions through a hierarchical graph representation. Specifically, molecules are first modeled based on their molecular graphs, which capture atomic-level structural information. They are then further refined using a chemical synthesis graph, where nodes represent reactant and product molecule sets, and edges encode chemical transformations between reactants and products (e.g. changes in molecular structures). CSGL optimizes molecular embeddings of reactant and product nodes in a fashion that ensures the embeddings conform to a chemical balance constraint. Results Experimental results show that our method CSGL achieves strong performance on a variety of tasks, including product prediction, reaction classification, and molecular property prediction. Availability and implementation https://github.com/li-2023/CSGL.

Biochemistry & Molecular Biology↗
DOE OSTI2022

The histone code of the fungal genus Aspergillus uncovered by evolutionary and proteomic analyses

Chemical modifications of DNA and histone proteins impact the organization of chromatin within the nucleus. Changes in these modifications, catalysed by different chromatin-modifying enzymes, influence chromatin organization, which in turn is thought to impact the spatial and temporal regulation of gene expression. While combinations of different histone modifications, the histone code, have been studied in several model species, we know very little about histone modifications in the fungal genus Aspergillus, whose members are generally well studied due to their importance as models in cell and molecular biology as well as their medical and biotechnological relevance. Here, we used phylogenetic analyses in 94 Aspergilli as well as other fungi to uncover the occurrence and evolutionary trajectories of enzymes and protein complexes with roles in chromatin modifications or regulation. We found that these enzymes and complexes are highly conserved in Aspergilli, pointing towards a complex repertoire of chromatin modifications. Nevertheless, we also observed few recent gene duplications or losses, highlighting Aspergillus species to further study the roles of specific chromatin modifications. SET7 (KMT6) and other components of PRC2 (Polycomb Repressive Complex 2), which is responsible for methylation on histone H3 at lysine 27 in many eukaryotes including fungi, are absent in Aspergilli as well as in closely related Penicillium species, suggesting that these lost the capacity for this histone modification. We corroborated our computational predictions by performing untargeted MS analysis of histone post-translational modifications in Aspergillus nidulans. This systematic analysis will pave the way for future research into the complexity of the histone code and its functional implications on genome architecture and gene regulation in fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Modeling carrier mobility in graphene as a sensitive probe of molecular magnets

Carrier mobility in graphene on a GaAs substrate and its change due to the adsorption of molecular magnets, Mn 12 and [Mn 3 ] 2 , on the surface of graphene, is calculated from first principles. Phonon limited mobility is also calculated for comparison. For Mn 12 adsorption on graphene, the mobility is compared for different organic ligands of Mn 12 (-H, -CH 3 , and - CHCl 2 ), while for [Mn 3 ] 2 dimers, the mobility is calculated for different linkers that yield the ferromagnetic and anti-ferromagnetic configurations of [Mn 3 ] 2 , as well as for different orientations of the molecule. Significant changes in carrier mobility due to the adsorption of the molecules and due to differences in linkers suggest mobility measurement as a possible sensitive probe of magnetic molecules.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantum many-body simulations of the two-dimensional Fermi-Hubbard model in ultracold optical lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. Finally, we observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Order fractionalization in a Kitaev-Kondo model

We describe a mechanism for order fractionalization in a two-dimensional Kondo lattice model, in which electrons interact with a gapless spin liquid of Majorana fermions described by the Yao-Lee (YL) model. When the Kondo coupling to the conduction electrons exceeds a critical value, the model develops a superconducting instability into a state with a a spinor order parameter with charge 4 and spin ( = 1 2. The broken symmetry state develops a gapless Majorana Dirac cone in the bulk. By including an appropriate gauge string, we can show that the charge 4, spinorial order develops off-diagonal long range order that allows electrons to coherently tunnel arbitrarily long distances through the spin liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling single-molecule stretching experiments using statistical thermodynamics

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. In conclusion, the asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Flux Hypothesis for Odd Transport Phenomena

Onsager's regression hypothesis makes a fundamental connection between macroscopic transport phenomena and the average relaxation of spontaneous microscopic fluctuations. This relaxation, however, is agnostic to odd transport phenomena, in which fluxes run orthogonal to the gradients driving them. To account for odd transport, we generalize the regression hypothesis, postulating that macroscopic linear constitutive laws are, on average, obeyed by microscopic fluctuations, whether they contribute to relaxation or not. From this "flux hypothesis," Green-Kubo and reciprocal relations follow, elucidating the separate roles of broken time-reversal and parity symmetries underlying various odd transport coefficients. As an application, we derive and verify the Green-Kubo relation for odd collective diffusion in chiral active matter, first in an analytically tractable model and subsequently through molecular dynamics simulations of concentrated active spinners.

Hargus, Cory↗