Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ionic forces”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Long range electrostatic forces in ionic liquids

Ionic liquids are pure salts that are liquid under ambient conditions. As liquids composed solely of ions, the scientific consensus has been that ionic liquids have exceedingly high ionic strengths and thus very short Debye screening lengths. However, several recent experiments from laboratories around the world have reported data for the approach of two surfaces separated by ionic liquids which revealed remarkable long range forces that appear to be electrostatic in origin. Evidence has accumulated demonstrating long range surface forces for several different combinations of ionic liquids and electrically charged surfaces, as well as for concentrated mixtures of inorganic salts in solvent. The original interpretation of these forces, that ionic liquids could be envisioned as “dilute electrolytes,” was controversial, and the origin of long range forces in ionic liquids remains the subject of discussion. Here we seek to collate and examine the evidence for long range surface forces in ionic liquids, identify key outstanding questions, and explore possible mechanisms underlying the origin of these long range forces. Long range surface forces in ionic liquids and other highly concentrated electrolytes hold diverse implications from designing ionic liquids for energy storage applications to rationalizing electrostatic correlations in biological self-assembly.

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

Stopping power beyond the adiabatic approximation

Energetic ions traveling in solids deposit energy in a variety of ways, being nuclear and electronic stopping the two avenues in which dissipation is usually treated. This separation between electrons and ions relies on the adiabatic approximation in which ions interact via forces derived from the instantaneous electronic ground state. In a more detailed view, in which non-adiabatic effects are explicitly considered, electronic excitations alter the atomic bonding, which translates into changes in the interatomic forces. In this work, we use time dependent density functional theory and forces derived from the equations of Ehrenfest dynamics that depend instantaneously on the time-dependent electronic density. With them we analyze how the inter-ionic forces are affected by electronic excitations in a model of a Ni projectile interacting with a Ni target, a metallic system with strong electronic stopping and shallow core level states. We find that the electronic excitations induce substantial modifications to the inter-ionic forces, which translate into nuclear stopping power well above the adiabatic prediction. Particularly, we observe that most of the alteration of the adiabatic potential in early times comes from the ionization of the core levels of the target ions, not readily screened by the valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Dissecting ion-specific dielectric spectra of sodium-halide solutions into solvation water and ionic contributions

Using extensive equilibrium molecular dynamics simulations we determine the dielectric spectra of aqueous solutions of NaF, NaCl, NaBr, and NaI. The ion-specific and concentration-dependent shifts of the static dielectric constants and the dielectric relaxation times match experimental results very well, which serves as a validation of the classical and non-polarizable ionic force fields used. The purely ionic contribution to the dielectric response is negligible, but determines the conductivity of the salt solutions. The ion-water cross correlation contribution is negative and reduces the total dielectric response by about 5%-10% for 1 M solutions. The dominating water dielectric response is decomposed into different water solvation shells and ion-pair configurations, by this the spectral blue shift and the dielectric decrement of salt solutions with increasing salt concentration is demonstrated to be primarily caused by first-solvation shell water. With rising salt concentration the simulated spectra show more pronounced deviations from a single-Debye form and can be well described by a Cole-Cole fit, in quantitative agreement with experiments. Our spectral decomposition into ionic and different water solvation shell contributions does not render the individual contributions more Debye-like, this suggests the non-Debye-like character of the dielectric spectra of salt solutions not to be due to the superposition of different elementary relaxation processes with different relaxation times. Rather, the non-Debye-like character is likely to be an inherent spectral signature of solvation water around ions.

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

DFT-FE – A massively parallel adaptive finite-element code for large-scale density functional theory calculations

We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching ~100,000 electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over 5-10 fold speedups in systems with more than 10,000 electrons. In conclusion, the benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗
DOE OSTI2019

Organophilic, Superparamagnetic, and Reversibly Thermoresponsive Silica-Polypeptide Core–Shell Particles

Particles with a superparamagnetic cobalt inner core, silica outer core, and covalently bound homopolypeptide shell were investigated under thermal and magnetic stimuli. The homopolypeptide was poly(ε-carbobenzyloxy- l -lysine), PCBL, which is known to exhibit a thermoreversible coil ⇔ helix transition when dissolved as a pure polymer in m- cresol. Tethering to a core particle did not prevent PCBL from undergoing this conformational transition, as confirmed by dynamic light scattering and optical rotation, but the transition was broadened compared to that of the untethered polymer. The Co@SiO 2 -PCBL hybrid particles retained the superparamagnetic properties of the cobalt inner nougat. Indeed, some response remains even after aging for >5 years. The aged PCBL shell also preserved its responsiveness to temperature, although differences in the shape of the size vs temperature transition curve were determined compared to the freshly made particles. A reversible coil ⇔ helix transition for a particle-bound polypeptide in a pure organic solvent is rare. Further to providing a convenient tool for characterizing coil ⇔ helix transitions for surface-bound polypeptides without interference from pH or the strong ionic forces that dominate behavior in aqueous systems, the Co@SiO 2 -PCBL/ m -cresol system may prove useful in studies of the effect of shell polymer conformation on colloid interactions. The stability of the magnetic core and polypeptide shell suggest a long shelf life for Co@SiO 2 -PCBL, which can, in principle, be deprotected to yield positively charged Co@SiO 2 -poly( l -lysine) particles for possible transfection or antimicrobial applications or chained magnetically to produce responsive poly(colloids).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Theory of electron-ion energy transfer applied to laser ablation

