Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “liquid dynamics”

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 163 records · Page 9

Two-component dynamics and the liquid-like to gas-like crossover in supercritical water

Molecular-scale dynamics in sub- to super-critical water is studied with inelastic X-ray scattering and molecular dynamics simulations. The obtained longitudinal current correlation spectra can be decomposed into two main components: a low-frequency (LF), gas-like component and a high-frequency (HF) component arising from the O--O stretching mode between hydrogen-bonded molecules, reminiscent of the longitudinal acoustic mode in ambient water. With increasing temperature, the hydrogen-bond network diminishes and the spectral weight shifts from HF to LF, leading to a transition from liquid-like to gas-like dynamics with rapid changes around the Widom line.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Isomerization reaction dynamics and equilibrium at the liquid-vapor interface of water. A molecular-dynamics study

The gauche-trans isomerization reaction of 1,2-dichloroethane at the liquid-vapor interface of water is studied using molecular-dynamics computer simulations. The solvent bulk and surface effects on the torsional potential of mean force and on barrier recrossing dynamics are computed. The isomerization reaction involves a large change in the electric dipole moment, and as a result the trans/gauche ratio is considerably affected by the transition from the bulk solvent to the surface. Reactive flux correlation function calculations of the reaction rate reveal that deviation from the transition-state theory due to barrier recrossing is greater at the surface than in the bulk water. This suggests that the system exhibits non-Rice-Ramsperger-Kassel-Marcus behavior due to the weak solvent-solute coupling at the water liquid-vapor interface.

Benjamin, Ilan↗

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature T o , the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of T o is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding–unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water–air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin–Wagner fluctuations in experiments and molecular simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Vibrational dynamics in the undercooled liquid of ultra-fragile metallic glasses

Here, the lattice vibrational dynamics in ultra-fragile Pt 57 Cu 23 P 20 and Pt 60 Cu 20 P 20 glasses were investigated using in situ inelastic neutron scattering to measure the evolution of the phonon density of states (DOS) through the glass transition, into the undercooled liquid phase, and through crystallization. The phonon DOS curves are remarkably similar for the undercooled liquids and the amorphous solids. The vibrational entropy was obtained with these phonon DOS curves, and changes in vibrational entropy caused by the glass transition are less than 0.05 k B /atom. For these highly fragile metal/metalloid glasses, the excess vibrational entropy of the undercooled liquid is about 0.5–1.5 Jmol -1 K -1 , compared to 0.1 Jmol -1 K -1 reported earlier for stronger Cu–Zr glasses. Our work suggests that, consistent with our previous report for less fragile metallic glasses, the fraction of the total excess entropy that is vibrational is small, and the discontinuity in excess entropy at the glass transition is mainly due to configurational entropy.

36 MATERIALS SCIENCE↗

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗

Atomic dynamics of electrified solid–liquid interfaces in liquid-cell TEM

Electrified solid–liquid interfaces (ESLIs) play a key role in various electrochemical processes relevant to energy, biology and geochemistry. The electron and mass transport at the electrified interfaces may result in structural modifications that markedly influence the reaction pathways. For example, electrocatalyst surface restructuring during reactions can substantially affect the catalysis mechanisms and reaction products. Despite its importance, direct probing the atomic dynamics of solid–liquid interfaces under electric biasing is challenging owing to the nature of being buried in liquid electrolytes and the limited spatial resolution of current techniques for in situ imaging through liquids. Here, with our development of advanced polymer electrochemical liquid cells for transmission electron microscopy (TEM), we are able to directly monitor the atomic dynamics of ESLIs during copper (Cu)-catalysed CO 2 electroreduction reactions (CO 2 ERs). Our observation reveals a fluctuating liquid-like amorphous interphase. It undergoes reversible crystalline–amorphous structural transformations and flows along the electrified Cu surface, thus mediating the crystalline Cu surface restructuring and mass loss through the interphase layer. Furthermore, the combination of real-time observation and theoretical calculations unveils an amorphization-mediated restructuring mechanism resulting from charge-activated surface reactions with the electrolyte. Our results open many opportunities to explore the atomic dynamics and its impact in broad systems involving ESLIs by taking advantage of the in situ imaging capability.

36 MATERIALS SCIENCE↗

Machine learning-driven descriptions of protein dynamics at solid-liquid interfaces

