Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,153 records · Page 64

Two-component dynamics in supercritical $\text {CO}_2$ from inelastic X-ray scattering

Supercritical fluids are characterized by unique thermodynamic properties. One of these properties is the existence of two-component dynamics that is associated with distinct low-frequency and high-frequency vibrational responses of the fluid. However, the origin of this behavior remains unknown. By combining inelastic X-ray scattering and molecular dynamics simulations, we show that this behavior can be connected to density heterogeneities arising from molecular clusters. Analyses of measurements and molecular trajectories suggest that the two-component dynamics emerges due to distinct momentum fluctuations of clustered and unbound molecules. This connection between clusters and two-component dynamics highlights the importance of molecular-structural heterogeneities in supercritical fluids, colloids, and condensed-matter systems.

biophysics↗

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Windowed Fourier transform and cross-correlation algorithms for molecular tagging velocimetry

In this study, simulated and experimental molecular tagging velocimetry (MTV) images have been analyzed with a technique commonly used to process grid images on surfaces, the windowed Fourier transform with local spectrum analysis (WFT-LSA). A systematic synthetic image study of the modulation transfer function (MTF) and error tendencies of the WFT-LSA was performed and compared with a PIV-style cross-correlation algorithm to see if advanced strategies such as iterative image deformation can improve analysis of gridded images with high noise levels. Testing of single-pass algorithms showed that in typical MTV images, the WFT-LSA yields significantly lower bias errors than cross-correlation (CC) at displacements greater than 1 pixel but slightly higher random error at all displacements and image conditions. Analysis of the MTF shows that CC provided better resolution of spatial fluctuations than the WFT-LSA in many combinations of grid size and interrogation window. Tests of image deformation algorithms showed that the gap in performance between CC and WFT-LSA is maintained even as both methods improve. Additionally, WFT-LSA and CC methods are applied to real MTV experiments in high-speed gas jet flows. A preliminary analysis of phosphorescence lifetime provided by acetone vapor excited at 266 nm is used for assessing the required gas speeds for making MTV application feasible. The application of WFT-LSA to real MTV images demonstrates the ability of the algorithm to handle further real-world effects that could not be considered in the synthetic image analysis, like reduced signal-to-noise ratio and non-uniform intensity of the tagging grid across the image introduced by the actual laser beam energy distribution. With experimental images, CC is more accurate with shear flows but less robust to high noise levels than WFT-LSA, as predicted by the synthetic image analysis.

42 ENGINEERING↗

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗

First principles molecular dynamics maps out complete mineral surface acidity landscape

In “Interfacial structures and acidity constants (pKa) of goethite from first principles molecular dynamics simulations,” authors Y. Zhang, X. Lui, J. Cheng, and X. Lu apply First Principles molecular dynamics (FPMD, also called Density Functional Theory MD, DFT/MD, or ab initio MD, AIMD), to evaluable the complete set of acidity constants (pK a ) of the hydroxyl groups on the most prominent goethite crystal facets. The pK a of these OH and OH$^+_2$ groups are compared with available data from the multisite complexation (MUSIC) model traditionally used to estimate pK a on mineral surfaces. The authors have presented eloquent rationale for the importance and implications of understanding goethite acidity constants in room temperature geochemistry settings. In this paper, I focus on the computational aspects, the strengths of FPMD, and its possibilities.

58 GEOSCIENCES↗

Proton Transport

The transport of protons across membranes is an essential process for both bioenergetics of modern cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATP), synthesized from adenosine diphosphate. ATP, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this Study was: how the same process can be accomplished with the aid of similar but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC) which is a good model of the biological membranes forming cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M$_2$ protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M$_2$ protein is 97 amino acids in length, but a fragment 25 amino acids long. which contains a transmembrane domain of 19 amino acids flanked by three amino acids on each side. is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules. and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small ions, such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel with the histidine residues in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel loses its structural integrity and one helix moves away from the remaining three. Additional information is contained within the original extended abstract.

Pohorille, Andrew↗

The M2 Proton Channel of Influenza Virus: How Does It Work?

