Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical studies”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Oxidation of ligand-protected aluminum clusters: An ab initio molecular dynamics study

We report Car-Parrinello molecular dynamics simulations of the oxidation of ligand-protected aluminum clusters that form a prototypical cluster-assembled material. These clusters contain a small aluminum core surrounded by a monolayer of organic ligand. The aromatic cyclopentadienyl ligands form a strong bond with surface Al atoms, giving rise to an organometallic cluster that crystallizes into a low-symmetry solid and is briefly stable in air before oxidizing. Our calculations of isolated aluminum/cyclopentadienyl clusters reacting with oxygen show minimal reaction between the ligand and O{sub 2} molecules at simulation temperatures of 500 and 1000 K. In all cases, the reaction pathway involves O{sub 2} diffusing through the ligand barrier, splitting into atomic oxygen upon contact with the aluminum, and forming an oxide cluster with aluminum/ligand bonds still largely intact. Loss of individual aluminum-ligand units, as expected from unimolecular decomposition calculations, is not observed except following significant oxidation. These calculations highlight the role of the ligand in providing a steric barrier against oxidizers and in maintaining the large aluminum surface area of the solid-state cluster material.

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

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2017

Can Excited State Electronic Coherence Be Tuned via Molecular Structural Modification? A First-Principles Quantum Electronic Dynamics Study of Pyrazolate-Bridged Pt(II) Dimers

Materials and molecular systems exhibiting long-lived electronic coherence can facilitate coherent transport, opening the door to efficient charge and energy transport beyond traditional methods. Recently, signatures of a possible coherent, recurrent electronic motion were identified in femtosecond pump-probe spectroscopy experiments on a binuclear platinum complex, where a persistent periodic beating in the transient absorption signal’s anisotropy was observed. In this study, we investigate the excitonic dynamics that underlie the suspected electronic coherence for a series of binuclear platinum complexes exhibiting a range of interplatinum distances. Results suggest that the long-lived coherence can only result when competitive electronic couplings are in balance. At longer Pt-Pt distances, the electronic couplings between the two halves of the binuclear system weaken, and exciton localization and recombination is favored on short time scales. For short Pt-Pt distances, electronic couplings between the states in the coherent superposition are stronger than the coupling with other excitonic states, leading to long-lived coherence.

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

Solute–solute correlations responsible for the prepeak in structure factors of undercooled Al-rich liquids: A molecular dynamics study

In this study, we have performed molecular dynamics simulations on a typical Al-based alloy Al 90 Sm 10 . The short-range and medium-range correlations of the system are reliably produced by ab initio calculations, whereas the long-range correlations are obtained with the assistance of a semi-empirical potential well-fitted to ab initio data. Our calculations show that a prepeak in the structure factor of this system emerges well above the melting temperature, and the intensity of the prepeak increases with increasing undercooling of the liquid. These results are in agreement with x-ray diffraction experiments. The interplay between the short-range order of the system originating from the large affinity between Al and Sm atoms, and the intrinsic repulsion between Sm atoms gives rise to a stronger correlation in the second peak than the first peak in the Sm–Sm partial pair correlation function (PPCF), which in turn produces the prepeak in the structure factor.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution

Hydrodynamics, oxygen transfer, and oxygen uptake in bubble-column and airlift bioreactors were studied using multiphase Euler-Euler computational fluid dynamics (CFD) simulations. Interphase mass transfer of oxygen was modeled with a liquid mass-transfer coefficient and sustained driving force provided by a phenomenological model for microbial oxygen uptake rate (OUR). Laboratory-scale reactor simulations of a bubble-column (0.15 m diameter and 0.75 m initial liquid height) showed reasonable agreement with gas holdup and mass transfer experiments reported in the literature. Commercial-scale bubble-column and airlift reactor simulations (5 m diameter and 25 m initial liquid height) were simulated using this validated model. Similar trends for the variation of overall gas holdup and oxygen concentrations with superficial velocity were observed for both of the commercial-scale reactors. However, the simulation results indicate differences in hydrodynamics, such as better recirculation in the airlift reactors that enables favorable distribution of oxygen in the reactor.

09 BIOMASS FUELS↗
DOE OSTI2018

Molecular dynamics studies of fundamental bulk properties of palladium hydrides for hydrogen storage

