Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural 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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Thermodynamic Stability of Molybdenum Oxycarbides Formed from Orthorhombic Mo 2 C in Oxygen-Rich Environments

Molybdenum carbide (Mo 2 C) nanoparticles and thin films are particularly suitable catalysts for catalytic fast pyrolysis (CFP) as they are effective for deoxygenation and can catalyze certain reactions that typically occur on noble metals. Oxygen deposited during deoxygenation reactions may alter the carbide structure, leading to the formation of oxycarbides, which can determine changes in catalytic activity or selectivity. Despite emerging spectroscopic evidence of bulk oxycarbides, so far there have been no reports of their precise atomic structure or their relative stability with respect to orthorhombic Mo 2 C. This knowledge is essential for assessing the catalytic properties of molybdenum (oxy)carbides for CFP. In this article, we use density functional theory (DFT) calculations to (a) describe the thermodynamic stability of surface and subsurface configurations of oxygen and carbon atoms for a commonly studied Mo-terminated surface of orthorhombic Mo 2 C and (b) determine atomic structures for oxycarbides with a Mo:C ratio of 2:1. The surface calculations suggest that oxygen atoms are not stable under the top Mo layer of the Mo 2 C(100) surface. Coupling DFT calculations with a polymorph sampling method, we determine (Mo 2 C) x O y oxycarbide structures for a wide range of oxygen compositions. Oxycarbides with lower oxygen content (y/x = 2) adopt layered structures reminiscent of the parent carbide phase, with flat Mo layers separated by layers of oxygen and carbon; for higher oxygen content, our results suggest the formation of amorphous phases, as the atomic layers lose their planarity with increasing oxygen content. We characterize the oxidation states of Mo in the oxycarbide structures determined computationally, and simulate their X-ray diffraction (XRD) patterns in order to facilitate comparisons with experiments. Our study may provide a platform for large-scale investigations of the catalytic properties of oxycarbides and their surfaces and for tailoring the catalytic properties for different desired reactions.

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

Automated Development of Molten Salt Machine Learning Potentials: Application to LiCl

The in silico modeling of molten salts is critical for emerging "carbon-free" energy applications but is inhibited by the cost of quantum mechanically treating the high polarizabilities of molten salts. Here, we integrate configurational sampling using classical force fields with active learning to automate and accelerate the generation of Gaussian approximation potentials (GAP) for molten salts. This methodology reduces the number of expensive ab initio evaluations required for training set generation to O(100), enabling the facile parametrization of a molten LiCl GAP model that exhibits a 19 000-fold speedup relative to AIMD. The developed molten LiCl GAP model is applied to sample extended spatiotemporal scales, permitting new physical insights into molten LiCl's coordination structure as well as experimentally validated predictions of structures, densities, self-diffusion constants, and ionic conductivities. The developed methodology significantly lowers the barrier to the in silico understanding and design of molten salts across the periodic table.

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

Hydride Accessibility and Reactivity in the Configurational and Stoichiometric Space of β-Ga 2 O 3 for CO 2 Hydrogenation

Understanding how surface species evolve under reaction conditions is essential for improving catalyst design for efficient CO 2 hydrogenation. This work combines systematic DFT calculations with grand canonical sampling to investigate the stability and reactivity of Ga–H species on β-Ga 2 O 3 across a range of reaction conditions. Initial DFT studies reveal that when Ga–H species are present, they facilitate formate formation via a low-barrier pathway, largely independent of the surface termination or hydrogen site. However, grand canonical sampling shows that under a broad range of reaction conditions─especially at high oxygen chemical potentials associated with high water content─Ga–H species are thermodynamically inaccessible. Furthermore, adsorbed water molecules can block reactive sites, inhibiting CO 2 activation even when hydrides are present. These findings suggest that the lack of accessible hydride species, rather than their intrinsic reactivity, could contribute to reduced catalytic performance of β-Ga 2 O 3 under more oxidizing, high-conversion conditions.

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

Solvents Polarity Governs Ion Interactions and Transport in a Solvated Room Temperature Ionic Liquid

We explore the influence of the solvent dipole moment on cation anion interactions and transport in 1-butyl-3-methyl-imidazolium bis-(trifluoromethylsulfonyl), [BMIM+][Tf2N ]. Free energy profiles derived from atomistic molecular dynamics (MD) simulations show a correlation of the cation anion separation and the equilibrium depth of the potential of mean force with the dipole moment of the solvent. Correlations of the ion diffusivity with the dipole moment and the concentration of the solvent were further demonstrated by classical MD simulations. Quasi-elastic neutron scattering experiments with deuterated solvents reveal a complex picture of nanophase separation into the ionic liquid-rich and solvent-rich phases. The experiment corroborates the trend of concentration- and dipole moment-dependent enhancement of ion mobility by the solvent, as suggested by the simulations. Despite the considerable structural complexity of ionic liquid solvent mixtures, we can rationalize and generalize the trends governing ionic transport in these complex electrolytes.

