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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Investigating Lignin Aggregation and Interactions with Solvents during γ-Valerolactone (GVL) Pretreatment: A Combined Small Angle Neutron Scattering and Molecular Simulations Study

The strong tendency of lignin to aggregate in solution, coupled with limited understanding of how its molecular structure governs this behavior, hinders its effective utilization in biorefineries. Here, in this study, we investigated the solution behavior of lignin extracted from poplar using γ-valerolactone/water (GVL/H2O, 9:1 wt/wt) through combined small-angle neutron scattering (SANS) and molecular dynamics (MD) simulations. Lignin samples obtained at 100 °C (L100) and 120 °C (L120) differed in β–O–4 content, hydroxyl distribution, and S/G ratio, enabling direct assessment of how molecular composition governs solvation and aggregation. SANS showed that L120 formed rigid and elongated cylindrical aggregates at 25 °C that transitioned to more flexible spheroidal structures by 50 °C and remained stable up to 80 °C, whereas L100 adopted globular aggregates that progressively collapsed with increasing temperature. MD simulations reinforced these observations by showing that S-rich (L120-like) oligomers had larger radii of gyration, stronger solvent coordination driven by methoxy groups, and fewer lignin–lignin contacts. In contrast, G-rich (L100-like) oligomers displayed persistent aggregation and lower solubility. Collectively, these results reveal that increased aromatic methoxylation enhances lignin–solvent interactions and suppresses self-association, whereas reduced methoxylation and higher β–O–4 content promote persistent aggregation into colloid-like structures with restricted solvent penetration into the aggregate interior.

biorefinery pretreatment↗
DOE OSTI2014

Communication: Molecular simulation study of kaolinite intercalation with realistic layer size

Intercalation phenomena of kaolinite in aqueous potassium acetate and in hexyl-amine solutions are studied by large scale molecular dynamics simulations. The simulated kaolinite particle is constructed from ∼6.5 × 10{sup 6} atoms, producing a particle size of ∼100 nm × 100 nm × 10 nm. The simulation with potassium acetate results in a stable kaolinite-potassium acetate complex, with a basal spacing that is in close agreement with experimental data. The simulation with hexyl-amine shows signs of the experimentally observed delamination of kaolinite (the initial phase of the formation of nanoscrolls from the external layers)

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

Development of ions-TIP4P-Ew force fields for molecular processes in bulk and at the aqueous interface using molecular simulations

Monovalent ions such as alkalis and halides play essential roles in aqueous and biological systems; thus a proper representation of these ions is significant in explicit molecular dynamics simulations. In this study, we re-parameterize ions interaction potentials based on TIP4PEw water at 298K with the additive non-polarizable models. The charges of the ions were kept as +1 for alkali metal ions (Li+, Na+, K+, Rb+, Cs+) and -1 for halides (F-, Cl-, Br-, I-). The experimental enthalpies of hydration were used to optimize the potential parameters. To validate the potential models, extensive molecular dynamics calculations were carried out to examine the bulk, interfacial, static structural and dynamical properties of the aqueous ionic solutions. These included radial distribution functions, angle distributions, velocity autocorrelation functions, diffusion coefficients, binding energies, mean residence time, surface potential, and potential of mean force. The computed results agreed with the experimental data and observations. LXD was funded by the Division of Chemical Sciences, Geosciences and Biosciences, Office of Basic Energy Sciences, US Department of Energy. Battelle operates the Pacific Northwest National Laboratory for the U.S. Department of Energy.

Peng, Tiefeng↗
DOE OSTI2015

Probing the folded state and mechanical unfolding pathways of T4 lysozyme using all-atom and coarse-grained molecular simulation

The Bacteriophage T4 Lysozyme (T4L) is a prototype modular protein comprised of an N-terminal and a C-domain domain, which was extensively studied to understand the folding/unfolding mechanism of modular proteins. To offer detailed structural and dynamic insights to the folded-state stability and the mechanical unfolding behaviors of T4L, we have performed extensive equilibrium and steered molecular dynamics simulations of both the wild-type (WT) and a circular permutation (CP) variant of T4L using all-atom and coarse-grained force fields. Our all-atom and coarse-grained simulations of the folded state have consistently found greater stability of the C-domain than the N-domain in isolation, which is in agreement with past thermostatic studies of T4L. While the all-atom simulation cannot fully explain the mechanical unfolding behaviors of the WT and the CP variant observed in an optical tweezers study, the coarse-grained simulations based on the Go model or a modified elastic network model (mENM) are in qualitative agreement with the experimental finding of greater unfolding cooperativity in the WT than the CP variant. Interestingly, the two coarse-grained models predict different structural mechanisms for the observed change in cooperativity between the WT and the CP variant—while the Go model predicts minor modification of the unfolding pathways by circular permutation (i.e., preserving the general order that the N-domain unfolds before the C-domain), the mENM predicts a dramatic change in unfolding pathways (e.g., different order of N/C-domain unfolding in the WT and the CP variant). Based on our simulations, we have analyzed the limitations of and the key differences between these models and offered testable predictions for future experiments to resolve the structural mechanism for cooperative folding/unfolding of T4L.

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

