Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RADIAL DISTRIBUTION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Cost-effective carbon fiber precursor selections of polyacrylonitrile-derived blend polymers: carbonization chemistry and structural characterizations

Blending polyacrylonitrile (PAN) with plastic wastes and bio-based polymers provides a convenient and inexpensive method to realize cost-effective carbon fiber (CF) precursors. In this work, PAN-based blend precursors are investigated using ReaxFF reactive molecular dynamics simulations with respect to the formation of all-carbon rings, the evolutions of oxygen-containing and nitrogen-containing species, and the migration of carbon atoms to form turbostratic graphene. From these simulations, we identify that PAN/cellulose (CL) blend manifests the highest carbon yield and the most substantial all-carbon ring formation. This ReaxFF-based finding is confirmed by Raman and TEM experiments indicating high crystallinity for PAN/CL-derived blend CFs. We trace the pathway of gasification and carbonization of PAN/CL to elaborate the mechanism of the formation of all-carbon ring networks. We discover that the acetals of CL can catalyze the cyclization of the blend precursor, allowing for the search for CL derivatives or the other kinds of bio-based polymers with similar functionalities as alternative blends. In addition, we examine the structural characteristics using the carbon–carbon (C–C) radial distribution functions, C–C bond length distributions, and sp2 C atom ratios for the four representative precursors, i.e., PAN, oxidized PAN, PAN/nylon 6,6, and PAN/CL. Our simulation results show the most extensive all-carbon ring cluster and graphitic structure growths for PAN/CL. Here, we propose PAN/CL as a cost-effective alternative CF precursor, since (a) CL is naturally abundant and eco-friendly for production, (b) the blend precursor PAN/CL does not require oxidation treatment, (c) PAN/CL has a high carbon yield with substantial all-carbon ring formation, and (d) PAN/CL based CFs potentially provide a mechanical property enhancement.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Feasibility study of Mg storage in a bilayer silicene anode via application of an external electric field

With the goal of developing a Si-based anode for Mg-ion batteries (MIBs) that is both efficient and compatible with the current semiconductor industry, the current research utilized classical Molecular Dynamics (MD) simulation in investigating the intercalation of a Mg 2+ ion under an external electric field (E-field) in a 2D bilayer silicene anode (BSA). First principles density functional theory calculations were used to validate the implemented EDIP potentials. Our simulation shows that there exists an optimum E-field value in the range of 0.2–0.4 V Å –1 for Mg 2+ intercalation in BSA. To study the effect of the E-field on Mg 2+ ions, an exhaustive spread of investigations was carried out under different boundary conditions, including calculations of mean square displacement (MSD), interaction energy, radial distribution function (RDF), and trajectory of ions. Our results show that the Mg 2+ ions form a stable bond with Si in BSA. The effects of E-field direction and operating temperature were also investigated. In the X–Y plane in the 0°–45° range, 15° from the X-direction was found to be the optimum direction for intercalation. The results of this work also suggest that BSA does not undergo drastic structural changes during the charging cycles with the highest operating temperature being ~300 K

25 ENERGY STORAGE↗
DOE OSTI2023

Computational insights into the interaction of water with the UiO-66 metal–organic framework and its functionalized derivatives

The UiO-66 metal–organic framework (MOF) has been identified as a promising hydrophilic material for water harvesting. Recent studies show that its water uptake ability at low relative humidity (RH) can be improved by incorporating hydrophilic functional groups into the framework. In this work, we provide computational insights into the adsorption of water in UiO-66 and its functionalized derivatives to reveal the role played by different adsorption sites and functional groups in the adsorption mechanism. We started by developing molecular models for UiO-66, UiO-66-NH 2 , UiO-66-OH, and UiO-66-(OH) 2 compatible with the MB-pol data-driven many-body potential of water. We then benchmarked these models against ab initio data. Here we used these models to perform molecular dynamics simulations and calculate radial distribution functions, IR spectra, and two-dimensional density distribution maps for water in the MOFs. These results consistently show that the μ 3 -OH sites are the preferential interaction sites for water in UiO-66 and all its variants, and the formation of localised water clusters inside the octahedral pores is responsible for the abrupt step in the experimental adsorption isotherms. Furthermore, the presence of functional groups in the framework allows water to cluster in the octahedral pores at lower RH, thus making the MOF a more efficient water harvester. Overall, this study provides molecular-level insights into the pore filling process of UiO-66 and its functionalized derivatives, which are needed for the design of efficient water harvesting materials based on MOFs.

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

Multiscale investigation of the mechanism of biomass deconstruction in the dimethyl isosorbide/water Co-solvent pretreatment system

