Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Partial Solvation of Lithium Ions Enhances Conductivity in a Nanophase-Separated Polymer Electrolyte

In this study, we demonstrate that a multiblock lithium-ion-conducting polymer can be swollen with ethylene carbonate solvent to increase the conductivity relative to the dry polymer material by nearly 4 orders of magnitude. This increase is due to the partial solvation of lithium ions by ethylene carbonate, which leads to Li + diffusion along the solvent–polymer interface. This differs from the vehicular transport mechanism for lithium ions in pure solvent. We use a combination of broadband dielectric spectroscopy, X-ray scattering, and all-atom molecular dynamics simulations to probe the effect of the solvent on the polymer morphology and to elucidate the mechanism of lithium ion transport.

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

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

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

Effect of Solvent on the Interaction of Lignin with a Zeolite Nanosheet in the Condensed Phase

Lignin is the essential building block of lignocellulosic biomass, an excellent renewable source of different aromatic monomers for the polymer and biofuel industry. The depolymerization of lignin into value-added chemicals and fuels through the catalytic process poses a significant challenge due to the complex structure of lignin. Understanding lignin’s conformational diversity and dynamics in the liquid phase is crucial to designing an effective depolymerization process. Here, we conducted all-atom molecular dynamics simulations to understand the conformation and dynamics of softwood lignin on the all-silica zeolite nanosheet based on the MFI topology in a binary mixture of water–methanol at three different molar compositions (0%, 50%, and 100% methanol). We observed that the methanol–surface interaction is stronger than the water–surface interaction, and methanol readily diffused into the MFI core. Lignin surface contacts decrease with increasing methanol composition due to higher solubility and dynamics. Lignin dynamics on the surface in neat water is an order of magnitude smaller than methanol. Here, we also found that lignin adopts a slightly extended conformation when it stays on the surface than in the bulk solution phase for the pure water case, whereas for pure methanol and the binary solution structures are statistically similar.

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

Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions

Here, we perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

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

Side-Chain and Ring-Size Effects on Permeability in Artificial Water Channels

Artificial water channels (AWCs) have emerged as a promising framework for stable water permeation, with water transport rates comparable to aquaporins (3.4–40.3 × 10 8 H 2 O/channel/s). In this study, we probe the influence of ring-size and side-chain length on the water permeability observed within a class of AWCs termed ligand-appended pillar[n]arenes (LAPs) that have an adjustable ring-size (m) and side-chain length (n). Through all-atom molecular dynamics simulations, we calculate the permeability of these channels using the collective diffusion model and find their permeabilities. We characterize the mechanistic influence of pillar[n]arene ring-size and side-chain length on the channel water permeability by analyzing the characteristics of the internal permeating water-wire and the surrounding channel structure. We observe that water permeability decreases as a function of increasing ring-size due to increases in hydrophilic contacts between the permeating water-wire and the oxygen groups on the channel wall. Further, we observe an increase in water permeability as a function of side-chain length due to increased partitioning of the channel terminal groups into the hydrophilic blocks of the surrounding bilayer. For the LAP6 channel, with increase in side-chain length, the distance between terminal groups increases and leads to an increase in pore size, thereby enhancing water permeability. In the case of LAP5, as side-chain length increases, the channel displays a compensatory effect between tilt and bend angle due to the flexible side-chains. Such flexibility leads to higher terminal group partitioning in the hydrophilic blocks of the bilayer and extends the permeating water-wire. Furthermore, this increase in water-wire length and hydrophilic block access overcomes the nonmonotonic pore size trend in pillar[5]arene channels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-Assembled Oligomers Facilitate Amino Acid-Driven CO 2 Capture at the Air–Aqueous Interface

Direct air capture of CO 2 using amino acid absorbents, such as glycine or sarcosine, is constrained by the relatively slow mass transfer of CO 2 through the air–aqueous interface. Our recent study showed a marked improvement in CO 2 capture by introducing CO 2 -permeable oligo-dimethylsiloxane (ODMS-MIM + ) oligomers with cationic (imidazolium, MIM + ) headgroups. Here, in this work, we have employed all-atom molecular dynamics simulations in combination with subensemble analysis using network theory to provide a detailed molecular picture of the behavior of CO 2 and the glycinate anions (Gly – ) at the ODMS-MIM + decorated air–aqueous interfaces. We show that the cationic head groups of the surfactants enhance the concentration and lifetime of Gly – in the interfacial region, while ODMS tails promote the physisorption of CO 2 in the interfacial region. Together, these two factors increase the effective region of contact and the probability of interactions between CO 2 and Gly – compared to that of the pure air–aqueous interface. The fundamental insights gained in this work establish essential foundations for developing hybrid systems with oligomer-decorated interfaces to maximize the overall CO 2 capture rates.

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

Subsurface Characteristics of Metal-Halide Perovskites Polished by an Argon Ion Beam