Molecular simulation data for 'Data-guided Multi-Map variables for ensemble refinement of molecular movies'

These trajectories, scripts, and analysis performed on Summit underly the work published as 'Data-guided Multi-Map variables for ensemble refinement of molecular movies'. The trajectories include equilibrium and non-equilibrium sampling of ADK, CODH, and FLPP3, the scripts used to build the systems, and the scripts used to analyze the output. The directory structure is explained further in an internal README file.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Effect of homopolymer matrix on diblock copolymer grafted nanoparticle conformation and potential of mean force: A molecular simulation study

ABSTRACT We study the effect of homopolymer A or B matrix on the conformations and effective interactions of AB diblock copolymer grafted particles using coarse‐grained molecular dynamics simulations. In an A homopolymer matrix we observe patchy conformations within the AB diblock copolymer grafted layer, where the number of B patches is controlled by the A‐A attractive interaction strength. In a B homopolymer matrix the grafted particle takes on a core‐corona conformation, where the inner A block aggregates near the particle surface and the outer B block forms a corona that interacts with the B matrix. The potential of mean force (PMF) between two particles in an A homopolymer matrix has a long‐ranged attractive well with a minima at intermediate distances corresponding to the location of the outer B block patches. The PMF between two particles in a B homopolymer matrix has an attractive well at short interparticle distances corresponding to the size of the inner A block. We isolate the contribution of the homopolymer matrix on the PMF between the two diblock copolymer grafted particles, by deducting the PMF in the absence of a matrix, assuming the contributions of the grafted particle and matrix to the PMF to be additive. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53 , 76–88

Estridge, Carla E.↗
DOE OSTI2020

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

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

Mesoscale Polymer Dissolution Probed by Raman Spectroscopy and Molecular Simulations

The diffusion of various solvents into a polystyrene (PS) matrix was probed experimentally by monitoring the temporal profiles of the Raman spectra and theoretically from molecular dynamics (MD) simulations of the binary system. The simulation results assist in providing a fundamental, molecular level connection between the mixing/dissolution processes and the difference = solvent – PS in the values of the Hildebrand parameter () between the two components of the binary systems: solvents having similar values of with PS (small ) exhibit fast diffusion into the polymer matrix, whereas the diffusion slows down considerably when the ’s are different (large ). To this end, the Hildebrand parameter was identified as a useful descriptor that governs the process of mixing in polymer – solvent binary systems. The experiments also provide insight into further refinements of the models specific to non-Fickian diffusion phenomena that need to be used in the simulations.

Chang, Tsun-Mei↗
DOE OSTI2014

Molecular simulation study of role of polymer–particle interactions in the strain-dependent viscoelasticity of elastomers (Payne effect)

The strain-amplitude dependence of viscoelastic behavior of model crosslinked elastomers containing various concentrations of spherical nanoparticles (NPs) was studied by non-equilibrium molecular dynamics simulation. All the filler NPs were in monodispersed state and the interactions between these particles were purely repulsive. The polymer–particle interactions were attractive and their interaction energies were tuned in a broad range. Through the computational study, many important features of the behavior of particle-reinforced elastomers observed in experiments, including the Payne effect, were successfully reproduced. It was shown that the magnitude of the Payne effect was found to depend on the polymer–particle interaction and the filler loading. By examining the microstructures of the simulation systems and their evolution during oscillatory shear, four different mechanisms for the role of the polymer–particle interactions in the Payne effect were revealed that consist of the debonding of polymer chains from NP surfaces, the breakage of polymer-shell-bridged NP network, the rearrangement of the NPs in the network into different layers and the shear-induced yielding of the rigid polymer shell in-between neighboring NPs.

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

Sulfonated polybenzimidazole membrane with graphene oxide additive for 2,3-butanediol/water separation: A molecular simulation