In the context of promoting a circular bioeconomy, the development of green and efficient lignocellulosic biomass pretreatment technologies so as to realize high value-added biomass utilization is of intense interest. We demonstrated the potential of the bio-based green solvent dimethyl isosorbide (DMI) for the fractionation of Eucalyptus biomass with excellent performance. Here, to investigate the mechanisms involved in biomass fractionation, microimaging and microspectroscopic techniques were employed together with molecular dynamics (MD) simulation and COSMO-RS quantum chemical calculations to derive multiscale information. Both the microstructure and regional chemistry of the cell wall vary significantly with the volume ratio of DMI/H 2 O. The strongest effects were found at DMI/H 2 O = 9:1 and showed visible cell wall tearing cracks and cell wall deformation and collapse as well as the lowest values of cell wall thickness and circularity. From the MD simulations, lignin exhibits collapsed-like structure in pure H 2 O with low solvent accessibility surface area (SASA) and radius of gyration (R g ). In contrast, lignin in DMI/H 2 O shows extended structure with high SASA and solvent interactions dominated by van der Waals forces, with maximal contact in the 9:1 (v/v) system. Further, the COSMO-RS calculated sigma (σ-) potential suggests the intermolecular interactions in DMI and DMI/H 2 O co-solvent are weak, leading to stronger interaction with lignin and correspondingly higher lignin dissolution. The radial distribution functions and σ-potential all show that again DMI/H 2 O at 9:1 is an optimal volume ratio for high lignin dissolution. Finally, this study provides a solvent-ratio dependent mechanism for the action of polar aprotic solvents in the deconstruction of biomass.

09 BIOMASS FUELS↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

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

Liquid state properties of SEI components in dimethoxyethane

The solid-electrolyte interphase (SEI) layer is a critical constituent of battery technology, which incorporates the use of lithium metals. Since the formation of the SEI is difficult to avoid, the engineering and harnessing of the SEI are absolutely critical to advancing energy storage. One problem is that much fundamental information about SEI properties is lacking due to the difficulty in probing a chemically complex interfacial system. One such property that is currently unknown is the dissolution of the SEI. This process can have significant effects on the stability of the SEI, which is critical to battery performance but is difficult to probe experimentally. Here, we report the use of ab initio computational chemistry simulations to probe the solution state properties of SEI components LiF, Li 2 O, LiOH, and Li 2 CO 3 in order to study their dissolution and other solution-based characteristics. Ab initio molecular dynamics was used to study the solvation structures of the SEI with a combination of radial distribution functions, discrete solvation structure maps, and vibrational density of states, which allows for the determination of free energies. From the change in free energy of dissolution, we determined that LiOH is the most likely component to dissolve in the electrolyte followed by LiF, Li 2 CO 3 , and Li 2 O although none were favored thermodynamically. Furthermore, this indicates that dissolution is not probable, but Li 2 O would make the most stable SEI with regard to dissolution in the electrolyte.

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

Structural heterogeneity in non-crystalline Te x Se1−x thin films

Rapid crystallization behavior of amorphous TexSe1−x thin films limits the use of these alloys as coatings and in optoelectronic devices. Understanding the short- and medium-range ordering of the amorphous structure and the fundamental physics governing the crystallization of the films is crucial. Although the lack of long range crystalline order restricts the characterization of the amorphous films, electron microscopy offers a way to extract information about the nanoscale ordering. In this paper, the local ordering of amorphous TexSe1−x thin films with x=0.22, 0.61, 0.70, 0.90, and 1 grown by thermal evaporation is investigated using radial distribution function (RDF) and fluctuation electron microscopy (FEM) analysis. RDF results show that the nearest-neighbor distances of selenium (Se) and tellurium (Te) in their crystalline structure are preserved, and their bond lengths increase with the addition of Te. Density functional theory (DFT) calculations predict structures with interatomic distances similar to those measured experimentally. Additionally, fluctuations in atomic coordination are analyzed. Medium range order (MRO) analysis obtained from FEM and DFT calculations suggests that there are at least two populations within the chain network structure, which are close to the Se–Se and Te–Te intrachain distances. For the binary alloy with x > 0.61, TexSe1−x, Te–Te like populations increase and Te fragments might form, suggesting that the glass forming ability decreases rapidly.

Sari, Bengisu (ORCID:0000000244217098)↗
DOE OSTI2024

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Chemistry↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

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

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

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

Elucidating ion capture and transport mechanisms of Preyssler anions in aqueous solutions using biased MACE-accelerated MD simulations