The transport of protons across membranes is an essential process for both bioenergetics of modem cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATR), synthesized from adenosine diphosphate. ATR, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this study was: how the same process can be accomplished with the aid of similar, but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC), which is a good model of the biological membranes focusing cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M2 protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M2 protein is 97 amino acids in length, but a fragment 25 amino acids long, which contains a transmembrane domain of 19 amino acids flanked by 3 amino acids on each side, is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules, and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional, important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel, in which the histidine residues were in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel looses its structural integrity and one helix moves away from the remaining three. This result indicates that such a mechanism of proton transport is unlikely. The simulations revealed that translocation along a network of water molecules in the channel and tautomerization of the histidine residues in the M2 proteins in the most likely mechanism of proton transport. The results not only explain how a remarkably simple protein system can efficiently aid in the formation of proton gradients across cell walls, but also suggest how this system can be genetically re-engineered to become a directional, reversible proton pump. Such a pump can provide energy to laboratory-built models of simple cellular systems. If they were successfully constructed it would greatly advance our understanding of the beginnings of life and find important applications in medicine and pharmacology.

Pohorille, Andrew↗

Electrokinetic comparison of CaF 2 and AgCl from a MDS perspective

CaF 2 has limited electrokinetic sensitivity to fluoride ions in solution and the positive CaF 2 surface charge is not reversed at high fluoride ion concentration. Similar insensitivity of surface charge to high fluoride ion concentration has been observed for other alkali and alkaline earth fluoride salts. On the contrary, chloride ions act as potential determining ions for AgCl and reverse the surface charge, as expected. FTIR-IRS results indicate that the AgCl/water interface has a water structure-breaking tendency, whereas the CaF 2 /water interface has a structure-making tendency. Here, it is speculated that the strong hydrogen bonding between fluoride ions and water prevents the accommodation of excess fluoride ions at surface lattice sites of fluorite, and therefore explains the low sensitivity of the fluorite surface charge sign and magnitude to fluoride ion concentration. However, this speculation has not been examined at the molecular level and is evaluated by molecular dynamics simulations (MDS) in the research presented in this paper. Higher hydrogen bond density at the CaF 2 mineral/water interface confirms that the water structure has stronger hydrogen bonding at both the CaF 2 mineral surface and with the F - ion in solution. In this way, the interaction of the F - ion with the CaF 2 surface is inhibited.

36 MATERIALS SCIENCE↗

Rational design of antimicrobial peptides targeting Gram-negative bacteria

Membrane-targeting host antimicrobial peptides (AMPs) can kill or inhibit the growth of Gram-negative bacteria. However, the evolution of resistance among microbes poses a substantial barrier to the long-term utility of the host AMPs. Combining experiment and molecular dynamics simulations, we demonstrate that terminal carboxyl capping enhances both membrane insertion and antibacterial activity of an AMP called P1. Furthermore, we show that a bacterial strain with evolved resistance to this peptide becomes susceptible to P1 variants with either backbone capping or lysine-to-arginine substitutions. Our results suggest that cocktails of closely related AMPs may be useful in overcoming evolved resistance.

59 BASIC BIOLOGICAL SCIENCES↗

Probing DNA structural heterogeneity by identifying conformational subensembles of a bicovalently bound cyanine dye

DNA is a re-configurable, biological information-storage unit, and much remains to be learned about its heterogeneous structural dynamics. For example, while it is known that molecular dyes templated onto DNA exhibit increased photostability, the mechanism by which the structural dynamics of DNA affect the dye photophysics remains unknown. Here, we use femtosecond, two-dimensional electronic spectroscopy measurements of a cyanine dye, Cy5, to probe local conformations in samples of single-stranded DNA (ssDNA–Cy5), double-stranded DNA (dsDNA–Cy5), and Holliday junction DNA (HJ–DNA–Cy5). A line shape analysis of the 2D spectra reveals a strong excitation–emission correlation present in only the dsDNA–Cy5 complex, which is a signature of inhomogeneous broadening. Molecular dynamics simulations support the conclusion that this inhomogeneous broadening arises from a nearly degenerate conformer found only in the dsDNA–Cy5 complex. These insights will support future studies on DNA’s structural heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Post-Marcus Theory and Simulation of Interfacial Charge Transfer Dynamics in Organic Semiconducting Materials (Final Report)

