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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

Si–Cl 2 –Ar + Atomic Layer Etching Window: A Fundamental Study Using Molecular Dynamics Simulations and a Reduced Order Model

Silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ) is studied by using molecular dynamics (MD) simulations and a reduced order model (ROM). Here, the purpose of this study is to elucidate the properties of the ALE window, a range of ion energies where the amount of Si etched over a series of cycles is nonzero and nearly independent of ion energy. Experimental studies of the Si–Cl 2 –Ar + ALE system report contradictory results related to the ALE window’s ion energy range. Both MD simulations and the ROM show that there is an ALE window present from approximately 15 to 20 eV for normal incidence argon ions. The Si–Cl 2 –Ar + system, therefore, exhibits a narrow ALE window. The amount of Si etched per cycle is less than one atomic layer because of the higher etch yield of Cl atoms relative to atomic Si and silicon chlorides. A modified version of the ROM with an artificially increased Si physical sputtering threshold energy expands the ALE window, illustrating the importance of the difference in chemical and physical sputtering threshold energies in the ALE window energy range. The ROM is also used to examine the dependence of the EPC on the Ar + ion fluence.

energy↗
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗
DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

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

Graphic contrastive learning analyses of discontinuous molecular dynamics simulations: Study of protein folding upon adsorption

A comprehensive understanding of the interfacial behaviors of biomolecules holds great significance in the development of biomaterials and biosensing technologies. In this work, we used discontinuous molecular dynamics (DMD) simulations and graphic contrastive learning analysis to study the adsorption of ubiquitin protein on a graphene surface. Our high-throughput DMD simulations can explore the whole protein adsorption process including the protein structural evolution with sufficient accuracy. Contrastive learning was employed to train a protein contact map feature extractor aiming at generating contact map feature vectors. Subsequently, these features were grouped using the k-means clustering algorithm to identify the protein structural transition stages throughout the adsorption process. The machine learning analysis can illustrate the dynamics of protein structural changes, including the pathway and the rate-limiting step. Our study indicated that the protein–graphene surface hydrophobic interactions and the π–π stacking were crucial to the seven-stage adsorption process. Upon adsorption, the secondary structure and tertiary structure of ubiquitin disintegrated. The unfolding stages obtained by contrastive learning-based algorithm were not only consistent with the detailed analyses of protein structures but also provided more hidden information about the transition states and pathway of protein adsorption process and structural dynamics. Our combination of efficient DMD simulations and machine learning analysis could be a valuable approach to studying the interfacial behaviors of biomolecules.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Investigations of water/oxide interfaces by molecular dynamics simulations

Water/oxide interfaces are ubiquitous on earth and show significant influence on many chemical processes. For example, understanding water and solute adsorption as well as catalytic water splitting can help build better fuel cells and solar cells to overcome our looming energy crisis; the interaction between biomolecules and water/oxide interfaces is one hypothesis to explain the origin of life. However, knowledge in this area is still limited due to the difficulty of studying water/solid interfaces. As a result, research using increasingly sophisticated experimental techniques and computational simulations has been carried out in recent years. Although it is difficult for experimental techniques to provide detailed microscopic structural information, molecular dynamics (MD) simulations have satisfactory performance. In this report, we discuss classical and ab initio MD simulations of water/oxide interfaces. Generally, we are interested in the following questions: How do solid surfaces perturb interfacial water structure? How do interfacial water molecules and adsorbed solutes affect solid surfaces and how do interfacial environments affect solvent and solute behavior? Finally, we discuss progress in the application of neural network potential based MD simulations, which offer a promising future because this approach has already enabled ab initio level accuracy for very large systems and long trajectories.

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

Modification of a force field for molecular dynamics simulations of silicon etching by chlorine atoms

Here, a modified classical molecular dynamics (MD) force field is presented for the interaction between silicon (Si) and chlorine (Cl). The original version of the force field is shown to significantly overestimate the probability of Si etching by thermal Cl atoms. However, the modified force field corrects this problem and results in generally good agreement with experimental data. Further, it is shown that while the modification of the force field improves the prediction of Si spontaneous etching with Cl atoms, it does not degrade predictions of atomic-layer etching of Si with Cl[Formula: see text] molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗
NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗
DOE OSTI2020

