Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Modeling”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reversing inverse Hall-Petch and direct computation of Hall-Petch coefficients

A 2D bicrystal atomistic model of dislocation transmission through a Σ11<101>{131} symmetric-tilt grain boundary reveals details of Hall-Petch breakdown in the single dislocation regime in Al, Ni, and Cu. Partly based on a previous study, this research determines the stress required for a single dislocation to be transmitted through the Σ11 boundary and finds that single dislocation transmission typically deviates from Hall-Petch behavior because the leading partial dislocation becomes trapped in the boundary and emits glissile grain boundary disconnections that cause boundary sliding deformation at stresses well below boundary transmission stresses. Thus, mechanisms of inverse Hall-Petch, namely grain boundary shear and sliding, are shown to operate in lieu of grain boundary transmission in the single dislocation regime. However, by controlling the applied stresses, inverse Hall-Petch can be reversed and typical Hall- Petch grain boundary transmission regained. This, in turn, allows the direct computation of Hall-Petch coefficients. Further, a second dislocation on the same slip plane preserves typical Hall-Petch behavior for Al and Ni and does not lead to boundary sliding events, thus implying a critical grain size for Hall-Petch breakdown based on the pileup model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Opportunities in multiscale modeling of mosquito-borne flaviviruses

Mosquito-borne flaviviruses, such as Zika, dengue, West Nile, and yellow fever virus, represent a growing public health concern due to their widespread distribution and the severe diseases they cause. These viruses are difficult to control as climate change and urbanization help mosquitoes expand into new areas, increasing the risk of outbreaks. Mathematical models play a key role in understanding their spread, providing insights at every level—from how the virus multiplies inside cells to how it circulates through entire populations. This review examines various approaches used in modeling arboviruses, including microscale models that focus on cellular and molecular dynamics, mesoscale models that address within-host processes, and macroscale models that capture population-level transmission. We briefly summarize the methodology used for models at each scale, which primarily consists of sets of differential equations with parameters that represent physical rates of change for different subprocesses. We particularly highlight how temperature affects virus transmission, which is key to understanding the impact of climate change. We also show how multiscale models can connect viral replication, immune response, and the spread of infection at a larger scale. This is essential for developing better vaccines and treatments, evaluating disease control measures, predicting the impact of climate change, and improving public health responses to outbreaks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

A fast computational framework for the design of solvent-based plastic recycling processes

Multicomponent plastics cannot be processed using mechanical recycling technologies, hindering efforts to deal with plastic waste. Multicomponent plastics include multilayer plastic films, which are widely used for food and healthcare packaging. Multilayer films combine several layers (potentially dozens) of different polymers to protect products from external factors (e.g., oxygen, water, temperature, shock, and light). Solvent-based separation processes have emerged as a promising alternative to recycle these complex materials. For instance, the Solvent-Targeted Recovery and Precipitation (STRAP TM ) process uses sequential solvent washes to selectively dissolve and separate constituent polymers from multicomponent plastic waste, including films. STRAP TM process design (separation sequence, type of solvents, and operating conditions) changes significantly depending on the design of the multilayer plastic film (e.g., number, types, and proportions of polymers). The ability to quickly quantify the economic and environmental benefits of diverse STRAP TM process designs is essential to accelerate the development of sustainable recycling processes and more recyclable multilayer film products. In this work, we present a fast computational framework that integrates molecular-scale models, process modeling, and techno-economic and life cycle analysis to quickly evaluate STRAP TM designs. The computational framework is general and can be used to study the processing of complex multilayer plastic waste streams that contain many layers. Furthermore, we highlight the different uses of the framework via targeted case studies.

Computational framework↗
DOE OSTI2022

Modeling of linear nanopores in a-SiO 2 tuning pore surface structure

Here, new strategies are presented for the generation of models for linear nanopores in amorphous silica (a-SiO 2 ) with surface structure tuned to match experimental observations. Specifically, the models successfully target not only the overall density of surface silanol groups, but also the proportion of geminal versus mono silanols for which additional experimental NMR data is reported. The latter quantity has not been appropriately described in previous modeling, and in fact has typically not been considered. Strategies include “pore drilling” of bulk a-SiO 2 , and “cylindrical resist” methodology forming a-SiO 2 around a cylindrical exclusion region, followed by dehydroxylation and hydroxylation processes, respectively. However, these latter processes must be judiciously tailored in order to tune the proportion of geminals, in addition to the overall silanol density, to achieve experimental values. Such tailoring has not been incorporated into previous modeling. Another approach considered tunes surface structure of pores obtained by “pore drilling” through mild annealing.

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

Enhancing molecular design efficiency: Uniting language models and generative networks with genetic algorithms

