Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular liquids”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Chiral Electrokinetic Phenomena in Single Nanopores

Abstract The arrangement of solvent molecules and ions at solid–liquid interfaces determines electrochemical properties that are important in separations platforms, sensing technologies, and energy‐storage systems. Here we show that single glass and polymer pores in contact with propylene carbonate (PC) solutions of LiClO 4 exhibit an effective surface potential that is modulated by the enantiomeric excess of the solvent. In particular, electrochemical and electrokinetic measurements of ionic transport through glass pipettes and polymer pores reveal that the effective surface potential is significantly lower in solutions prepared using enantiomerically pure PC than in solutions prepared using racemic PC. Both pore systems became positively charged in all racemic solutions examined in the range of LiClO 4 concentrations between 1 mM and 100 mM, whereas solutions in ( R )‐(+)‐PC induced a positive surface potential only at concentrations above ~5 mM. The effective surface potential is quantified through asymmetry in current–voltage curves and zeta‐potential measurements. Vibrational sum‐frequency‐generation experiments on LiClO 4 solutions in racemic and enantiomerically pure PC indicate that the surface lipid‐bilayer‐like region in the former is more strongly organized than in the latter, dictating the favorable positions for lithium and perchlorate ions in each case. The more ordered molecular packing in the racemic liquid leads to accumulation of lithium ions on the outside of the bilayer, creating a higher effective positive charge. Our results highlight the extreme sensitivity of the interfacial potential on molecular organization of the solvent, and the relatively unexplored role that chirality can play in electrokinetic phenomena.

Alanis, Kristen↗
DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic dynamics in fluids: Normal mode analysis revisited

Developing microscopic understanding of the thermal properties of liquids is challenging due to their strong dynamic disorder, which prevents characterization of the atomic degrees of freedom. There have been significant research interests in the past few decades to extend the normal mode analysis for solids to instantaneous structures of liquids. However, the nature of normal modes that arise from these unstable structures is still elusive. Here, in this paper, we explore the instantaneous eigenmodes of dynamical matrices of various Lennard-Jones argon liquid and gas systems at high temperatures and show that the normal modes can be interpreted as an interpolation of T →∞ (gas) and T=0 (solid) mode descriptions. We find that normal modes become increasingly collisional and translational, recovering atomistic gaslike behavior rather than vibrational with increase in temperature, suggesting that normal modes in liquids may be described by both solidlike and gaslike modes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Thermoosmotic Transport in Nanochannels Grafted with pH-responsive Polyelectrolyte Brushes Modelled Using Augmented Strong Stretching Theory

In this paper, we develop a theory to establish that the thermoosmotic (TOS) effects, induced by the application of an axial temperature gradient, leads to a massive enhancement in liquid transport in nanochannels grafted with charged polyelectrolyte (PE) brushes. Here, we quantify the TOS transport by quantifying the induced electric field and the induced TOS flow field. The different components of the electric field, namely the ionic component, the thermal component, and the osmotic component, as well as the contributions of different ions to these components, are quantified. Furthermore, we express the TOS velocity as a combination of chemiosmotic (COS), thermal, and electroosmotic (EOS) components. The COS and the thermal components augment each other and the overall strength and direction of the TOS flow are dictated by the direction and the relative strength of the EOS component. Most importantly, we compare the cases of brush-grafted nanochannels with those of the brush-free nanochannels of identical surface charge densities: the TOS transport is massively augmented in the brush-grafted nanochannels attributed to the combination of the localization of the EDL (and hence any body force that depends on the EDL charge density) away from the nanochannel wall (i.e., the location of the maximum drag force) and the presence of a possible molecular slip (experienced by the liquid) along the brush surface.

42 ENGINEERING↗
DOE OSTI2023

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Toward sustainable production of N-containing products via nonthermal plasma-enhanced conversion of natural gas resources

