Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ATOMIC MODELS”

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,099 records · Page 61

SciDAC - Center for Simulation of Wave Interactions with MHD -- General Atomics Support of ORNL Collaboration

The Center for Simulation of Wave Interactions with Magnetohydrodynamics (SWIM) project is dedicated to conduct research on integrated multi-physics simulations. The Integrated Plasma Simulator (IPS) is a framework that was created by the SWIM team. It provides an integration infrastructure for loosely coupled component-based simulations by facilitating services for code execution coordination, computational resource management, data management, and inter-component communication. The IPS framework features improving resource utilization, implementing application-level fault tolerance, and support of the concurrent multi-tasking execution model. The General Atomics (GA) team worked closely with other team members on this contract, and conducted research in the areas of computational code monitoring, meta-data management, interactive visualization, and user interfaces. The original website to monitor SWIM activity was developed in the beginning of the project. Due to the amended requirements, the software was redesigned and a revision of the website was deployed into production in April of 2010. Throughout the duration of this project, the SWIM Monitoring Portal (http://swim.gat.com:8080/) has been a critical production tool for supporting the project's physics goals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Revisiting pressure-induced phase transition in silicon clathrates using Ge substitution

Ba 8 Si 39 Ge 7 and Ba 8 Si 29 Ge 17 have been studied at high pressure using x-ray diffraction and x-ray absorption spectroscopy (XAS) at the Ge K edge. In Ba 8 Si 39 Ge 7 , a transition is observed similar to the one in Ba 8 Si 46 , apparently isostructural. However, the XAS data analysis shows that the transformation is related to the off-centering of the Ba atoms. A theoretical model based on a Landau potential suggests that this transition is second order, with a symmetry-breaking mechanism related to the Ba displacement probably initiated by the vacancy creation or local distortion predicted theoretically. Lastly, this analysis gives a coherent picture of the phase transition mechanism. In the case of Ba 8 Si 29 Ge 17 , such phase transition is not observed as the Ba atoms appear already off-center at ambient pressure.

36 MATERIALS SCIENCE↗

The Bohr–Kalckar correspondence principle and a new construction of partitions in number theory

The author attempts to change and supplement the standard scheme of partitions of integers in number theory to make it completely concur with the Bohr–Kalckar correspondence principle. In order to make the analogy between the the atomic nucleus and the theory of partitions of natural numbers more complete, to the notion of defect of mass author assigns the “defect” ({a})-bar = [a + 1] − a of any real number a (i.e., the fractional value that must be added to a in order to obtain the nearest larger integer). This allows to carry over the Einstein relation between mass and energy to a relation between the whole numbers M and N, where N is the number of summands in the partition ofM into positive summands, as well as to define the forbidding factor for the number M, and apply this to the Bohr–Kalckar model of heavy atomic nuclei and to the calculation of the maximal number of nucleons in the nucleus.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolutionary problems of Cepheids and other giants investigated with new radiative opacities

Comparison of evolutionary tracks, pulsation constants, and linearized pulsational-stability coefficients for stellar models applicable to the problems of classical Cepheids, whose structures were calculated using the Cox-Stewart (1965) opacities and a recently computed set of opacities. The latter are based on the hot 'Thomas-Fermi' statistical model of the atom for all elements heavier than hydrogen and helium; they contain larger helium and metals contributions, but a smaller hydrogen contribution than the former ones for the same chemical composition. The difference in metals contribution affects mainly the location and shape of the evolutionary tracks on the H-R diagram, while the difference in hydrogen and helium contributions has its greatest effect on the pulsational properties of the Cepheid models. From the comparison of evolutionary tracks it is concluded that: (1) the theoretical M/L relation for evolved giants is changed very little by using the second set of opacities; (2) Q-values for the fundamental mode of radial pulsation in Cepheid envelope models increase if the second set is used, but the classical mass discrepancy remains; and (3) the second set leads to pulsational-instability.

Carson, T. R.↗

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗

An evaluation of equilibrium conditions and temperature-dependent speciation in a laser-produced air plasma

A laser-induced air breakdown plasma is a dynamical system and its fundamental parameters change significantly during its lifetime. The temporal evolution of the spectral features from an air plasma shows emission from bremsstrahlung followed by discrete electronic transitions from ions, atoms, and rotational and vibrations bands from molecules. The molecular band emission from air plasma typically appears at times ≳5 μs and persists for 100s of μs. Here, an evaluation is made on the temperature evolution and speciation in a laser-produced air plasma. The air plasmas are generated by focusing 1064 nm, 6 ns pulses from an Nd:YAG laser. High-resolution emission spectra of various molecules such as N 2 , N 2 + , CN, OH, NH, NO are acquired for fits to spectral models. Fitting of atomic and molecular emission features permits tracking of the air spark temperature evolution from 1 to 200 µs. Though the excitation and molecular temperatures show good overlap at times ≲10 μs, the molecular temperatures obtained from different species show a discontinuity at times ~ 30 μs and this is related to shock collapse and subsequent changes in hydrodynamics and chemistry of the plume. The fitting of multiple species in broadband spectra has permitted calculation of the relative concentrations of various molecules as a function of temperature that provides insight into the air spark speciation. Finally, the measured mole fractions from laser-induced air plasma are within a factor of 2 - 6 of the values estimated by speciation model under the assumption of thermal and chemical equilibrium across the temperature ranges studied and probable reasons for the observed deviation are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic properties of CrX 3 (X = Cl, Br, I) monolayers in excited states

The presence of intrinsic ferromagnetism in a new class of magnetic semiconductors CrX 3 (X = Cl, Br, I) in a reduced dimension is associated with magnetic anisotropy energy (MAE). MAE is relevant for the magnetic order of CrX 3 . The control of magnetism in such 2D materials through external stimuli is gaining importance in view of potential technological applications. Here, we use all-electron density functional theory to investigate the magnetism of monolayers of chromium halides in excited states using a cluster model. We study the magnetization reversal barrier for two types of electronic singly excited states: (a) spin-preserving transitions, and (b) spin-flip transitions. While all the three halides show enhancement of MAE in excited states, the MAEs are significant for CrI 3 due to the large spin–orbit coupling in the iodines. The changes in the MAE in the excited state are influenced by the local charge distribution surrounding the iodine atoms. Using a model system, we also show that the MAE in chromium tri-iodides is strongly influenced by the charge transfer to the iodines in an excited state. Furthermore, controlling the charge transfer to the iodines can potentially be used to control the magnetization reversal barrier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exposing hidden alternative backbone conformations in X-ray crystallography using qFit

Proteins must move between different conformations of their native ensemble to perform their functions. Crystal structures obtained from high-resolution X-ray diffraction data reflect this heterogeneity as a spatial and temporal conformational average. Although movement between natively populated alternative conformations can be critical for characterizing molecular mechanisms, it is challenging to identify these conformations within electron density maps. Alternative side chain conformations are generally well separated into distinct rotameric conformations, but alternative backbone conformations can overlap at several atomic positions. Our model building program qFit uses mixed integer quadratic programming (MIQP) to evaluate an extremely large number of combinations of sidechain conformers and backbone fragments to locally explain the electron density. Here, we describe two major modeling enhancements to qFit: peptide flips and alternative glycine conformations. We find that peptide flips fall into four stereotypical clusters and are enriched in glycine residues at the n+1 position. The potential for insights uncovered by new peptide flips and glycine conformations is exemplified by HIV protease, where different inhibitors are associated with peptide flips in the “flap” regions adjacent to the inhibitor binding site. Our results paint a picture of peptide flips as conformational switches, often enabled by glycine flexibility, that result in dramatic local rearrangements. Our results furthermore demonstrate the power of large-scale computational analysis to provide new insights into conformational heterogeneity. Furthermore, improved modeling of backbone heterogeneity with high-resolution X-ray data will connect dynamics to the structure-function relationship and help drive new design strategies for inhibitors of biomedically important systems.

59 BASIC BIOLOGICAL SCIENCES↗

Advanced EXAFS analysis techniques applied to the L -edges of the lanthanide oxides

The unique properties of the lanthanide (Ln) elements make them critical components of modern technologies, such as lasers, anti-corrosive films and catalysts. Thus, there is significant interest in establishing structure–property relationships for Ln-containing materials to advance these technologies. Extended X-ray absorption fine structure (EXAFS) is an excellent technique for this task considering its ability to determine the average local structure around the Ln atoms for both crystalline and amorphous materials. However, the limited availability of EXAFS reference spectra of the Ln oxides and challenges in the EXAFS analysis have hindered the application of this technique to these elements. The challenges include the limited k-range available for the analysis due to the superposition of L-edges on the EXAFS, multielectron excitations (MEEs) creating erroneous peaks in the EXAFS and the presence of inequivalent absorption sites. Herein, we removed MEEs to model the local atomic environment more accurately for light Ln oxides. Further, we investigated the use of cubic and non-cubic lattice expansion to minimize the fitting parameters needed and connect the fitting parameters to physically meaningful crystal parameters. The cubic expansion reduced the number of fitting parameters but resulted in a statistically worse fit. The non-cubic expansion resulted in a similar quality fit and showed non-isotropic expansion in the crystal lattice of Nd 2 O 3 . In total, the EXAFS spectra and the fits for the entire set of Ln oxides (excluding promethium) are included. The knowledge developed here can assist in the structural determination of a wide variety of Ln compounds and can further studies on their structure–property relationships.

36 MATERIALS SCIENCE↗

Energetic Neutral Atom Emissions From Venus: VEX Observations and Theoretical Modeling

Venus has almost no intrinsic magnetic field to shield itself from its surrounding environment. The solar wind thus directly interacts with the planetary ionosphere and atmosphere. One of the by-products of this close encounter is the production of energetic neutral atom (ENA) emissions. Theoretical studies have shown that significant amount of ENAs are emanated from the planet. The launch of the Venus Express (VEX) in 2005 provided the first light ever of the Venus ENA emissions. The observed ENA flux level and structure are in pretty good agreement with the theoretical studies. In this paper, we present VEX ENA data and the comparison with numerical simulations. We seek to understand the solar wind interaction with the planet and the impacts on its atmospheres.

Fok, M.-C.↗

Interrelated structures of the transport shock and collisional relaxation layer in a multitemperature, multilevel ionized gas

The gas dynamic structures of the transport shock and the downstream collisional relaxation layer are evaluated for partially ionized monatomic gases. Elastic and inelastic collisional nonequilibrium effects are taken into consideration. Three electronic levels are accounted for in the microscopic model of the atom. Nonequilibrium processes with respect to population of levels and species plus temperature are considered. By using an asymptotic technique the shock morphology is found on a continuum flow basis. The asymptotic procedure gives two distinct layers in which the nonequilibrium effects to be considered are different. A transport shock appears as the inner solution to an outer collisional relaxation layer in which the gas reaches local equilibrium. A family of numerical examples is displayed for different flow regimes. Argon and helium models are used in these examples.

Vinolo, A. R.↗

Machine Learning for Quantum Mechanical Properties of Atoms in Molecules

Here we introduce machine learning models of quantum mechanical observables of atoms in molecules. Instant out-of-sample predictions for proton and carbon nuclear chemical shifts, atomic core level excitations, and forces on atoms reach accuracies on par with density functional theory reference. Locality is exploited within nonlinear regression via local atom-centered coordinate systems. The approach is validated on a diverse set of 9 k small organic molecules. Linear scaling of computational cost in system size is demonstrated for saturated polymers with up to submesoscale lengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of $\mathrm{AMOEBA}$ Polarizable Force Field for Rare-Earth La 3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. In this report we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La 3+ ) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La 3+ –ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La 3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion–water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion–protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Acceleration of isolated atoms by radiation pressure

