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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Formation of a single quasicrystal upon collision of multiple grains

Abstract Quasicrystals exhibit long-range order but lack translational symmetry. When grown as single crystals, they possess distinctive and unusual properties owing to the absence of grain boundaries. Unfortunately, conventional methods such as bulk crystal growth or thin film deposition only allow us to synthesize either polycrystalline quasicrystals or quasicrystals that are at most a few centimeters in size. Here, we reveal through real-time and 3D imaging the formation of a single decagonal quasicrystal arising from a hard collision between multiple growing quasicrystals in an Al-Co-Ni liquid. Through corresponding molecular dynamics simulations, we examine the underlying kinetics of quasicrystal coalescence and investigate the effects of initial misorientation between the growing quasicrystalline grains on the formation of grain boundaries. At small misorientation, coalescence occurs following rigid rotation that is facilitated by phasons. Our joint experimental-computational discovery paves the way toward fabrication of single, large-scale quasicrystals for novel applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

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

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

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

Solid Electrolyte Interphase Formation at the Ionic Liquid Electrolyte—Lithium-Metal Interface Using an Ab Initio Molecular Dynamics Approach

An interfacial study is performed using ab initio molecular dynamics (AIMD) simulations to elucidate the electrochemical interfacial stability of an ionic liquid (IL) as component of an electrolyte in contact with a Li-metal anode when an additive, vinylene carbonate (VC), is added to the electrolyte. The IL electrolyte is composed of 1-butyl-1-methylpyrrolidinium bis(fluoro-sulfonyl)-imide (PYR14 + FSI − ) with lithium bis(trifluoromethane-sulfonyl)imide salt (Li + TFSI − ) in a 8:2 concentration ratio. Reactions observed at the interface electrolyte/Li-metal during the AIMD simulations are verified by electronic structure calculations using density functional theory calculations considering a continuum solvent model. We found that the addition of VC to the IL prevents further reduction of the counter-ion TFSI − and retards reduction of FSI − in comparison to the case when the additive VC is absent. The PYR14 + reduction is never observed with or without the additive. The 1st, 2nd, and 3rd electron affinities of the PYR14 + are −0.08, −0.36 and −3.65 eV, indicating little to large repulsive behavior against electrons. The opposite is observed for the two negative ions, FSI − and TFSI − , which yield 1st electron affinities of 1.73 and 1.83 eV, respectively. Reactions obtained under open circuit conditions. Dissociation mechanisms are presented for all IL components reduced at the IL/Li-metal interface.

Electrochemistry↗
DOE OSTI2019

Liquid-like and rigid-body motions in molecular-dynamics simulations of a crystalline protein

To gain insight into crystalline protein dynamics, we performed molecular-dynamics (MD) simulations of a periodic 2 × 2 × 2 supercell of staphylococcal nuclease. We used the resulting MD trajectories to simulate X-ray diffraction and to study collective motions. The agreement of simulated X-ray diffraction with the data is comparable to previous MD simulation studies. We studied collective motions by analyzing statistically the covariance of alpha-carbon position displacements. The covariance decreases exponentially with the distance between atoms, which is consistent with a liquidlike motions (LLM) model, in which the protein behaves like a soft material. To gain finer insight into the collective motions, we examined the covariance behavior within a protein molecule (intraprotein) and between different protein molecules (interprotein). The interprotein atom pairs, which dominate the overall statistics, exhibit LLM behavior; however, the intraprotein pairs exhibit behavior that is consistent with a superposition of LLM and rigid-body motions (RBM). Our results indicate that LLM behavior of global dynamics is present in MD simulations of a protein crystal. They also show that RBM behavior is detectable in the simulations but that it is subsumed by the LLM behavior. Finally, the results provide clues about how correlated motions of atom pairs both within and across proteins might manifest in diffraction data. Overall, our findings increase our understanding of the connection between molecular motions and diffraction data and therefore advance efforts to extract information about functionally important motions from crystallography experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Atomistic simulations of the thermal conductivity of liquids

We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Polyampholyte physics: Liquid–liquid phase separation and biological condensates

Here we review the current understanding of molecular interactions that govern liquid-liquid phase separation (LLPS) of biological condensates. The connection between sequence, chain conformation, and phase separation of intrinsically disordered proteins (IDPs) and their model polyampholytes is discussed. In particular, we highlight how the charge pattern influences the conformation and phase behavior of natural IDPs. We then describe recent results from theoretical treatments of polyampholytes implementing random phase approximation, field-theoretic simulations, and transfer matrix theory that show an increase in charge segregation results in an increased tendency to phase separate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Room Temperature to 150 °C Lithium Metal Batteries Enabled by a Rigid Molecular Ionic Composite Electrolyte