Focused ion beam (FIB) techniques have been frequently used to section metal-halide perovskites for microstructural investigations. However, the ion beams directly irradiating the sample surface may alter its properties far different from those of pristine, potentially leading to modified deterioration mechanisms under aging stressors. Here, we combine complementary approaches to measure the subsurface characteristics of polished perovskites and identify the chemical species responsible for the measured properties. Analysis of the experimental results in conjunction with Monte Carlo simulations indicates that atomic displacements and local heating occur in the first ≈15 nm of the subsurface of methylammonium lead iodide (MAPbI 3 ) by glancing-angle Ar + beam irradiation (4 kV at 3°). The lead-rich, iodine-deficient surface promotes rapid phase segregation under thermal aging conditions. On the other hand, despite the subsurface modification, our experiments confirm that the rest of the MAPbI 3 bulk retains the material integrity. Our observation supports that polished perovskites could serve in studying the properties of bulk or buried junctions far away from the altered subsurface with care.

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

A Hotspot’s Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material’s microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. Furthermore, the origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elucidating the Role of Hydrogen Bonding in the Optical Spectroscopy of the Solvated Green Fluorescent Protein Chromophore: Using Machine Learning to Establish the Importance of High-Level Electronic Structure

Hydrogen bonding interactions with chromophores in chemical and biological environments play a key role in determining their electronic absorption and relaxation processes, which are manifested in their linear and multidimensional optical spectra. For chromophores in the condensed phase, the large number of atoms needed to simulate the environment has traditionally prohibited the use of high-level excited-state electronic structure methods. By leveraging transfer learning, we show how to construct machine-learned models to accurately predict the high-level excitation energies of a chromophore in solution from only 400 high-level calculations. Here, we show that when the electronic excitations of the green fluorescent protein chromophore in water are treated using EOM-CCSD embedded in a DFT description of the solvent the optical spectrum is correctly captured and that this improvement arises from correctly treating the coupling of the electronic transition to electric fields, which leads to a larger response upon hydrogen bonding between the chromophore and water.

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

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Specific Ion Effects in Lanthanide–Amphiphile Structures at the Air–Water Interface and Their Implications for Selective Separation

The use of surfactants to attract dissolved ions to water surfaces and interfaces is an essential step in both solvent-based and solvent-free separation processes. Here, we have studied the interactions of lanthanide ions in the aqueous subphase with monolayers of dihexadecyl phosphate at air–water interfaces. With heavier lanthanides (atomic number Z ≥ 65) in the subphase, the floating layer can be compressed to an area/molecule of about half the molecular cross section, indicating bilayer formation. X-ray fluorescence and reflectivity data support this conclusion. In the presence of lighter lanthanides (Z < 65), only monolayers are observed. Subphase-concentration-dependent studies using Er 3+ (heavier) and Nd 3+ (lighter) lanthanides show a stepwise progression, with ions attaching to the monolayer only when the solution concentration is >3 × 10 –7 M. Above ~10 –5 M, bilayers form but only in the presence of the heavier lanthanide. Grazing incidence X-ray diffraction shows evidence of lateral ion–ion correlations in the bilayer structure but not in monolayers. Explicit solvent all-atom molecular dynamics simulations confirm the elevated ion–ion correlation in the bilayer system. This bilayer structure isolates heavier lanthanides but not lighter lanthanides from an aqueous solution and is therefore a potential mechanism for the selective separation of heavier lanthanides.

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

Melatonin-Induced Modulation of Cholesterol-Enriched Model Neuronal Membranes

Melatonin, a hormone primarily produced by the brain’s pineal gland, not only regulates circadian rhythms, but also influences the structural and biophysical properties of neuronal membranes. Its amphiphilic nature enables direct incorporation into lipid bilayers and preferential interactions with cholesterol-rich lipid rafts, critical hubs for cellular signaling and membrane organization. Despite increasing recognition of its membrane activity, the molecular basis of melatonin’s interactions with coexisting liquid-ordered (L o ) and liquid-disordered (L d ) phases remains unclear. Here, in this study, we combine small-angle neutron scattering (SANS) and all-atom molecular dynamics simulations to examine model neuronal membranes composed of DSPC, DOPC, POPC, and cholesterol. Our results show that melatonin preserves domain morphology while adopting distinct orientations within the bilayer and at the membrane interface, allowing both lateral and transmembrane bridging across lipid phases. These findings establish the molecular underpinnings of melatonin’s modulation of membrane heterogeneity and provide strong support for its receptor-independent actions.

atomistic simulations↗
DOE OSTI2021

Overscreening, Coion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na + ) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl – ) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na + ) ensures a water transport in the direction of motion of the counterions. In this work, the net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.

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

Structural insights into the modulation of coronavirus spike tilting and infectivity by hinge glycans

Abstract Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectrometry. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a role of hinge glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays implicated involvement of N1242-glyan in virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds’ intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

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

Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABABR1a

A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABA B R1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABA B R1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = -1.91 ± 0.66 kcal mol -1 ), indicating a potentially enhanced GABA B R1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABA B R1a-specific therapeutics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