Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 541 records · Page 30

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↗

Monte Carlo Simulation of Borosilicate Glass Dissolution using Molecular Dynamics-Generated Glass Structures

Atomic-scale insights into the dissolution mechanisms of borosilicate glasses are critical to the development of the waste form corrosion models that inform performance assessments of nuclear waste disposal facilities. Monte Carlo (MC) simulations have been instrumental in advancing the state of knowledge of borosilicate glass dissolution; however, all such simulations to date have relied on lattices to represent the amorphous glass structure, thereby potentially introducing artifacts. Therefore, a new approach, referred to as the amorphous MC approach, has been developed whereby glass structures generated from molecular dynamics simulations are used as starting points for MC simulations of glass dissolution. MC simulations of sodium borosilicate glasses that covered a wide compositional range were performed and revealed correlations between dissolution rate and structural features of the glasses. Notably, the lattice MC approach predicted the fraction of tetrahedral boron to have a lesser influence on the dissolution rate than computed with the amorphous MC approach. The lattice MC simulations also overestimated the dissolution rate for a given glass composition. These findings highlight the importance of using truly amorphous structures in MC simulations of glass corrosion.

Kerisit, Sebastien N.↗

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗

Modeling error PDF optimization based wavelet neural network modeling of dynamic system and its application in blast furnace ironmaking

In general, the modeling errors of dynamic system model are a set of random variables. The traditional performance index of modeling such as means square error (MSE) and root means square error (RMSE) can not fully express the connotation of modeling errors with stochastic characteristics both in the dimension of time domain and space domain. Therefore, the probability density function (PDF) is introduced to completely describe the modeling errors in both time scales and space scales. Based on it, a novel wavelet neural network (WNN) modeling method is proposed by minimizing the two-dimensional (2D) PDF shaping of modeling errors. First, the modeling error PDF by the tradional WNN is estimated using data-driven kernel density estimation (KDE) technique. Then, the quadratic sum of 2D deviation between the modeling error PDF and the target PDF is utilized as performance index to optimize the WNN model parameters by gradient descent method. Since the WNN has strong nonlinear approximation and adaptive capability, and all the parameters are well optimized by the proposed method, the developed WNN model can make the modeling error PDF track the target PDF, eventually. Simulation example and application in a blast furnace ironmaking process show that the proposed method has a higher modeling precision and better generalization ability compared with the conventional WNN modeling based on MSE criteria. Furthermore, the proposed method has more desirable estimation for modeling error PDF that approximates to a Gaussian distribution whose shape is high and narrow.

Modelling↗

Dynamical coupled-channel models for hadron dynamics

Dynamical coupled-channel (DCC) approaches parametrize the interactions and dynamics of two and more hadrons and their response to different electroweak probes. The inclusion of unitarity, three-body channels, and other properties from scattering theory allows for a reliable extraction of resonance spectra and their properties from data. Here, we review the formalism and application of the ANL-Osaka, the Juelich-Bonn-Washington, and other DCC approaches in the context of light baryon resonances from meson, (virtual) photon, and neutrino-induced reactions, as well as production reactions, strange baryons, light mesons, heavy meson systems, exotics, and baryon-baryon interactions. Finally, we also provide a connection of the formalism to study finite-volume spectra obtained in Lattice QCD, and review applications involving modern statistical and machine learning tools.

Amplitude analysis↗

Heterogeneous CPU+GPU-Enabled Simulations for DFTB Molecular Dynamics of Large Chemical and Biological Systems

We introduce a new heterogeneous CPU+GPU-enhanced DFTB approach for the routine and efficient simulation of large chemical and biological systems. Here, compared to homogenous computing with conventional CPUs, heterogeneous computing approaches exhibit substantial performance with only a modest increase in power consumption, both of which are essential to upcoming exascale computing initiatives. We show that DFTB-based molecular dynamics is a natural candidate for heterogeneous computing since the computational bottleneck in these simulations is the diagonalization of the Hamiltonian matrix, which is performed several times during a single molecular dynamics trajectory. To thoroughly test and understand the performance of our heterogeneous CPU+GPU approach, we examine a variety of algorithmic implementations, benchmarks of different hardware configurations, and applications of this methodology on several large chemical and biological systems. Finally, to demonstrate the capability of our implementation, we conclude with a large-scale DFTB MD simulation of explicitly solvated HIV protease (3,974 atoms total) as a proof-of-concept example of an extremely large/complex system which, to the best of our knowledge, is the first time that an entire explicitly-solvated protein has been treated at a quantum-based MD level of detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photo-Induced Charge Transfer of Fullerene and Non-Fullerene Conjugated Polymer Blends via Ab Initio Excited-State Dynamics