Membrane separation for 2,3-butanediol (2,3-BDO) recovery from fermentation broth is highly valued for sustainable and renewable processes, but it requires efficient membrane materials. Here, this work evaluates the sulfonated polybenzimidazole (sPBI) and its graphene oxide (GO) doped composite membrane for separating 2,3-BDO and water via atomistic simulations. Density functional theory calculations are applied to identify various forms of sPBI structures and quantify their binding interactions with 2,3-BDO and water. Classical molecular dynamic simulations are used to evaluate the structural changes, diffusivity, and selectivity of 2,3-BDO and water in different sPBI models, GO surfaces, and GO-doped sPBI composite models. Our results suggest that sPBI slightly increases the crystallinity of the membrane structures, enhances the adsorption strength for both 2,3-BDO and water, and improves the water/2,3-BDO selectivity by 2–3 times. The GO surfaces display a maximum selectivity at a surface coverage of 0.1–0.15 for both hydroxyl and epoxy surface groups. The addition of GO flakes to sPBI creates new interaction sites for 2,3-BDO and water at the interface of sPBI and GO, and the water/2,3-BDO selectivity of GO-doped sPBI models is further increased up to 3 times. This work illustrates how the integrated addition of sPBI and GO flakes offers a promising approach to selective separation of 2,3-BDO and water, providing theoretical guidance for polybenzimidazole-based membranes in the potential application of 2,3-BDO recovery.

2,3-butanediol↗
DOE OSTI2014

Rupture mechanism of liquid crystal thin films realized by large-scale molecular simulations

The ability of liquid crystal (LC) molecules to respond to changes in their environment makes them an interesting candidate for thin film applications, particularly in bio-sensing, bio-mimicking devices, and optics. Yet the understanding of the (in)stability of this family of thin films has been limited by the inherent challenges encountered by experiment and continuum models. Using unprecedented largescale molecular dynamics (MD) simulations, we address the rupture origin of LC thin films wetting a solid substrate at length scales similar to those in experiment. Our simulations show the key signatures of spinodal instability in isotropic and nematic films on top of thermal nucleation, and importantly, for the first time, evidence of a common rupture mechanism independent of initial thickness and LC orientational ordering. We further demonstrate that the primary driving force for rupture is closely related to the tendency of the LC mesogens to recover their local environment in the bulk state. Our study not only provides new insights into the rupture mechanism of liquid crystal films, but also sets the stage for future investigations of thin film systems using peta-scale molecular dynamics simulations.

Nguyen, Trung D↗
DOE OSTI2015

Molecular simulation of fluids with non-identical intermolecular potentials: Thermodynamic properties of 10-5 + 12-6 Mie potential binary mixtures

General methods for combining interactions between particles characterised by non-identical intermolecular potentials are investigated. The combination methods are tested by performing molecular dynamics simulations to determine the pressure, energy, isochoric and isobaric heat capacities, thermal expansion coefficient, isothermal compressibility, Joule-Thomson coefficient, and speed of sound of 10-5 + 12-6 Mie potential binary mixtures. In addition to the two non-identical Mie potentials, mixtures are also studied with non-identical intermolecular parameters. The combination methods are compared with results obtained by simply averaging the Mie exponents. When either the energy or size parameters are non-identical, very significant differences emerge in the thermodynamic properties predicted by the alternative combination methods. The isobaric heat capacity is the thermodynamic property that is most affected by the relative magnitude of the intermolecular potential parameters and the method for combining non-identical potentials. Either the arithmetic or geometric combination of potentials provides a simple and effective way of performing simulations involving mixtures of components characterised by non-identical intermolecular potentials, which is independent of their functional form.

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

Understanding Cellulose Through Molecular Simulation and Electron Tomography

High-resolution cellulose crystal structures have been determined from diffraction experiments using large diameter microfibrils as the sample material. However, cellulose microfibrils in plants are much smaller in diameter, and are more difficult to directly examine experimentally. Molecular dynamics simulation combined with quantum chemical calculations can help to elucidate the structure and dynamics of small diameter cellulose microfibrils. These simulation techniques also aid in the interpretation of electron tomography volumetric structural data from maize cell walls, where pretreatment with dilute acid or ammonia reveals microfibril geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Inelastic neutron scattering and molecular simulation of the dynamics of interlayer water in smectite clay minerals

The study of mineral–water interfaces is of great importance to a variety of applications including oil and gas extraction, gas subsurface storage, environmental contaminant treatment, and nuclear waste repositories. Understanding the fundamentals of that interface is key to the success of those applications. Confinement of water in the interlayer of smectite clay minerals provides a unique environment to examine the interactions among water molecules, interlayer cations, and clay mineral surfaces. Smectite minerals are characterized by a relatively low layer charge that allows the clay to swell with increasing water content. Montmorillonite and beidellite varieties of smectite were investigated to compare the impact of the location of layer charge on the interlayer structure and dynamics. Inelastic neutron scattering of hydrated and dehydrated cation-exchanged smectites was used to probe the dynamics of the interlayer water (200–900 cm –1 spectral region) and identify the shift in the librational edge as a function of the interlayer cation. Molecular dynamics simulations of equivalent phases and power spectra, derived from the resulting molecular trajectories, indicate a general shift in the librational behavior with interlayer cation that is generally consistent with the neutron scattering results for the monolayer hydrates. Both neutron scattering and power spectra exhibit librational structures affected by the location of layer charge and by the charge of the interlayer cation. Furthermore, divalent cations (Ba 2+ and Mg 2+ ) characterized by large hydration enthalpies typically exhibit multiple broad librational peaks compared to monovalent cations (Cs + and Na + ), which have relatively small hydration enthalpies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular simulation of the structural and thermodynamic properties of n-alkane/brine interfacial systems with nonionic surfactants