Electron-ion temperature equilibration will be considered for states of matter typically occurring during laser ablation that is for strongly heated electrons and densities similar and lower than the one in the solid. The focus will be on fluid systems after the melt of the lattice as these dominate the hydrodynamic phase of the ablation process. It is shown that simple expressions derived for dilute gases and plasmas become highly questionable for solid-density fluids whereas the electron-phonon picture is inapplicable as well. Instead, one has to consider the collective behavior of electrons and ions in the strongly coupled fluid. The strong inter-ionic forces also strongly affect the heat capacities applied in a two-temperature model. The results presented here are based on a quantum-statistical approach applying nonequilibrium Green's functions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

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

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion–ion and ion–solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K + ion. Diffusion-based Nernst–Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

25 ENERGY STORAGE↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2016

Examination of ionic wind and cathode sheath effects in a E-field premixed flame with ion density measurements

The effect of the ionic wind on a premixed methane-air flame under a DC electric field is studied via mapping of the ion density with Langmuir probes. Ion densities were observed to increase near the burner with increasing electrode voltage up to 6 kV. Past this electrode supply voltage, ion densities ceased increasing and began to decline in some locations within the premixed flame. The increased ion density is caused by an increase in ionic wind force and cathode sheath thickness. The plateau in density is due to the cathode sheath fully encompassing the flame front which is the ion source, thereby collecting all ions in the flame. The spatial density data support the ionic wind hypothesis and provide further explanation of its limits based on the plasma sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Carbon dioxide in an ionic liquid: Structural and rotational dynamics

Ionic liquids (ILs), which have widely tunable structural motifs and intermolecular interactions with solutes, have been proposed as possible carbon capture media. In order to inform the choice of an optimal ionic liquid system, it can be useful to understand the details of dynamics and interactions on fundamental time scales (femtoseconds to picoseconds) of dissolved gases, particularly carbon dioxide (CO 2 ), within the complex solvation structures present in these uniquely organized materials. The rotational and local structural fluctuation dynamics of CO 2 in the room temperature ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EmimNTf2) were investigated by using ultrafast infrared spectroscopy to interrogate the CO 2 asymmetric stretch. Polarization-selective pump probe measurements yielded the orientational correlation function of the CO 2 vibrational transition dipole. It was found that reorientation of the carbon dioxide occurs on 3 time scales: 0.91 ± 0.03, 8.3 ± 0.1, 54 ± 1 ps. The initial two are attributed to restricted wobbling motions originating from a gating of CO 2 motions by the IL cations and anions. The final (slowest) decay corresponds to complete orientational randomization. Two-dimensional infrared vibrational echo (2D IR) spectroscopy provided information on structural rearrangements, which cause spectral diffusion, through the time dependence of the 2D line shape. Analysis of the time-dependent 2D IR spectra yields the frequency-frequency correlation function (FFCF). Polarization-selective 2D IR experiments conducted on the CO 2 asymmetric stretch in the parallel- and perpendicular-pumped geometries yield significantly different FFCFs due to a phenomenon known as reorientation-induced spectral diffusion (RISD), revealing strong vector interactions with the liquid structures that evolve slowly on the (independently measured) rotation time scales. To separate the RISD contribution to the FFCF from the structural spectral diffusion contribution, the previously developed first order Stark effect RISD model is reformulated to describe the second order (quadratic) Stark effect—the first order Stark effect vanishes because CO 2 does not have a permanent dipole moment. Through this analysis, we characterize the structural fluctuations of CO 2 in the ionic liquid solvation environment, which separate into magnitude-only and combined magnitude and directional correlations of the liquid’s time dependent electric field. This new methodology will enable highly incisive comparisons between CO 2 dynamics in a variety of ionic liquid systems.

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

Clean surface processing of rubrene single crystal immersed in ionic liquid by using frequency modulation atomic force microscopy

Surface processing of a rubrene single crystal immersed in ionic liquids is valuable for further development of low voltage transistors operated by an electric double layer. We performed a precise and clean surface processing based on the tip-induced dissolution of rubrene molecules at the ionic liquid/rubrene single crystal interfaces by using frequency modulation atomic force microscopy. Molecular resolution imaging revealed that the tip-induced dissolution proceeded via metastable low density states derived from the anisotropic intermolecular interactions within the crystal structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