A model is considered of a point source of radiation with a power-law spectrum in which atoms are both accelerated and photoionized by the radiation pressure, and the resulting spectrum of escape velocities is calculated. It is found that the velocity spectrum will depend only on the power-law index and will be independent of the radiation intensity. It is shown that relativistic velocities are not possible for power-law indices typical of astrophysical objects.

Gordon, R. H.↗

Observational constraints on the multiphase nature of outflows using large spectroscopic surveys at z ∼ 0

ABSTRACT Mass outflow rates and loading factors are typically used to infer the quenching potential of galactic-scale outflows. However, these generally rely on observations of a single gas phase that can severely underestimate the total ejected gas mass. To address this, we use observations of high mass (≥1010 M⊙), normal star-forming galaxies at z ∼ 0 from the MaNGA, xCOLD GASS, xGASS, and ALFALFA surveys and a stacking of Na d, Hα, CO(1–0), and H i 21 cm tracers with the aim of placing constraints on an average, total mass outflow rate, and loading factor. We find detections of outflows in both neutral and ionized gas tracers, with no detections in stacks of molecular or atomic gas emission. Modelling of the outflow components reveals velocities of |vNa d| = 131 km s−1 and |vHα| = 439 km s−1 and outflow rates of $\dot{M}_{\rm {Na\,\small{D}}}$ = 7.55 M⊙ yr−1 and $\dot{M}_{\text{H}\alpha }$ = 0.10 M⊙ yr−1 for neutral and ionized gas, respectively. Assuming a molecular/atomic outflow velocity of 200 km s−1, we derive upper limits of $\dot{M}_{\text{CO}}\lt 19.43$ M⊙ yr−1 and $\dot{M}_{\rm {H\,\small {I}}}\lt $ 26.72 M⊙ yr−1 for the molecular and atomic gas, respectively. Combining the detections and upper limits, we find average total outflow rates of $\dot{M}_{\text{tot}}\lesssim$27 M⊙ yr−1 and a loading factor of ηtot ≲ 6.39, with molecular gas likely contributing ≲72 per cent of the total mass outflow rate, and neutral and ionized gas contributing ∼28 and <1 per cent, respectively. Our results suggest that, to first order, a degree of quenching via ejective feedback could occur in normal galaxies when considering all gas phases, even in the absence of an active galactic nucleus.