Modelling dislocation-graphene interactions in a BCC Fe matrix by molecular dynamics simulations and gradient plasticity theory

We report graphene nanosheets (GNS) can enhance the strength and ductility of metal-based composites as they can obstruct the propagation of dislocations. The present article employs Molecular Dynamics (MD) simulations to investigate dislocation-GNS interaction mechanisms and possible influencing factors, including the number of GNS layers, the thickness of the metallic amorphous layer and the C - C bond strength. The results indicated that the shear strength of the metal/GNS interface and the bending stiffness of GNS determined the ability of GNS to block dislocation transmission. A physically based phenomenological parameter that can capture such dislocation-GNS interactions is the mechanical interface energy that has been put forth within gradient plasticity. By fitting the theoretical expressions to the simulation data, it was possible to obtain estimates for the mechanical interface energy for the GNS. It was found that increasing the GNS layers and adding an amorphous layer resulted in a strengthening in the stress–strain response and increased the value of this interfacial parameter. This indicates that the mechanical interfacial energy can be a unified measure for capturing and tuning the strength of various interfaces such as grain boundaries, GNS, amorphous-crystalline interface and bimetal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structure of Water Adsorbed on Nanocrystalline Calcium Silicate Hydrate Determined from Neutron Scattering and Molecular Dynamics Simulations

Calcium silicate hydrate (C-S-H) is a disordered, nanocrystalline material that acts as a primary binding phase in Portland cement. Thin films of water are present on the surfaces and in nanopores of C-S-H, impacting many of its chemical and mechanical properties, such as ion transport, creep, or thermal behavior. Despite decades of research, a full understanding of the structural details of adsorbed, confined, and bulk water in C-S-H remains elusive. In this work, we applied a multitechnique study involving molecular dynamics (MD) simulations validated by neutron diffraction with isotopic substitution (NDIS) and X-ray scattering methods to investigate the structure of water in C-S-H and C-A-S-H (an Al-bearing, low-CO 2 C-S-H substitute). Direct comparison of NDIS data with the MD results reveals that the structure of confined and interfacial water differs significantly from the bulk water and exhibits a larger degree of mesoscale ordering for more hydrated C-S-H structures. As a result, this observation suggests an important role of water as a stabilizer of the atomistic-level structure of C-S-H.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Dynamics Simulation of the Structures, Dynamics, and Aggregation of Dissolved Organic Matter

Dissolved organic matter (DOM) plays a significant role in the transport and transformation of pollutants in the aquatic environment. However, the experimental characterization of DOM has been limited mainly to bulk properties, and the molecular-level interactions among various components of DOM remain to be fully characterized. In this study, we use molecular dynamics (MD) simulations to probe the structural properties of model DOM systems at atomic detail. The 200 ns simulations, validated by available experimental data, reveal processes and mechanisms by which chemical species (cations, peptides, lipids, lignin, carbohydrates, and some low-molecular-weight aliphatic and aromatic compounds) aggregate to form complex DOM. The DOM aggregates are dynamic, consisting of a hydrophobic core and amphiphilic exterior. The lipid tails and other hydrophobic fragments form the core, with hydrophilic and amphiphilic groups exposed to water, making DOM accessible to both polar and nonpolar species. Thus, the lipid component acts as a nucleator, whereas cations (especially Ca 2+ ) connect the molecular fragments on the surface by coordinating with the O-containing functional groups of DOM. The structural details revealed here provide new insights including surface accessible atoms, overall assemblage, and interactions among the molecules of DOM for understanding the kinetics and mechanisms through which DOM interacts with metal and other contaminants.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Cardiovascular Responses to Simulated Spaceflight: Molecular Signatures and Surrogate Outputs to Measure CVD Risk