Osti, Naresh C↗
DOE OSTI2024

Vacancy-Mediated Increases in Brine–Salt Surface Energies

Salt formations have been explored for the permanent isolation of spent nuclear fuel based on their high thermal conductivity, self-healing nature, and low hydraulic permeability to brine flow. Vacancy defect concentrations in salt complicate fracture mechanics not driven by dislocation dynamics and can influence the resulting surface structure. Classical molecular dynamic simulations were used to simulate tensile testing of salt crystals (halite) with vacancy defect concentrations of up to 0.5 defects/nm 3 . Increasing defect concentrations resulted in a decrease in ultimate tensile strength and fracture surface energies, driven by increased surface roughness rather than changes in the amount of surface area. Brine–salt surface energies of the fractured surfaces were 0.22 to 0.26 J/m 2 , significantly higher than values reported for atomically flat (100) surfaces at the same brine composition. This change in surface energy increased the brine–salt dihedral angle by ~27°. The dihedral angle threshold for percolation in salt is 60°, and a 27° increase due to rough fracture surfaces identifies a reduction in porosity percolation and a decrease in salt permeability. Therefore, bedded salt and salt domes may be even more stable than those previously predicted from dihedral angle calculations.

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

Stability of Ligands on Nanoparticles Regulating the Integrity of Biological Membranes at the Nano–Lipid Interface

When nanoparticles interact with cellular or organelle membranes, the coating ligands are known to affect the integrity of the membranes, which regulate cell death and inflammation. However, the molecular mechanisms of this modulation remain unresolved. Here, we use synchrotron X-ray liquid surface scattering and molecular dynamics simulations to study interface structures between phospholipids and gold nanorods (AuNRs) coated by surfactant and polyelectrolyte. These ligands are two types of widely used surface modification with different self-assembled structures and stabilities on the surface of nanoparticles. We reveal distinct mechanisms of the ligand stability in disrupting membrane integrity. We find that cationic surfactant ligand cetyltrimethylammonium bromide (CTAB) detaches from AuNRs and inserts into phospholipids, resulting in reduced membrane thickness by compressing the phospholipids to align with the shorter ligand. Conversely, the cationic polyelectrolyte ligand poly(diallyldimethylammonium chloride) (PDC) is more stable on AuNRs; although it adsorbs onto the membrane, it does not cause much impairment. The distinct coating ligand interactions with phospholipids are further verified by cellular responses including impaired lysosomal membranes and triggered inflammatory effects in macrophages. Together, the quantitative analysis of interface structures elucidates key bio-nano interactions and highlights the importance of surface ligand stability for safety and rational design of nanoparticles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

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

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cyclodextrin-Derived Porous Liquids Enabled by In Situ Solvation Shell Formation

Porous liquids (PLs) represent a unique platform for molecular separations by combining permanent porosity with liquid-phase mobility. However, it remains a formidable challenge to construct and stabilize PLs with sub-5 Å pores using readily available porous host and liquid media. Here, we report the construction of cyclodextrin (CD)-derived PLs enabled by in situ solvation shell formation. The acid–base neutralization reaction between CD and an organic base was leveraged to generate a thin ionic solvation shell around the CD host, effectively liquefying CD and preventing its segregation in the liquid base medium while preserving accessible molecular-scale cavities. Spectroscopic analysis, neutron scattering, density functional theory calculations, and molecular dynamics simulations collectively confirm the structural evolution and existence of abundant internal porosity in PLs. The unique architectures of CD-derived PLs enable highly selective encapsulation of fluorinated alkanes and significantly enhanced uptake of inert gases. This facile and generalizable strategy enables construction of high-quality PLs with engineered ultramicroporosity to facilitate molecular separations.

Cavities↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗
DOE OSTI2017

Crystal structure of a multi-domain human smoothened receptor in complex with a super stabilizing ligand

Here, the Smoothened receptor (SMO) belongs to the Class Frizzled of the G protein-coupled receptor (GPCR) superfamily, constituting a key component of the Hedgehog signalling pathway. Here we report the crystal structure of the multi-domain human SMO, bound and stabilized by a designed tool ligand TC114, using an X-ray free-electron laser source at 2.9 Å. The structure reveals a precise arrangement of three distinct domains: a seven-transmembrane helices domain (TMD), a hinge domain (HD) and an intact extracellular cysteine-rich domain (CRD). This architecture enables allosteric interactions between the domains that are important for ligand recognition and receptor activation. By combining the structural data, molecular dynamics simulation, and hydrogen-deuterium-exchange analysis, we demonstrate that transmembrane helix VI, extracellular loop 3 and the HD play a central role in transmitting the signal employing a unique GPCR activation mechanism, distinct from other multi-domain GPCRs.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Non-conventional mechanism of ferroelectric fatigue via cation migration

