Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics simulation”

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 1,783 records · Page 99

Structural and Na-ion conduction characteristics of Na 3 PS x Se 4–x

The recent discovery of the isostructrual cubic Na 3 PS 4 and Na 3 PSe 4 as fast Na-ion conductors provided a general structural framework for the exploration of new sodium superionic conductors. In this work, we systematically investigated the structures and ionic conduction characteristics of a series of compounds with the general chemical formula of Na 3 PS x Se 4–x . Synthesis of Na 3 PS 4 under different conditions (e.g., temperature, reaction vessel, mass of the precursors) reveals the reactivity of the precursors with the reaction tubes, producing different polymorphs. X-ray diffraction studies on the solid solution phases Na 3 PS x Se 4–x identified a tetragonal-to-cubic phase transition with increasing Se concentration. This observation is consistent with the computed stability of the tetragonal and cubic polymorphs, where the energy difference between the two polymorphs becomes very close to zero in Se-rich compositions. Furthermore, ab initio molecular dynamic simulations suggest that the fast Na-ion conduction in Na 3 PS x Se 4–x may not be causally related with the symmetry or the composition of these phases. The formation of defects, instead, enables fast Na-ion conduction in this class of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nanoscale view of assisted ion transport across the liquid–liquid interface

During solvent extraction, amphiphilic extractants assist the transport of metal ions across the liquid–liquid interface between an aqueous ionic solution and an organic solvent. Investigations of the role of the interface in ion transport challenge our ability to probe fast molecular processes at liquid–liquid interfaces on nanometer-length scales. Recent development of a thermal switch for solvent extraction has addressed this challenge, which has led to the characterization by X-ray surface scattering of interfacial intermediate states in the extraction process. Here, we review and extend these earlier results. We find that trivalent rare earth ions, Y(III) and Er(III), combine with bis(hexadecyl) phosphoric acid (DHDP) extractants to form inverted bilayer structures at the interface; these appear to be condensed phases of small ion–extractant complexes. The stability of this unconventional interfacial structure is verified by molecular dynamics simulations. The ion–extractant complexes at the interface are an intermediate state in the extraction process, characterizing the moment at which ions have been transported across the aqueous–organic interface, but have not yet been dispersed in the organic phase. In contrast, divalent Sr(II) forms an ion–extractant complex with DHDP that leaves it exposed to the water phase; this result implies that a second process that transports Sr(II) across the interface has yet to be observed. Calculations demonstrate that the budding of reverse micelles formed from interfacial Sr(II) ion–extractant complexes could transport Sr(II) across the interface. Our results suggest a connection between the observed interfacial structures and the extraction mechanism, which ultimately affects the extraction selectivity and kinetics.

Liang, Zhu↗

In Situ and Ex Situ TEM Study of Lithiation Behaviours of Porous Silicon Nanostructures

Here in this work, we study the lithiation behaviours of both porous silicon (Si) nanoparticles and porous Si nanowires by in situ and ex situ transmission electron microscopy (TEM) and compare them with solid Si nanoparticles and nanowires. The in situ TEM observation reveals that the critical fracture diameter of porous Si particles reaches up to 1.52 μm, which is much larger than the previously reported 150 nm for crystalline Si nanoparticles and 870 nm for amorphous Si nanoparticles. After full lithiation, solid Si nanoparticles and nanowires transform to crystalline Li 15 Si 4 phase while porous Si nanoparticles and nanowires transform to amorphous Li x Si phase, which is due to the effect of domain size on the stability of Li 15 Si 4 as revealed by the first-principle molecular dynamic simulation. Ex situ TEM characterization is conducted to further investigate the structural evolution of porous and solid Si nanoparticles during the cycling process, which confirms that the porous Si nanoparticles exhibit better capability to suppress pore evolution than solid Si nanoparticles. The investigation of structural evolution and phase transition of porous Si nanoparticles and nanowires during the lithiation process reveal that they are more desirable as lithium-ion battery anode materials than solid Si nanoparticles and nanowires.

25 ENERGY STORAGE↗

Polymer Conformations and Diffusion Through a Monolayer of Confining Nanoparticles