Nonthermal plasmas can directly activate and cleave the strong chemical bonds in molecular nitrogen and methane to facilitate the transformation of these inherently stable molecules through the application of electrical energy. Here, we report a low temperature, atmospheric pressure, nonthermal plasma for the “one-pot” synthesis of olefins, alkynes, higher molecular weight hydrocarbons, ammonia, and nitrogen-containing liquids from a representative shale gas feed enriched with nitrogen. Furthermore, we reproducibly observe a wide range of valuable and synthetically challenging gas-phase and liquid-phase products containing C-N, C-C, and N-H bonding by controlling the N 2 concentration in the inlet feed stream. In nitrogen-lean regimes, hydrocarbon products dominate (e.g., ethylene, acetylene, etc.), while nitrogen-rich regimes promote incorporation of nitrogen into the products, leading to the formation of ammonia and liquid products containing a variety of functionalities (e.g., nitriles, amines, heterocycles). High resolution electrospray ionization mass spectrometry was used to measure molecular weights and identify the chemical formulas of the liquid products. Van Krevelen diagrams were created and showed many products with compositions around H/C = 2 and N/C = 0.5, indicating the potential importance of intermediate species with these ratios for liquid formation (e.g., CH 3 CN + H).

03 NATURAL GAS↗
DOE OSTI2022

Deep eutectic solvent-based polymer electrolyte for solid-state lithium metal batteries

Poly(ethylene) oxide (PEO)-based electrolytes have been widely studied for solid-state lithium batteries while their ionic conductivity and lithium-ion transference number still need to be further improved. Herein, using the combined experimental and theoretical approach, we demonstrate a novel, solid-state PEO-deep eutectic solvent (DES) electrolyte for the first time. We found that the in situ formation of DES can reduce the crystallinity of PEO matrix and more Li + ions can move freely owing to the weakened coordination between ether oxygens and Li-ions. Besides, we show that more Li + ions can be dissociated from Li salts in PEO-DES electrolyte using the molecular dynamics simulations. Such liquid-free PEO-DES electrolytes showed good ionic conductivity (2.1 × 10 −4 S cm −1 ) which is 160% higher than that of conventional PEO-LiTFSI (8.1 × 10 −5 S cm −1 ) electrolyte at 60 °C. Additionally, the PEO-DES electrolyte showed 136% increase of Li-ion transference number (0.33) compared with ionic liquid-doped PEO-LiTFSI (0.14) at 60 °C. Moreover, the PEO-DES exhibited good compatibility with Li metal and stable Li plating/stripping behavior with little morphology change of Li metal. In conclusion, this research also provides new insights into the enhancement mechanisms of novel polymer electrolytes, improving our fundamental understanding of critical challenges that have impeded the adoption of solid-state lithium metal batteries.

Deep eutectic solvent↗
DOE OSTI2020

Solvation Energy of Ions in a Stockmayer Fluid

We calculate the solvation energy of monovalent and divalent ions in various liquids with coarse-grained molecular dynamics simulations. Our theory treats the solvent as a Stockmayer fluid, which accounts for the intrinsic dipole moment of molecules and the rotational dynamics of the dipoles. Despite the simplicity of the model, we obtain qualitative agreement between the simulations and experimental data for the free energy and enthalpy of ion solvation, which indicates that the primary contribution to the solvation energy arises mainly from the first and possibly second solvation shells near the ions. Our results suggest that a Stockmayer fluid can serve as a reference model that enables direct comparison between theory and experiment and may be invoked to scale up electrostatic interactions from the atomic to the molecular length scale.

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

Tuning aminopolycarboxylate chelators for efficient complexation of trivalent actinides

The complexation of trivalent lanthanides and minor actinides (Am 3+ , Cm 3+ , and Cf 3+ ) by the acyclic aminopolycarboxylate chelators 6,6'-((ethane-1,2-diylbis–((carboxymethyl)azanediyl))bis–(methylene))dipicolinic acid (H 4 octapa) and 6,6'-((((4-(1-(2-(2-(2-hydroxyethoxy)ethoxy)ethyl)-1H-1,2,3-triazol-4-yl)pyridine-2,6-diyl)bis–(methylene))bis–((carboxymethyl)azanediyl))bis–(methylene)) dipicolinic acid (H 4 pypa-peg) were studied using potentiometry, spectroscopy, competitive complexation liquid–liquid extraction, and ab initio molecular dynamics simulations. Two studied reagents are strong multidentate chelators, well-suited for applications seeking radiometal coordination for in-vivo delivery and f-element isolation. The previously reported H 4 octapa forms a compact coordination packet, while H 4 pypa-peg is less sterically constrained due to the presence of central pyridine ring. The solubility of H 4 octapa is limited in a non-complexing high ionic strength perchlorate media. However, the introduction of a polyethylene glycol group in H 4 pypa-peg increased the solubility without influencing its ability to complex the lanthanides and minor actinides in solution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Hydration structure of flat and stepped MgO surfaces

