Engineering Papers⌕ Search

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.

At least 541 records · Page 30

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↗

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↗

Nanoparticle generation and transport resulting from femtosecond laser ablation of ultrathin metal films: Time-resolved measurements and molecular dynamics simulations

The synthesis of metal nanoparticles by ultrafast laser ablation of nanometers-thick metal films has been studied experimentally and computationally. Near-threshold backside laser ablation of 2–20 nm-thick Pt films deposited on fused silica substrates was found to produce nanoparticles with size distributions that were bimodal for the thicker films, but collapsed into a single mode distribution for the thinnest film. Time-resolved imaging of blackbody emission from the Pt nanoparticles was used to reveal the nanoparticle propagation dynamics and estimate their temperatures. The observed nanoparticle plume was compact and highly forward-directed with a well-defined collective velocity that permitted multiple rebounds with substrates to be revealed. Large-scale molecular dynamics simulations were used to understand the evolution of compressive and tensile stresses in the thicker melted liquid films that lead to their breakup and ejection of two groups of nanoparticles with different velocity and size distributions. Ultrafast laser irradiation of ultrathin (few nm) metal films avoids the splitting of the film and appears to be a method well-suited to cleanly synthesize and deposit nanoparticles from semitransparent thin film targets in highly directed beams.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Constant pressure and temperature discrete-time Langevin molecular dynamics

We present a new and improved method for simultaneous control of temperature and pressure in molecular dynamics simulations with periodic boundary conditions. The thermostat-barostat equations are built on our previously developed stochastic thermostat, which has been shown to provide correct statistical configurational sampling for any time step that yields stable trajectories. Here, we extend the method and develop a set of discrete-time equations of motion for both particle dynamics and system volume in order to seek pressure control that is insensitive to the choice of the numerical time step. The resulting method is simple, practical, and efficient. The method is demonstrated through direct numerical simulations of two characteristic model systems—a one-dimensional particle chain for which exact statistical results can be obtained and used as benchmarks, and a three-dimensional system of Lennard-Jones interacting particles simulated in both solid and liquid phases. The results, which are compared against the method of Kolb and Dünweg [J. Chem. Phys. 111, 4453 (1999)], show that the new method behaves according to the objective, namely that acquired statistical averages and fluctuations of configurational measures are accurate and robust against the chosen time step applied to the simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interplay between the structure and dynamics in liquid and undercooled boron: An ab initio molecular dynamics simulation study

In the present work, the structural and dynamic properties of liquid and undercooled boron are investigated by means of ab initio molecular dynamics simulation. Our results show that both liquid and undercooled states present a well pronounced short-range order (SRO) mainly due to the formation of inverted umbrella structural units. Moreover, we observe the development of a medium-range order (MRO) in the undercooling regime related to the increase of inverted umbrella structural units and of their interconnection as the temperature decreases. We also evidence that this MRO leads to a partial crystallization in the β-rhombohedral crystal below T = 1900 K. Finally, we discuss the role played by the SRO and MRO in the nearly Arrhenius evolution of the diffusion and the non-Arrhenius temperature dependence of the shear viscosity, in agreement with the experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-resolved terahertz dynamics in thin films of the topological insulator Bi 2 Se 3

We use optical pump–THz probe spectroscopy at low temperatures to study the hot carrier response in thin Bi 2 Se 3 films of several thicknesses, allowing us to separate the bulk from the surface transient response. We find that for thinner films the photoexcitation changes the transport scattering rate and reduces the THz conductivity, which relaxes within 10 picoseconds (ps). For thicker films, the conductivity increases upon photoexcitation and scales with increasing both the film thickness and the optical fluence, with a decay time of approximately 5 ps as well as a much higher scattering rate. Furthermore, these different dynamics are attributed to the surface and bulk electrons, respectively, and demonstrate that long-lived mobile surface photo-carriers can be accessed independently below certain film thicknesses for possible optoelectronic applications.

36 MATERIALS SCIENCE↗

Mechanical and dynamical behaviors of ZrSi and ZrSi{sub 2} bulk metallic glasses: A molecular dynamics study