We present coarse-grained molecular dynamics simulations to probe chain-scale polymer conformations and diffusion between confining nanoparticles (NPs). By constructing a monolayer of hexagonally-packed NPs in a polymer melt with athermal interactions, we observe the magnitude and length-scale over which homogeneously confining NPs impact the polymer behavior, which provides fundamental insights into more complex polymer nanocomposites. We show that polymer conformations are more perturbed under strong confinement (i.e. when the interparticle distance, ID, is less than twice the polymer radius of gyration, 2R g ) as compared to around an isolated NP, and the effect depends on the ratio of R NP /R g rather than either independently. In fact, these conformations can be quantitatively replicated by executing a simple random walk in a similarly confining environment. We also show that polymer diffusion is slowed by the presence of NPs and the slowing persists far beyond the length-scale over which polymer conformations are perturbed. Although the slowing is strongest ~R g from the NPs, the diffusion coefficient is slower even beyond ~5R g from the NPs. Furthermore, by analyzing the directional van Hove distributions, we show polymer diffusion away from the NP monolayer is bulk-like while diffusion through the monolayer is slower with increasing NP confinement. Furthermore, these molecular dynamics simulations provide fundamental insights into the temporal and spatial effect of confining, athermal NPs on chain-scale polymer conformations and diffusion.

36 MATERIALS SCIENCE↗

Forging Fast Ion Conducting Nanochannels with Swift Heavy Ions: The Correlated Role of Local Electronic and Atomic Structure

Atomically disordered oxides have attracted significant attention in recent years due to the possibility of enhanced ionic conductivity. However, the correlation between atomic disorder, corresponding electronic structure, and the resulting oxygen diffusivity is not well understood. The disordered variants of the ordered pyrochlore structure in gadolinium titanate (Gd 2 Ti 2 O 7 ) are seen as a particularly interesting prospect due to intrinsic presence of a vacant oxygen site in the unit atomic structure, which could provide a channel for fast oxygen conduction. As such, in the present work, we provide insights into the subangstrom scale on the disordering-induced variations in the local atomic environment and its effect on the electronic structure in high-energy ion irradiation-induced disordered nanochannels, which can be utilized as pathways for fast oxygen ion transport. With the help of an atomic plane-by-plane-resolved analyses, the work shows how the presence of various types of TiO$_x$ polyhedral that exist in the amorphous and disordered crystalline phase modify the electronic structures relative to the ordered pyrochlore phase in Gd 2 Ti 2 O 7 . The correlated molecular dynamics simulations on the disordered structures show a remarkable enhancement in oxygen diffusivity as compared with ordered pyrochlore lattice and make that a suitable candidate for applications requiring fast oxygen conduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Life-cycle implications of hydrogen fuel cell electric vehicle technology for medium- and heavy-duty trucks

This study provides a comprehensive and up-to-date life-cycle comparison of hydrogen fuel cell electric trucks (FCETs) and their conventional diesel counterparts in terms of energy use and air emissions, based on the ensemble of well-established methods, high-fidelity vehicle dynamic simulations, and real-world vehicle test data. For the centralized steam methane reforming (SMR) pathway, hydrogen FCETs reduce life-cycle or well-to-wheel (WTW) petroleum energy use by more than 98% compared to their diesel counterparts. The reduction in WTW air emissions for gaseous hydrogen (G.H2) FCETs ranges from 20-45% for greenhouse gases, 37-65% for VOC, 49-77% for CO, 62-83% for NOx, 19-43% for PM10, and 27-44% for PM2.5, depending on vehicle weight classes and truck types. With the current U.S. average electricity generation mix, FCETs tend to create more WTW SOx emissions than their diesel counterparts, mainly because of the upstream emissions related to electricity use for hydrogen compression/liquefaction. Compared to G.H2, liquid hydrogen (L.H2) FCETs generally provide smaller WTW emissions reductions. For both G.H2 and L.H2 pathways for FCETs, because of electricity consumption for compression and liquefaction, spatio-temporal variations of electricity generation can affect the WTW results. FCETs retain the WTW emission reduction benefits, even when considering aggressive diesel engine efficiency improvement.