Solid-state hydrogen storage materials undergo complex phase transformations whose behavior are collectively determined by thermodynamic (e.g., Gibbs free energy), mechanical (e.g., lattice and elastic constants), and mass transport (e.g., diffusivity) properties. These properties depend on the reaction conditions and evolve continuously during (de)hydrogenation. Thus, they are difficult to measure in experiments. Because of this, past progress to improve solid-state hydrogen storage materials has been prolonged. Using PdH x as a representative example for interstitial metal hydride, we have recently applied molecular dynamics simulations to quantify hydrogen diffusion in the entire reaction space of temperature and composition. In this paper, we have further applied molecular dynamics simulations to obtain well-converged expressions for lattice constants, Gibbs free energies, and elastic constants of PdH x at various stages of the reaction. Finally, our studies confirm significant dependence of elastic constants on temperature and composition. Specifically, a new dynamic effect of hydrogen diffusion on elastic constants is discovered and discussed.

08 HYDROGEN↗
NASA NTRS1992

A molecular dynamics study of freezing in a confined geometry

The dynamics of freezing of a Lennard-Jones liquid in narrow channels bounded by molecular walls is studied by computer simulation. The time development of ordering is quantified and a novel freezing mechanism is observed. The liquid forms layers and subsequent in-plane ordering within a layer is accompanied by a sharpening of the layer in the transverse direction. The effects of channel size, the methods of quench, the liquid-wall interaction and the roughness of walls on the freezing mechanism are elucidated. Comparison with recent experiments on freezing in confined geometries is presented.

Ma, Wen-Jong↗
DOE OSTI2022

Simulations of graphite boronization: A molecular dynamics study of amorphization resulting from bombardment

The molecular dynamics code LAMMPS was used to simulate the bombardment of a graphite structure by atomic boron with impact energies ranging from 50–250 eV. The transient structural evolution, penetration depth, and amorphous layer thickness were analyzed. Simulations show that larger impact energies lead to a greater volume of amorphization and penetration of boron, but that the growth rate of the amorphous layer decreases with increasing fluence. Furthermore, the change in surface chemistry of the amorphized structures was studied using the ReaxFF formalism, which found that the amorphization process introduces dangling bonds thus increasing reactivity in the amorphous region.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1981

Cluster adsorption on amorphous and crystalline surfaces - A molecular dynamics study of model Pt on Cu and model Pd on Pt

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Garofalini, S. H.↗
DOE OSTI2017

Molecular dynamics studies of electron-ion temperature equilibration in hydrogen plasmas within the coupled-mode regime

Here, we use classical molecular dynamics (MD) to study electron-ion temperature equilibration in two-component plasmas in regimes for which the presence of coupled collective modes has been predicted to substantively reduce the equilibration rate. Guided by previous kinetic theory work, we examine hydrogen plasmas at a density of n = 10 26 cm –3 , T i = 10 5 K, and 10 7 K < Te < 10 9 K. The nonequilibrium classical MD simulations are performed with interparticle interactions modeled by quantum statistical potentials (QSPs). Our MD results indicate (i) a large effect from time-varying potential energy, which we quantify by appealing to an adiabatic two-temperature equation of state, and (ii) a notable deviation in the energy equilibration rate when compared to calculations from classical Lenard-Balescu theory including the QSPs. In particular, it is shown that the energy equilibration rates from MD are more similar to those of the theory when coupled modes are neglected. We suggest possible reasons for this surprising result and propose directions of further research along these lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Molecular Dynamics Study of the Proposed Proton Transport Pathways in [FeFe]-Hydrogenase

