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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Stokes–Einstein relation and excess entropy in Al-rich Al-Cu melts

We investigate the conditions for the validity of the Stokes-Einstein relation that connects diffusivity to viscosity in melts using entropy-scaling relationships developed by Rosenfeld. Employing ab initio molecular dynamics simulations to determine transport and structural properties of liquid Al{sub 1−x}Cu{sub x} alloys (with composition x ≤ 0.4), we first show that reduced self-diffusion coefficients and viscosities, according to Rosenfeld's formulation, scale with the two-body approximation of the excess entropy except the reduced viscosity for x = 0.4. Then, we use our findings to evidence that the Stokes-Einstein relation using effective atomic radii is not valid in these alloys while its validity can be related to the temperature dependence of the partial pair-excess entropies of both components. Finally, we derive a relation between the ratio of the self-diffusivities of the components and the ratio of their pair excess entropies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Molecular dynamics simulation studies of structural and dynamical properties of rapidly quenched Al

The structural and dynamical properties of rapidly quenched Al are studied by molecular dynamics simulations. The pair-correlation function of high temperature liquid Al agrees well with the experimental results. Different cooling rates are applied with high cooling rates leading to glass formation, while low cooling rates leading to crystallization. The local structures are characterized by Honeycutt–Andersen indices and Voronoi tessellation analysis. The results show that for high cooling rates, the local structures of the liquid and glassy Al are predominated by icosahedral clusters, together with considerable amount of face-centered cubic and hexagonal close packed short-range orders. These short-range order results are further confirmed using the recently developed atomic cluster alignment method. Moreover, the atomic cluster alignment clearly shows the crystal nucleation process in supercooled liquid of Al. Finally, the mean square displacement for the liquid is also analyzed, and the corresponding diffusion coefficient as a function of temperature is calculated.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Local structure origin of ultrafast crystallization driven by high-fidelity octahedral clusters in amorphous Sc 0.2 Sb 2 Te 3

Phase-change material Sc 0.2 Sb 2 Te 3 (SST) can remarkably boost the writing speed of memory devices due to the extremely fast crystallization. It was reflected that the fast crystallization is because Sc stabilizes the 4-fold rings which act as precursors of nuclei in the amorphous phase. Here, by using first-principles molecular dynamics simulations, we studied the local structures in the liquid and amorphous SST at various temperatures. The findings reveal that Sc-centered configurations are almost in the form of stable octahedral clusters, which enhances the local order of the amorphous phase. Different from Sb- and Te-centered clusters which have lower coordination numbers, Sc-centered clusters are mainly in high-coordinated octahedral structures. These Sc-centered octahedral clusters present a high stability in supercooled liquid and amorphous states, remarkably reducing the incubation time of nucleation and speeding up the crystallization. Our study reveals the role of Sc atoms in the liquid and amorphous structure, paving the way for the application of Sc-based phase-change memory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Structural Anomalies and Electronic Properties of an Ionic Liquid under Nanoscale Confinement

Ionic liquids (ILs) promise far greater electrochemical performance compared to aqueous systems, yet key physicochemical properties governing their assembly at interfaces within commonly used graphitic nanopores remain poorly understood. Here, we combine synchrotron X-ray scattering with first-principles molecular dynamics simulations to unravel key structural characteristics of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([TFSI] – ) ionic liquids confined in carbon slit pores. X-ray scattering reveals selective pore filling due to size exclusion, while filled pores exhibit disruption in the IL intermolecular structure, the extent of which increases for narrower slit pores. First-principles simulations corroborate this finding and quantitatively describe how perturbations in the local IL structure, particularly the hydrogen-bond network, depend strongly on the degree of confinement. Despite significant deviations in structure under confinement, electrochemical stability remains intact, which is important for energy storage based on nanoporous carbon electrodes (e.g., supercapacitors).

36 MATERIALS SCIENCE↗
DOE OSTI2016

Structure and dynamics of aqueous solutions from PBE-based first-principles molecular dynamics simulations

Establishing an accurate and predictive computational framework for the description of complex aqueous solutions is an ongoing challenge for density functional theory based first-principles molecular dynamics (FPMD) simulations. In this context, important advances have been made in recent years, including the development of sophisticated exchange-correlation functionals. On the other hand, simulations based on simple generalized gradient approximation (GGA) functionals remain an active field, particularly in the study of complex aqueous solutions due to a good balance between the accuracy, computational expense, and the applicability to a wide range of systems. In such simulations we often perform them at elevated temperatures to artificially “correct” for GGA inaccuracies in the description of liquid water; however, a detailed understanding of how the choice of temperature affects the structure and dynamics of other components, such as solvated ions, is largely unknown. In order to address this question, we carried out a series of FPMD simulations at temperatures ranging from 300 to 460 K for liquid water and three representative aqueous solutions containing solvated Na + , K + , and Cl - ions. We show that simulations at 390–400 K with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional yield water structure and dynamics in good agreement with experiments at ambient conditions. Simultaneously, this computational setup provides ion solvation structures and ion effects on water dynamics consistent with experiments. These results suggest that an elevated temperature around 390–400 K with the PBE functional can be used for the description of structural and dynamical properties of liquid water and complex solutions with solvated ions at ambient conditions.

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

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

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

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

