Search NASA⌕ Search

SEARCH · Search NASA

Results for “Liquid structure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Molecular dynamics simulation of the solid-liquid interface migration in terbium

In this work, we developed a Tb embedded atom method potential which properly reproduces the liquid structure obtained from the ab initio molecular dynamics simulation, the hexagonal close packed (hcp)-body-centered cubic (bcc) phase transformation, and melting temperatures. At least three crystal phases [hcp, face-centered cubic (fcc), and bcc] described by this potential can coexist with the liquid phase. Thus, the developed potential provides an excellent test bed for studies of the completive phase nucleation and growth in a single component system. The molecular dynamics simulation showed that all crystal phases can grow from the liquid phase close to their melting temperatures. However, in the cases of the hcp and fcc growth from the liquid phase at very large supercoolings, the bcc phase forms at the solid-liquid interface in the close packed orientations in spite of the fact that both hcp and fcc phases are more stable than the bcc phase at these temperatures. Finally, this bcc phase closes the hcp and fcc phase from the liquid such that the remaining liquid solidifies into the bcc phase. The initial hcp phase then slowly continues growing in expense of the bcc phase.

36 MATERIALS SCIENCE↗

Reconfigurable structure and tunable transport in synchronized active spinner materials

Ensembles of actuated colloids are excellent model systems to explore emergent out-of-equilibrium structures, complex collective dynamics, and design rules for the next generation materials. Here, we demonstrate that ferromagnetic microparticles suspended at an air-water interface and energized by an external rotating magnetic field spontaneously form dynamic ensembles of synchronized spinners in a certain range of the excitation field parameters. Each spinner generates strong hydrodynamic flows, and collective interactions of the multiple spinners promote a formation of dynamic lattices. On the basis of experiments and simulations, we reveal structural transitions from liquid to nearly crystalline states in this novel active spinner material and demonstrate that dynamic spinner lattices are reconfigurable, capable of self-healing behavior and that the transport of embedded inert cargo particles can be remotely tuned by the parameters of the external excitation field. Our findings provide insights into the behavior of active spinner materials with reconfigurable structural order and tunable functionalities.

36 MATERIALS SCIENCE↗

X-ray Absorption Spectral Signature of Quantum Nuclei in Liquid Water

Based on electron-hole excitation theory, we investigate the X-ray absorption spectral signature of nuclear quantum effect in liquid water, whose molecular structure is simulated by path-integral molecular dynamics using the MB-pol model. Compared to spectra generated from classically modeled water structure, quantum nuclei has important effect on spectra in terms of both the spectral energies and their line shapes. At the short-range ordering of H-bond network, the delocalized protons increase the fluctuations on the intramolecular covalency and broaden the pre-edge of the spectra. For intermediate-range and long-range orderings, the observed red and blue shifts of the main-edge and post-edge are attributed to the so-called competing quantum effects, under which both the weak and well-formed H-bonds are promoted. The theoretical spectra are in nearly quantitative agreement with the available experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A molecular dynamics study of the effect of thermal boundary conductance on thermal transport of ideal crystal of n-alkanes with different number of carbon atoms

Phase change materials such as n-alkanes that exhibit desirable characteristics such as high latent heat, chemical stability, and negligible supercooling are widely used in thermal energy storage applications. However, n-alkanes have the drawback of low thermal conductivity values. The low thermal conductivity of n-alkanes is linked to formation of randomly oriented nano-domains of molecules in their solid structure that is responsible for excessive phonon scattering at the grain boundaries. Thus, understanding the thermal boundary conductance at the grain boundaries can be crucial for improving the effectiveness of thermal storage systems. The concept of the ideal crystal is proposed in this paper, which describes a simplified model such that all the nano-domains of long-chain n-alkanes are artificially aligned perfectly in one direction. In order to study thermal transport of the ideal crystal of long-chain n-alkanes, four (4) systems (C{sub 20}H{sub 42}, C{sub 24}H{sub 50}, C{sub 26}H{sub 54}, and C{sub 30}H{sub 62}) are investigated by the molecular dynamics simulations. Thermal boundary conductance between the layers of ideal crystals is determined using both non-equilibrium molecular dynamics (NEMD) and equilibrium molecular dynamics (EMD) simulations. Both NEMD and EMD simulations exhibit no significant change in thermal conductance with the molecular length. However, the values obtained from the EMD simulations are less than the values from NEMD simulations with the ratio being nearly three (3) in most cases. This difference is due to the nature of EMD simulations where all the phonons are assumed to be in equilibrium at the interface. Thermal conductivity of the n-alkanes in three structures including liquid, solid, and ideal crystal is investigated utilizing NEMD simulations. Our results exhibit a very slight rise in thermal conductivity values as the number of carbon atoms of the chain increases. The key understanding is that thermal transport can be significantly altered by how the molecules and the nano-domains are oriented in the structure rather than by the length of the n-alkane molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A method of orbital analysis for large-scale first-principles simulations

