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 685 records · Page 38

Measurement Error in Atomic-Scale Scanning Transmission Electron Microscopy—Energy-Dispersive X-Ray Spectroscopy (STEM-EDS) Mapping of a Model Oxide Interface

Abstract With the development of affordable aberration correctors, analytical scanning transmission electron microscopy (STEM) studies of complex interfaces can now be conducted at high spatial resolution at laboratories worldwide. Energy-dispersive X-ray spectroscopy (EDS) in particular has grown in popularity, as it enables elemental mapping over a wide range of ionization energies. However, the interpretation of atomically resolved data is greatly complicated by beam–sample interactions that are often overlooked by novice users. Here we describe the practical factors—namely, sample thickness and the choice of ionization edge—that affect the quantification of a model perovskite oxide interface. Our measurements of the same sample, in regions of different thickness, indicate that interface profiles can vary by as much as 2–5 unit cells, depending on the spectral feature. This finding is supported by multislice simulations, which reveal that on-axis maps of even perfectly abrupt interfaces exhibit significant delocalization. Quantification of thicker samples is further complicated by channeling to heavier sites across the interface, as well as an increased signal background. We show that extreme care must be taken to prepare samples to minimize channeling effects and argue that it may not be possible to extract atomically resolved information from many chemical maps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A New Approach to Calculating the Potential Energy of Interaction between Two Atoms

A Fourier component of the potential energy of interaction between two atoms has been presented as a polynomial of the biquadratic atomic form factor. A numerical calculation has been performed in the screened Coulomb potential approximation. It is demonstrated that the account for the Pauli principle leads to the occurrence of a potential barrier and an additional region of attraction of two atoms. This model is shown to agree qualitatively with the results of the density-functional-theory calculation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Atomization data requirements for rocket combustor modeling

The complex computer codes, which model liquid rocket combustors, require information about the distribution and atomization of these liquid reactants. The available information is, in general, of questionable validity and applicability. Authors and users of combustion codes are often unaware of, or underestimate the importance of, these deficiencies in atomization data. These deficiencies and their importance are examined. Results of analyses performed with a state-of-the-art rocket combustion code are presented which demonstrate the important effects of such atomization information as initial droplet sizes and size distribution on vaporization rate and losses. Also, the questionable aspects and inapplicability of the available atomization data are discussed. One important and often neglected or misunderstood aspect of atomization data is the differences between spatial (concentration) and flux (often called temporal) droplet size distributions. These are described, and a computer model constructed to assess the difference between concentration and flux droplet size distributions is described and results presented. Experimental data are also given to demonstrate this difference. Finally, experimental results are presented that demonstrate the very great, and often neglected effect, of the local gas velocity field on atomization.

Ferrenberg, A. J.↗

Complex Density Wave Orders and Quantum Phase Transitions in a Model of Square-Lattice Rydberg Atom Arrays

In this work, we describe the zero-temperature phase diagram of a model of a two-dimensional square-lattice array of neutral atoms, excited into Rydberg states and interacting via strong van der Waals interactions. Using the density-matrix renormalization group algorithm, we map out the phase diagram and obtain a rich variety of phases featuring complex density wave orderings, upon varying lattice spacing and laser detuning. While some of these phases result from the classical optimization of the van der Waals energy, we also find intrinsically quantum-ordered phases stabilized by quantum fluctuations. These phases are surrounded by novel quantum phase transitions, which we analyze by finite-size scaling numerics and Landau theories. Our work highlights Rydberg quantum simulators in higher dimensions as promising platforms to realize exotic many-body phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Local Structures of High-Entropy Alloys (HEAs) on Atomic Scales: An Overview

The high-entropy alloys (HEAs), containing several elements mixed in equimolar or near-equimolar ratios, have shown exceptional engineering properties. Local structures on atomic level are essential to understand the mechanical behaviors and related mechanisms. In this paper, the local structure and stress on the atomic level are reviewed by the pair-distribution function (PDF) of neutron-diffraction data, ab-initio-molecular-dynamics (AIMD) simulations, and atomic-probe microscopy (APT).

36 MATERIALS SCIENCE↗

Measuring the oscillator strength of intercombination lines of helium-like V ions in a laser-produced-plasma

In this work, we present results of measurements of the oscillator strength of intercombination lines of He-like Vanadium ions in high energy density (HED) laser-produced-plasmas and compare them with the simulations from commonly used codes and data from the NIST database. Whilst not yet sufficiently accurate to constrain different trusted atomic-physics models for the particular system studied, our results are in agreement with the available data within experimental error bars, yet differ from cruder approximations of the oscillator strength used in certain atomic-kinetics packages, suggesting that this general method could be further extended to be used as a measurement of the oscillator strength of additional atomic transitions under the extreme conditions that are achieved in HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toward Spatial Control of Reaction Selectivity on Photocatalysts Using Area-Selective Atomic Layer Deposition on the Model Dual Site Electrocatalyst Platform

