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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1993

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗
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↗
DOE OSTI2024

Process Intensification by One-Step, Plasma-Assisted Catalytic Synthesis of Liquid Chemicals from Light Hydrocarbons

In this project, we present a plasma-assisted, catalytic process that can handle variations in production rates and gas compositions of producing wells and reliably converts the hydrocarbons to gas phase olefins, high molecular weight alkanes, and liquid chemicals. The low-temperature plasma serves as a reactive chemical environment that activates and converts the light hydrocarbons to these valuable products. We present results from the integration of a catalyst into the low-temperature plasma zone to improve reaction efficiency and product selectivity. We show that the gas composition, bulk gas temperature, and input power of the plasma, with and without catalysts, influence the production rates of liquids from natural gas feeds.

03 NATURAL GAS↗
DOE OSTI2023

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Wetting and spreading at the molecular scale

We have studied the microscopic aspects of the spreading of liquid drops on a solid surface by molecular dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor and solid. We consider both spherically symmetric atoms and chain-like molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observed a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with qualitative behavior resembling recent experimental findings, but with interesting differences in the spreading rate.

Koplik, Joel↗
DOE OSTI2025

Liquid-liquid phase transition of hydrogen and its critical point: Analysis from ab initio simulation and a machine-learned potential

We simulate high-pressure hydrogen in its liquid phase close to molecular dissociation using a machine-learned interatomic potential. The model is trained with density functional theory (DFT) forces and energies, with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. We show that an accurate NequIP model, an E(3)-equivariant neural network potential, accurately reproduces the phase transition present in PBE. Moreover, the computational efficiency of this model allows for substantially longer molecular dynamics trajectories, enabling us to perform a finite-size scaling (FSS) analysis to distinguish between a crossover and a true first-order phase transition. Here, we locate the critical point of this transition, the liquid-liquid phase transition (LLPT), at 1200-1300 K and 155-160 GPa, a temperature lower than most previous estimates and close to the melting transition.

08 HYDROGEN↗
DOE OSTI2023

Liquid Redoxmers for Nonaqueous Redox Flow Batteries

Redoxmers are organic active molecules storing electrochemical energy in nonaqueous redox flow batteries (NRFBs). Increasing the solubility of redoxmers is an important approach for increasing energy density of NRFBs as effective redoxmer concentration determines how much electricity can be stored in a given volume. Molecular engineering redoxmers towards liquid forms is regarded as one promising strategy as liquid redoxmers represent an extreme scenario where fluidity is maintained at maximum concentration using a minimum amount of supporting solvents. In our Perspective, recent examples of liquid redoxmers as well as their development strategy will be discussed.

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

Mechanistic Pathways for the Molecular Step Growth of Calcium Oxalate Monohydrate Crystal Revealed by In Situ Liquid-Phase Atomic Force Microscopy

Calcium oxalate monohydrate (COM) crystal is the most common crystalline component of human kidney stones. The molecular-scale inhibitory mechanisms of COM crystal growth by urinary biomolecules such as citrate and osteopontin adsorbed onto the crystal surface are now well understood. However, the pathways by which dissolved calcium and oxalate ions are incorporated into the molecular step of the COM crystal surface, leading to COM crystal growth—a prerequisite to be elucidated for developing effective therapeutics to inhibit COM stones—remain unknown. Here, using in situ liquid-phase atomic microscopy along with a step kinetic model, we reveal the pathways of the calcium and oxalate ions into the COM molecular step via the growth speed analysis of the molecular steps with respect to their step width at the nanoscale. Our results show that, primarily, the ions are adsorbed onto the terrace of the crystal surface from the solution—the rate-controlling stage for the molecular step growth, i.e., COM crystal growth—and then diffuse over it and are eventually incorporated into the steps. This primary pathway of the ions is unaffected by the model peptide D-Asp6 adsorbed on the COM crystal surface, suggesting that urinary biomolecules will not alter the pathway. These new findings rendering an essential understanding of the fundamental growth mechanism of COM crystal at the nanoscale provide crucial insights beneficial to the development of effective therapeutics for COM kidney stones.

59 BASIC BIOLOGICAL SCIENCES↗