All-atom molecular dynamics (MD) simulations are used to study the structural and thermodynamic properties of water + n-heptane interfacial systems in the presence of salinity (NaCl, NaI, KCl, and KI) and two different nonionic surfactants. Excellent qualitative and quantitative agreement with interfacial tension experimental data is obtained, but this requires applying a charge scaling factor to the ionic species, suggesting the likely role of polarizability effects. Tensoactive behavior is observed for the NaI and KI systems, coinciding with layering of the iodide species near the interface. Further, there are significant interaction differences between the surfactants and the different ionic species, but the structural and thermodynamic behavior of the surfactant molecules is relatively unaffected by the specific salt species. Our analysis suggest that the surfactant-surfactant interactions play a critical role in determining the interfacial behavior, with a much smaller impact associated with the composition of the brine phase.

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

Predicting peak spectral sensitivities of vertebrate cone visual pigments using atomistic molecular simulations

Vision is the dominant sensory modality in many organisms for foraging, predator avoidance, and social behaviors including mate selection. Vertebrate visual perception is initiated when light strikes rod and cone photoreceptors within the neural retina of the eye. Sensitivity to individual colors, i.e., peak spectral sensitivities (λ max ) of visual pigments, are a function of the type of chromophore and the amino acid sequence of the associated opsin protein in the photoreceptors. Large differences in peak spectral sensitivities can result from minor differences in amino acid sequence of cone opsins. To determine how minor sequence differences could result in large spectral shifts we selected a spectrally-diverse group of 14 teleost Rh2 cone opsins for which sequences and λ max are experimentally known. Classical molecular dynamics simulations were carried out after embedding chromophore-associated homology structures within explicit bilayers and water. These simulations revealed structural features of visual pigments, particularly within the chromophore, that contributed to diverged spectral sensitivities. Statistical tests performed on all the observed structural parameters associated with the chromophore revealed that a two-term, first-order regression model was sufficient to accurately predict λ max over a range of 452–528 nm. The approach was accurate, efficient and simple in that site-by-site molecular modifications or complex quantum mechanics models were not required to predict λ max . These studies identify structural features associated with the chromophore that may explain diverged spectral sensitivities, and provide a platform for future, functionally predictive opsin modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Molecular Simulation of Cesium Adsorption at the Basal Surface of Phyllosilicate Minerals

A better understanding of the thermodynamics of radioactive cesium uptake at the surfaces of phyllosilicate minerals is needed to understand mechanisms of its selective adsorption and help guide the development of practical and inexpensive decontamination techniques. In this work, molecular dynamics simulations were carried out to determine the thermodynamics of adsorption of Cs + at the basal surface of six 2:1 phyllosilicate minerals, namely pyrophyllite, illite, muscovite, phlogopite, celadonite, and margarite. These minerals were selected to isolate the effects of the magnitude of the permanent layer charge (≤ 2), its location (tetrahedral versus octahedral sheet), and the structure of the octahedral sheet (dioctahedral versus trioctahedral). Good agreement was obtained with experiment in terms of the hydration free energy of Cs + and the structure and thermodynamics of Cs + adsorption at the muscovite basal surface, for which published data were available for comparison. With the exception of pyrophyllite, which did not exhibit an inner-sphere free energy minimum, all phyllosilicate minerals showed similar behavior with respect to Cs + adsorption; notably, Cs + adsorption was predominantly inner-sphere whereas outer-sphere adsorption was very weak with the simulations predicting the formation of an extended outer-sphere complex. For a given location of the layer charge, the free energy of adsorption as an inner-sphere complex was found to vary linearly with the magnitude of the layer charge. For a given location and magnitude of the layer charge, adsorption at phlogopite (trioctahedral sheet structure) was much less favorable than at muscovite (dioctahedral sheet structure) due to the electrostatic repulsion between the adsorbed Cs + and the hydrogen atom of the hydroxyl group directly below the six-membered siloxane ring cavity. For a given magnitude of the layer charge and structure of the octahedral sheet, adsorption at celadonite (layer charge located in the octahedral sheet) was favored over muscovite (layer charge located in the tetrahedral sheet) due to the increased distance with surface potassium ions.

(001) surface↗
DOE OSTI2024

Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons

Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