Photocatalytic water splitting is a promising route to low-cost, green H 2 . However, this approach is currently limited in its solar-to-hydrogen conversion efficiency. One major source of efficiency loss is attributed to the high rates of undesired side and back reactions, which are exacerbated by the proximity of neighboring oxidation and reduction sites. Nanoscopic oxide coatings have previously been used to selectively block undesired reactants from reaching active sites; however, a coating encapsulating the entire photocatalyst particle limits activity as it cannot facilitate both half-reactions. In this work, area selective atomic layer deposition (AS-ALD) was used to selectively deposit semipermeable TiO 2 films onto model metallic cocatalysts for enhancing reaction selectivity while maintaining a high overall activity. Pt and Au were used as exemplary reduction and oxidation cocatalyst sites, respectively, where Au was deactivated toward ALD growth through self-assembled thiol monolayers while TiO 2 was coated onto Pt sites. Electroanalytical measurements of monometallic thin film electrodes showed that the TiO 2 -encapsulated Pt effectively suppressed undesired H 2 oxidation and Fe(II)/Fe(III) redox reactions while still permitting the desired hydrogen evolution reaction (HER). A planar model photocatalyst platform containing patterned interdigitated arrays of Au and Pt microelectrodes was further assessed using scanning electrochemical microscopy (SECM), demonstrating the successful use of AS-ALD to enable local reaction selectivity in a dual-reaction-site (photo)electrocatalytic system. Finally, interdigitated microelectrodes having independent potential control were used to show that selectively deposited TiO 2 coatings can suppress the rate of back reactions on neighboring active sites by an order of magnitude compared with uncoated control samples.

08 HYDROGEN↗

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

We have extended here a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017)], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the thermodynamic efficiency of a nickel-based multiferroic thermomagnetic generator: From bulk to atomic scale

A model is developed to correlate the effects of size on the thermodynamic efficiency for a nickel-based multiferroic thermomagnetic generator device. Three existing models are combined in order to estimate this correlation, they are (1) thermodynamic efficiency relations, (2) a model of ferromagnetic transition behavior, and (3) the bond-order length strength correlation. At the smallest size considered, a monolayer of nickel atoms shows a reduction in Curie temperature from its bulk value of T{sub c,Bulk}=630 K to T{sub c,ML}=240 K. This difference is analytically shown to affect the thermodynamic efficiency values when compared to bulk. Various nickel nanofilms are considered as a working body, such that the combined model predicts relative efficiency values that are comparable to the bulk scale, but operating closer to room-temperature when compared to bulk form. This result is unexpected since the absolute efficiency is shown to increase as a function of decreasing size, this discrepancy is explained as a consequence of Curie point suppression. The combined model is also applied to a hypothetical composite made of separated layers of nickel with distinct thicknesses. This composite material is predicted to spread the ferromagnetic transition across a much larger temperature range as compared to bulk nickel, such that this material may be better suited for different applications; for example, as a sensor or thermal switch. Moreover, this combined model is also shown to give a lower-bound estimate for thermodynamic efficiency, since the actual performance depends on material characterizations that have yet to be determined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating the magnetic and atomic interface configuration for a model Fe/CrN bilayer system

A bilayer of iron on chromium nitride (Fe/CrN) is an interesting system for exchange biasing and sensing applications as the Néel temperature of CrN is 280 K and the Curie temperature of Fe is 1043 K. In this paper, we study the crystal and magnetic structures of the Fe/CrN interface at the atomic level. High quality epitaxial Fe/CrN bilayers prepared by molecular beam epitaxy grow in 001 orientation on MgO(001) substrates with uniform layer thicknesses and sharp interfaces. Our data reveal the epitaxial correlation between Fe and CrN crystals and their magnetic structures at the interface. The magnetic anisotropy directions of Fe and CrN are found parallel to [110] MgO . We studied the electronic and magnetic properties of the interface by performing the first-principles total-energy calculations. Here, we present a model that combines the crystal and magnetic structures of the Fe/CrN bilayer and fully explains all results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fast and Accurate Predictions of Total Energy for Solid Solution Alloys with Graph Convolutional Neural Networks

We use graph convolutional neural networks (GCNNs) to produce fast and accurate predictions of the total energy of solid solution binary alloys. GCNNs allow us to abstract the lattice structure of a solid material as a graph, whereby atoms are modeled as nodes and metallic bonds as edges. This representation naturally incorporates information about the structure of the material, thereby eliminating the need for computationally expensive data pre-processing which would be required with standard neural network (NN) approaches. We train GCNNs on ab-initio density functional theory (DFT) for copper-gold (CuAu) and iron-platinum (FePt) data that has been generated by running the LSMS-3 code, which implements a locally self-consistent multiple scattering method, on OLCF supercomputers Titan and Summit. GCNN outperforms the ab-initio DFT simulation by orders of magnitude in terms of computational time to produce the estimate of the total energy for a given atomic configuration of the lattice structure. We compare the predictive performance of GCNN models against a standard NN such as dense feedforward multi-layer perceptron (MLP) by using the root-mean-squared errors to quantify the predictive quality of the deep learning (DL) models. We find that the attainable accuracy of GCNNs is at least an order of magnitude better than that of the MLP.

