Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Formation, chemical evolution and solidification of the dense liquid phase of calcium (bi)carbonate

Metal carbonates, which are ubiquitous in the near-surface mineral record, are a major product of biomineralizing organisms and serve as important targets for capturing anthropogenic CO 2 emissions. However, pathways of carbonate mineralization typically diverge from classical predictions due to the involvement of disordered precursors, such as the dense liquid phase (DLP), yet little is known about DLP formation or solidification processes. Using in situ methods we report that a highly hydrated bicarbonate DLP forms via liquid–liquid phase separation and transforms into hollow hydrated amorphous CaCO 3 particles. Acidic proteins and polymers extend DLP lifetimes while leaving the pathway and chemistry unchanged. Molecular simulations suggest that the DLP forms via direct condensation of solvated Ca 2+ •(HCO 3 – ) 2 complexes that react due to proximity effects in the confined DLP droplets. Furthermore, our findings provide insight into CaCO 3 nucleation that is mediated by liquid–liquid phase separation, advancing the ability to direct carbonate mineralization and elucidating an often-proposed complex pathway of biomineralization.

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

Impact of surface hydrophilicity on the ordering and transport properties of bicontinuous microemulsions

Microemulsions (MEs) have many industrial applications, where recent developments have shown that MEs can be utilized for electrochemical applications, including potentially in redox flow batteries. However, understanding the structure and dynamics of these systems, including at a surface, is needed to direct and rationally control their electrochemical behavior. While bulk solution measurements have provided insight into their structure, their assembly at an interface also impacts the electron (to the electrode) and ion (across the surfactant) charge transfer processes in the system. To address this shortcoming, neutron reflectivity experiments and molecular simulations have been completed that document the near surface structure of a series of deuterated water (D 2 O)/polysorbate-20/toluene MEs on hydrophilic and amphiphilic surfaces. These results show that the microemulsions form complex layered structures near a hard electrode surface, where most layers are not purely one component. Decreasing the D 2 O concentration in the ME increases the number of and purity of the layers established on the solid surface. These lamellar-type layers transition from the surface to the bulk microemulsion as a series of mixed layers (i.e., containing oil, water, and surfactant) that are consistent with perforated lamellae. Additionally, these mixed lamellae appear to become more perforated with oil and water pathways on an amphiphilic surface. Furthermore, the purity and thickness of these layers will influence the accessibility of an electrode by redox active species, as well as ion transport required to satisfy the electroneutrality condition.

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

Near-cryogenic direct air capture using adsorbents

Direct air capture (DAC) of CO 2 is a key component in the portfolio of negative emissions technologies for mitigating global warming. However, even with the most potent amine sorbents, large-scale DAC deployment remains limited by high energy and capital costs. Recently, adsorbents relying on weak interactions with CO 2 have emerged as a potential alternative, thanks to their rapid adsorption kinetics and superior long-term stability, particularly under sub-ambient conditions (∼253 K). Despite these advantages, their use is hindered by the need for a water-removal process, location-specific constraints, and insufficient working capacity even in cold climates. In this study, we hypothesized that further reducing the adsorption temperature to a near-cryogenic range (160–220 K) could enable cost-effective DAC by utilizing the full potential of physisorbents. We primarily consider integrating DAC with a relatively untapped source of cold energy—liquified natural gas (LNG) regasification—to perform near-cryogenic DAC. From large-scale molecular simulations, Zeolite 13X and CALF-20 were identified as promising candidates. These materials were subsequently examined through experiments, including breakthrough analyses at 195 K. Their high CO 2 sorption capacity (4.5–5.5 mmol g −1 ), combined with a low desorption enthalpy and robust long-term stability, led to a threefold reduction in the levelized cost of capture (down to 68.2 USD per tonne CO 2 ). Estimates of the global LNG regasification resource suggest that LNG–DAC coupling could potentially enable the capture of 103–142 megatonnes of CO 2 annually as of 2050.