The unique properties of ferroelectric materials enable a plethora of applications, which are hindered by the phenomenon known as ferroelectric fatigue that leads to the degradation of ferroelectric properties with polarization cycling. Multiple microscopic models explaining fatigue have been suggested; however, the chemical origins remain poorly understood. Here, we utilize multimodal chemical imaging that combines atomic force microscopy with time-of-flight secondary mass spectrometry to explore the chemical phenomena associated with fatigue in PbZr 0.2 Ti 0.8 O 3 (PZT) thin films. Investigations reveal that the degradation of ferroelectric properties is correlated with a local chemical change and migration of electrode ions into the PZT structure. Density functional theory simulations support the experimental results and demonstrate stable doping of the thin surface PZT layer with copper ions, leading to a decrease in the spontaneous polarization. Overall, the performed research allows for the observation and understanding of the chemical phenomena associated with polarization cycling and their effects on ferroelectric functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2021

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation

Abstract Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized 1 B 2u (ππ*) (S 2 ) and 1 B 3u (nπ*) (S 1 ) states, the participation of the optically dark 1 A u (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite 1 A u (nπ*) and 1 B 3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The 1 A u (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

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

Conformational restriction shapes the inhibition of a multidrug efflux adaptor protein

Abstract Membrane efflux pumps play a major role in bacterial multidrug resistance. The tripartite multidrug efflux pump system from Escherichia coli , AcrAB-TolC, is a target for inhibition to lessen resistance development and restore antibiotic efficacy, with homologs in other ESKAPE pathogens. Here, we rationalize a mechanism of inhibition against the periplasmic adaptor protein, AcrA, using a combination of hydrogen/deuterium exchange mass spectrometry, cellular efflux assays, and molecular dynamics simulations. We define the structural dynamics of AcrA and find that an inhibitor can inflict long-range stabilisation across all four of its domains, whereas an interacting efflux substrate has minimal effect. Our results support a model where an inhibitor forms a molecular wedge within a cleft between the lipoyl and αβ barrel domains of AcrA, diminishing its conformational transmission of drug-evoked signals from AcrB to TolC. This work provides molecular insights into multidrug adaptor protein function which could be valuable for developing antimicrobial therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structural determinants for pH-dependent activation of a plant metacaspase

Arabidopsis thaliana metacaspase 9 (AtMC9) plays roles in clearing dead cells, forming xylem vessels, and regulating immunity and programmed cell death in plants. The protease's activation is controlled by pH levels, but the exact structural mechanism behind this has not been elucidated. In this work, high-resolution crystal structures for AtMC9 at active (pH 5.5 and pH 4.2) and inactive (pH 7.5) conditions are reported. The three structures are similar except for local conformations where their hydrogen bonding interactions with solvents are mediated through the protonation of specific titratable amino acid residues' side chains. By combining structural analysis, molecular dynamics simulations under constant pHs, and biochemical assays coupled with site-directed mutagenesis, we show that the regulation of AtMC9 activation involves multiple titratable glutamate and histidine residues across the three domains of p20, linker, and p10. Specifically, deprotonated Glu112, His193, and His208 can suppress AtMC9 proteolytic activity, while protonation of Glu255 and His307 at acidic pH may promote it. This study provides valuable insights into the pH-dependent activation of AtMC9 and could potentially lead to improving crops with enhanced immunity and controlled cell death, ultimately increasing agricultural productivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Phosphorylation toggles the SARS-CoV-2 nucleocapsid protein between two membrane-associated condensate states

Abstract The Nucleocapsid protein (N) of SARS-CoV-2 plays a critical role in the viral lifecycle by regulating RNA replication and by packaging the viral genome. N and RNA phase separate to form condensates that may be important for these functions. Both functions occur at membrane surfaces, but how N toggles between these two membrane-associated functional states is unclear. Here, we reveal that phosphorylation switches how N condensates interact with membranes, in part by modulating condensate material properties. Our studies also show that phosphorylation alters N’s interaction with viral membrane proteins. We gain mechanistic insight through structural analysis and molecular simulations, which suggest phosphorylation induces a conformational change in N that softens condensate material properties. Together, our findings identify membrane association as a key feature of N condensates and provide mechanistic insights into the regulatory role of phosphorylation. Understanding this mechanism suggests potential therapeutic targets for COVID infection.

Science & Technology - Other Topics↗
DOE OSTI2021

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Abstract Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Computational modeling predicts the stability of both Pd + and Pd2 + ion-exchanged into H-CHA

Passive NO x adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NO x adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework before the introduction of NO x . Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next–next–nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd + instead of Pd 2+ . The stability of Pd + in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-Vis, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