air emissions↗

Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation

Tight oil formations such as those in the Bakken petroleum system are known to hold hundreds of billions of barrels of oil in place; however, the primary recovery factor for these plays is typically less than 10%. Tight oil formations, including the Bakken Formation, therefore, may be attractive candidates for enhanced oil recovery (EOR) using CO 2 . Multiphase fluid behavior and flow in fluid-rich shales can vary substantially depending on the size of pore throats, and properties such as fluid viscosity and density are much different in nanoscale pores than in macroscale pores. Thus it is critical to understand the nature and distribution of nano-, micro-, and macroscale pores and fracture networks. To address these issues, the Energy & Environmental Research Center (EERC) has been conducting a research program entitled “Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation.” The objectives of the project are 1) the use of advanced characterization methods to better understand and quantify the petrophysical and geomechanical factors that control CO 2 and oil mobility within tight oil formation samples, 2) the determination of CO 2 permeation and oil extraction rates in tight reservoir rocks and organic-rich shales of the Bakken, and 3) the integration of the laboratory-based CO 2 permeation and oil extraction data and the characterization data into geologic models and dynamic simulations to develop predictions of CO 2 storage resource and EOR in the Bakken tight oil formation. A combination of standard and advanced petrophysical characterization techniques were applied to characterize samples of Bakken Formation tight reservoir rock and shales from multiple wells. Techniques included advanced computer tomography (CT) imaging, scanning electron microscopy (SEM) techniques, whole-core and micro x-ray CT imaging, field emission (FE) SEM, and focused ion beam (FIB) SEM. Selected samples were also analyzed for geomechanical properties. X-ray CT imaging yielded information on the occurrence of fractures, bedding planes, fossils, and bioturbation in core, as well as data on bulk density and photoelectric factor logs, which were used to interpret porosity, organic content, and mineralogy. FESEM was used for characterization of nano- and microscale features, including nanoscale pore visualization and micropore and pore throat mineralogy. FIBSEM yielded micro- to nanoscale visualization of fracture networks, porosity and pore-size distribution, connected versus isolated porosity, and distribution of organics. Results from the characterization activities provide insight on nanoscale fracture properties, pore throat mineralogy and connectivity, rock matrix characteristics, mineralogy, and organic content. Laboratory experiments demonstrated that CO 2 can permeate the tight matrix of Bakken shale and nonshale reservoir samples and mobilize oil from those samples. Geologic models were created at scales ranging from the core plug to the reservoir, and dynamic simulations were conducted. The data from the characterization and laboratory-based activities were integrated into modeling research activities to determine the fundamental mechanisms controlling fluid transport in the Bakken, which support EOR scheme design and estimation of CO 2 storage potential in tight oil formations. Simulation results suggest a CO 2 storage resource estimate range of 169 million to 1.5 billion tonnes for the Bakken in North Dakota, possibly resulting in 1.8 billion to 16 billion barrels of incremental oil.

04 OIL SHALES AND TAR SANDS↗

A Redox-Active, Compact Molecule for Cross-Linking Amyloidogenic Peptides into Nontoxic, Off-Pathway Aggregates: In Vitro and In Vivo Efficacy and Molecular Mechanisms

Chemical reagents targeting and controlling amyloidogenic peptides have received much attention for helping identify their roles in the pathogenesis of protein-misfolding disorders. Herein, we report a novel strategy for redirecting amyloidogenic peptides into nontoxic, off-pathway aggregates, which utilizes redox properties of a small molecule (DMPD, N,N-dimethyl-p-phenylenediamine) to trigger covalent adduct formation with the peptide. In addition, for the first time, biochemical, biophysical, and molecular dynamics simulation studies have been performed to demonstrate a mechanistic understanding for such an interaction between a small molecule (DMPD) and amyloid-β (Aβ) and its subsequent anti-amyloidogenic activity, which, upon its transformation, generates ligand–peptide adducts via primary amine-dependent intramolecular cross-linking correlated with structural compaction. Furthermore, in vivo efficacy of DMPD toward amyloid pathology and cognitive impairment was evaluated employing 5xFAD mice of Alzheimer’s disease (AD). Such a small molecule (DMPD) is indicated to noticeably reduce the overall cerebral amyloid load of soluble Aβ forms and amyloid deposits as well as significantly improve cognitive defects in the AD mouse model. Overall, our in vitro and in vivo studies of DMPD toward Aβ with the first molecular-level mechanistic investigations present the feasibility of developing new, innovative approaches that employ redox-active compounds without the structural complexity as next-generation chemical tools for amyloid management.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of surface initial condition on aqueous corrosion of glass-A study with a model nuclear waste glass