Molecular ionic composites (MICs), made from ionic liquids and a rigid-rod polymer poly(2,2'-disulfonyl-4,4'-benzidine terephthalamide) (PBDT), are a new type of rigid gel electrolyte that combine fast ion transport with high thermal stability and mechanical strength. In this work, a MIC electrolyte membrane is prepared that is composed of PBDT, lithium bis(trifluoromethylsulfonyl)imide (LiTFSI), and 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (Pyr 14 TFSI) in a mass ratio of 10:10:80. The ionic conductivity at 25 °C is 0.56 mS cm –1 with no added flammable/volatile components. Although the polymer content is only 10 wt%, this MIC membrane is rigid with a tensile modulus of 410 MPa at room temperature. The MIC membrane remains stable and rigid at 200 °C with the shear storage modulus (G') only slightly decreasing by 35%. Li/MIC/LiFePO 4 cells demonstrate stable cycling performance over a wide temperature range from 23 to 150 °C. The specific discharge capacity at 100 and 150 °C at 1 C rate exceeds 160 mAh g –1 . The discharge capacity retention is 99% after 50 cycles at 150 °C. Furthermore, this stable battery performance shows that this low polymer content MIC membrane qualifies for use as a solid electrolyte in lithium metal batteries operating over a wide temperature range.

25 ENERGY STORAGE↗
DOE OSTI2026

Cyclodextrin-Derived Porous Liquids Enabled by In Situ Solvation Shell Formation

Porous liquids (PLs) represent a unique platform for molecular separations by combining permanent porosity with liquid-phase mobility. However, it remains a formidable challenge to construct and stabilize PLs with sub-5 Å pores using readily available porous host and liquid media. Here, we report the construction of cyclodextrin (CD)-derived PLs enabled by in situ solvation shell formation. The acid–base neutralization reaction between CD and an organic base was leveraged to generate a thin ionic solvation shell around the CD host, effectively liquefying CD and preventing its segregation in the liquid base medium while preserving accessible molecular-scale cavities. Spectroscopic analysis, neutron scattering, density functional theory calculations, and molecular dynamics simulations collectively confirm the structural evolution and existence of abundant internal porosity in PLs. The unique architectures of CD-derived PLs enable highly selective encapsulation of fluorinated alkanes and significantly enhanced uptake of inert gases. This facile and generalizable strategy enables construction of high-quality PLs with engineered ultramicroporosity to facilitate molecular separations.

Cavities↗
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

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

Thermophysical Properties of Liquid Tritium: A Path Integral Monte Carlo Study

Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

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

Entropy-Driven Porous Liquids Allowing Gas Solubility in Solvent-Filled Imine-Based Porous Organic Cages

Porous liquids offer a promising platform for gas separation by combining fluid processability with intrinsic molecular porosity. Traditional Type II porous liquids are formed by dissolving porous molecular cages in size-excluded solvents, limiting solvent options and practical applications. In this work, we introduce a novel method of creating Type II porous liquids using common small solvents, where intrinsic porosity is achieved at elevated pressures due to the selective displacement of solvent molecules by gas molecules within the cage structures. Using molecular simulations, we investigate the behavior of CO 2 in solutions of the imine-based porous organic cage CC13 dissolved in small molecular solvents such as chloroform and 1,2-dimethoxyethane (DME). Grand canonical Monte Carlo simulations of solid-state CC13 reveal that selectivity reversal, where smaller CO 2 molecules displace larger solvent molecules inside the cage, is achievable at sufficiently high pressures. Temperature quench molecular dynamics simulations confirm that while CO 2 displacement within chloroform-filled cages is limited, DME enables entropy-driven cage CO 2 occupancy at pressures as low as ∼23 bar, setting up the foundation of an alternative way of forming Type II porous liquids.

adsorption↗
DOE OSTI2023

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density–density and hydrostatic pressure–pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. Furthermore, these results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

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

Origin of short- and medium-range order in supercooled liquid Ge 3 Sb 2 Te 6 from ab initio molecular dynamics simulations

Phase-change materials such as Ge–Sb–Te compounds have attracted much attention due to their potential value in electrical data storage. In contrast to the amorphous and crystalline phases, supercooled liquids are far from being deeply understood despite their inevitable role in both amorphization and crystallization processes. To this end, we have studied the dynamics properties and structural characteristics of liquid and supercooled liquid Ge 3 Sb 2 Te 6 during the fast cooling process. As the temperature decreases, chemical bonds become more homogeneous, but coordination numbers of Ge, Sb and Te atoms change very little. Meanwhile, the structural order of short-range configuration is obviously enhanced. Further studies suggest that Ge-centered, Sb-centered and Te-centered configurations change to the more ordered defective octahedrons mainly by adjusting the bond-angle relationship and bond length, rather than just by changing the coordination environment. It is the more ordered octahedrons that promote the formation of medium-range order. Our findings provide a deep insight into the origin of local structural order in supercooled liquid Ge 3 Sb 2 Te 6 , which is of great importance for the comprehensive understanding of amorphization and crystallization processes.

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

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗
DOE OSTI2023

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multiscale molecular simulations for the solvation of lignin in ionic liquids

Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (H E ) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and H E ) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet–Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA] + -based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA] + -based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.

09 BIOMASS FUELS↗