Organic conjugated polymers (CPs) are promising candidates for organic photovoltaic (OPV) devices due to their unique tunable mechanical and optoelectronic performance. Over the last decade, optoelectronic properties of narrow band gap CPs as a blend with acceptor units are largely optimized, which leads to noticeable progress in OPV technology. However, their power conversion efficiency is still lower than their organic counterparts (i.e., silicon), limiting their practical usage. In this study, we employ ab initio molecular dynamics to explore photo-induced charge transfer (CT) of the diketopyrrolopyrrole-based polymer as a blend with non-fullerene (i.e., ITIC) and fullerene (i.e., PCBM) acceptor units. The results of charge carrier dynamics induced by selected photoexcitation show that hole density redistribution in space is much faster than electron relaxation. We track the relaxation rates of charge carriers over time, where the derivative of the difference between the rate of electron and hole implies the current density at zero voltage. Furthermore, this can be utilized to characterize the CT performance of CPs blended with different acceptor units. Relaxation rate results indicate that CP blend with ITIC promises a better PV performance, illustrating that the current computational approach opens the door to determine bulk heterojunctions’ electronic performance for OPV devices and narrowing down the list of potential donor–acceptor candidates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

4D-STEM Mapping of Nanocrystal Reaction Dynamics and Heterogeneity in a Graphene Liquid Cell

Chemical reaction kinetics at the nanoscale are intertwined with heterogeneity in structure and composition. However, mapping such heterogeneity in a liquid environment is extremely challenging. Here, in this work, we integrate graphene liquid cell (GLC) transmission electron microscopy and four-dimensional scanning transmission electron microscopy to image the etching dynamics of gold nanorods in the reaction media. Critical to our experiment is the small liquid thickness in a GLC that allows the collection of high-quality electron diffraction patterns at low dose conditions. Machine learning-based data-mining of the diffraction patterns maps the three-dimensional nanocrystal orientation, groups spatial domains of various species in the GLC, and identifies newly generated nanocrystallites during reaction, offering a comprehensive understanding on the reaction mechanism inside a nanoenvironment. This work opens opportunities in probing the interplay of structural properties such as phase and strain with solution-phase reaction dynamics, which is important for applications in catalysis, energy storage, and self-assembly.

four-dimensional scanning transmission electron mi↗

Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations

Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.

Agarwal, Garvit↗

Photocarrier Dynamics in MoTe 2 Nanofilms with 2 H and Distorted 1 T Lattice Structures

Molybdenum telluride (MoTe 2 ), an emerging layered two-dimensional (2D) material, possesses excellent phase-changing properties. Previous studies revealed its reversible transition between 2H and 1T' phases with a transition energy as small as 35 meV. Since 1T'-MoTe 2 is metallic, it can serve as an electrical contact for semiconducting 2H-MoTe 2 -based optoelectronic devices. Here, the photocarrier dynamics in MoTe 2 nanofilms synthesized by a one-step method and with coexisting multiple phases are investigated by transient absorption measurements. Both the energy relaxation time and the recombination lifetime of the excitons are shorter in the 1T'-MoTe 2 compared to its 2H phase. Furthermore, these results provide information on the different photocarrier dynamical properties of these two phases, which is important for future 2D optoelectronic and phase-change electronic devices based on MoTe 2 .

36 MATERIALS SCIENCE↗

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗

Dynamics of Back Electron Transfer in Dye-Sensitized Solar Cells Featuring 4-tert-Butyl-Pyridine and Atomic-Layer-Deposition Alumina as Surface Modifiers.