Equilibrium and biased multi-atomic cluster expansion (MACE) accelerated molecular dynamics (MD) simulations in aqueous solutions are performed to investigate the ion capture and transport mechanisms of the {P 5 W 30 } Preyssler anion (PA) as the smallest representative member of the extended polyoxometalate (POM) family with an internal cavity. The unique interatomic interactions present in the internal cavity vs the exterior of PA are carefully investigated using equilibrium MACE MD simulations for two representative Na(H 2 O)@PA and Na@PA complexes in aqueous solutions. Our careful analyses of radial distribution functions and coordination numbers show that the presence of confined water in Na(H 2 O)@PA has profound modulating effects on the nature of the interactions of the encapsulated ion with the oxygens of the PA cavity. Using well-converged MACE-accelerated multiple walker well-tempered metadynamics simulations with nanosecond timescales, two different associative (ion exchange) and dissociative (ion ejection) ion transport mechanisms were carefully investigated for Na + as one of the most abundant and representative ions present in seawater and saline solutions. By comparing systems with and without confined water, it was found that the presence of only one pre-encapsulated confined water in Na(H 2 O)@PA dramatically changes the free energy landscape of ion transport processes. It was also found that the contraction and dilation of the two windows present in PA directly influence the Na + and H 2 O transport. Furthermore, the results from this work are helpful, as they show a viable path toward tuning the ion exchange and transport phenomena in aqueous solutions of POM molecular clusters and frameworks.

25 ENERGY STORAGE↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2023

In-situ study of rules of nanostructure evolution, severe plastic deformations, and friction under high pressure

Severe plastic deformations under high pressure are used to produce nanostructured materials but were studied ex-situ. Rough diamond anvils are introduced to reach maximum friction equal to yield strength in shear and the first in-situ study of the evolution of the pressure-dependent yield strength and radial distribution of nanostructural parameters are performed for severely pre-deformed Zr. ω-Zr behaves like perfectly plastic, isotropic, and strain-path-independent and reaches steady values of the crystallite size and dislocation density, which are pressure-, strain- and strain-path-independent. However, steady states for a-Zr obtained with smooth and rough anvils are different, causing major challenge in plasticity theory. In-situ study of severe plastic deformation of α-Zr with rough diamond anvils revealed that pressure-dependent yield strength, crystallite size, and dislocation density are getting steady and plastic strain- and strain-path-independent.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Towards Mitigating Electron Beam-Induced Damage in MOFs Using Low Dose Electron Microscopy for In-Situ Room-Temperature Measurements

Metal-Organic Frameworks (MOFs) have emerged as a versatile class of materials with applications in gas storage and separation, catalysis, and drug delivery [1]. Understanding their damage mechanisms in under various environmental conditions is crucial for optimizing performance and stability [2]. This study employs a combination of Length Structural Coherence (LSC) analysis in electron scattering and Electron Energy Loss Spectroscopy (EELS) to elucidate the damage mechanisms in two representative MOFs. MIL101(Fe), and MIL101(Cr). Here, the primary objective is to reveal structural and chemical changes that occur under electron beam irradiation in Scanning Transmission Electron Microscopy (STEM), both in the pristine state and during in-situ CO and CO 2 adsorption. By combining LSC in Radial Distribution Function (RDF) analysis with EELS, comprehensive insights into the nanoscale alterations are obtained and correlated to the electron dose rate.

dos Santos, Gabriel T. [Northwestern University, E↗
DOE OSTI2026

Machine learning-based interatomic potential development and phase transition analysis of ferroelectric hafnium dioxide

The ferroelectric phase (𝑃⁢𝑐⁢𝑎⁢2 1 , which is in orthorhombic symmetry) of hafnium dioxide (HfO 2 ) has gained much attention due to its potential applications in nanoelectronics and advanced memory devices. However, its complex phase behavior under external stimuli, such as pressure and temperature, remains a subject of intense investigation. This study focuses on developing a machine learning-based interatomic potential (MLIP) that is trained with data from density-functional theory (DFT) calculations to simulate phase transitions and mechanical properties of HfO 2 . The developed MLIP predicts lattice parameters, equations of state, bulk and shear moduli, and elastic constants that closely align with DFT predictions for several phases and at various pressures. Once validated, the MLIP is used to investigate the phase transitions of ferroelectric HfO 2 (𝑃⁢𝑐⁢𝑎⁢2 1 ) under both isobaric and constant stress conditions at elevated temperatures ranging from 200 to 2500 K. We used several complementary methods, including local symmetry identification, radial distribution function, and x-ray diffraction characterization, to identify interesting phase transitions among several competitive hafnia phases predicted from our simulations. The suggested methods uniformly reveal that under pure deviatoric condition, the system favors a transition from the orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase to a tetragonal (𝑃⁢4 2 /𝑛⁢𝑚⁢𝑐) phase, whereas a zero stress condition drives the system from the 𝑃⁢𝑐⁢𝑎⁢2 1 phase to another orthorhombic (𝑃⁢𝑏⁢𝑐⁢𝑛) phase. These findings provide crucial insights into stress and temperature-induced phase behavior of hafnia, guiding future experimental and theoretical studies for optimizing hafnia-based ferroelectric devices.

Ferroelectric HfO2↗