Being a nonequilibrium material, the structure of glass varies with the sample history. Thus, the initial surface condition of a glass can vary with the preparation condition and have a large impact on its reactivity. This paper shows that the aqueous corrosion behavior of international simple glass (ISG) varies depending on the initial surface state. The ISG glass samples were prepared as-polished-only and polished-then-annealed and they were immersed in aqueous solution saturated with soluble SiO 2 at 30°C (modeling a mild condition) and at 90°C (modeling a severe condition). Molecular dynamics simulations were performed to obtain coordination numbers of network formers of ISG to assist oxygen speciation calculations. The surface structures of as-prepared and corroded ISG samples were analyzed using various imaging and spectroscopic techniques. Among these analyses, only the oxygen speciation with x-ray photoelectron spectroscopy showed discernable differences between two uncorroded surfaces with different preparation histories; all other methods could not differentiate the surface preparation history before aqueous corrosion. Such minor difference in chemical structures was found to have a profound impact on corrosion behaviors in the mild condition. In the harsh condition, the surface history dependence was not as drastic as the corrosion in the mild condition. The analysis results of the corroded surfaces suggested that the thickness and structure of the alteration layer formed on ISG in aqueous corrosion can vary with the initial surface state.

nuclear, defects, corrosion↗

A Minimal Membrane Metal Transport System: Dynamics and Energetics of mer Proteins

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg 2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg 2+ -bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg 2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free-energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg 2+ -bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse-grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Finally, our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Combined Experimental and Theoretical Study on the Activity and Selectivity of the Electrocatalytic Hydrogenation of Aldehydes

A detailed mechanistic study of the electrochemical hydrogenation of aldehydes is presented toward the goal of identifying how organic molecules in solution behave at the interface with charged surfaces and what is the best manner to convert them. Specifically, this study focuses on designing an electrocatalytic route for ambient temperature post-pyrolysis treatment of bio-oil. Aldehyde reductions are needed to convert biomass into fuels or chemicals. A combined experimental and computational approach is taken toward catalyst design to provide testable hypotheses regarding catalyst composition, activity, and selectivity. Electrochemical hydrogenation mechanisms for benzaldehyde and pentanal reduction are found to proceed by a coupled proton-electron transfer process. Initial results show that Au, Ag, Cu, and C catalysts exhibit the highest conversion to alcohol products. These catalysts are suitable because they show high cathodic onset potentials for H2 formation and low cathodic onset potentials for organic reduction. Conversion of aromatic aldehydes is found to be appreciably higher than aliphatic ones. Classical molecular dynamics simulations of solvent and substrate mixtures in an electrolytic cell were performed to assess how species concentrations vary at the solid/liquid interface and in the bulk as a function of applied voltage. Results show that an increase in voltage in the electrolytic cell decreases organic and increases water mole fractions at the solid/liquid interface. In this current study, charged cathodic surfaces result in carbonyl orientations at the surface that do not favor electron transfer. Repulsion of organic substrates to the bulk must be compensated by strong adhesion to the electrode surface. Implications on catalyst choice and process design are discussed.

Pyrolysis Oil upgrading, Atomistic Simulations, El↗

Cellulose–hemicellulose interactions at elevated temperatures increase cellulose recalcitrance to biological conversion