This study examines the effectiveness of generative models in drug discovery, material science, and polymer science, aiming to overcome constraints associated with traditional inverse design methods relying on heuristic rules. Generative models generate synthetic data resembling real data, enabling deep learning model training without extensive labeled datasets. They prove valuable in creating virtual libraries of molecules for material science and facilitating drug discovery by generating molecules with specific properties. While generative adversarial networks (GANs) are explored for these purposes, mode collapse restricts their efficacy, limiting novel structure variability. To address this, we introduce a masked language model (LM) inspired by natural language processing. Although LMs alone can have inherent limitations, we propose a hybrid architecture combining LMs and GANs to efficiently generate new molecules, demonstrating superior performance over standalone masked LMs, particularly for smaller population sizes. This hybrid LM-GAN architecture enhances efficiency in optimizing properties and generating novel samples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Thermoplasticity of metallic glasses: Processing and applications

Owing to their extremely low shear viscosity and fragility, bulk metallic glasses (BMGs) are envisioned as new net-shaping material replacing conventionally used crystalline metals for various applications. The first part of this review describes the general technique and kinetics of thermoplastic forming. The second part elaborates on the thermal processing capability of metallic glasses from atomic- to meter-scale. The micro-/nano-structure obtained after various quenching rates followed by temperature, time and heating rate dependent thermoplastic forming is described in the third section. The deformation behavior and flow kinetics of BMGs in terms of composition, mold features, and applied conditions are elaborated in the fourth section. The variation of volume and enthalpy at the supercooled liquid region, strength retention of metallic glasses compared to other conventional metals and alloys, and kinetics of liquid fragility are given in the fifth section. Finite element modeling and molecular dynamics simulations of high-temperature deformation in BMGs are presented in the sixth part. Finally, thermally-formed BMGs used for different applications, including energy, biomedical and micro-optics, are presented in the final part. Altogether, this review provides an overview of shaping capabilities and modifications in the macro-scale properties and short-to-medium range order of BMGs upon thermoplastic forming.

36 MATERIALS SCIENCE↗
DOE OSTI2024

From biowaste to BioPave: Biological pathways for sequestration of anthropogenic CO 2 and enhancing durability of roadway infrastructures

Biomass-derived bio-oils are emerging as sustainable, low-carbon alternatives for construction materials, particularly in innovative pavement applications, 'BioPave'. Here, this study evaluates bio-oils from various sources, including algae and wastewater biomass, revealing variations in carbon, nitrogen, and sulfur contents, along with viscosity. We focused on their application in outdoor construction, assessing their resistance to thermal and solar radiation, and moisture. The study also examines their interaction with siliceous surfaces and their influence on asphalt adhesion in BioPave systems under moisture conditions. Laboratory experiments and computational modeling demonstrated molecular composition significantly influenced bio-oils' responses to thermal and UV exposure. Bio-oils rich in polar groups showed hardening upon thermal treatment, while those with higher concentrations of saturated aliphatics remained more stable. Furthermore, our findings highlight the significance of dosage control in maximizing beneficial effects of bio-oils in intermolecular interactions at bitumen-aggregate interface in BioPave applications. These findings offer insights into the potential of bio-oils in sustainable construction and the need for dosage control to optimize pavement performance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Correlating structure, self-assembly chemistry and conductivity of trithiocyanuric acid on Au(111)

The majority of candidates for simple model molecular-electronic components consist of a conductive π-conjugated hydrocarbon linker attached to at least two anchoring groups, such as thiols or isocyanides. It has been found that select molecules self-assemble on gold surfaces, creating one-dimensional conductive structures that act as “molecular wires”. Furthermore, these oligomers can form molecular bridges between gold nanoparticles, leading to the creation of simple molecular-electronic devices. This raises the question whether other π-conjugated molecular linkers could exhibit similar behavior that might offer a broader range of candidates for fabricating electronic devices. Trithiocyanuric acid (1,3,5-triazine-2,4,6-trithiol, TTCA) provides a possible candidate. TTCA (C3N3(SH)3) can exist as a trithiol or as a trithione in which hydrogens transfer to the sulfurs so that they are present with three C=N groups within the ring. TTCA exists naturally in the trithione form but converts into a trithiol when adsorbed onto an Ag(111) where it is vertically oriented. The structure of TTCA adsorbed on Au(111) is studied here using reflection-absorption infrared spectroscopy (RAIRS) where it is found to remain as the trithione isomer, but changes orientation as the coverage increases. Scanning-tunneling microscopy (STM) reveals that TTCA oligomerizes on Au(111) to form chains and triangular structures. The influence on molecular conductivity due to the differences in the adsorbate's isomeric structure was investigated using devices comprising either silver or gold nanoparticles deposited in the gap between gold nanoelectrodes. Both devices were found to conduct when dosed with TTCA, but the devices fabricated using silver were about 13 time more conductive than those made from gold nanoparticles, consistent with the π-conjugated structure formed on silver but not on gold. This implies that oligomers form both on silver and on gold and potentially increases the range of molecule-metal combinations that might be used to fabricate molecular-electronic devices.