An efficient method of calculating the natural bond orbitals (NBOs) based on a truncation of the entire density matrix of a whole system is presented for large-scale density functional theory calculations. The method recovers an orbital picture for O(N) electronic structure methods which directly evaluate the density matrix without using Kohn-Sham orbitals, thus enabling quantitative analysis of chemical reactions in large-scale systems in the language of localized Lewis-type chemical bonds. With the density matrix calculated by either an exact diagonalization or O(N) method, the computational cost is O(1) for the calculation of NBOs associated with a local region where a chemical reaction takes place. As an illustration of the method, we demonstrate how an electronic structure in a local region of interest can be analyzed by NBOs in a large-scale first-principles molecular dynamics simulation for a liquid electrolyte bulk model (propylene carbonate + LiBF{sub 4})

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of interatomic potentials appropriate for simulation of devitrification of Al 90 Sm 10 alloy

In this study, a semi-empirical potential for the Al 90 Sm 10 alloy is presented. The potential provides satisfactory reproduction of pure Al properties, the formation energies of a set of Al–Sm crystal phases with Sm content about 10%, and the structure of the liquid Al 90 Sm 10 alloy. During molecular dynamics simulation in which the liquid alloy is cooled at a rate of 10 10 K/s, the developed potential produces a glass structure with lower ab initio energy than that produced by ab initio molecular dynamics (AIMD) itself using a typical AIMD cooling rate of 8 ∙10 13 K/s. Based on these facts the developed potential should be suitable for simulations of phase transformations in the Al 90 Sm 10 alloy.

36 MATERIALS SCIENCE↗

Weak links between fast mobility and local structure in molecular and atomic liquids

We investigate by molecular-dynamics simulations, the fast mobility—the rattling amplitude of the particles temporarily trapped by the cage of the neighbors—in mildly supercooled states of dense molecular (linear trimers) and atomic (binary mixtures) liquids. The mixture particles interact by the Lennard-Jones potential. The non-bonded particles of the molecular system are coupled by the more general Mie potential with variable repulsive and attractive exponents in a range which is a characteristic of small n-alkanes and n-alcohols. Possible links between the fast mobility and the geometry of the cage (size and shape) are searched. The correlations on a per-particle basis are rather weak. Instead, if one groups either the particles in fast-mobility subsets or the cages in geometric subsets, the increase of the fast mobility with both the size and the asphericity of the cage is revealed. The observed correlations are weak and differ in states with equal relaxation time. Local forces between a tagged particle and the first-neighbour shell do not correlate with the fast mobility in the molecular liquid. It is concluded that the cage geometry alone is unable to provide a microscopic interpretation of the known, universal link between the fast mobility and the slow structural relaxation. We suggest that the particle fast dynamics is affected by regions beyond the first neighbours, thus supporting the presence of collective, extended fast modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polyamorphic Transformations in Fe-Ni-C Liquids: Implications for Chemical Evolution of Terrestrial Planets

During the formation of the Earth's core, the segregation of metallic liquids from silicate mantle should have left behind evident geochemical imprints on both the mantle and the core. Some distinctive geochemical signatures of the mantle-derived rocks likely own their origin to the metal-silicate differentiation of the primitive Earth, setting our planet apart from undifferentiated meteorites as well as terrestrial planets or moons isotopically and compositionally. Understanding the chemical evolution of terrestrial planetary bodies requires knowledge on properties of both liquid iron alloys and silicates equilibrating under physicochemical conditions pertinent to the deep magma ocean. Here we report experimental and computational results on the pressure-induced structural evolution of iron-nickel liquids alloyed with carbon. Our X-ray diffraction experiments up to 7.3 gigapascals (GPa) demonstrate that Fe-Ni (Fe 90 Ni 10 ) liquids alloyed with 3 and 5 wt % carbon undergo a polyamorphic liquid structure transition at approximately 5 GPa. Corroborating the experimental observations, our first-principles molecular dynamic calculations reveal that the structural transitions result from the marked prevalence of three-atom face-sharing polyhedral connections in the liquids at >5 GPa. The structure and polyamorphic transitions of liquid iron-nickel-carbon alloys govern their physical and chemical properties and may thus cast fresh light on the chemical evolution of terrestrial planets and moons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing the structural and dynamical properties of liquid water with models including non-local electron correlation

