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 217 records · Page 12

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗

A structural coarse-grained model for clays using simple iterative Boltzmann inversion

Cesium-137 is a crucial byproduct of nuclear energy generation and is environmentally threatening due to its long half-life and affinity for naturally occurring micaceous clays. Current experimental observations of illite and phlogopite mica indicate that Cs+ is capable of exchanging with K+ bound in the anhydrous interlayers of layered silicates, forming sharp exchange fronts, leading to interstratification of Cs- and K-illite. We introduce here a coarse-grained (CG) model of the anhydrous illite interlayer developed using iterative Boltzmann inversion that qualitatively and quantitatively reproduces features of a previously proposed feedback mechanism of ion exchange. The CG model represents a 70-fold speedup over all-atom models of clay systems and predicts interlayer expansion for K-illite near ion exchange fronts. Contrary to the longstanding theory that ion exchange in a neighboring layer increases the binding of K in lattice counterion sites leading to interstratification, we find that the presence of neighboring exchanged layers leads to short-range structural relaxations that increase basal spacing and decrease cohesion of the neighboring K-illite layers. We also provide evidence that the formation of alternating Cs- and K-illite interlayers (i.e., ordered interstratification) is both thermodynamically and mechanically favorable compared to exchange in adjacent interlayers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray spectroscopy of planar laser-plasma interaction experiments at the National Ignition Facility

X-ray spectroscopy has been newly used to diagnose electron temperatures in planar-geometry experiments at the National Ignition Facility (NIF) designed to study laser-plasma interactions at plasma conditions relevant to direct-drive ignition. These experiments used a buried co-mixed Mn/Co microstrip in a CH ablator in experiments that generated long scale-length plasmas susceptible to stimulated Raman scattering (SRS). Time-resolved Mn and Co K-shell spectra, diagnosed using the NIF x-ray spectrometer, were analyzed by fi tting to synthetic spectra based on a detailed atomic model of emission from the microstrip. The electron temperature at the time when the microstrip passes through the quarter-critical density surface, the key region for the development of SRS, was inferred to be around 2 to 3 keV. These measurements constrain 2-D DRACO radiation-hydrodynamics modeling of the planar experiments, important for determining plasma conditions pertinent to SRS, and demonstrate that this platform approaches direct-drive ignition-relevant conditions. The modeling is also assessed by a direct comparison of measured spectra to modeled spectra generated by DRACO in conjunction with the atomic physics postprocessor code SPECT3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the physical and chemical details of alumina atomic layer deposition: A combined experimental and numerical approach

Alumina thin film is typically studied as a model atomic layer deposition (ALD) process due to its high dielectric constant, high thermal stability, and good adhesion on various wafer surfaces. Despite extensive applications of alumina ALD in microelectronics industries, details on the physical and chemical processes are not yet well understood. ALD experiments are not able to shed adequate light on the detailed information regarding the transient ALD process. Most of current numerical approaches lack detailed surface reaction mechanisms, and their results are not well correlated with experimental observations. In this paper, the authors present a combined experimental and numerical study on the details of flow and surface reactions in alumina ALD using trimethylaluminum and water as precursors. Results obtained from experiments and simulations are compared and correlated. By experiments, growth rate on five samples under different deposition conditions is characterized. The deposition rate from numerical simulation agrees well with the experimental results. Details of precursor distributions in a full cycle of ALD are studied numerically to bridge between experimental observations and simulations. The 3D transient numerical model adopts surface reaction kinetics and mechanisms based on atomic-level studies to investigate the surface deposition process. Surface deposition is shown as a strictly self-limited process in our numerical studies. ALD is a complex strong-coupled fluid, thermal and chemical process, which is not only heavily dependent on the chemical kinetics and surface conditions but also on the flow and material distributions.

42 ENGINEERING↗

Computed profiles of the solar C I multiplets at 1561 and 1657 A

A nine-level atomic model is presented for C I lines synthesis in the sun. All the atomic parameters are based on recent experimental and theoretical data. Profiles of the multiplets at 1561 and 1657 A are computed from a model of the solar chromosphere. The sensitivity of the results to changes in the atomic parameters is examined. Partial frequency redistribution must be included in these lines to obtain agreement between calculated and observed profiles.

Mauas, Pablo J.↗

Numerical modeling for primary atomization of liquid jets