The mechanical and dynamical properties of ZrSi and ZrSi{sub 2} bulk metallic glasses (BMGs) have been investigated by molecular dynamics simulation. The Honeycutt-Anderson (HA) index analysis indicates that the major indexes in ZrSi and ZrSi{sub 2} bulk metallic glasses are 1551, 1541, and 1431, which refers to the liquid structure. For uniaxial tension, the results show that the ZrSi and ZrSi{sub 2} BMGs are more ductile than their crystal counterparts. The evolution of the distribution of atomic local shear strain clearly shows the initialization of shear transformation zones (STZs), the extension of STZs, and the formation of shear bands along a direction 45° from the tensile direction when the tensile strain gradually increases. The self-diffusion coefficients of ZrSi and ZrSi{sub 2} BMGs at temperatures near their melting points were calculated by the Einstein equation according to the slopes of the MSD profiles at the long-time limit. Because the HA fraction summation of icosahedral-like structures of ZrSi BMG is higher than that of ZrSi{sub 2} BMG, and these local structures are more dense, the self-diffusion coefficients of the total, Zr, and Si atoms of ZrSi{sub 2} BMG are larger than those of ZrSi BMG. This can be attributed to the cage effect, where a denser local structure has a higher possibility of atoms jumping back to form a backflow and then suppress atomic diffusivity. For ZrSi{sub 2} BMG, the self-diffusion coefficient of Si increases with temperature more significantly than does that of Zr, because more open packing rhombohedra structures are formed by the Si-Si pair.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations

Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into a few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation, which can be rigorously derived by means of a linear projection technique. Within this framework, a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here, we discuss that this task is more naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD{sub 3}CN treated with a parallel multi-state EVB model

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD{sub 3}CN → DF + CD{sub 2}CN reaction in CD{sub 3}CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD{sub 3}CN solvent, equilibrium power spectra of DF in CD{sub 3}CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol{sup −1} localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD{sub 3}CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-equilibrium dynamics in disordered materials: Ab initio molecular dynamics simulations

The dynamic properties of liquid B{sub 2}O{sub 3} under pressure and highly-charged bromophenol molecule are studied by using molecular dynamics (MD) simulations based on density functional theory (DFT). Diffusion properties of covalent liquids under high pressure are very interesting in the sense that they show unexpected pressure dependence. It is found from our simulation that the magnitude relation of diffusion coefficients for boron and oxygen in liquid B{sub 2}O{sub 3} shows the anomalous pressure dependence. The simulation clarified the microscopic origin of the anomalous diffusion properties. Our simulation also reveals the dissociation mechanism in the coulomb explosion of the highly-charged bromophenol molecule. When the charge state n is 6, hydrogen atom in the hydroxyl group dissociates at times shorter than 20 fs while all hydrogen atoms dissociate when n is 8. After the hydrogen dissociation, the carbon ring breaks at about 100 fs. There is also a difference on the mechanism of the ring breaking depending on charge states, in which the ring breaks with expanding (n = 6) or shrink (n = 8)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of configurational changes on molecular dynamics in polyvinylidene fluoride and poly(vinylidene fluoride-trifluoroethylene) ferroelectric polymers

Here, we present a comparative study of proton dynamics in unpoled non-ferroelectric polymer polyvinylidene fluoride (PVDF) and in its trifluoroethylene containing ferroelectric copolymer (with 70/30 molar proportion), using quasi-elastic neutron scattering. The neutron data reveal the existence of two distinct types of molecular motions in the temperature range investigated. Moreover, the slower motion, which is characterized in details here, is ascribed to protons jump diffusion along the polymeric carbon chains, while the faster motion could be attributed to localized rotational motion of methylene groups. At temperatures below the Curie point (T-c similar to 385 K) of the composite polymer, the slower diffusive mode experiences longer relaxation times in the ferroelectric blend than in the bare PVDF, although the net corresponding diffusion coefficient remains comparatively the same in both polymers with characteristic activation energy of E-A approximate to 27-33 kJ/mol. This arises because of a temperature dependent jump length r(0), which we observe to be effectively longer in the copolymer, possibly due to the formation of ordered ferroelectric domains below Tc. Above Tc, there is no appreciable difference in r(0) between the two systems. Our observation directly relates the known dependence of Tc on molar ratio to changes in r(0), providing fundamental insight into the ferroelectric properties of PVDF-based copolymers.

36 MATERIALS SCIENCE↗