Roberts-Borsani, G. W.↗

Li-atoms-induced structure changes of Guinier–Preston–Bagaryatsky zones in AlCuLiMg alloys

Guinier–Preston–Bagaryatsky (GPB) zones are the well-known strengthening precipitates of AlCuMg alloys formed upon thermal ageing. Here we report that when formed in AlCuLiMg alloys the GPB zones can change significantly in morphology and structure. It is shown that though they do still consist of Al, Cu and Mg elements fundamentally, the GPB zones in AlCuLiMg alloys have a rather different structure due to a featured Li-segregation at their interfaces with the matrix and possible Li-replacement of partial Mg atoms in the structure. As such the Li-containing GPB zones often develop from one-dimensional to quasi-two-dimensional precipitates. - Highlights: • We observe Guinier–Preston–Bagaryatsky zone variants in AlCuLiMg alloys. • We obtain atomic-resolution images of the precipitates and model their structures. • Li-atoms play a key role in modifying the structure of these precipitate variants.

36 MATERIALS SCIENCE↗

Atomic Layer Deposition for Improved Biomass Conversion Catalysts

Heterogeneous catalysts are a key enabler of the transition towards a sustainable, bio-based economy for fuels and chemicals. However, the harsh conditions in many biomass conversion processes lead to nanoparticle sintering, support collapse, and metal leaching in conventional PGM catalysts. Next-generation catalysts must be developed to address these stability challenges. This presentation will discuss ongoing work between the Catalytic Carbon Transformation and Scale-Up Center at NREL and various industrial partners to develop scalable and cost-effective atomic layer deposition (ALD) coatings for improving the performance of biomass conversion catalysts. The substantial focus will be devoted to the benefits of TiO2 ALD coatings on supported Pd hydrogenation catalysts. Ten cycles of TiO2 ALD were found to dramatically improve the activity of a conventional Pd/Al2O3 catalyst towards aromatic hydrogenation, despite partial coverage of the Pd sites by the ALD layer. Subsequent advanced characterization and atomic-scale computational modeling revealed that the ALD coating weakened the adsorption strength of hydrogenation surface intermediates, leading to higher reaction rates. Reaction testing after exposure to sulfur impurities, high temperature oxidation, and hydrothermal treatment demonstrated the improved stability of the ALD-coated catalyst. Additionally, the ALD synthesis process was found to be scalable over two orders of magnitude with minimal deviation in synthesized catalyst properties. These results were contextualized with cost models of industrial ALD coating and aromatic hydrogenation processes to further refine the value proposition of ALD coatings. Given the demonstrated improvements in hydrogenation, TiO2 ALD coatings have also been applied to supported Pt catalysts for use in other biomass conversion reactions, such as hydrodeoxygenation for the production of sustainable aviation fuel (SAF) and the oxidation of glucose to bio-derived gluconic acid.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Accelerating Ab Initio Simulation via Nested Monte Carlo and Machine Learned Reference Potentials

As a corollary of the rapid advances in computing, ab initio simulation is playing an increasingly important role in modeling materials at the atomic scale. Two strategies are possible, ab initio Monte Carlo (AIMC) and molecular dynamics (AIMD) simulation. The former benefits from exact sampling from the correct thermodynamic distribution, while the latter is typically more efficient with its collective all-atom coordinate updates. In this study, using a relatively simple test model comprised of helium and argon, we show that AIMC can be brought up to, and even above, the performance levels of AIMD via a hybrid nested sampling/machine learning (ML) strategy. Here, ML provides an accurate classical reference potential (up to three-body explicit interactions) that can pilot long collective Monte Carlo moves that are accepted or rejected in toto à la nested Monte Carlo (NMC); this is in contrast to the single move nature of a naive implementation. Our proposed method only requires a small up front expense from evaluating the ab initio energies and forces of (100) random configurations for training. Importantly, our method does not totally rely on the trained, assuredly imperfect, interaction. We show that high performance and exact sampling at the desired level of theory can be realized even when the trained interaction has appreciable differences from the ab initio potential. Remarkably, at the highest levels of performance realized via our approach, a pair of statistically uncorrelated atomic configurations can be generated with (1) ab initio calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