Possible proton channels in Clostridium pasteurianum [FeFe]-hydrogenase were investigated with molecular dynamics simulations. This study was undertaken to discern proposed channels, compare their properties, evaluate the functional channel, and to provide insight into the features of an active proton channel. Our simulations suggest that protons are not transported through water wires. Instead, a five-residue motif (E282, S319, E279, HOH, C299) was found to be the likely channel, consistent with experimental observations. This channel connects the surface of the enzyme and the di-thiomethylamine bridge of the catalytic H-cluster, permitting the transport of protons. The channel was found to have a persistent hydrogen bonded core (residues E279 to S319), with less persistent hydrogen bonds at the ends of the channel. The hydrogen bond occupancy in this network was found to be sensitive to the protonation state of the residues in the channel, with different protonation states enhancing or stabilizing hydrogen bonding in different regions of the network. Single site mutations to non-hydrogen bonding residues break the hydrogen bonding network at that residue, consistent with experimental observations showing catalyst inactivation. In many cases, these mutations alter the hydrogen bonding in other regions of the channel which may be equally important in catalytic failure. A correlation between the protein dynamics near the proton channel and the redox partner binding regions was also found as a function of protonation state. The likely mechanism of proton movement in [FeFe]-hydrogenases is discussed based on the structural analysis presented here. This work was funded by the DOE Office of Science Early Career Research Program through the Office of Basic Energy Sciences. Computational resources were provided at W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at Pacific Northwest National Laboratory, and a portion of the research was performed using PNNL Institutional Computing at Pacific Northwest National Laboratory. Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy.

biomimetics↗
DOE OSTI2015

Molecular-dynamics study of amorphous SiO{sub 2} relaxation

Using Molecular-Dynamics simulation we observed the generation of amorphous SiO{sub 2} target from a randomly distributed Si and O atoms. We applied a sequence of annealing of the target with various temperature and quenching to room temperature. The relaxation time required by the system to form SiO{sub 4} tetrahedral mesh after a relatively long simulation time, is studied. The final amorphous target was analyzed using the radial distribution function method, which can be compared with the available theoretical and experimental data. We found that up to 70% of the target atoms form the tetrahedral SiO{sub 4} molecules. The number of formed tetrahedral increases following the growth function and the rate of SiO{sub 4} formation follows Arrhenius law, depends on the annealing temperature. The local structure of amorphous SiO{sub 2} after this treatment agrees well with those reported in some literatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Ab initio molecular dynamics study of SiO2 lithiation

Li-ion batteries (LIBs) are sought to meet the demand for high energy storage applications. Due to its high specific charge capacity and low discharge potential, SiO 2 is a promising candidate material for LIBs anodes. To design high performance anodes and coating materials using SiO 2 , understanding the structural transformation and Li-ion kinetics in different structural forms of SiO 2 is essential. Here, we performed ab initio molecular dynamics to study the lithiation mechanism for crystalline and amorphous SiO 2 and the effect of surface termination to elucidate the lithiation process of conversion oxides and contribute to the development of future LIBs.

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

Mixed Molecular and Dissociative Water Adsorption on Hydroxylated TiO 2 (110): An Infrared Spectroscopy and Ab Initio Molecular Dynamics Study

For this work, we have investigated the structure and dynamics of water (D 2 O) adsorbed on TiO 2 (110) for coverages between 0 and 1 monolayer with infrared reflection absorption spectroscopy (IRAS) and ab initio molecular dynamics (AIMD) simulations. For D 2 O coverages as low as 0.4 monolayers (ML) on a hydroxylated surface, IR spectra typical of hydrogen-bonded chains of water molecules are observed. However, for D 2 O coverages ≥ 0.3 ML, a sharp, high frequency peak is also observed in the p-polarized spectra that is red-shifted relative to the bridging hydroxyl peak. This new peak is not observed for water adsorbed on an oxidized surface. Based on the AIMD simulations and comparisons with previous IR spectra for TiO 2 nanoparticles, the new peak is assigned to terminal hydroxyl groups produced by dissociative adsorption of some of the water on TiO 2 (110). The simulations indicate that water dissociation is related to the presence of defect electrons in the system, but not due to direct interactions between adsorbed water and bridging hydroxyls.

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

Effects of sub-Tg annealing on Cu64.5Zr35.5 glasses: A molecular dynamics study

Creating metallic glasses by cooling liquid melts in molecular dynamics simulations faces a well-known challenge that the cooling rate is too fast compared with experiments. Taking the prototypical Cu 64.5Zr35.5 glasses as an example, we propose an efficient cooling strategy in which most of the computer time is spent on a prolonged isothermal process slightly below the glass-transition temperature, Tg . The glassy sample prepared in this way demonstrates significant energetic stability, slow dynamics, and well-developed short-range icosahedral order. By conventional uniform cooling, similar properties can only be obtained using a cooling rate more than 15 times slower.

36 MATERIALS SCIENCE↗