Kim, Seo-Yul [Georgia Institute of Technology, Atl↗
DOE OSTI2025

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗
DOE OSTI2024

On the engulfment of antifreeze proteins by ice

Antifreeze proteins (AFPs) are remarkable biomolecules that suppress ice formation at trace concentrations. To inhibit ice growth, AFPs must not only bind to ice crystals, but also resist engulfment by ice. The highest supercooling, ΔT*, for which AFPs are able to resist engulfment is widely believed to scale as the inverse of the separation, L, between bound AFPs, whereas its dependence on the molecular characteristics of the AFP remains poorly understood. By using specialized molecular simulations and interfacial thermodynamics, here, we show that in contrast with conventional wisdom, ΔT* scales as L –2 and not as L –1 . We further show that ΔT* is proportional to AFP size and that diverse naturally occurring AFPs are optimal at resisting engulfment by ice. By facilitating the development of AFP structure–function relationships, we hope that our findings will pave the way for the rational design of AFPs.

Science & Technology - Other Topics↗
DOE OSTI2024

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

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

Why Seeding Works When Nucleation Barriers Vanish

This paper explains why seeds can speed up crystallization even when the formation of new crystal nuclei is not the slow step. In zeolite synthesis, local crystalline order can form relatively early, but the material may still appear poorly crystalline by X-ray diffraction because many small crystallites are misoriented and lack long-range coherence. Using coarse-grained molecular simulations, the study showed that seeds help solve this coherence problem. Rather than primarily lowering a nucleation barrier, seeds impose a common orientation on nearby growing crystallites, allowing them to merge into larger, better aligned domains. This accelerates the development of X-ray-detectable crystallinity even if the total amount of locally ordered material has not changed much. The work changes how induction times in zeolite synthesis should be interpreted: the delay before X-ray crystallinity appears can reflect the time needed to build long-range structural coherence, not simply the time needed to form crystalline material. This provides a new explanation for why seeding is effective in growth-limited crystallization and offers design principles for controlling crystallization rate and crystal quality.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Derivation of scaled surface reflectances from AVIRIS data

A method for retrieving 'scaled surface reflectances' assuming horizontal surfaces having Lambertian reflectances from spectral data collected by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is presented here. In this method, the integrated water vapor amount on a pixel by pixel basis is derived from the 0.94 micron and 1.14 micron water vapor absorption features. The transmission spectra of H2O, CO2, O3, N2O, CO, CH4, and O2 in the 0.4-2.5 micron region are simulated. The scattering effect due to atmospheric molecules and aerosols is modeled with the 5S computer code. The AVIRIS radiances are divided by solar irradiances above the atmosphere to obtain the apparent reflectances. The scaled surface reflectances are derived from the apparent reflectances using the simulated atmospheric gaseous transmittances and the simulated molecular and aerosol scattering data. The scaled surface reflectances differ from the real surface reflectances by a multiplicative factor. In order to convert the scaled surface reflectances into real surface reflectances, the slopes and aspects of the surfaces must be known.

Gao, Bo-Cai↗
NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗
NASA NTRS1997

Toroidal Single Wall Carbon Nanotubes in Fullerene Crop Circles

We investigate energetics and structure of circular and polygonal single wall carbon nanotubes (SWNTs) using large scale molecular simulations on NAS SP2, motivated by their unusual electronic and magnetic properties. The circular tori are formed by bending tube (no net whereas the polygonal tori are constructed by turning the joint of two tubes of (n, n), (n+1, n-1) and (n+2, n-2) with topological pentagon-heptagon defect, in which n =5, 8 and 10. The strain energy of circular tori relative to straight tube decreases by I/D(sup 2) where D is torus diameter. As D increases, these tori change from buckling to an energetically stable state. The stable tori are perfect circular in both toroidal and tubular geometry with strain less than 0. 03 eV/atom when D greater than 10, 20 and 40 nm for torus (5,5), (8,8) and (10, 10). Polygonal tori, whose strain is proportional to the number of defects and I/D are energetically stable even for D less than 10 nm. However, their strain is higher than that of perfect circular tori. In addition, the local maximum strain of polygonal tori is much higher than that of perfect circular tori. It is approx. 0.03 eV/atom or less for perfect circular torus (5,5), but 0.13 and 0.21 eV/atom for polygonal tori (6,4)/(5,5) and (7,3)/(5,5). Therefore, we conclude that the circular tori with no topological defects are more energetically stable and kinetically accessible than the polygonal tori containing the pentagon-heptagon defects for the laser-grown SWNTs and Fullerene crop circles.

Han, Jie↗
NASA NTRS2003

CO2 Acquisition Membrane (CAM)

The objective of CAM is to develop, test, and analyze thin film membrane materials for separation and purification of carbon dioxide (CO2) from mixtures of gases, such as those found in the Martian atmosphere. The membranes are targeted toward In Situ Resource Utilization (ISRU) applications that will operate in extraterrestrial environments and support future unmanned and human space missions. A primary application is the Sabatier Electrolysis process that uses Mars atmosphere CO2 as raw material for producing water, oxygen, and methane for rocket fuel and habitat support. Other applications include use as an inlet filter to collect and concentrate Mars atmospheric argon and nitrogen gases for habitat pressurization, and to remove CO2 from breathing gases in Closed Environment Life Support Systems (CELSS). CAM membrane materials include crystalline faujasite (FAU) zeolite and rubbery polymers such as silicone rubber (PDMS) that have been shown in the literature and via molecular simulation to favor adsorption and permeation of CO2 over nitrogen and argon. Pure gas permeation tests using commercial PDMS membranes have shown that both CO2 permeance and the separation factor relative to other gases increase as the temperature decreases, and low (Delta)P(Sub CO2) favors higher separation factors. The ideal CO2/N2 separation factor increases from 7.5 to 17.5 as temperature decreases from 22 C to -30 C. For gas mixtures containing CO2, N2, and Ar, plasticization decreased the separation factors from 4.5 to 6 over the same temperature range. We currently synthesize and test our own Na(+) FAU zeolite membranes using standard formulations and secondary growth methods on porous alumina. Preliminary tests with a Na(+) FAU membrane at 22 C show a He/SF6 ideal separation factor of 62, exceeding the Knudsen diffusion selectivity by an order of magnitude. This shows that the membrane is relatively free from large defects and associated non-selective (viscous flow) transport mechanisms. The Membrane Test Facility (MTF) has been developed to measure membrane permeance over a wide range of temperature and pressure. The facility uses two volume compartments separated by the membrane that are instrumented to measure temperature, delta pressure across the membrane, and gas composition. A thermal shroud supports and encloses the membrane, and provides temperature control. Methods were developed to determine membrane permeance using the first order decay of the pressure difference between the sealed compartments, using the total pressure for pure gases, and partial pressure of each species in gas mixtures. The technique provides an end-to-end measurement of gas permeance that includes concentration polarization effects. Experiments have shown that in addition to membrane permeance properties, the geometry and design of associated structures play an important role in how membrane systems will function on Mars.

Mason, Larry W.↗
NASA NTRS2004

Novel Microstructures for Polymer-Liquid Crystal Composite Materials

There are a number of interface-dominated composite materials that contain a liquid crystalline (LC) phase in intimate contact with an isotropic phase. For example, polymer- dispersed liquid crystals, used in the fabrication of windows with switchable transparency, consist of micron size LC droplets dispersed in an isotropic polymer matrix. Many other types of liquid crystal composite materials can be envisioned that might have outstanding optical properties that could be exploited in novel chemical sensors, optical switches, and computer displays. This research project was based on the premise that many of these potentially useful LC composite materials can only be fabricated under microgravity conditions where gravity driven flows are absent. In the ground-based research described below, we have focused on a new class of LC composites that we call thermotropic- lyotropic liquid crystal systems (TLLCs). TLLCs consist of nanosize droplets of water dispersed in an LC matrix, with surfactants at the interface that stabilize the structure. By varying the type of surfactant one can access almost an infinite variety of unusual LC composite microstructures. Due to the importance of the interface in these types of systems, we have also developed molecular simulation models for liquid crystals at interfaces, and made some of the first measurements of the interfacial tension between liquid crystals and water.

Magda, Jules J.↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspect of out understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to to simulate these process in numerical models of chemistry and transport. These models depend on laboratory-measured kinetic reaction rates and photlysis cross section to simulate molecular interactions. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluation are applied in random combinations. We determine the key reaction and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspects of our understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to simulate these processes in numerical models of chemistry and transport. The fidelity of the models is assessed in comparison with a wide range of observations. These models depend on laboratory-measured kinetic reaction rates and photolysis cross sections to simulate molecular interactions. A typical stratospheric chemistry mechanism has on the order of 50- 100 species undergoing over a hundred intermolecular reactions and several tens of photolysis reactions. The rates of all of these reactions are subject to uncertainty, some substantial. Given the complexity of the models, however, it is difficult to quantify uncertainties in many aspects of system. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluations are applied in random combinations. We determine the key reactions and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2018

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗
NASA NTRSDate not supplied

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