In the proposed numerical model for primary atomization, surface-wave dispersion equations are solved in conjunction with the jet-embedding technique of solving mean flow equations of a liquid jet. Linear and approximate nonlinear models have been considered. In each case, the dispersion equation is solved over the whole wavelength spectrum to predict drop sizes, frequency, and liquid-mass breakup rates without using any empirical constants. The present model has been applied to several low-speed and high-speed jets. For the high-speed case (the LOX/H2 coaxial injector of the Space Shuttle Main Engine Preburner), predicted drop sizes and liquid breakup rates are in good agreement with the results of the CICM code, which have been calibrated against measured data.

Przekwas, A. J.↗

X-Ray Lines from Dark Matter Annihilation at the keV Scale

In 2014, several groups reported hints for a yet unidentified line in astrophysical x-ray signals from galaxies and galaxy clusters at an energy of 3.5 keV. While it is not unlikely that this line is simply a reflection of imperfectly modeled atomic transitions, it has renewed the community’s interest in models of keV-scale dark matter, whose decay would lead to such a line. The alternative possibility of dark matter annihilation into monochromatic photons is far less explored, a lapse that we strive to amend in this Letter. More precisely, we introduce a novel model of fermionic dark matter χ with O(keV) mass, annihilating to a scalar state Φ which in turn decays to photons, for instance via loops of heavy vectorlike fermions. The resulting photon spectrum is box shaped, but if χ and Φ are nearly degenerate in mass, it can also resemble a narrow line.We discuss dark matter production via two different mechanisms—misalignment and freeze-in —which both turn out to be viable in vast regions of parameter space. Furthermore, we constrain the model using astrophysical x-ray data, and we demonstrate that, thanks to the velocity dependence of the annihilation cross section, it has the potential to reconcile the various observations of the 3.5 keV line. We finally argue that the model can easily avoid structure formation constraints on keV-scale dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding the Cu-Zn brass alloys using a short-range-order cluster model: Significance of specific compositions of industrial alloys

Metallic alloys show complex chemistries that are not yet understood so far. It has been widely accepted that behind the composition selection lies a short-range-order mechanism for solid solutions. The present paper addresses this fundamental question by examining the face-centered-cubic Cu-Zn a-brasses. A new structural approach, the cluster-plus-glue-atom model, is introduced, which suits specifically for the description of short-range-order structures in disordered systems. Two types of formulas are pointed out, [Zn-Cu 12 ]Zn 1,6 and [Zn-Cu 12 ](Zn,Cu) 6 , which explain the a-brasses listed in the American Society for Testing and Materials (ASTM) specifications. In these formulas, the bracketed parts represent the 1 st -neighbor cluster, and each cluster is matched with one to six 2 nd -neighbor Zn atoms or with six mixed (Zn,Cu) atoms. Such a cluster-based formulism describes the 1 st - and 2 nd -neighbor local atomic units where the solute and solvent interactions are ideally satisfied. The Cu-Ni industrial alloys are also explained, thus proving the universality of the cluster-formula approach in understanding the alloy selections. The revelation of the composition formulas for the Cu-(Zn,Ni) industrial alloys points to the common existence of simple composition rules behind seemingly complex chemistries of industrial alloys, therefore offering a fundamental and practical method towards composition interpretations of all kinds of alloys.

36 MATERIALS SCIENCE↗

Ionic Transport Coefficients of Dense Plasmas without Molecular Dynamics

In this paper, we present a theoretical model that allows a fast and accurate evaluation of ionic transport properties of realistic plasmas spanning from warm and dense to hot and dilute conditions, including mixtures. This is achieved by combining a recent kinetic theory based on effective interaction potentials with a model for the equilibrium radial density distribution based on an average atom model and the integral equations theory of fluids. Lastly, the model should find broad use in applications where nonideal plasma conditions are traversed, including inertial confinement fusion, compact astrophysical objects, solar and extrasolar planets, and numerous present-day high energy density laboratory experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

AToM (Advanced Tokamak Modeling)

This final report for the AToM project documents the contributions from General Atomics over the period 1 September 2014 to 31 August 2017. This includes significant progress on the development of tokamak integrated modeling frameworks (OMFIT and IPS), and development of new workflows for computing plasma profile (temperature, density, etc) evolution. Significant effort was also carried out in the high-performance computing arena, by adding the COGENT edge-simulation code to the AToM portfolio, and performing optimization of the CGYRO plasma turbulence code. Issues of data management and user outreach are also documented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisional-radiative modeling of tungsten at temperatures of 1200–2400 eV

