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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2024

Effect of reaction media on hydrogenolysis of polyethylene plastic waste: Polymer-surface interactions in small alkane/polymer blends

The polymer reaction media and its properties can be altered by recycling a fraction of liquid products or adding alkane solvents. Less clear is whether this strategy affects hydrogenolysis. Herein, we investigated the effect of short-chain alkanes C n consisting of n carbons (n=8, 16, and 32) on the upcycling of high-density polyethylene (HDPE) plastic waste to lubricant-range products over Ru/TiO 2 catalysts by multiscale simulations and experiments. First, we trained a force field for polymer/surface interactions on a Ru 22 nanoparticle (NP) supported on TiO 2 . Using replica exchange molecular dynamics simulations, we studied the effect of small hydrocarbons on the adsorption of a surrogate polymer, C 142 , on the catalyst. We found segregation of long chains (C 142 ) at the catalyst surface due to the enthalpy gained by adsorbing more C-C bonds of the long chains, compensating for entropic losses upon adsorption. Short-chain molecules decrease the adsorbed carbons of long chains on the Ru NP due to blocking Ru active sites. Compared to the bulk chains, competitive adsorption results in a broader, heavy-tailed distribution of end-to-end distance of adsorbed chains. Our experiments demonstrated that catalyst activity declines significantly beyond simple dilution due to changes in polymer adsorption, and tuning the reaction media by creating suitable blends impacts hydrogenolysis. Density distributions for a 50:50%wt mixture of PP and PE show that PE chains are segregated at the surface, so they are prone to C-C bond breaking much faster than PP chains. H/D exchange experiments show preferential deuteration of PE, while CH 3 groups of PP remain undeuterated. Furthermore, this may be explained by the preferential sorption of PE over PP, leading to specific distribution in the polymer blend.

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

Polymorphism in Self-Assembly of Short Peptoid Sequences

Due to various applications enabled by diverse morphologies of self-assembled sequence-defined polymers, controlling the self-assembly of synthetic peptidomimetics into designed morphologies has emerged as a promising route for the development of bioinspired functional materials. Herein, we report morphological control over the assembly of a series of short peptoids, or poly-N-substituted glycines, that contain asymmetric hydrophobic domains. We demonstrate that the inherent flexibility of amphiphilic peptoid bilayers drives assembly polymorphism, resulting in the coexistence of nanosheets, twisted ribbons, and nanofibers three distinct morphologies. By tuning peptoid molecular interactions through variations in sequence design, solution pH, and temperature, we demonstrate precise control over the twisting and folding of peptoid bilayers, enabling the formation of well-defined nanosheets and nanohelices. Molecular dynamics simulations further unravel how the introduction of asymmetric hydrophobic domains enables the flexibility of peptoid bilayers and results in peptoid assembly polymorphism. By tuning peptoid molecular interactions through heating, we further demonstrate the transformation of nanosheets into nanohelices. We envision that our mechanistic investigation of peptoid assembly polymorphism provides a strong foundation for leveraging peptoid sequences and chemistries to achieve controlled molecular interactions, driving the creation of biomimetic materials with tailored morphologies and functionalities.

assembly polymorphism↗
DOE OSTI2025

Molecular-resolution imaging of ice crystallized from liquid water by cryogenic liquid-cell TEM

Despite the ubiquity of ice, a molecular-resolution image of nanoscopic defects or microstructures in ice crystallized from liquid water has never been obtained. This is mainly due to the difficulties in preparing and preserving crystalline ice samples that can survive under high-resolution imaging conditions. Here, we report the stabilization and Å-resolution electron imaging of ice Ih crystallized from liquid water by developing cryogenic liquid-cell transmission electron microscopy (CRYOLIC-TEM). We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defects such as misoriented subdomains and trapped gas bubbles, which are stabilized by molecular-scale structural motifs. Importantly, bubble surfaces adopt low-energy nanofacets and create negligible strain fields in the surrounding crystal. These bubbles can dynamically nucleate, grow, migrate, dissolve, and coalesce under electron irradiation and be monitored in situ near a steady state. This work improves our understanding of water crystallization behaviors at a molecular spatial resolution.

Materials science↗
DOE OSTI2020

Concentration and velocity profiles in a polymeric lithium-ion battery electrolyte

Predictive knowledge of ion transport in electrolytes which bridges microscopic and macroscopic length scales is imperative to design new ion conductors and to simulate device performance. Here, we employed a novel approach combining operando X-ray photon correlation spectroscopy, X-ray absorption microscopy, continuum modelling, and molecular dynamics simulations to probe the ion transport in a baseline polymeric lithium-ion battery electrolyte. In a Li/PEO–LiTFSI/Li symmetric cell under polarization, we determined and rationalized microscopic properties including local electrolyte velocities and ion correlations and connected this insight to measured and simulated macroscopic ion concentration gradients. By relating our results across length scales, we suggest a fairly concentration-independent transference number of about 0.2. Our study shows the broad applicability of operando X-ray photon correlation spectroscopy to the understanding of dynamic phenomena.

25 ENERGY STORAGE↗
DOE OSTI2023

