Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, Molecular”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Modeling Carbon and Hydrocarbon Molecular Structures in EZTB

A software module that models the electronic and mechanical aspects of hydrocarbon molecules and carbon molecular structures on the basis of first principles has been written for incorporation into, and execution within, the Easy (Modular) Tight-Binding (EZTB) software infrastructure, which is summarized briefly in the immediately preceding article. Of particular interest, this module can model carbon crystals and nanotubes characterized by various coordinates and containing defects, without need to adjust parameters of the physical model. The module has been used to study the changes in electronic properties of carbon nanotubes, caused by bending of the nanotubes, for potential utility as the basis of a nonvolatile, electriccharge- free memory devices. For example, in one application of the module, it was found that an initially 50-nmlong carbon, (10,10)-chirality nanotube, which is a metallic conductor when straight, becomes a semiconductor with an energy gap of .3 meV when bent to a lateral displacement of 4 nm at the middle.

Lee, Seungwon↗
DOE OSTI2025

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗
DOE OSTI2025

RLMolLM: Reinforcement Learning-Enhanced Language Model Framework for Inverse Molecular Design

Inverse molecular design faces significant challenges due to vast chemical space and complex property requirements. While language models show promise for molecular generation, they struggle with validity, multi-property optimization, and structural constraints. This work presents RLMolLM, a reinforcement learning framework combining Proximal Policy Optimization (PPO) with genetic algorithms to address these limitations. Our approach optimizes multiple user-specified properties including quantitative estimates of drug-likeness (QED), synthetic accessibility (SA), and ADMET (absorption, distribution, metabolism, excretion, and toxicity) endpoints without requiring complete model retraining, while maintaining capability for scaffold-constrained generation where specific substructures must be preserved. We outperform state-of-the-art methods for molecular optimization, achieving best QED scores across GDB13, Moses, and Zinc datasets with up to 31% improvement over previous methods while maintaining excellent validity, uniqueness, and novelty metrics. For simultaneous multi-property optimization, our framework achieves substantial improvements in ADMET properties including 4.5-fold reduction in hERG toxicity and enhanced Caco-2 permeability compared to Moses dataset. Under structural constraints, the framework significantly improves molecular validity while preserving scaffolds and effectively optimizing properties. In conclusion, this versatile solution advances pharmaceutical and materials molecular design through effective integration of reinforcement learning and genetic algorithms with multi-property optimization and scaffold preservation.

Genetic algorithms↗
NASA NTRS1974

Model calculations for diffuse molecular clouds

A local steady-state isobaric model of the interstellar clouds is developed and applied to the OAO-3 observations. The model allows for cooling by C II fine-structure transitions and for inhomogeneities resulting from attenuations in the external radiations which heat the cloud constituents. Factors affecting the thermal properties of the clouds are emphasized; these include H and H2 balance, interstellar grain heating, electron density, and calcium equilibrium. The thermodynamic consequences of the attenuation of the ultraviolet interstellar radiation field are studied. The effect of varying the most important parameters and their allowed ranges is considered. Results indicate that the model is capable of reproducing the H2 fractional abundances and kinetic temperatures observed by OAO-3.

Glassgold, A. E.↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
NASA NTRS2016

A DSMC Surface Chemistry Model for Carbon-Based Ablators

A detailed molecular surface chemistry model for the DSMC (Direct Simulation Monte Carlo) method is proposed and implemented into the SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) DSMC solver. Molchanova et al. constructed a molecular model for surface recombination in DSMC that includes different surface processes (adsorption, desoprtion, Eley-Rideal and Langmuir-Hinshelwood). All surface processes can be divided into two groups: surface mechanisms, which involve only the particle adsorbed by the surface (desorption and Langmuir-Hinshelwood), and impact mechanisms, which also involve gas-phase particles (adsorption, Eley-Rideal). Using a similar approach, the 14-reaction kinetic model of oxygen-carbon interaction suggested by Zhlukhtov and Abe, as well as more recent models by Alba et al., Poovathinghal et al., and a new model developed in the scope of this work, are implemented in SPARTA. The computational results for the different oxidation models are compared with experimental results from Murray et al. (oxidation of a vitreous carbon surface due to a hyperthermal beam of O and O2), with a particular focus on fluxes, angular and Time-Of-Flight distributions of scattered particles.

oxidation↗
NASA NTRS2006

Modeling Ellipsometry Measurements of Molecular Thin-Film Contamination on Genesis Array Samples

