SEARCH · Engineering Papers
Results for “dynamic”
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.
Massively parallel simulations of relativistic fluid dynamics on graphics processing units with CUDA
Not provided.
A generalized parallel replica dynamics
Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming–Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated with various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.
Nitrate bioreduction dynamics in hyporheic zone sediments under cyclic changes of chemical compositions
Not Available
Revertant mutation V48G alters conformational dynamics of highly drug resistant HIV protease PRS17
Explore the source record for details and available documents.
Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories
We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.
Imaging solid–electrolyte interphase dynamics using operando reflection interference microscopy
The quality of the solid–electrolyte interphase is crucial for the performance of most battery chemistries, but its formation dynamics during operation are not well understood due to a lack of reliable operando characterization techniques. Herein, we report a dynamic, non-invasive, operando reflection interference microscope to enable the real-time imaging of the solid–electrolyte interphase during its formation and evolution processes with high sensitivity. The stratified structure of the solid–electrolyte interphase formed during four distinct steps includes the emergence of a permanent inner inorganic layer enriched in LiF, a transient assembly of an interfacial electrified double layer and a consequent emergence of a temporary outer organic-rich layer whose presence is reversible with electrochemical cycling. Reflection interference microscope imaging reveals an inverse correlation between the thicknesses of two interphasial subcomponents, implying that the permanent inorganic-rich inner layer dictates the organic-rich outer layer formation and lithium nucleation. In conclusion, the real-time visualization of solid–electrolyte interphase dynamics provides a powerful tool for the rational design of battery interphases.
Dynamics of water-alcohol mixtures: Insights from nuclear magnetic resonance, broadband dielectric spectroscopy, and triplet solvation dynamics
We combine {sup 2}H nuclear magnetic resonance (NMR), broadband dielectric spectroscopy (BDS), and triplet solvation dynamics (TSD) to investigate molecular dynamics in glass-forming mixtures of water and propylene glycol in very broad time and temperature ranges. All methods yield consistent results for the α process of the studied mixtures, which hardly depends on the composition and shows Vogel-Fulcher temperature dependence as well as Cole-Davidson spectral shape. The good agreement between BDS and TDS data reveals that preferential solvation of dye molecules in microheterogeneous mixtures does not play an important role. Below the glass transition temperature T{sub g}, NMR and BDS studies reveal that the β process of the mixtures shows correlation times, which depend on the water concentration, but exhibit a common temperature dependence, obeying an Arrhenius law with an activation energy of E{sub a} = 0.54 eV, as previously reported for mixtures of water with various molecular species. Detailed comparison of NMR and BDS correlation functions for the β process unravels that the former decay faster and more stretched than the latter. Moreover, the present NMR data imply that propylene glycol participates in the β process and, hence, it is not a pure water process, and that the mechanism for molecular dynamics underlying the β process differs in mixtures of water with small and large molecules.
Coherence penalty functional: A simple method for adding decoherence in Ehrenfest dynamics
We present a new semiclassical approach for description of decoherence in electronically non-adiabatic molecular dynamics. The method is formulated on the grounds of the Ehrenfest dynamics and the Meyer-Miller-Thoss-Stock mapping of the time-dependent Schrödinger equation onto a fully classical Hamiltonian representation. We introduce a coherence penalty functional (CPF) that accounts for decoherence effects by randomizing the wavefunction phase and penalizing development of coherences in regions of strong non-adiabatic coupling. The performance of the method is demonstrated with several model and realistic systems. Compared to other semiclassical methods tested, the CPF method eliminates artificial interference and improves agreement with the fully quantum calculations on the models. When applied to study electron transfer dynamics in the nanoscale systems, the method shows an improved accuracy of the predicted time scales. The simplicity and high computational efficiency of the CPF approach make it a perfect practical candidate for applications in realistic systems.
Structural, dynamic, electronic, and vibrational properties of flexible, intermediate, and stressed rigid As-Se glasses and liquids from first principles molecular dynamics
The structural, vibrational, electronic, and dynamic properties of amorphous and liquid As{sub x}Se{sub 1-x} (0.10 <x < 0.45) are studied by First Principles Molecular Dynamics. Within the above range of compositions, thresholds, and anomalies are found in the behavior of reciprocal and real space properties that can be correlated to the experimental location of the Boolchand intermediate phase in these glassy networks, observed at 0.27 <x < 0.37. These findings are associated with diffusion anomalies for the parent liquid phase, thereby linking structural and dynamical atomic-scale fingerprints for the onset of rigidity within the network, while also providing a much more complex picture than the one derived from mean-field approaches of stiffness transitions.
Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity
In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.
Time lagged ordinal partition networks for capturing dynamics of continuous dynamical systems
We investigate a generalised version of the recently proposed ordinal partition time series to network transformation algorithm. First, we introduce a fixed time lag for the elements of each partition that is selected using techniques from traditional time delay embedding. The resulting partitions define regions in the embedding phase space that are mapped to nodes in the network space. Edges are allocated between nodes based on temporal succession thus creating a Markov chain representation of the time series. We then apply this new transformation algorithm to time series generated by the Rössler system and find that periodic dynamics translate to ring structures whereas chaotic time series translate to band or tube-like structures—thereby indicating that our algorithm generates networks whose structure is sensitive to system dynamics. Furthermore, we demonstrate that simple network measures including the mean out degree and variance of out degrees can track changes in the dynamical behaviour in a manner comparable to the largest Lyapunov exponent. We also apply the same analysis to experimental time series generated by a diode resonator circuit and show that the network size, mean shortest path length, and network diameter are highly sensitive to the interior crisis captured in this particular data set.
Effects of dynamic heterogeneity and density scaling of molecular dynamics on the relationship among thermodynamic coefficients at the glass transition
In this paper, we define and experimentally verify thermodynamic characteristics of the liquid-glass transition, taking into account a kinetic origin of the process. Using the density scaling law and the four-point measure of the dynamic heterogeneity of molecular dynamics of glass forming liquids, we investigate contributions of enthalpy, temperature, and density fluctuations to spatially heterogeneous molecular dynamics at the liquid-glass transition, finding an equation for the pressure coefficient of the glass transition temperature, dTg/dp. This equation combined with our previous formula for dTg/dp, derived solely from the density scaling criterion, implies a relationship among thermodynamic coefficients at Tg. Since this relationship and both the equations for dTg/dp are very well validated using experimental data at Tg, they are promising alternatives to the classical Prigogine-Defay ratio and both the Ehrenfest equations in case of the liquid-glass transition.
Effect of water on structure and dynamics of [BMIM][PF{sub 6}] ionic liquid: An all-atom molecular dynamics simulation investigation
Composition dependent structural and dynamical properties of aqueous hydrophobic 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF{sub 6}]) ionic liquid (IL) have been investigated by using all-atom molecular dynamics simulation. We observe that addition of water does not increase significant number of dissociated ions in the solution over the pure state. As a consequence, self-diffusion coefficient of the cation and anion is comparable to each other at all water concentration similar to that is observed for the pure state. Voronoi polyhedra analysis exhibits strong dependence on the local environment of IL concentration. Void and neck distributions in Voronoi tessellation are approximately Gaussian for pure IL but upon subsequent addition of water, we observe deviation from the Gaussian behaviour with an asymmetric broadening with long tail of exponential decay at large void radius, particularly at higher water concentrations. The increase in void space and neck size at higher water concentration facilitates ionic motion, thus, decreasing dynamical heterogeneity and IL reorientation time and increases self-diffusion coefficient significantly.
Structural evolution and atomic dynamics in Ni–Nb metallic glasses: A molecular dynamics study
In this work, the composition and temperature dependence of static and dynamic structures in Ni x Nb 1– x ( x = 50–70 at. %) were systematically studied using molecular dynamics with a new-released semi-empirical embedded atom method potential by Mendelev. The calculated pair correlation functions and the structure factor match well with the experimental data, demonstrating the reliability of the potential within relatively wide composition and temperature ranges. The local atomic structures were then characterized by bond angle distributions and Voronoi tessellation methods, demonstrating that the icosahedral $\langle$0,0,12,0$\rangle$ is only a small fraction in the liquid state but increases significantly during cooling and becomes dominant at 300 K. The most abundant clusters are identified as $\langle$0,0,12,0$\rangle$ and distorted icosahedron$\langle$0,2,8,2$\rangle$. The large fraction of these two clusters hints that the relatively good glass forming ability is near the eutectic point. Unlike Cu–Zr alloys, both the self-diffusion coefficient and shear viscosity are insensitive to compositions upon cooling in Ni–Nb alloys. The breakdown of the Stokes-Einstein relation happens at around 1.6 T g ( T g : glass transition temperature). In the amorphous state, the solid and liquid-like atoms can be distinguished based on the Debye-Waller factor $\langle$ u 2 $\rangle$. The insensitivity of the dynamic properties of Ni–Nb alloys to compositions may result from the relatively simple solidification process in the phase diagram, in which only one eutectic point exists in the studied composition range.
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.