Local cation ordering in compositionally complex Ruddlesden–Popper n = 1 oxides

The Ruddlesden–Popper (RP) layered perovskite structure is of great interest due to its inherent tunability, and the emergence and growth of the compositionally complex oxide (CCO) concept endows the RP family with further possibilities. Here, a comprehensive assessment of thermodynamic stabilization, local order/disorder, and lattice distortion was performed in the first two reported examples of lanthanum-deficient La n+1 B n O 3n+1 (n = 1, B = Mg, Co, Ni, Cu, Zn) obtained via various processing conditions. Chemical short-range order (CSRO) at the B-site and the controllable excess interstitial oxygen (δ) in RP-CCOs are uncovered by neutron pair distribution function analysis. Reverse Monte Carlo analysis of the data, Metropolis Monte Carlo simulations, and extended x-ray absorption fine structure analysis implies a modest degree of magnetic element segregation on the local scale. Further, ab initio molecular dynamics simulations results obtained from special quasirandom structure disagree with experimentally observed CSRO but confirm Jahn–Teller distortion of CuO 6 octahedra. These findings highlight potential opportunities to control local order/disorder and excess interstitial oxygen in layered RP-CCOs and demonstrate a high degree of freedom for tailoring application-specific properties. They also suggest a need for expansion of theoretical and data modeling approaches in order to meet the innate challenges of CCO and related high-entropy phases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems

The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.

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

Influence of acid leaching surface treatment on indentation cracking of soda lime silicate glass

Past work has shown that water or acid soaking treatments can increase the mechanical strength of soda lime silicate (SLS) glasses. In this work, we show that acid leaching treatments result in an increase in the apparent crack resistance of the acid-leached surface of SLS glass during indentation. Vickers indentation tests in controlled environments show a humidity dependence of radial cracking, suggesting that the transport of water through the leached layer plays a critical role in the propagation of cracks to the glass surface. Here, molecular dynamics simulations with reactive force fields indicate that the leached surface layer can undergo pressure-induced mechanochemical reactions during indentation, which increases the bridging oxygen connectivity in the silica network of the leached layer. Such structural changes can hinder transport of water molecules from the environment to the subsurface crack tip. Based on experimental observations and simulation results, a new hypothesis is proposed that mechanochemical restructuring in the leached layer in response to the applied load may lower the transport kinetics of molecular water to critical flaws, resulting in an apparent enhancement in the crack resistance of the acid-leached surface of SLS glass.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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

A Robotic High-Throughput Grid-Search Platform for Mapping Phase Behavior in Triblock Copolymer–Homopolymer Blends

We present a high-throughput experimental investigation of the phase behavior in triblock copolymers (PS-b-PB-b-PS and PS-b-PI-b-PS) and polystyrene (PS) homopolymer blends as a function of homopolymer molecular weight (MW) and blend ratio. Using a robotic thin-film processing platform (NOVA) integrated with Grazing Incidence Small-Angle X-ray Scattering (GISAXS) and Atomic Force Microscopy (AFM), we systematically mapped the order–disorder transition (ODT) boundaries and domain spacing evolution across a broad MW range (4.0–101.3 kDa) with varying homopolymer loadings (10% to 90%). The results reveal three distinct regimes: low-MW homopolymers, corresponding to the wet-brush regime produced only gradual domain swelling before disordering at high blend ratios (weight fraction); medium-MW homopolymers, corresponding to thedry-brush regime induced significant domain spacing increase up to 80% followed by earlier disordering, while high-MW homopolymers led to macrophase separation with minimal changes in domain spacing. Additionally, coarse-grained molecular dynamics simulations confirmed our experimental finding that in the low-MW region, the PS homopolymer uniformly distributed in the PS domain. These findings demonstrate that homopolymer molecular weight critically governs both the extent of domain swelling and the onset of disorder in triblock copolymer systems. This high-throughput platform enables the rapid mapping of composition–morphology relationships and can be integrated with AI/ML tools for designing next-generation nanostructured polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nature of the bonded-to-atomic transition in liquid silica to TPa pressures

First-principles calculations and analysis of the thermodynamic, structural, and electronic properties of liquid SiO 2 characterize the bonded-to-atomic transition at 0.1–1.6 TPa and 10 4 –10 5 K (1–9 eV), the high-energy-density regime relevant to understanding planetary interiors. We find strong ionic bonds that become short-lived due to high kinetics during the transition, with sensitivity of the transition temperature to pressure, and our calculated Hugoniots agree with past experimental data. Furthermore, these results reconcile previous experimental and theoretical findings by clarifying the nature of the bond dissociation process in early Earth and “rocky” (oxide) constituents of large planets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Atomistic simulation of brittle-to-ductile transition in silicon carbide embedded with nano-sized helium bubbles

