Ionization wave dynamics of a plasma jet in contact with liquid water.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
The nonlinear dynamics of the free surface of an ideal dielectric liquid that is exposed to an external oblique electric field has been studied theoretically. In the framework of the Hamiltonian formalism, a system of nonlinear integro-differential equations has been derived that describes the dynamics of nonlinear waves in the small-angle approximation. It is established that for a liquid with high dielectric permittivity, these equations have a solution in the form of plane waves of arbitrary shape that propagate without distortion in the direction of the horizontal component of the external field.
Phase-change material such as Ge 2 Sb 2 Te 5 is usually utilized to store data due to the pronounced contrast in optical and electrical properties between crystalline and amorphous phases. As the density differs in the two phases, it is necessary to explore the influence of pressure on the structures of Ge 2 Sb 2 Te 5 , especially for the supercooled liquid which is an inevitable state in the formation of the two phases. The short-range structures in supercooled liquid Ge 2 Sb 2 Te 5 under compression have been investigated by using ab initio molecular dynamics simulation. The supercooled liquid eventually changes to a solid with an increase in pressure. During the process, tetrahedrons decrease slightly, revealing that tetrahedral structures are insensitive to the pressure. Octahedrons increase as the pressure is less than 3.7 GPa and then decrease, suggesting that a moderate pressure can promote the formation of octahedrons. The body-centered-cubic, face-centered-cubic and hexagonal closed-packed structures are observed at 9.9 GPa and then increase gradually. Additionally, it is noticed that Sb- and Te-centered clusters prefer Ge-centered clusters to forming the high-coordinated short-range structures during the compression process. Here our research makes us aware of the effects of pressure on short-range structures in supercooled liquid, which is of great importance in the application of phase-change materials.
Most liquid cell transmission electron microscopy (LC TEM) studies focus on nanoparticles or nanowires, in large part because the preparation and study of materials in this size range is straightforward. By contrast, this is not true for samples in the micrometre size range, in large part because of the difficulties associated with sample preparation starting from a ‘bulk’ material. There are also many advantages inherent to the study of micrometre-sized samples compared to their nanometre-sized counterparts. Here, we present a liquid cell transmission electron study that employed an innovative sample preparation technique using focused ion beam (FIB) milling to fabricate micrometre-sized electron transparent lamellae that were then welded to the liquid cell substrate. This technique, for which we have described in detail all of the fabrication steps, allows for samples having dimensions of several square micrometres to be observed by TEM in situ in a liquid. We applied this technique to test whether we could observe and measure in situ dissolution of a crystalline material called wollastonite, a calcium silicate mineral. More specifically, this study was used to observe and record surface dynamics associated with step and terrace edge movement, which are ultimately linked to the overall rate of dissolution. The wollastonite lamella underwent chemical reactions in pure deionized water at ambient temperature in a liquid cell with a 5-math formulam-spacer thickness. The movement of surface steps and terraces was measured periodically over a period of almost 5 h. Quite unexpectedly, the one-dimensional rates of retreat of these surface features were not constant, but changed over time. In addition, there were noticeable quantitative differences in retreat rates as a function crystallographic orientation, indicating that surface retreat is anisotropic. Several bulk rates of dissolution were also determined (1.6–4.2 • 10 -7 mol m -2 s -1 ) using the rates of retreat of representative terraces and steps, and were found to be within one order of magnitude of dissolution rates in the literature based on aqueous chemistry data.
Relativistic Mott insulators known as “Kitaev materials” potentially realize spin liquids hosting non-Abelian anyons. Motivated by fault-tolerant quantum-computing applications in this setting, we introduce a dynamical anyon-generation protocol that exploits universal edge physics. The setup features holes in the spin liquid, which define energetically cheap locations for non-Abelian anyons, connected by a narrow bridge that can be tuned between spin liquid and topologically trivial phases. We show that modulating the bridge from trivial to spin liquid over intermediate time scales—quantified by analytics and extensive simulations—deposits non-Abelian anyons into the holes with O(1) probability. The required bridge manipulations can be implemented by integrating the Kitaev material into magnetic tunnel junction arrays that engender locally tunable exchange fields. Combined with existing readout strategies, our protocol reveals a path to topological qubit experiments in Kitaev materials at zero applied magnetic field.
Heterogeneous nucleation of water droplet on surfaces with different solid-liquid interaction intensities is investigated by molecular dynamics simulation. The interaction potentials between surface atoms and vapor molecules are adjusted to obtain various surface free energies, and the nucleation process and wetting state of nuclei on surfaces are investigated. The results indicate that near-constant contact angles are already established for nano-scale nuclei on various surfaces, with the contact angle decreasing with solid-liquid interaction intensities linearly. Meanwhile, noticeable fluctuation of vapor-liquid interfaces can be observed for the nuclei that deposited on surfaces, which is caused by the asymmetric forces from vapor molecules. The formation and growth rate of nuclei are increasing with the solid-liquid interaction intensities. For low energy surface, the attraction of surface atoms to water molecules is comparably weak, and the pre-existing clusters can depart from the surface and enter into the bulk vapor phase. The distribution of clusters within the bulk vapor phase becomes competitive as compared with that absorbed on surface. For moderate energy surfaces, heterogeneous nucleation predominates and the formation of clusters within bulk vapor phase is suppressed. The effect of high energy particles that embedded in low energy surface is also discussed under the same simulation system. The nucleation preferably initiates on the high energy particles, and the clusters that formed on the heterogeneous particles are trapped around their original positions instead of migrating around as that observed on smooth surfaces. This feature makes it possible for the heterogeneous particles to act as fixed nucleation sites, and simulation results also suggest that the number of nuclei increases monotonously with the number of high energy particles. The growth of nuclei on high energy particles can be divided into three sub-stages, beginning with the formation of a wet-spot, increase of contact angle with near-constant contact line, and finally growth with constant contact angle. The growth rate of nuclei also increases with the size of high energy particles.
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.
Unlike phonons in crystals, the collective excitations in liquids cannot be treated as propagation of harmonic displacements of atoms around stable local energy minima. The viscoelasticity of liquids, reflected in transition from the adiabatic to elastic high-frequency speed of sound and in absence of the long-wavelength transverse excitations, results in dispersions of longitudinal (L) and transverse (T) collective excitations essentially different from the typical phonon ones. Practically, nothing is known about the effect of high pressure on the dispersion of collective excitations in liquids, which causes strong changes in liquid structure. Here dispersions of L and T collective excitations in liquid Li in the range of pressures up to 186 GPa were studied by ab initio simulations. Two methodologies for dispersion calculations were used: direct estimation from the peak positions of the L/T current spectral functions and simulation-based calculations of wavenumber-dependent collective eigenmodes. It is found that at ambient pressure, the longitudinal and transverse dynamics are well separated, while at high pressures, the transverse current spectral functions, density of vibrational states, and dispersions of collective excitations yield evidence of two types of propagating modes that contribute strongly to transverse dynamics. Emergence of the unusually high-frequency transverse modes gives evidence of the breakdown of a regular viscoelastic theory of transverse dynamics, which is based on coupling of a single transverse propagating mode with shear relaxation. The explanation of the observed high-frequency shift above the viscoelastic value is given by the presence of another branch of collective excitations. With the pressure increasing, coupling between the two types of collective excitations is rationalized within a proposed extended viscoelastic model of transverse dynamics.
Electron transport through a nanostructure can be characterized in part using concepts from classical fluid dynamics. Hence, it is natural to ask how far the analogy can be taken and whether the electron liquid can exhibit nonlinear dynamical effects such as turbulence. Here we present an ab initio study of the electron dynamics in nanojunctions which reveals that the latter indeed exhibits behavior quite similar to that of a classical fluid. In particular, we find that a transition from laminar to turbulent flow occurs with increasing current, corresponding to increasing Reynolds numbers. These findings reveal unexpected features of electron dynamics and shed new light on our understanding of transport properties of nanoscale systems.
The structure and dynamics of a task-specific ionic liquid (TSIL), trihexyl(tetradecyl)phosphonium imidazolate, before and after absorbing CO 2 were studied with a molecular dynamics (MD) simulation. This particular ionic liquid is one of several newly discovered azole-based TSILs for equimolar CO 2 capture. Unlike other TSILs whose viscosity increases drastically upon reaction with CO 2 , its viscosity decreases after CO 2 absorption. This unique behavior was confirmed in our MD simulation. In this study, we find that after CO 2 absorption the translational dynamics of the whole system is accelerated, accompanied by an accelerated rotational dynamics of the cations. Radial distribution function and spatial distribution function analyses show that the anions become asymmetric after reaction with CO 2 , and this causes the imbalance of the interaction between the positive and negative regions of the ions. The interaction between the phosphorus atom of the cation and oxygen atoms of the carboxyl group on the anion is enhanced, while that between the phosphorus atom and the naked nitrogen atom of the anion is weakened. The ion-pair correlation functions further support that the weakened interaction leads to faster dissociation of cation–anion pairs, thereby causing an accelerated dynamics. Hence, the asymmetry of anions influences the dynamics of the system and affects the viscosity. In conclusion, this insight may help design better TSILs with decreased viscosity for CO 2 capture.
First principles molecular dynamics simulation protocol is established using revised functional of Perdew-Burke-Ernzerhof (revPBE) in conjunction with Grimme's third generation of dispersion (D3) correction to describe properties of water at ambient conditions. This study also demonstrates the consistency of the structure of water across both isobaric (NpT) and isothermal (NVT) ensembles. Going beyond the standard structural benchmarks for liquid water, we compute properties that are connected to both local structure and mass density fluctuations that are related to concepts of solvation and hydrophobicity. Finally, we directly compare our revPBE results to the Becke-Lee-Yang-Parr (BLYP) plus Grimme dispersion corrections (D2) and both the empirical fixed charged model (SPC/E) and many body interaction potential model (MB-pol) to further our understanding of how the computed properties herein depend on the form of the interaction potential.
Using molecular dynamics simulations, we investigate the temperature-dependent evolution of the first peak position/shape in pair distribution functions of liquids. For metallic liquids, the peak skews towards the left (shorter distance side) with increasing temperature, similar to the previously reported anomalous peak shift. Making use of constant-volume simulations in the absence of thermal expansion and change in inherent structure, we demonstrate that the apparent shift of the peak maximum can be a result of the asymmetric shape of the peak, as the asymmetry increases with temperature-induced spreading of neighboring atoms to shorter and longer distances due to the anharmonic nature of the interatomic interaction potential. These findings shed light on the first-shell expansion/contraction paradox for metallic liquids, aside from possible changes in local topological or chemical short-range ordering. The melts of covalent materials are found to exhibit an opposite trend of peak shift, which is attributed to an effect of the directionality of the interatomic bonds.
A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM - ) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM - with the local environment and dynamics of TCM - have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM - in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM - in all ILs exhibit a nonexponential behavior with a short time component of ~2 ps and a long time component spanning ~9 ps to 14 ps. The TCM - vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM - in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM - and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM - based ILs compared to commonly used ILs. In addition, this work demonstrates that TCM - can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.
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.
We compare predictions from two familiar models of the metastable supercooled liquid, respectively, constructed with thermodynamic and dynamic approaches. In the so called density functional theory the free energy F[ρ] of the liquid is a functional of the inhomogeneous density ρ(r). The metastable state is identified as a local minimum of F[ρ]. The sharp density profile characterizing ρ(r) is identified as a single particle oscillator, whose frequency is obtained from the parameters of the optimum density function. On the other hand, a dynamic approach to supercooled liquids is taken in the mode coupling theory (MCT) which predict a sharp ergodicity-non-ergodicity transition at a critical density. The single particle dynamics in the non-ergodic state, treated approximately, represents a propagating mode whose characteristic frequency is computed from the corresponding memory function of the MCT. The mass localization parameters in the above two models (treated in their simplest forms) are obtained, respectively, in terms of the corresponding natural frequencies depicted and are shown to have comparable magnitudes.
Liquid sodium tunnel for magneto-fluid-dynamic flows
We present multiscale molecular dynamics simulation results comparing the conductivity and transference numbers in lithium salt-doped polymeric ionic liquids (PolyILs) and the lamellar phases of block copolymeric ionic liquids (coPolyILs). In both systems, the anion mobilities decreased with salt loading. Lithium ions exhibited negative mobilities in both systems, but the magnitudes decreased with an increase in salt concentrations. More interestingly, the anion mobilities were lower in the lamellar systems compared to homopolymers in magnitude, but the lithium ion mobilities and transference numbers were less negative in such systems. We examine the anion–cation and lithium–anion interactions in terms of radial distribution functions, coordination characteristics, and ion-pair relaxation timescales. Based on such analyses, we rationalize the salt concentration dependencies as a result of the interfacial interactions in lamellar systems and the competition between anion–cation and lithium–anion interactions in both PolyILs and coPolyILs. Altogether, the findings presented in this study demonstrate that the modified anion–cation and lithium–anion interactions in the microphase-separated coPolyILs may provide a strategy for realizing higher lithium ion transference numbers relative to the homopolymeric counterparts.
The unsteady evaporation of a thin planar liquid film is studied by molecular dynamics simulations of Lennard-Jones fluid. The obtained results are compared with the predictions of a diffuse interface model in which capillary Korteweg contributions are added to hydrodynamic equations, in order to obtain a unified description of the liquid bulk, liquid-vapor interface and vapor region. Particular care has been taken in constructing a diffuse interface model matching the thermodynamic and transport properties of the Lennard-Jones fluid. The comparison of diffuse interface model and molecular dynamics results shows that, although good agreement is obtained in equilibrium conditions, remarkable deviations of diffuse interface model predictions from the reference molecular dynamics results are observed in the simulation of liquid film evaporation. It is also observed that molecular dynamics results are in good agreement with preliminary results obtained from a composite model which describes the liquid film by a standard hydrodynamic model and the vapor by the Boltzmann equation. The two mathematical model models are connected by kinetic boundary conditions assuming unit evaporation coefficient.