Conformation and interactions of dopamine hydrochloride in solution

The aqueous solution of dopamine hydrochloride has been investigated using neutron and X-ray total scattering data together with Monte-Carlo based modelling using Empirical Potential Structure Refinement. The conformation of the protonated dopamine molecule is presented and the results compared to the conformations found in crystal structures, dopamine-complexed protein crystal structures and predicted from theoretical calculations and pharmacophoric models. It is found that protonated dopamine adopts a range of conformations in solution, highlighting the low rotational energy barrier between different conformations, with the preferred conformation being trans-perpendicular. The interactions between each of the species present (protonated dopamine molecules, water molecules, and chloride anions) have been determined and are discussed with reference to interactions observed in similar systems both in the liquid and crystalline state, and predicted from theoretical calculations. The expected strong hydrogen bonds between the strong hydrogen bond donors and acceptors are observed, together with evidence of weaker CH hydrogen bonds and π interactions also playing a significant role in determining the arrangement of adjacent molecules.

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

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗
DOE OSTI2014

Structure and dynamics in liquid bismuth and Bi{sub n} clusters: A density functional study

Density functional/molecular dynamics simulations with more than 500 atoms have been performed on liquid bismuth at 573, 773, 923, and 1023 K and on neutral Bi clusters with up to 14 atoms. There are similar structural patterns (coordination numbers, bond angles, and ring patterns) in the liquid and the clusters, with significant differences from the rhombohedral crystalline form. We study the details of the structure (structure factor, pair, and cavity distribution functions) and dynamical properties (vibration frequencies, diffusion constants, power spectra), and compare with experimental results where available. While the three short covalent bonds typical to pnictogens are characteristic in both liquid and clusters, the number of large voids and the total cavity volume is much larger in the liquid at 1023 K, with larger local concentration variations. The inclusion of spin-orbit coupling results in a lowering of the cohesive energies in Bi{sub n} clusters of 0.3–0.5 eV/atom.

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

Structure and dynamics of POPC bilayers in water solutions of room temperature ionic liquids

Molecular dynamics simulations in the NPT ensemble have been carried out to investigate the effect of two room temperature ionic liquids (RTILs), on stacks of phospholipid bilayers in water. We consider RTIL compounds consisting of chloride ([bmim][Cl]) and hexafluorophosphate ([bmim][PF{sub 6}]) salts of the 1-buthyl-3-methylimidazolium ([bmim]{sup +}) cation, while the phospholipid bilayer is made of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Our investigations focus on structural and dynamical properties of phospholipid and water molecules that could be probed by inelastic and quasi-elastic neutron scattering measurements. The results confirm the fast incorporation of [bmim]{sup +} into the lipid phase already observed in previous simulations, driven by the Coulomb attraction of the cation for the most electronegative oxygens in the POPC head group and by sizeable dispersion forces binding the neutral hydrocarbon tails of [bmim]{sup +} and of POPC. The [bmim]{sup +} absorption into the bilayer favours the penetration of water into POPC, causes a slight but systematic thinning of the bilayer, and further stabilises hydrogen bonds at the lipid/water interface that already in pure samples (no RTIL) display a lifetime much longer than in bulk water. On the other hand, the effect of RTILs on the diffusion constant of POPC (D{sub POPC}) does not reveal a clearly identifiable trend, since D{sub POPC} increases upon addition of [bmim][Cl] and decreases in the [bmim][PF{sub 6}] case. Moreover, because of screening, the electrostatic signature of each bilayer is only moderately affected by the addition of RTIL ions in solution. The analysis of long wavelength fluctuations of the bilayers shows that RTIL sorption causes a general decrease of the lipid/water interfacial tension and bending rigidity, pointing to the destabilizing effect of RTILs on lipid bilayers.

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

Structural connection between gallium crystals and near- T m liquids under ambient pressure