We discuss new collisional-radiative modeling calculations of tungsten at moderate temperatures of 1200 to 2400 eV. Such plasma conditions are relevant to ongoing experimental work at ASDEX Upgrade and are expected to be relevant for ITER. Our calculations are made using the Los Alamos National Laboratory (LANL) collisional-radiative modeling ATOMIC code. These calculations formed part of a submission to the recent NLTE-8 workshop that was held in November 2013. This series of workshops provides a forum for detailed comparison of plasma and spectral quantities from NLTE collisional-radiative modeling codes. We focus on the LANL ATOMIC calculations for tungsten that were submitted to the NLTE-8 workshop and discuss different models that were constructed to predict the tungsten emission. In particular, we discuss comparisons between semi-relativistic configuration-average and fully relativistic configuration-average calculations. We also present semi-relativistic calculations that include fine-structure detail, and discuss the difficult problem of ensuring completeness with respect to the number of configurations included in a CR calculation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Adsorption and recombination of hydrogen atoms on a model graphite surface

The adsorption and recombination of atomic hydrogen on a model graphite grain have been examined in a series of calculations in which a modified, iterative, extended Hueckel program was used. The hydrogen atom is found to be chemisorbed at a site with a zero-point binding energy of 0.7 eV and at an equilibrium distance of 2.25 A above the site. Despite a barrier of about 0.4 eV between adjacent sites, calculations suggest that at temperatures as low as 10 K, an H atom will tunnel through to adjacent sites in less than one nanosecond. However, a potential barrier to the recombination of two hydrogen atoms has been found which displays high sensitivity to the mutual arrangement of the two hydrogen atoms with respect to the graphite surface. Results show that at very low temperatures, recombinations can occur only by tunneling. Consistent with experiment, the region in which H2 begins to form exhibits a repulsive potential with respect to possible chemisorption of the incipient H2 entity.

Aronowitz, S.↗

Coupled Experimental–Theoretical Characterization of a Carbon Electrode in Vanadium Redox Flow Batteries using X-ray Absorption Spectroscopy

Vanadium redox flow batteries (VRFBs) have emerged as promising solutions for stationary grid energy storage due to their high efficiency, scalability, safety, near room-temperature operation conditions, and the ability to independently size power and energy capacities. The performance of VRFBs heavily relies on the redox couple reactions of V 2+ /V 3+ and VO 2+ /VO 2 + on carbon electrodes. Therefore, a thorough understanding of the surface functionality of carbon electrodes and their propensity for degradation during electrochemical cycles is crucial for designing VRFBs with extended lifespans. In this study, we present a coupled experimental–theoretical approach based on carbon K edge X-ray absorption spectroscopy (XAS) to characterize carbon electrodes prepared under different conditions and identify relevant functional groups that contribute to unique spectroscopic features. Atomic models were created to represent functional groups, such as hydroxyl, carboxyl, methyl, and aldehyde, bonded to carbon atoms in either sp 2 or sp 3 environments. The interactions between functionalized carbon and various solvated vanadium complexes were modeled using density functional theory. A library of carbon K-edge XAS spectra was generated for distinct carbon atoms in different functional groups, both before and after interacting with solvated vanadium complexes. Here we demonstrate how these simulated spectra can be used to deconvolve ex situ experimental spectra measured from carbon electrodes and to track changes in the electrode composition following immersion in different electrolytes or extended cycling within a functional VRFB. By doing so, we identify the active species present on the carbon electrodes, which play a crucial role in determining their electrochemical performance.

25 ENERGY STORAGE↗

Improved chemistry restraints for crystallographic refinement by integrating the Amber force field into Phenix

The refinement of biomolecular crystallographic models relies on geometric restraints to help to address the paucity of experimental data typical in these experiments. Limitations in these restraints can degrade the quality of the resulting atomic models. Here, an integration of the full all-atom Amber molecular-dynamics force field into Phenix crystallographic refinement is presented, which enables more complete modeling of biomolecular chemistry. The advantages of the force field include a carefully derived set of torsion-angle potentials, an extensive and flexible set of atom types, Lennard–Jones treatment of nonbonded interactions and a full treatment of crystalline electrostatics. The new combined method was tested against conventional geometry restraints for over 22 000 protein structures. Structures refined with the new method show substantially improved model quality. On average, Ramachandran and rotamer scores are somewhat better, clashscores and MolProbity scores are significantly improved, and the modeling of electrostatics leads to structures that exhibit more, and more correct, hydrogen bonds than those refined using traditional geometry restraints. In general it is found that model improvements are greatest at lower resolutions, prompting plans to add the Amber target function to real-space refinement for use in electron cryo-microscopy. This work opens the door to the future development of more advanced applications such as Amber -based ensemble refinement, quantum-mechanical representation of active sites and improved geometric restraints for simulated annealing.