The research program developed, validated and applies predictive computational protocols for calculating charge transfer (CT) rates in complex molecular systems, including molecular dyads and triads in liquid solution and solid-state organic semiconducting (OSC) materials. We have established a transformative computational scheme that goes beyond widely used simplifications, to achieve realistic descriptions of CT processes. The approach properly addresses the contribution of molecular environment at ambient conditions to CT processes. Our approach achieves unique insight on investigated CT processes in relevant experimental efforts. The collaborative team included three principal investigators (PIs), with complimentary expertise in classical molecular dynamics simulations and data science (Cheung), state-of-the-art electronic structure calculations (Dunietz), and cutting-edge theory and simulation techniques for modeling energy, charge and coherence transfer dynamics in molecular systems (Geva).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Environmentally Benign Electrolyte for High Energy Lithium Metal Batteries

In this work, a hybrid electrolyte comprising a high content of H 2 O for a lithium metal cell is reported. At high LiFSI salt concentration, the N-methyl-N-propyl-piperidinium bis(fluorosulfonyl) imide (PMpipFSI) electrolyte can tolerate up to 1 M H2O addition without sacrificing its redox stability on both lithium nickel manganese cobalt oxide (NMC) cathode and lithium metal anode. Molecular dynamics simulations revealed the underpinned mechanism that, at high salt concentrations, H 2 O molecules are embedded in the Li + , PMpip + , and FSI - bulk as a structural material with a strong solvation with Li + and are orderly distributed at the surface of both electrodes. This electrolyte eliminates the critical moisture controls required for the state-of-the-art (SOA) carbonate/LiPF 6 electrolyte, electrode, separator and cell assembly, thus significantly reducing the cost of the mass production of the batteries.

molecular dynamics simulation↗

Molecular Dynamics Investigation of Lipid-Specific Interactions with a Fusion Peptide

The HIV-1 fusion peptide, which is a short hydrophobic peptide from the gp41 coat glycoprotein that participates in the infection of a cell, interacts with model lipid bilayer membranes in a concentration-dependent manner. The interaction of the peptide with the bilayer also strongly depends on the lipid composition. Here, molecular dynamics simulations were performed to investigate lipid-specific interactions that arise shortly after the binding of a less-fusogenic variant of the HIV-1 fusion peptide to a lipid bilayer composed of a mixture of dimyristoyl phosphatidylcholine and dimyristoyl phosphatidylglycerol. The impact of peptide concentration was also studied. An improved understanding was gained of the lipid-specific interactions experienced by the FP. New insight was also gained into how the peptide concentration changes these interactions.

60 APPLIED LIFE SCIENCES↗

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN–nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

60 APPLIED LIFE SCIENCES↗

Direct Correlation of the Salt-Reduced Diffusivities of Organic Solvents with the Solvent’s Mole Fraction

We report that lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in organic solvents (especially propylene carbonate) has demonstrated extraordinary pseudocapacitive performance as an electrolyte in the supercapacitor configuration. However, the influence of the solvated ions on the diffusivity of the solvent molecules is yet to be understood. We examine the impact of LiTFSI on the diffusivity in five organic solvents: acetonitrile (ACN), tetrahydrofuran (THF), methanol (MeOH), dimethyl sulfoxide (DMSO), and propylene carbonate (PC) using a combination of neutron scattering, conductivity measurements, and molecular dynamics simulations. The extent of the diffusivity reduction in the concentration regime of ≤1 M directly correlates with the solvent mole fraction at which the solvation shells around Li + ions are of similar size in all the solvents, resulting in a universal ~50% reduction in the solvent diffusivity. These results provide guidance for formulation of the new electrolytes to enhance the performance of energy storage devices.

conductivity↗

Homologous mutations in human β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman–Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin’s highly allosteric nature.

59 BASIC BIOLOGICAL SCIENCES↗