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 OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2023

Intergranular Hotspots: A Molecular Dynamics Study on the Influence of Compressive and Shear Work

Numerous crystal- and microstructural-level mechanisms are at play in the formation of hotspots, which are known to govern high explosives initiation behavior. Most of these mechanisms, including pore collapse, interfacial friction, and shear banding, involve both compressive and shear work done within the material and have thus far remained difficult to separate. Here, we assess hotspots formed at shocked crystal–crystal interfaces using quasi-1D molecular dynamics simulations that isolate effects due to compression and shear. Two high explosive materials are considered (TATB and PETN) that exhibit distinctly different levels of molecular conformational flexibility and crystal packing anisotropy. Temperature and intramolecular strain energy localization in the hotspot are assessed through parametric variation of the crystal orientation and two velocity components that respectively modulate compression and shear work. The resulting hotspots are found to be highly localized to a region within 5–20 nm of the crystal–crystal interface. Compressive work plays a considerably larger role in localizing temperature and intramolecular strain energy for both materials and all crystal orientations considered. Shear induces a moderate increase in energy localization relative to unsheared cases only for relatively weak compressive shock pressures of approximately 10 GPa. These results help isolate and rank the relative importance of hotspot generation mechanisms and are anticipated to guide the treatment of crystal–crystal interfaces in coarse-grained models of polycrystalline high-explosive materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