59 BASIC BIOLOGICAL SCIENCES↗

Simulation of Heterogeneous Atom Probe Tip Shape Evolution during Field Evaporation Using a Level Set Method and Different Evaporation Models

Accurate reconstruction of spatial positions of field evaporated ions from detector patterns depends on a correct understanding of the tip shape evolution and evaporation laws of component atoms. Reconstruction artifacts introduced in the heterogeneous tip can be contributed by the assumption on homogeneous evaporation behavior. A level set method based specimen shape evolution model is developed in this study to simulate the evaporation of synthetic layered-structured atom probe tip. The simulation results of the shape evolution by the present level set model agree well with the literature data using other numerical approaches, i.e. finite difference method. The asymmetric evolving shape predicted by the level set model demonstrates that the complex evaporation behavior of heterogeneous tips and interface curvature can potentially lead to the tip reconstruction distortion. Compared with other shape evaporation methods the new method provides considerably higher resolution and computational efficiency with the aid of the intrinsic subgrid scale accuracy.

atom probe tomography, field evaporation, tip shap↗

Using Beta Cen to Probe the Photosphere-Wind Connection to B Stars

The goal of the project was to probe the connection between the photospheric pulsations and the X-ray production in the stellar winds of B stars. The B1 III star Beta Centauri is a Beta Cephei variable with a high X-ray count rate, and was observed by XMM. Analysis of XMM X-ray spectra showed characteristic plasma temperatures of 0.1, 0.2 and 0.6 keV. Analysis of He-like forbidden-to-intercombination lines was performed to determine the dominant locations of the X-ray emission. At Prism, work focused on developing and upgrading spectral models to calculate forbidden-to-intercombination line ratios using more detailed atomic models than the simple models typically used in the analysis of X-ray spectral lines, thereby allowing us to assess the effects of a broader variety of populating and depopulating channels for energy levels contributing to the observed X-ray line emission. Based on Ne IX lines, the source location was determined to be no farther than 4 stellar radii. Examination of the time variability of the XMM data for Beta Cen showed no significant correlation with the known optical pulsation period for the star, and also did not find any significant periods of variation at X-ray wavelengths.

MacFarlane, Joseph J.↗

Improved chemistry restraints for crystallographic refinement by integrating the Amber force field into Phenix

The refinement of biomolecular crystallographic models relies on geometric restraints to help address the paucity of experimental data typical in these experiments. Limitations in these restraints can degrade the quality of the resulting atomic models. Here we present an integration of the full all-atom Amber molecular dynamics force field into Phenix crystallographic refinement, which enables a more complete modeling of biomolecular chemistry. The advantages of the force field include a carefully derived set of torsion angle potentials, an extensive and flexible set of atom types, Lennard-Jones treatment of non-bonded interactions and a full treatment of crystalline electrostatics. The new combined method was tested against conventional geometry restraints for over twenty-two thousand protein structures. Structures refined with the new method show substantially improved model quality. On average, Ramachandran and rotamer scores are somewhat better; clash scores and MolProbity scores are significantly improved; and the modelling of electrostatics leads to structures that exhibit more, and more correct, hydrogen bonds than those refined with traditional geometry restraints. We find in general that model improvements are greatest at lower resolutions, prompting plans to add the Amber target function to real-space refinement for use in electron cryo-microscopy. This work opens the door to the future development of more advanced applications such as Amber-based ensemble refinement, quantum mechanical representation of active sites and improved geometric restraints for simulated annealing. This document contains embedded data - to preserve data integrity, please ensure where possible that the IUCr Word tools (available from http://journals.iucr.org/services/docxtemplate/ ) are installed when editing this document. The full Amber force field has been integrated into Phenix as an alternative refinement target. With a slight loss in speed, it achieves improved stereochemistry, fewer steric clashes and better hydrogen bonds.

59 BASIC BIOLOGICAL SCIENCES↗