During extended space missions beyond low Earth orbit, astronauts will encounter prolonged periods of weightlessness and low dose space radiation. Previous studies have shown that exposure to small doses of high LET radiation (< 50 cGy) can lead to both short-term and long-term alterations in heart function, structure and underlying molecular mechanisms. In this study, we aim to identify the molecular signature associated with the cardiovascular response to simulated galactic cosmic radiation (5-ion GCR) alone or in combination with simulated weightlessness at time intervals relevant to mission length and recovery. Additionally, we aim to determine whether sex impacts cardiovascular responses to these spaceflight factors. Our overarching goal is to enhance our understanding of the cardiovascular risks associated with extended space missions and the clinical endpoints they suggest. We hypothesize that exposure to simulated space radiation leads to enduring alterations in the transcriptome, redox signaling and cytokine environment of cardiovascular tissue, some which have known links with reduced cardiovascular performance, aging, and increased risk of cardiovascular disease (CVD). Furthermore, we posit that simulated space radiation exposure in combination with simulated microgravity exacerbates cardiovascular deficits compared to single factor exposure. Female and male C57BL/6J mice, aged 23-24 weeks, were exposed to a single dose of 5, 15, or 50 cGy of 5-ion GCR, or sham-treated (0 cGy). Euthanasia was performed at 14 days and ~4 months post-irradiation. Hearts, aorta and blood plasma were collected shortly thereafter. RNA-sequencing of left ventricles at ~4 months post-GCR exposure revealed sex differences in the heart transcriptome with a few genes showing radiation-dependent changes in expression levels. Notably, some of the differentially expressed genes in 15 and 50 cGy GCR groups are known to play roles in the development of CVD. Analysis of protein levels of a subset of inflammatory cytokines in the heart indicated sex differences but no differences between sham and 50 cGy groups. Results also showed correlations among differentially expressed genes and a subset of inflammatory cytokines, with some correlations altered by GCR exposure. These findings suggest that GCR exposure can modify protein and gene networks linked to inflammation and CVD progression. In the aorta, telomere lengths were comparable across treatment groups sexes. Mitochondrial copy number is a biomarker for mitochondrial function with decreased copy numbers associated with cardiometabolic disease traits. Mitochondrial copy numbers of aorta also showed no sex nor dose differences. In a second study, mice underwent one week of simulated microgravity by hindlimb unloading (HU) and then exposed to a single dose of 15 cGy of 5-ion GCR. HU was conducted for an additional two weeks following GCR exposure. Single factor exposure groups (HU or GCR only) also were included in the study. Euthanasia was then performed and the same tissues were collected. Protein levels of select inflammatory cytokines in the heart showed sex-dependent differences in expression. In the aorta, telomere lengths and mitochondrial copy number also showed sex differences. In summary, our results indicate differences between sexes in biomarkers related to cardiovascular health. Exposure to 5-ion GCR or HU, alone or in combination, did not result in changes in most of the cardiovascular biomarkers that were examined. However, in the heart, simulated space radiation at doses of 15 and 50 cGy led to long-term alterations in the expression levels of a small group of genes known to be associated with the progression of CVD. The long-term transcriptomic changes resulting from exposure to simulated space radiation should be carefully investigated to mitigate adverse cardiovascular events during and after deep space missions. Our results also highlight the importance of sex-specific strategies in monitoring and maintaining cardiovascular health during and after deep space missions.

cardiovascular↗
DOE OSTI2020

Orientation and Conformation of Proteins at the Air–Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. Here in this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

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

Molecular Dynamics Simulation Study of the Protonation State Dependence of Glutamic Acid Transport through a Cyclic Peptide Nanotube