The discovery of a molecular thin-film contamination on Genesis flown array samples changed the course of preliminary assessment strategies. Analytical techniques developed to measure solar wind elemental abundances must now compensate for a thin-film contamination. Currently, this is done either by experimental cleaning before analyses or by depth-profiling techniques that bypass the surface contamination. Inside Johnson Space Center s Genesis dedicated ISO Class 4 (Class 10) cleanroom laboratory, the selection of collector array fragments allocated for solar wind analyses are based on the documentation of overall surface quality, visible surface particle contamination greater than 1 m, and the amount of thin film contamination measured by spectroscopic ellipsometry. Documenting the exact thickness, surface topography, and chemical composition of these contaminates is also critical for developing accurate cleaning methods. However, the first step in characterization of the molecular film is to develop accurate ellipsometry models that will determine an accurate thickness measurement of the contamination film.

Calaway, Michael J.↗
DOE OSTI2024

Multinucleon transfer in 48 Ti + 48 Ti collisions at 11.5 MeV/nucleon

In this work, experimental data of projectile fragment distributions from a preliminary experimental study of the reaction of 11.5 MeV/nucleon 48 Ti on 48 Ti analyzed with the MARS recoil separator at the Cyclotron Institute of Texas A&M University are presented. Production cross sections, momentum distributions and excitation-energy distributions are extracted and compared with calculations with the Deep Inelastic Transfer (DIT) model and the Constrained Molecular Dynamics model (CoMD). Despite the limited extent of the present dataset, an overall agreement of the models with the data is found pointing at the prevailing nucleon-exchange character of the collisions at this energy.

Souliotis, Georgios↗
NASA NTRS1992

A source model for the L134N molecular cloud

The dark molecular cloud L134N is observed at millimeter wavelengths in the CS (J = 3-2), OCS (J = 7-6), and HDO (1 sub 11-1 sub 10) transitions. The CS (J = 3-2) transition was observed at four positions within the cloud, while the other two transitions were observed at one position each. Fractional abundances in the LTE approximation are calculated for each emission line detected. L134N appears to have a high-density core characterized by NH3, C3H2, and H(C-13)O(+) emission maps. A lower density envelope characterized by C(0-18), CS (J = 2-1), and SO emission surrounds the core. There appears to be a gas-phase oxygen abundance gradient in L134N with atomic oxygen depleted in the high-density core. Observed molecular distributions within L134N can be explained by a model in which chemical and physical processes in icy-dust-grain mantles influence the gas-phase molecular abundances.

Swade, Daryl A.↗
NASA NTRS2006

Isotope Effects in Collisional VT Relaxation of Molecular Hydrogen

A simple exponential-potential model of molecular collisions leads to a two-parameter analytic expression for rates of collisionally induced vibrational-translation (VT) energy exchange that has been shown to be accurate over variations of orders of magnitude as a function of temperature in a variety of systems. This includes excellent agreement with reported experimental and theoretical results for the fundamental self-relaxation rate of molecular hydrogen H2(v = 1) + H2 yields H2(v = 0) + H2. The analytic rate successfully follows the five-orders-of-magnitude change in experimental values for the temperature range 50-2000 K. This approach is now applied to isotope effects in the vibrational relaxation rates of excited HD and D2 in collision with H2: HD(v = 1)+H2 yields HD(v = 0)+H2 and D2(v = 1)+H2 yields D2(v = 0)+H2. The simplicity of the analytic expression for the thermal rate lends itself to convenient application in modeling the evolving vibrational populations of molecular hydrogen in shocked astrophysical environments.

Bieniek, R. J.↗
NASA NTRS1985

Spectral identification/elimination of molecular species in spacecraft glow

Computer models of molecular electronic and vibrational emission intensities were developed. Known radiative emission rates (Einstein coefficients) permit the determination of relative excited state densities from spectral intensities. These codes were applied to the published spectra of glow above shuttle surface and to the Spacelab 1 results of Torr and Torr. The theoretical high-resolution spectra were convolved with the appropriate instrumental slit functions to allow accurate comparison with data. The published spacelab spectrum is complex but N2+ Meinel emission can be clearly identified in the ram spectrum. M2 First Positive emission does not correlate well with observed features, nor does the CN Red System. Spectral overlay comparisons are presented. The spectrum of glow above shuttle surfaces, in contrast to the ISO data, is not highly structured. Diatomic molecular emission was matched to the observed spectral shape. Source excitation mechanisms such as (oxygen atom)-(surface species) reaction product chemiluminescence, surface recombination, or resonance fluorescent re-emission will be discussed for each tentative assignment. These assignments are the necessary first analytical step toward mechanism identification. Different glow mechanisms will occur above surfaces under different orbital conditions.

Green, B. D.↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Geometry-complete diffusion for 3D molecule generation and optimization