A series of dye-sensitized solar cells (DSCs) was constructed with TiO2 nanoparticles and N719 dye. The standard I-3(-)/I- redox shuttle and the Co(1,10-phenanthroline)(3)(3+/2+) shuttle were employed. DSCs were modified with atomic-layered-deposited (ALD) coatings of Al2O3 and/or with the surface-adsorbing additive 4-tert-butyl-pyridine. Current-voltage data were collected to ascertain the influence of each modification upon the back electron transfer (ET) dynamics of the DSCs. The primary effect of the additives alone or in tandem is to increase the open-circuit voltage. A second is to alter the short-circuit current density, J(SC). With dependence on the specifics of the system examined, any of a myriad of dynamics-related effects were observed to come into play, in both favorable (efficiency boosting) and unfavorable (efficiency damaging) ways. These effects include modulation of (a) charge-injection yields, (b) rates of interception of injected electrons by redox shuttles, and (c) rates of recombination of injected electrons with holes on surface-bound dyes. In turn, these influence charge-collection lengths, charge-collection yields, and onset potentials for undesired dark current. The microscopic origins of the effects appear to be related mainly to changes in driving force and/or electronic coupling for underlying component redox reactions. Perhaps surprisingly, only a minor role for modifier-induced shifts in conduction-band-edge energy was found. The combination of DSC-efficiency-relevant effects engendered by the modifiers was found to vary substantially as a function of the chemical identity of the redox shuttle employed. While types of modifiers are effective, a challenge going forward will be to construct systems in ways in which the benefits of organic and inorganic modifiers can be exploited in fully additive, or even synergistic, fashion.

Katz, Michael J.↗

Charge Hopping Dynamics along a Disordered Chain in Quantum Environments: Comparative Study of Different Rate Kernels

This work presents a computational study of charge hopping dynamics along a one-dimensional chain with Gaussian site energy disorder and linearly coupled quantum bath. Time dependent square displacements are calculated directly from numerical solutions of Pauli master equations, for five different rate kernels: exact Fermi golden rule (FGR) rate expression, stationary phase interpolation (SPI) approximation, semiclassical (SC) approximation, classical Marcus rate expression, and Miller-Abrahams expression. All results demonstrate diffusive behavior in the steady state limit. The results based on the FGR rate expression show that the charge transport in the quantum bath can be much more sensitive to the disorder than the prediction from the classical Marcus expression. While the SPI approximation captures this general trend reasonably well, the SC approximation tends to be unreliable at both quantitative and qualitative levels and becomes even worse than the classical Marcus expression under certain conditions. Furthermore, these results offer useful guidance in the choice of approximate rate kernels for larger-scale simulations and also demonstrate significant but fragile positive effects of quantum environments on the charge hopping dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of micelle-nanoparticle systems undergoing shear. A coarse-grained molecular dynamics approach

Recent experimental work has shown that polymeric micelles can template nanoparticles via interstitial sites in shear-ordered micelle solutions. In the current study, we report simulation results based on a coarse-grained molecular dynamics (CGMD) model of a solvent/polymer/nanoparticle system. Our results demonstrate the importance of polymer concentration and the micelle corona length in 2D shear-ordering of neat block copolymer solutions. Although our results do not show strong 3D ordering during shear, we find that cessation of shear allows the system to relax into a 3D configuration of greater order than without shear. It is further shown that this post-shear relaxation is strongly dependent on the length of the micelle corona. For the first time, we demonstrate the presence and importance of a flow disturbance surrounding micelles in simple shear flow at moderate Péclet numbers. This disturbance is similar to what is observed around simulated star polymers and ellipsoids. The extent of the flow disturbance increases as expected with a longer micelle corona length. It is further suggested that without proper consideration of these dynamics, a stable nanoparticle configuration would be difficult to obtain.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Inter-cage dynamics in structure I, II, and H fluoromethane hydrates as studied by NMR and molecular dynamics simulations

Prospective industrial applications of clathrate hydrates as materials for gas separation require further knowledge of cavity distortion, cavity selectivity, and defects induction by guest-host interactions. The results presented in this contribution show that under certain temperature conditions the guest combination of CH{sub 3}F and a large polar molecule induces defects on the clathrate hydrate framework that allow intercage guest dynamics. {sup 13}C NMR chemical shifts of a CH{sub 3}F/CH{sub 4}/TBME sH hydrate and a temperature analysis of the {sup 2}H NMR powder lineshapes of a CD{sub 3}F/THF sII and CD{sub 3}F/TBME sH hydrate, displayed evidence that the populations of CH{sub 4} and CH{sub 3}F in the D and D{sup ′} cages were in a state of rapid exchange. A hydrogen bonding analysis using molecular dynamics simulations on the TBME/CH{sub 3}F and TBME/CH{sub 4} sH hydrates showed that the presence of CH{sub 3}F enhances the hydrogen bonding probability of the TBME molecule with the water molecules of the cavity. Similar results were obtained for THF/CH{sub 3}F and THF/CH{sub 4} sII hydrates. The enhanced hydrogen bond formation leads to the formation of defects in the water hydrogen bonding lattice and this can enhance the migration of CH{sub 3}F molecules between adjacent small cages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