In this study, the tensile response of cubic silicon carbide (SiC) bulk containing cavities (voids and He bubbles) has been investigated using molecular dynamic simulations. The formation of cavities in SiC leads to a significant degradation in the mechanical properties of SiC with more influence on material fracture than initial elastic deformation. The brittle-to-ductile transition occurs in cavity-embedded SiC as the pressure in He bubbles increases. This is associated with the deformation mechanism that bond breaking at a low He bubble pressure transfers to extensive dislocation activities at a higher He bubble pressure. The cavities can effectively concentrate stress around them in the direction perpendicular to the tension, which leads to preferred cracking in the region with a higher tensile stress. The failure mechanism as revealed by this study improves understanding of property degradation in SiC that may be useful for applications of SiC in advanced nuclear energy systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Scaling Field-Theoretic Simulation for Multicomponent Mixtures with Neural Operators

Multicomponent polymer mixtures are ubiquitous in biological self-organization but are notoriously difficult to study computationally. Plagued by both slow single molecule relaxation times and slow equilibration within dense mixtures, molecular dynamics simulations are typically infeasible at the spatial scales required to study the stability of mesophase structure. Polymer field theories offer an attractive alternative, but analytical calculations are only tractable for mean-field theories and nearby perturbations, constraints that become especially problematic for fluctuation-induced effects such as coacervation. Here, we show that a recently developed technique for obtaining numerical solutions to partial differential equations based on operator learning, neural operators, lends itself to a highly scalable training strategy by parallelizing per-species operator maps. We illustrate the efficacy of our approach on six-component mixtures with randomly selected compositions and that it significantly outperforms the state-of-the-art pseudospectral integrators for field-theoretic simulations, especially as polymer lengths become long.

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

DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations

DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman’s path integral formalism of quantum mechanics. Here, the flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

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

Hydrodynamically Controlled Self‐Organization in Mixtures of Active and Passive Colloids

Abstract Active particles are known to exhibit collective behavior and induce structure in a variety of soft‐matter systems. However, many naturally occurring complex fluids are mixtures of active and passive components. The authors examine how activity induces organization in such multi‐component systems. Mixtures of passive colloids and colloidal micromotors are investigated and it is observed that even a small fraction of active particles induces reorganization of the passive components in an intriguing series of phenomena. Experimental observations are combined with large‐scale simulations that explicitly resolve the near‐ and far‐field effects of the hydrodynamic flow and simultaneously accurately treat the fluid–colloid interfaces. It is demonstrated that neither conventional molecular dynamics simulations nor the reduction of hydrodynamic effects to phoretic attractions can explain the observed phenomena, which originate from the flow field that is generated by the active colloids and subsequently modified by the aggregating passive units. These findings not only offer insight into the organization of biological or synthetic active–passive mixtures, but also open avenues to controlling the behavior of passive building blocks by means of small amounts of active particles.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mapping the Interfacial Chemistry and Structure of Partially Fluorinated Bottlebrush Polymers and Their Linear Analogues

Polymer interfaces are key to a range of applications including membranes for chemical separations, hydrophobic coatings, and passivating layers for antifouling. While important, challenges remain in probing the interfacial monolayer where the molecular ordering and orientation can change depending on the chemical makeup or processing conditions. In this work, we leverage surface specific vibrational sum frequency generation (SFG) and the associated dependence on molecular symmetry to elucidate the ordering and orientations of key functional groups for poly(2,2,2-trifluoroethyl methacrylate) bottlebrush polymers and their linear polymer analogues. These measurements were framed by atomistic molecular dynamic simulations to provide a complementary physical picture of the gas–polymer interface. Simulations and SFG measurements show that methacrylate backbones are buried beneath a layer of trifluoroethyl containing side groups that result in structurally similar interfaces regardless of the polymer molecular weight or architecture. The average orientational angles of the trifluoroethyl containing side groups differ depending on polymer linear and bottlebrush architectures, suggesting that the surface groups can reorient via available rotational degrees of freedom. Results show that the surfaces of the bottlebrush and linear polymer samples do not strongly depend on molecular weight or architecture. As such, one cannot rely on increasing the molecular weight or altering the architecture to tune surface properties. This insight into the polymer interfacial structure is expected to advance the design of new material interfaces with tailored chemical/functional properties.

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

Evaluating Cryo–TEM Reconstruction Accuracy of Self–Assembled Polymer Nanostructures

Cryogenic transmission electron microscopy (cryo–TEM) combined with single particle analysis (SPA) is an emerging imaging approach for soft materials. However, the accuracy of SPA–reconstructed nanostructures, particularly those formed by synthetic polymers, remains uncertain due to potential packing heterogeneity of the nanostructures. In this study, the combination of molecular dynamics (MD) simulations and image simulations is utilized to validate the accuracy of cryo–TEM 3D reconstructions of self–assembled polypeptoid fibril nanostructures. Using CryoSPARC software, image simulations, 2D classifications, ab initio reconstructions, and homogenous refinements are performed. By comparing the results with atomic models, the recovery of molecular details is assessed, heterogeneous structures are identified, and the influence of extraction location on the reconstructions is evaluated. In conclusion, these findings confirm the fidelity of single particle analysis in accurately resolving complex structural characteristics and heterogeneous structures, exhibiting its potential as a valuable tool for detailed structural analysis of synthetic polymers and soft materials.

36 MATERIALS SCIENCE↗