We study the short-range structural order in liquid gallium prepared by ab initio molecular dynamics simulations. We found that at 400 K, which is close to the melting point (Tm) of Ga, the dominant motif in the liquid phase is identical to that in the stable solid phase Ga-I. Strong directional bonding in this motif prevents it from nucleation in the liquid phase. Meanwhile, a newly identified motif, topologically distinct from yet related to those in β-Ga and Ga-III phases, is also abundant in Ga liquid. This new motif could serve as precursors for β-Ga or Ga-III during crystallization.

Zhang, F.↗
DOE OSTI2022

Structural and Electronic Transitions in Liquid FeO Under High Pressure

Abstract FeO represents an important end‐member for planetary interiors mineralogy. However, its properties in the liquid state under high pressure are poorly constrained. Here, in situ high‐pressure and high‐temperature X‐ray diffraction experiments, ab initio simulations, and thermodynamic calculations are combined to study the local structure and density evolution of liquid FeO under extreme conditions. Our results highlight a strong shortening of the Fe‐Fe distance, particularly pronounced between ambient pressure and ∼40 GPa, possibly related with the insulator to metal transition occurring in solid FeO over a similar pressure range. Liquid density is smoothly evolving between 60 and 150 GPa from values calculated for magnetic liquid to those calculated for non‐magnetic liquid, compatibly with a continuous spin crossover in liquid FeO. The present findings support the potential decorrelation between insulator/metal transition and the high‐spin to low‐spin continuous transition, and relate the changes in the microscopic structure with macroscopic properties, such as the closure of the Fe‐FeO miscibility gap. Finally, these results are used to construct a parameterized thermal equation of state for liquid FeO providing densities up to pressure and temperature conditions expected at the Earth's core‐mantle boundary.

Morard, Guillaume↗
DOE OSTI2015

Molecular structure of the discotic liquid crystalline phase of hexa-peri-hexabenzocoronene/oligothiophene hybrid and their charge transport properties

Using atomistic molecular dynamics simulation, we study the discotic columnar liquid crystalline (LC) phases formed by a new organic compound having hexa-peri-Hexabenzocoronene (HBC) core with six pendant oligothiophene units recently synthesized by Nan Hu et al. [Adv. Mater. 26, 2066 (2014)]. This HBC core based LC phase was shown to have electric field responsive behavior and has important applications in organic electronics. Our simulation results confirm the hexagonal arrangement of columnar LC phase with a lattice spacing consistent with that obtained from small angle X-ray diffraction data. We have also calculated various positional and orientational correlation functions to characterize the ordering of the molecules in the columnar arrangement. The molecules in a column are arranged with an average twist of 25° having an average inter-molecular separation of ∼5 Å. Interestingly, we find an overall tilt angle of 43° between the columnar axis and HBC core. We also simulate the charge transport through this columnar phase and report the numerical value of charge carrier mobility for this liquid crystal phase. The charge carrier mobility is strongly influenced by the twist angle and average spacing of the molecules in the column.

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

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

The effect of hydroxyl spacing in diols on the solvation structure, dynamics, and transport properties of choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are a class of liquids that offer great potential in alleviating some of the challenges present in today's long-term energy storage methods because they have physical properties that are favorable for storable electrolyte solutions. In this work, a series of glycols (ethylene glycol, 1,3-propanediol, 1,4-butanediol, and 1,5-pentanediol) were studied as potential hydrogen bond donors (HBD) with a common choline chloride (ChCl) as the hydrogen bond acceptor (HBA). The solvation dynamics of the prepared systems were studied by measuring the solvent reorganization response using femtosecond transient absorption spectroscopy (fs-TA). Conductivity, viscosity, density, ET(30) polarity, and dynamics of the prepared DESs were analyzed, with a particular interest in determining the effect of HBD chain length on these parameters. Here, classical molecular dynamics simulations were employed to investigate how the local liquid structure, solvent dynamics, and bulk solvent properties vary with changes in glycol chain length.

Conductivity↗
DOE OSTI2020

Modelling mixed-phase clouds with the large-eddy model UCLALES–SALSA

The large-eddy model UCLALES–SALSA, with an exceptionally detailed aerosol description for both aerosol number and chemical composition, has been extended for ice and mixed-phase clouds. Comparison to a previous mixed-phase cloud model intercomparison study confirmed the accuracy of newly implemented ice microphysics. A further simulation with a heterogeneous ice nucleation scheme, in which ice-nucleating particles (INPs) are also a prognostic variable, captured the typical layered structure of Arctic mid-altitude mixed-phase cloud: a liquid layer near cloud top and ice within and below the liquid layer. In addition, the simulation showed a realistic freezing rate of droplets within the vertical cloud structure. The represented detailed sectional ice microphysics with prognostic aerosols is crucially important in reproducing mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