This chapter has described how ML has enabled quantitative analysis of HS-AFM data to discover the physical phenomena governing protein dynamics and ordering at solid-liquid interfaces. The research detailed in this chapter modeled the rotation models of protein nanorods, the discovery of which would otherwise not be possible. By tracking the trajectories of individual protein rods from frame to frame, it was possible to model Brownian type motion and behaviors and Levy-flight dynamics that had not previously been shown. We also described the application of the Python package AtomAI, which has been developed specifically to analyze and extract physical phenomena, providing exemplar code for training an ensemble of deep neural networks to produce the semantic segmentation of AFM data and functions for encoding and decoding local environments. We last described a combinatorial approach to analyze very noisy data with a densely covered substrate where the emergence of order for the protein liquid crystals could be elucidated. By combining the methods from Case 1 and 2, it was possible to obtain the center of mass and angle for each rod in the images and track the assembly of the rods over time into a 2D liquid crystal array on the surface of mica.

protein dynamics, solid-liquid interfaces, atomic ↗

Nonlinear dielectric relaxation of polar liquids

Molecular dynamics of two water models, SPC/E and TIP3P, at a number of temperatures are used to test the Kivelson-Madden equation connecting single-particle and collective dielectric relaxation times through the Kirkwood factor. The relation is confirmed by simulations and used to estimate the nonlinear effect of the electric field on the dielectric relaxation time. We show that the main effect of the field comes through slowing down of the single-particle rotational dynamics and the relative contribution of the field-induced alteration of the Kirkwood factor is insignificant for water. Theories of nonlinear dielectric relaxation need to mostly account for the effect of the field on rotations of a single dipole in a polar liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics of the water liquid-vapor interface

The results of molecular dynamics calculations on the equilibrium interface between liquid water and its vapor at 325 K are presented. For the TIP4P model of water intermolecular pair potentials, the average surface dipole density points from the vapor to the liquid. The most common orientations of water molecules have the C2 nu molecular axis roughly parallel to the interface. The distributions are quite broad and therefore compatible with the intermolecular correlations characteristic of bulk liquid water. All near-neighbor pairs in the outermost interfacial layers are hydrogen bonded according to the common definition adopted here. The orientational preferences of water molecules near a free surface differ from those near rigidly planar walls which can be interpreted in terms of patterns found in hexagonal ice 1. The mean electric field in the interfacial region is parallel to the mean polarization which indicates that attention cannot be limited to dipolar charge distributions in macroscopic descriptions of the electrical properties of this interface. The value of the surface tension obtained is 132 +/- 46 dyn/cm, significantly different from the value for experimental water of 68 dyn/cm at 325 K.

NASA Center ARC↗

Rotating Liquid Drops: Plateaus Experiment Revisited

The dynamics of liquid drops rotating in another liquid were studied experimentally with an oil drop suspended in a neutral buoyancy tank. New stable shapes not predicted by the theory were observed.

Tagg, R.↗

Electrochemical dynamics of imidazolium ionic liquids at graphene electrodes for energy storage applications

Electric double layer capacitors (EDLCs) are prominent energy storage systems that constitute the foundation of more reliable and sustainable energy infrastructures. Modern EDLCs often incorporate ionic liquids (ILs) as a key component in their electrolytes, leveraging the high electrochemical stability of ILs to enhance device performance. The performance and functionality of these capacitors also hinge on the interfacial behavior of ILs at electrode surfaces, which remain insufficiently understood. Here, we performed synchrotron infrared nanospectroscopy (SINS) in combination with density functional theory (DFT) calculations to investigate the electric double layers (EDLs) of three imidazolium-based ILs in a custom-designed graphene liquid cell. This approach revealed new insights into the dynamics of IL EDLs and the underpinning factors originating from the IL structures. Variations in anion size and structure were found to tune the ILs’ ability to form EDLs with compact and closely-correlated ion arrangements, which are critical for enhancing their capacitive performance. These findings highlight the intricate interactions between ions governed by their structures and charge behaviors, which underscores the opportunities for targeted design of IL-based electrolytes to optimize EDLC functionality.

Electric double layers↗

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

Abstract A computational approach by an implementation of the principle component analysis (PCA) with K -means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers ( W C N s ‾ ) of particles calculated from particles’ positions are used as an order parameter to build a low-dimensional representation of feature (structural) space for K -means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states’ particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid–liquid phase separation after quenching to a supercooled state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