We investigate the solvation structure of flat and stepped MgO(001) in neutral liquid water using ab initio molecular dynamics based on a hybrid density functional with dispersion corrections. Our simulations show that the MgO surface is covered by a densely packed layer of mixed intact and dissociated adsorbed water molecules in a planar arrangement with strong intermolecular H-bonds. The water dissociation fractions in this layer are >20% and >30% on the flat and stepped surfaces, respectively. Slightly above the first water layer, we observe metastable OH groups perpendicular to the interface, similar to those reported in low temperature studies of water monolayers on MgO. Here, these species receive hydrogen bonds from four nearby water molecules in the first layer and have their hydrophobic H end directed toward bulk water, while their associated protons are bound to surface oxygens. The formation of these OH species is attributed to the strong basicity of the MgO surface and can be relevant for understanding various phenomena from morphology evolution and growth of (nano)crystalline MgO particles to heterogeneous catalysis.

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

Efficient real space formalism for hybrid density functionals

We present an efficient real space formalism for hybrid exchange-correlation functionals in generalized Kohn–Sham density functional theory (DFT). In particular, we develop an efficient representation for any function of the real space finite-difference Laplacian matrix by leveraging its Kronecker product structure, thereby enabling the time to solution of associated linear systems to be highly competitive with the fast Fourier transform scheme while not imposing any restrictions on the boundary conditions. We implement this formalism for both the unscreened and range-separated variants of hybrid functionals. We verify its accuracy and efficiency through comparisons with established planewave codes for isolated as well as bulk systems. In particular, we demonstrate up to an order-of-magnitude speedup in time to solution for the real space method. We also apply the framework to study the structure of liquid water using ab initio molecular dynamics, where we find good agreement with the literature. Overall, the current formalism provides an avenue for efficient real-space DFT calculations with hybrid density functionals.

Chemistry↗
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2022

Multiphase equation of state for Ta 2 O 5

A new equation of state for Ta 2 O 5 is presented. The EOS is constructed using the OpenSesame software and is referred to as SESAME 3530. The EOS uses a combination of density functional theory (DFT) calculations and experimental data. DFT calculations include cold curves and phonons of the solid phases, as well as DFT-based molecular dynamics simulations of the liquid phase. Experimental data includes isobaric, diamond anvil cell, and porous shock Hugoniot data. To fit the data, we create a multiphase EOS consisting of two solid phases and the liquid. Overall agreement with experimental data is shown, and we provide some suggestions for future experiments that could improve our knowledge of the phase diagram.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Biomolecular Films for Direct Air Capture of CO 2

Efficient utilization of CO 2 is amongst the most critical cost drivers in algal biomass production in open pond systems. CO 2 delivery costs represent approximately 20% of the final biomass selling price in algal mass cultivation systems. Technologies that enable direct air capture (DAC) of atmospheric CO 2 to decouple algae cultivation from CO 2 point sources thus present an opportunity to improve the economics and resource potential of algal biomass. Current DAC technologies typically employ amine- or caustic-based absorption, demanding significant water and/or energy inputs and incurring substantial capital expenditures. Conversely, bio-based approaches to DAC offer a means to bypass conventional technoeconomic and sustainability hurdles. We integrate recent advances in computational metabolic modeling, algal genetic engineering, algal cultivation, and algal biomass upgrading to enable directed localization and self-assembly of carbonic anhydrase molecular films to gas-liquid interfaces for enhanced CO 2 capture and conversion.

09 BIOMASS FUELS↗
DOE OSTI2020

The Terahertz Dynamics of an Aqueous Nanoparticle Suspension: An Inelastic X-ray Scattering Study

We used the high-resolution Inelastic X-ray Scattering beamline of the Advanced Photon Source at Argonne National Laboratory to measure the terahertz spectrum of pure water and a dilute aqueous suspension of 15 nm diameter spherical Au nanoparticles (Au-NPs). We observe that, despite their sparse volume concentration of about 0.5%, the immersed NPs strongly influence the collective molecular dynamics of the hosting liquid. We investigate this effect through a Bayesian inference analysis of the spectral lineshape, which elucidates how terahertz transport properties of water change upon Au-NP immersion. In particular, we observe a nearly complete disappearance of the longitudinal acoustic mode and a mildly decreased ability to support shear wave propagation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