Abstract Generative deep learning methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a denoising diffusion framework. However, such methods are unable to learn important geometric properties of 3D molecules, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which notably hinders their ability to generate valid large 3D molecules. In this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset and the larger GEOM-Drugs dataset, respectively. Importantly, we demonstrate that GCDM’s generative denoising process enables the model to generate a significant proportion of valid and energetically-stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but can be repurposed to consistently optimize the geometry and chemical composition of existing 3D molecules for molecular stability and property specificity, demonstrating new versatility of molecular diffusion models. Code and data are freely available on GitHub .

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

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
NASA NTRS2026

Gateway Element and Payload Materials Outgassing Analyses: HALO, HERMES, and ERSA

Gateway was intended to be humanity’s first space station around the Moon, but its development has been paused as the National Aeronautics and Space Administration (NASA) shifts focus to achieving the United States’ National Space Policy goals. Instead of an orbiting lunar outpost, NASA will now pursue the development of a lunar surface base to support a sustained human presence on the Moon. Before the program’s pause, Gateway’s Induced Environments team worked to ensure payloads and elements (i.e., modules) complied with induced environment requirements. Methods developed and insights gained from this work will have applicability to NASA’s Moon Base and the potential repurposing of Gateway elements and payloads, as well as to induced environments modeling for future space stations. The Gateway program’s induced environment included molecular contamination, electric thruster plume sputter and redeposition, and lunar dust transfer from the Human Landing System (HLS). Primary sources of external molecular contamination included materials outgassing, chemical thruster plume contamination, and vacuum venting. The focus of this paper will be on element- and payload-level materials outgassing analyses performed for Gateway Configuration 1, extending the previously-developed framework for Gateway system-level external molecular contamination modeling. Gateway Configuration 1 consisted of the Power and Propulsion Element (PPE) and the Habitation and Logistics Outpost (HALO). It also included payloads like the European Radiation Sensor Array (ERSA) attached to PPE and the Heliophysics Environmental and Radiation Measurement Experiment Suite (HERMES) attached to HALO. The element- and payload-level analyses to be introduced in this paper for HALO, HERMES, and ERSA enabled high-fidelity descriptions of Gateway’s external molecular contamination environment. Approaches to geometric modeling, meshing, outgassing rate assignment, molecular transport modeling, and analysis methodology will be presented. Element and payload contaminant deposition onto sensitive Gateway receiver surfaces will be summarized and results compared to induced environment requirements. While these results incorporate refinements made over the course of the program, they were not intended to be final. Therefore, modeling assumptions and inputs, potential improvements, and lessons-learned will be documented to inform future work on Moon Base, repurposed elements and payloads, and other space stations.

Gateway↗
NASA NTRS1980

Laboratory and modeling studies of chemistry in dense molecular clouds

A chemical evolutionary model with a large number of species and a large chemical library is used to examine the principal chemical processes in interstellar clouds. Simple chemical equilibrium arguments show the potential for synthesis of very complex organic species by ion-molecule radiative association reactions.

Huntress, W. T., Jr.↗
NASA NTRS1988

A molecular dynamics study of chloride binding by the cryptand SC24

The capture of chloride from water by the tetraprotonated form of the spherical macrotricyclic molecule SC24 was studied using molecular dynamics simulation methods. This model ionophore represents a broad class of molecules which remove ions from water. Two binding sites for the chloride were found, one inside and one outside the ligand. These sites are separated by a potential energy barrier of approximately 20 kcal mol-1. The major contribution to this barrier comes from dehydration of the chloride. The large, unfavorable dehydration effect is compensated for by an increase in electrostatic attraction between the oppositely charged chloride and cryptand, and by energetically favorable rearrangements of water structure. Additional assistance in crossing the barrier and completing the dehydration of the ion is provided by the shift of three positively charged hydrogen atoms of the cryptand towards the chloride. This structural rigidity is partially responsible for its selectivity.

NASA Discipline Exobiology↗
NASA NTRS1994

Modeling radiative transfer in molecular clouds. 1: HCO(+) in the star-forming region W49A North

A new general multilevel, non-Local Thermodynamic Equilibrium (LTE) radiative transfer code, valid for any velocity field, is applied to HCO(+) observations of W49A North. Three classes of collapse models are considered: free-fall collapse (v proportional to 1/sq. root of r), rho proportional to r(exp -3/2) throughout the molecular cloud, successfully reproduces the features of the observations and gives the best fit to the J = 1-0 and J = 3-2 profiles both toward the prominent H II component G of W49A North and off the center. In addition to a slow radial fall-off of density, the theoretical modeling implies the following for the molecular cloud: the large line widths result from motions occurring within the inner 1 pc, and there are probably one or more fragments with peculiar velocities within this same region.

Dickel, Helene R.↗