The effect of the protonation state of glutamic acid on its translocation through cyclic peptide nanotubes (CPNs) was assessed by using molecular dynamics (MD) simulations. Anionic (GLU–), neutral zwitterionic (GLU0), and cationic (GLU+) forms of glutamic acid were selected as three different protonation states for an analysis of energetics and diffusivity for acid transport across a cyclic decapeptide nanotube. Based on the solubility-diffusion model, permeability coefficients for the three protonation states of the acid were calculated and compared with experimental results for CPN-mediated glutamate transport through CPNs. Potential of mean force (PMF) calculations reveal that, due to the cation-selective nature of the lumen of CPNs, GLU–, so-called glutamate, shows significantly high free energy barriers, while GLU+ displays deep energy wells and GLU0 has mild free energy barriers and wells inside the CPN. The considerable energy barriers for GLU– inside CPNs are mainly attributed to unfavorable interactions with DMPC bilayers and CPNs and are reduced by favorable interactions with channel water molecules through attractive electrostatic interactions and hydrogen bonding. Unlike the distinct PMF curves, position-dependent diffusion coefficient profiles exhibit comparable frictional behaviors regardless of the charge status of three protonation states due to similar confined environments imposed by the lumen of the CPN. The calculated permeability coefficients for the three protonation states clearly demonstrate that glutamic acid has a strong protonation state dependence for its transport through CPNs, as determined by the energetics rather than the diffusivity of the protonation state. In addition, the permeability coefficients also imply that GLU– is unlikely to pass through a CPN due to the high energy barriers inside the CPN, which is in disagreement with experimental measurements, where a considerable amount of glutamate permeating through the CPN was detected. To resolve the discrepancy between this work and the experimental observations, several possibilities are proposed, including a large concentration gradient of glutamate between the inside and outside of lipid vesicles and bilayers in the experiments, the glutamate activity difference between our MD simulations and experiments, an overestimation of energy barriers due to the artifacts imposed in MD simulations, and/or finally a transformation of the protonation state from GLU– to GLU0 to reduce the energy barriers. Altogether, our study demonstrates that the protonation state of glutamic acid has a strong effect on the transport of the acid and suggests a possible protonation state change for glutamate permeating through CPNs.

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

Correlation between time and angular alignment in molecular dynamics simulations of heavy ion collisions

Neutron-proton equilibration is a process which has been used to study the density dependence of the symmetry energy term in the nuclear equation-of-state. Here, this study utilizes constrained molecular dynamics (CoMD) simulations of 70 Zn + 70 Zn with collision energies of 35 and 45 MeV/nucleon. An algorithm is used which searches through CoMD events and identifies the PLF* after it separates from the target and determines its lifetime, Δt. It also determines the fragments that the PLF* breaks apart into and determines their angular alignment. This technique gives an opportunity to explore how the average alignment of dynamically produced fragments, $\langle$α $\rangle$ dyn , evolves with PLF* lifetime. An approximately linear relationship was determined with d$\langle$α $\rangle$ dyn /dΔt = 0.98 ± 0.08 rad / zs and 1.06 ± 0.09 rad / zs for the 35 and 45 MeV / nucleon, respectively, indicating a correlation with magnitude consistent with classically determined values which were used for prior experimental studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Solvent-induced membrane stress in biofuel production: molecular insights from small-angle scattering and all-atom molecular dynamics simulations

The disruptive effect of organic solvents on microbial membranes represents a significant challenge to the economical production of green fuels and value-added chemicals from lignocellulosic feedstocks. One route to overcoming this challenge is to engineer microbes with membranes capable of resisting organic solvent stresses. In this regard, it is useful to understand the mechanisms by which organic solvents disrupt typical biomembranes. In this study, molecular dynamics (MD) simulation, complemented by small-angle X-ray and neutron scattering (SANS/SAXS), provide a molecular-scale view of the disruption of a microbial model membrane by 1-butanol and tetrahydrofuran (THF), two common water–organic cosolvent mixtures of importance in biofuel production. Solvent interactions at the interface between the head-group and fatty acid tail regions lead to more dramatic membrane changes than interactions solely at the head-groups or tails. Although both organic solvents are found to partition into the membrane, the depth of solvent penetration into the membrane is quite different. Specifically, 1-butanol localizes near the interface between the lipid heads and tails at low concentrations, but partitions into both the head and tail regions at high concentrations. In contrast, THF, overall, partitions less than 1-butanol and prefers the lipid tail regions. Importantly, the presence of 1-butanol near the head/tail interface introduces drastic membrane changes not seen with THF. The organic solvent interactions with the lipids lead to membrane thinning and fluidization, but more so for 1-butanol than for THF. These results suggest that an aim for the future engineering of robust membranes could be to design lipid head groups that reduce the accumulation of organic solvents at the head–tail interface and that rational designs need also be cognizant of the different solvent-specific mechanisms responsible for membrane disruption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