Lupo Pasini, Massimiliano↗

Measurement of atom resolvability in cryo-EM maps with Q-scores

Cryogenic electron microscopy (cryo-EM) maps are now at the point where resolvability of individual atoms can be achieved. However, resolvability is not necessarily uniform throughout the map. We introduce a quantitative parameter to characterize the resolvability of individual atoms in cryo-EM maps, the map Q-score. Q-scores can be calculated for atoms in proteins, nucleic acids, water, ligands and other solvent atoms, using models fitted to or derived from cryo-EM maps. Q-scores can also be averaged to represent larger features such as entire residues and nucleotides. Averaged over entire models, Q-scores correlate very well with the estimated resolution of cryo-EM maps for both protein and RNA. Here, assuming the models they are calculated from are well fitted to the map, Q-scores can be used as a measure of resolvability in cryo-EM maps at various scales, from entire macromolecules down to individual atoms. Q-score analysis of multiple cryo-EM maps of the same proteins derived from different laboratories confirms the reproducibility of structural features from side chains down to water and ion atoms.

59 BASIC BIOLOGICAL SCIENCES↗

Coherent Reaction between Molecular and Atomic Bose-Einstein Condensates: Integrable Model

Here, we solve a model that describes a stimulated conversion between ultracold bosonic atoms and molecules. The reaction is triggered by a linearly time-dependent transition throughout the Feshbach resonance. Our solution predicts a dependence, with a dynamic phase transition, of the reaction efficiency on the transition rate for both atoms-to-molecule pairing and molecular dissociation processes. We find that for the latter process with a linear energy dispersion of atomic modes, the emerging phase can have a thermalized energy distribution of noninteracting bosons with the temperature defined by the rate of the transition. This provides a simple interpretation of the phase transition in terms of the creation of equilibrium Bose-Einstein condensate.

74 ATOMIC AND MOLECULAR PHYSICS↗

Multi-qubit entanglement and algorithms on a neutral-atom quantum computer

Gate model quantum computers promise to solve currently intractable computational problems if they can be operated at scale with long coherence times and high fidelity logic. Neutral atom hyperfine qubits provide inherent scalability due to their identical characteristics, long coherence times, and ability to be trapped in dense multi-dimensional arrays. Combined with the strong entangling interactions provided by Rydberg states, all the necessary characteristics for quantum computation are available. Here we demonstrate several quantum algorithms on a programmable gate model neutral atom quantum computer in an architecture based on individual addressing of single atoms with tightly focused optical beams scanned across a two-dimensional array of qubits. Preparation of entangled Greenberger-Horne-Zeilinger (GHZ) states with up to 6 qubits, quantum phase estimation for a chemistry problem, and the Quantum Approximate Optimization Algorithm (QAOA) for the MaxCut graph problem are demonstrated. These results highlight the emergent capability of neutral atom qubit arrays for universal, programmable quantum computation, as well as preparation of non-classical states of use for quantum enhanced sensing.

97 MATHEMATICS AND COMPUTING↗

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗

VeloxChem: A Python-driven density-functional theory program for spectroscopy simulations in high-performance computing environments

An open-source program named VeloxChem has been developed for the calculation of electronic real and complex linear response functions at the levels of Hartree–Fock and Kohn–Sham density functional theories. With an object-oriented program structure written in a Python/C++ layered fashion, VeloxChem enables time-efficient prototyping of novel scientific approaches without sacrificing computational efficiency, so that molecular systems involving up to and beyond 500 second-row atoms (or some 10,000 contracted and in part diffuse Gaussian basis functions) can be routinely addressed. In addition, VeloxChem is equipped with a polarizable embedding scheme for the treatment of the classical electrostatic interactions with an environment that in turn is modeled by atomic site charges and polarizabilities. The underlying hybrid message passing interface (MPI)/open multiprocessing (OpenMP) parallelization scheme makes VeloxChem suitable for execution in high-performance computing cluster environments, showing even slightly beyond linear scaling for the Fock matrix construction with use of up to 16,384 central processing unit (CPU) cores. An efficient—with respect to convergence rate and overall computational cost—multifrequency/gradient complex linear response equation solver enables calculations not only of conventional spectra, such as visible/ultraviolet/X-ray electronic absorption and circular dichroism spectra, but also time-resolved linear response signals as due to ultra-short weak laser pulses.

97 MATHEMATICS AND COMPUTING↗