Water is a ubiquitous liquid that displays a wide range of anomalous properties and has a delicate structure that challenges experiment and simulation alike. The various intermolecular interactions that play an important role, such as repulsion, polarization, hydrogen bonding, and van der Waals interactions, are often difficult to reproduce faithfully in atomistic models. Here, electronic structure theories including all these interactions at equal footing, which requires the inclusion of non-local electron correlation, are used to describe structure and dynamics of bulk liquid water. Isobaric-isothermal (NpT) ensemble simulations based on the Random Phase Approximation (RPA) yield excellent density (0.994 g/ml) and fair radial distribution functions, while various other density functional approximations produce scattered results (0.8-1.2 g/ml). Molecular dynamics simulation in the microcanonical (NVE) ensemble based on Møller-Plesset perturbation theory (MP2) yields dynamical properties in the condensed phase, namely, the infrared spectrum and diffusion constant. At the MP2 and RPA levels of theory, ice is correctly predicted to float on water, resolving one of the anomalies as resulting from a delicate balance between van der Waals and hydrogen bonding interactions. For several properties, obtaining quantitative agreement with experiment requires correction for nuclear quantum effects (NQEs), highlighting their importance, for structure, dynamics, and electronic properties. A computed NQE shift of 0.6 eV for the band gap and absorption spectrum illustrates the latter. Giving access to both structure and dynamics of condensed phase systems, non-local electron correlation will increasingly be used to study systems where weak interactions are of paramount importance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Hydrogen Bonding on the Dynamics and Structure of Protic Ionic Liquid/Water Binary Mixtures

In this work, the orientational dynamics and microscopic liquid structure of a protic ionic liquid, 1- ethylimidazolium bis(trifluoromethylsulfonyl)imide (EhimNTf 2 ), and its aprotic analogue, 1- ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EmimNTf 2 ), were studied at various water concentrations using optical heterodyne-detected optical Kerr Effect (OHD-OKE) spectroscopy, linear infrared spectroscopy, and atomistic simulations. The OHD-OKE experiments essentially measure the orientational relaxation of the Ehim + and Emim + cations. The experiments and simulations show a significant dynamical and structural change in EhimNTf 2 between the 2:1 ion pair:water and the 1:1 ion pair:water concentrations. The OHDOKE data show that EmimNTf 2 /water mixtures exhibit hydrodynamic behavior at all water concentrations up to saturation. In contrast, EhimNTf 2 /water mixtures deviate from hydrodynamic dynamic behavior at water concentrations above 2:1. At the 1:1 concentration, the orientational randomization of the Ehim + cation is slower than that predicted using viscosity data. Atomistic simulation results reveal the microscopic ionic structures of dry liquids and the preferential hydrogen bonding of water to the H atom of the N-H of Ehim + over other sites on the Ehim + and Emim + cations. Atomistic simulation results demonstrate that in EhimNTf 2 RTIL/water mixtures there is a substantial jump in the formation of water-water hydrogen bonds in addition to N-H-water hydrogen bonds upon increasing the water concentration from 2:1 to 1:1. Water-water hydrogen bonding strengthens the spatial coordination of the H atom of the N-H moiety of Ehim + to neighboring water molecules through preferential hydrogen bonding. The jump in the concentration of water-water hydrogen bonds occurs at the Ehim + /water concentration at which the orientational relaxation deviates from hydrodynamic behavior. This structural observation is confirmed with FT-IR spectra that show asymmetry in the peak for the O-D stretch that is indicative of water clusters. The formation of water clusters and the strengthening of the N-H … OH 2 hydrogen bonds slow the orientational relaxation of Ehim + cations as observed by the OHD-OKE experiments.

14 SOLAR ENERGY↗

Interfaces of dicationic ionic liquids and graphene: a molecular dynamics simulation study

Molecular dynamics simulations were performed to investigate the interfacial structure and capacitance of electrical double layers (EDLs) in dicationic ionic liquids (DILs) 1-alkyl-3-dimethylimidazolium tetrafluoroborate [Cn(mim)2](BF4)2 (n = 3, 6, 9), with respect to a baseline of a monocationic ionic liquid [C3mim][BF4], near planar carbon electrodes consisting of graphene sheets. The simulation results show that an adsorbed layer with double peaks is exclusively found for [C3(mim)2](BF4)2, while a single peak of the other three cations is observed at the neutral electrode, due to the difference in ion–wall interaction and cation–anion association. As the electrode becomes negatively charged, the second peak of [C3(mim)2]2+ is dramatically reduced, whereas those of [C6(mim)2]2+ and [C9(mim)2]2+ become non-trivial. The capacitance–potential curve of EDLs in DILs manifests a transition from camel shape to bell shape as the cation chain length increases, which is attributed to the enlargement of ion adsorption (per unit charge) on the electrode and the decrease of attractive interaction between ions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Solvents Polarity Governs Ion Interactions and Transport in a Solvated Room Temperature Ionic Liquid