Molecular electronics↗
DOE OSTI2025

Structural and Mechanistic Advances in the Chemistry of Methyl-Coenzyme M Reductase (MCR)

Methane represents 34% of U.S. energy consumption and is a major greenhouse gas related to the global carbon cycle and energy production. However, current industrial practices significantly increase atmospheric methane levels, necessitating a deeper understanding of its biosynthesis and oxidation. Methyl-coenzyme M reductase (MCR) is central to biological methane metabolism, catalyzing the final step of methanogenesis and the first step in anaerobic methane oxidation. It is also a key target for strategies to capture and transform methane into value-added chemicals.The active site of MCR is a buried Ni-based cofactor only accessible by the substrates via a 50 Å long tunnel. Although the Ni(I) state is required to initiate catalysis, capturing this state remains a challenge for the current structural techniques. Recent advances in structural biology using X-ray Free-Electron Laser serial crystallography have provided insights into MCR's inactive Ni(II) state at room temperature and show promise for capturing its active Ni(I) form.Our team has established several critical aspects of the MCR mechanism using a combination of experimental and computational studies. MCR uses CH 3 -SCoM and CoBSH as substrates, producing methane and a disulfide product CoMSSCoB. Kinetic analysis showed that productive substrate binding requires CH 3 -SCoM to bind first, inducing conformational changes that optimize the active site for subsequent CoBSH binding. Following substrate binding, four proposed methane production/oxidation mechanisms were examined, establishing whether the reaction proceeds through an organometallic methyl-nickel(III), methyl anion ion, or methyl radical intermediate. Experimental measurements using CoBSH analogs successfully slowed the reaction, allowing for mechanistic insight that demonstrated the methyl radical pathway, where the initial interactions involve homolytic cleavage of the methyl-sulfur bond, generating a methyl radical that quickly abstracts the thiol hydrogen atom of CoBSH to form methane. Computational studies further confirmed that, compared to other mechanisms, the methyl radical mechanism is thermodynamically more favorable and accessible under physiological conditions.Spectroscopic and computational studies challenged the conventional understanding of substrate binding in MCR by proposing an alternative positioning of CH 3 -SCoM and CoMSSCoB in the active site pocket. The research suggested that CH 3 -SCoM (substrate) and CoMSSCoB (product) bind via their sulfonate groups to the Ni(I) center of cofactor F 430 . This binding allows for the reaction without substrate reorganization in the pocket but would require a long-range electron transfer. Altogether, the work summarized in this review reflects our current understanding of the enzyme's catalytic mechanism and structural dynamics. This is essential for developing efficient methane conversion technologies that could mitigate its environmental impact while harnessing energy-storage properties.

03 NATURAL GAS↗
DOE OSTI2025

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Invariant Molecular Representations for Heterogeneous Catalysis

Catalyst screening is a critical step in the discovery and development of heterogeneous catalysts, which are vital for a wide range of chemical processes. In recent years, computational catalyst screening, primarily through density functional theory (DFT), has gained significant attention as a method for identifying promising catalysts. However, the computation of adsorption energies for all likely chemical intermediates present in complex surface chemistries is computationally intensive and costly due to the expensive nature of these calculations and the intrinsic idiosyncrasies of the methods or data sets used. This study introduces a novel machine learning (ML) method to learn adsorption energies from multiple DFT functionals by using invariant molecular representations (IMRs). To do this, we first extract molecular fingerprints for the reaction intermediates and later use a Siamese-neural-network-based training strategy to learn invariant molecular representations or the IMR across all available functionals. Our Siamese network-based representations demonstrate superior performance in predicting adsorption energies compared with other molecular representations. Notably, when considering mean absolute values of adsorption energies as 0.43 eV (PBE-D3), 0.46 eV (BEEF-vdW), 0.81 eV (RPBE), and 0.37 eV (scan+rVV10), our IMR method has achieved the lowest mean absolute errors (MAEs) of 0.18 0.10, 0.16, and 0.18 eV, respectively. These results emphasize the superior predictive capacity of our Siamese network-based representations. The empirical findings in this study illuminate the efficacy, robustness, and dependability of our proposed ML paradigm in predicting adsorption energies, specifically for propane dehydrogenation on a platinum catalyst surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