It has been previously shown that cellulose-lignin droplets’ strong interactions, resulting from lignin coalescence and redisposition on cellulose surface during thermochemical pretreatments, increase cellulose recalcitrance to biological conversion, especially at commercially viable low enzyme loadings. However, information on the impact of cellulose–hemicellulose interactions on cellulose recalcitrance following relevant pretreatment conditions are scarce. Here, to investigate the effects of plausible hemicellulose precipitation and re-association with cellulose on cellulose conversion, different pretreatments were applied to pure Avicel® PH101 cellulose alone and Avicel mixed with model hemicellulose compounds followed by enzymatic hydrolysis of resulting solids at both low and high enzyme loadings. Solids produced by pretreatment of Avicel mixed with hemicelluloses (AMH) were found to contain about 2 to 14.6% of exogenous, precipitated hemicelluloses and showed a remarkably much lower digestibility (up to 60%) than their respective controls. However, the exogenous hemicellulosic residues that associated with Avicel following high temperature pretreatments resulted in greater losses in cellulose conversion than those formed at low temperatures, suggesting that temperature plays a strong role in the strength of cellulose–hemicellulose association. Molecular dynamics simulations of hemicellulosic xylan and cellulose were found to further support this temperature effect as the xylan–cellulose interactions were found to substantially increase at elevated temperatures. Furthermore, exogenous, precipitated hemicelluloses in pretreated AMH solids resulted in a larger drop in cellulose conversion than the delignified lignocellulosic biomass containing comparably much higher natural hemicellulose amounts. In conclusion, increased cellulase loadings or supplementation of cellulase with xylanases enhanced cellulose conversion for most pretreated AMH solids; however, this approach was less effective for solids containing mannan polysaccharides, suggesting stronger association of cellulose with (hetero) mannans or lack of enzymes in the mixture required to hydrolyze such polysaccharides.

59 BASIC BIOLOGICAL SCIENCES↗

Size- and temperature-dependent magnetization of iron nanoclusters

In this paper, the magnetic behavior of bcc iron nanoclusters, with diameters between 2 and 8 nm, is investigated by means of spin dynamics simulations coupled to molecular dynamics, using a distance-dependent exchange interaction. Finite-size effects in the total magnetization as well as the influence of the free surface and the surface/core proportion of the nanoclusters are analyzed in detail for a wide temperature range, going beyond the cluster and bulk Curie temperatures. Comparison is made with experimental data and with theoretical models based on the mean-field Ising model adapted to small clusters, and taking into account the influence of low coordinated spins at free surfaces. Our results for the temperature dependence of the average magnetization per atom M (T), including the thermalization of the transnational lattice degrees of freedom, are in very good agreement with available experimental measurements on small Fe nanoclusters. In contrast, significant discrepancies with experiment are observed if the translational degrees of freedom are artificially frozen. The finite-size effects on M (T) are found to be particularly important near the cluster Curie temperature. Simulated magnetization above the Curie temperature scales with cluster size as predicted by models assuming short-range magnetic ordering. Analytical approximations to the magnetization as a function of temperature and size are proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Differences in the dynamics of the tandem-SH2 modules of the Syk and ZAP-70 tyrosine kinases

The catalytic activity of Syk-family tyrosine kinases is regulated by a tandem Src homology 2 module (tSH2 module). In the autoinhibited state, this module adopts a conformation that stabilizes an inactive conformation of the kinase domain. The binding of the tSH2 module to phosphorylated immunoreceptor tyrosine-based activation motifs necessitates a conformational change, thereby relieving kinase inhibition and promoting activation. We determined the crystal structure of the isolated tSH2 module of Syk and find, in contrast to ZAP-70, that its conformation more closely resembles that of the peptide-bound state, rather than the autoinhibited state. Hydrogen–deuterium exchange by mass spectrometry, as well as molecular dynamics simulations, reveal that the dynamics of the tSH2 modules of Syk and ZAP-70 differ, with most of these differences occurring in the C-terminal SH2 domain. Our data suggest that the conformational landscapes of the tSH2 modules in Syk and ZAP-70 have been tuned differently, such that the autoinhibited conformation of the Syk tSH2 module is less stable. This feature of Syk likely contributes to its ability to more readily escape autoinhibition when compared to ZAP-70, consistent with tighter control of downstream signaling pathways in T cells.

59 BASIC BIOLOGICAL SCIENCES↗