We explore the influence of the solvent dipole moment on cation anion interactions and transport in 1-butyl-3-methyl-imidazolium bis-(trifluoromethylsulfonyl), [BMIM+][Tf2N ]. Free energy profiles derived from atomistic molecular dynamics (MD) simulations show a correlation of the cation anion separation and the equilibrium depth of the potential of mean force with the dipole moment of the solvent. Correlations of the ion diffusivity with the dipole moment and the concentration of the solvent were further demonstrated by classical MD simulations. Quasi-elastic neutron scattering experiments with deuterated solvents reveal a complex picture of nanophase separation into the ionic liquid-rich and solvent-rich phases. The experiment corroborates the trend of concentration- and dipole moment-dependent enhancement of ion mobility by the solvent, as suggested by the simulations. Despite the considerable structural complexity of ionic liquid solvent mixtures, we can rationalize and generalize the trends governing ionic transport in these complex electrolytes.

Osti, Naresh C↗

High-Performance First-Principles Molecular Dynamics for Predictive Theory and Modeling

This project focused on developing high-performance software tools for First-Principles Molecular Dynamics (FPMD) simulations, and applying them in investigations of materials relevant to energy conversion processes. FPMD is an atomistic simulation method that combines a quantum-mechanical description of electronic structure with the statistical description provided by molecular dynamics (MD) simulations. This reliance on fundamental principles allows FPMD simulations to provide a consistent description of structural, dynamical and electronic properties of a material. This is particularly useful in systems for which reliable empirical models are lacking. FPMD simulations are increasingly used as a predictive tool for applications such as batteries, solar energy conversion, light-emitting devices, electro-chemical energy conversion devices and other materials. During the course of the project, several new features were developed and added to the open-source Qbox FPMD code. The code was further optimized for scalable operation of large-scale, Leadership-Class DOE computers. When combined with Many-Body Perturbation Theory (MBPT) calculations, this infrastructure was used to investigate structural and electronic properties of liquid water, ice, aqueous solutions, nanoparticles and solid-liquid interfaces. Computing both ionic trajectories and electronic structure in a consistent manner enabled the simulation of several spectroscopic properties, such as Raman spectra, infrared spectra, and sum-frequency generation spectra. The accuracy of the approximations used allowed for direct comparisons of results with experimental data such as optical spectra, X-ray and neutron diffraction spectra. The software infrastructure developed in this project, as applied to various investigations of solids, liquids and interfaces, demonstrates that FPMD simulations can provide a detailed, atomic-scale picture of structural, vibrational and electronic properties of complex systems relevant to energy conversion devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics simulations of the structure and single-particle dynamics of mixtures of divalent salts and ionic liquids

We report a molecular dynamics study of the structure and single-particle dynamics of mixtures of a protic (ethylammonium nitrate) and an aprotic (1-butyl-3-methylimidazolium hexaflurophosphate [BMIM][PF{sub 6}]) room-temperature ionic liquids doped with magnesium and calcium salts with a common anion at 298.15 K and 1 atm. The solvation of these divalent cations in dense ionic environments is analyzed by means of apparent molar volumes of the mixtures, radial distribution functions, and coordination numbers. For the protic mixtures, the effect of salt concentration on the network of hydrogen bonds is also considered. Moreover, single-particle dynamics of the salt cations is studied by means of their velocity autocorrelation functions and vibrational densities of states, explicitly analyzing the influence of salt concentration, and cation charge and mass on these magnitudes. The effect of the valency of the salt cation on these properties is considered comparing the results with those for the corresponding mixtures with lithium salts. We found that the main structural and dynamic features of the local solvation of divalent cations in ionic liquids are similar to those of monovalent salts, with cations being localized in the polar nanoregions of the bulk mixture coordinated in monodentate and bidentate coordination modes by the [NO{sub 3}]{sup −} and [PF{sub 6}]{sup −} anions. However, stronger electrostatic correlations of these polar nanoregions than in mixtures with salts with monovalent cations are found. The vibrational modes of the ionic liquid (IL) are seen to be scarcely affected by the addition of the salt, and the effect of mass and charge on the vibrational densities of states of the dissolved cations is reported. Cation mass is seen to exert a deeper influence than charge on the low-frequency vibrational spectra, giving a red shift of the vibrational modes and a virtual suppression of the higher energy vibrational modes for the heavier Ca{sup 2+} cations. No qualitative difference with monovalent cations was found in what solvation is concerned, which suggests that no enhanced reduction of the mobility of these cations and their complexes in ILs respective to those of monovalent cations is to be expected.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