Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomic energy”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Application of Simulated Annealing and Related Algorithms to TWTA Design

Simulated Annealing (SA) is a stochastic optimization algorithm used to search for global minima in complex design surfaces where exhaustive searches are not computationally feasible. The algorithm is derived by simulating the annealing process, whereby a solid is heated to a liquid state and then cooled slowly to reach thermodynamic equilibrium at each temperature. The idea is that atoms in the solid continually bond and re-bond at various quantum energy levels, and with sufficient cooling time they will rearrange at the minimum energy state to form a perfect crystal. The distribution of energy levels is given by the Boltzmann distribution: as temperature drops, the probability of the presence of high-energy bonds decreases. In searching for an optimal design, local minima and discontinuities are often present in a design surface. SA presents a distinct advantage over other optimization algorithms in its ability to escape from these local minima. Just as high-energy atomic configurations are visited in the actual annealing process in order to eventually reach the minimum energy state, in SA highly non-optimal configurations are visited in order to find otherwise inaccessible global minima. The SA algorithm produces a Markov chain of points in the design space at each temperature, with a monotonically decreasing temperature. A random point is started upon, and the objective function is evaluated at that point. A stochastic perturbation is then made to the parameters of the point to arrive at a proposed new point in the design space, at which the objection function is evaluated as well. If the change in objective function values (Delta)E is negative, the proposed new point is accepted. If (Delta)E is positive, the proposed new point is accepted according to the Metropolis criterion: rho((Delta)f) = exp((-Delta)E/T), where T is the temperature for the current Markov chain. The process then repeats for the remainder of the Markov chain, after which the temperature is decremented and the process repeats. Eventually (and hopefully), a near-globally optimal solution is attained as T approaches zero. Several exciting variants of SA have recently emerged, including Discrete-State Simulated Annealing (DSSA) and Simulated Tempering (ST). The DSSA algorithm takes the thermodynamic analogy one step further by categorizing objective function evaluations into discrete states. In doing so, many of the case-specific problems associated with fine-tuning the SA algorithm can be avoided; for example, theoretical approximations for the initial and final temperature can be derived independently of the case. In this manner, DSSA provides a scheme that is more robust with respect to widely differing design surfaces. ST differs from SA in that the temperature T becomes an additional random variable in the optimization. The system is also kept in equilibrium as the temperature changes, as opposed to the system being driven out of equilibrium as temperature changes in SA. ST is designed to overcome obstacles in design surfaces where numerous local minima are separated by high barriers. These algorithms are incorporated into the optimal design of the traveling-wave tube amplifier (TWTA). The area under scrutiny is the collector, in which it would be ideal to use negative potential to decelerate the spent electron beam to zero kinetic energy just as it reaches the collector surface. In reality this is not plausible due to a number of physical limitations, including repulsion and differing levels of kinetic energy among individual electrons. Instead, the collector is designed with multiple stages depressed below ground potential. The design of this multiple-stage collector is the optimization problem of interest. One remaining problem in SA and DSSA is the difficulty in determining when equilibrium has been reached so that the current Markov chain can be terminated. It has been suggested in recent literature that simulating the thermodynamic properties opecific heat, entropy, and internal energy from the Boltzmann distribution can provide good indicators of having reached equilibrium at a certain temperature. These properties are tested for their efficacy and implemented in SA and DSSA code with respect to TWTA collector optimization.

Radke, Eric M.↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2026

Consistent inclusion of triple substitutions within a coupled cluster based static quantum embedding theory

We have previously proposed the MPCC static embedding framework for quantum chemistry that self-consistently couples a high-level coupled cluster (CC) treatment of the fragment (active region) with a lower level, Møller–Plesset perturbation treatment of the environment. Our initial implementation was limited to single and double (SD) substitutions, with CCSD for the fragment and first-order perturbative SD amplitudes for the environment. Here, we extend the MPCC embedding treatment to triple substitutions, which is essential for achieving chemical accuracy in energy differences. To this end, we employ a CCSDT solver for the fragment subsystem. For the environment subsystem, we construct a perturbative estimate of the triples amplitudes, explicitly accounting for feedback from all fragment amplitudes. The resulting approach is denoted MPCCSDT(pt). We further introduce a more complete formulation in which feedback from the environment amplitudes to the fragment amplitudes is also included. This scheme involves an iterative treatment of the environment triples amplitudes and is denoted MPCCSDT(it). In addition, we assess the accuracy of the previously proposed low-level method by introducing a modified low-level approach that incorporates a lowest-order treatment of selected long-range effects, including spin fluctuations and charge polarization. All resulting approaches may be viewed as post-CCSD(T) methods. We therefore consider test cases for which CCSD(T) exhibits substantial deviations from CCSDT. These include (i) single- and triple-bond stretching in F 2 and N 2 , (ii) bond dissociation energies of selected molecules from the W4-11 dataset, and (iii) total atomization energies of transition metal hydrides. Our results demonstrate that inclusion of triples amplitudes at the fragment level alone is insufficient; a perturbative treatment of the environment triples amplitudes is required. For many energy-difference applications, feedback from the environment triples amplitudes to the fragment amplitudes is not essential, but it does play a role in the very challenging CoH and FeH molecules. A very interesting finding from our study is that in some challenging cases, we need an improved (second-order) perturbative method for the SD amplitudes, going beyond the first-order one used in our earlier work. In conclusion, considering both cost and accuracy, the MP2CCSDT(pt) model is the most promising for future applications among the candidates considered here.

Shee, Avijit [University of California, Berkeley, ↗
DOE OSTI2026

Particle-hole symmetric slave-boson method for the mixed valence problem

We introduce an analytic slave-boson method for treating the finite-𝑈 Anderson impurity model. Our approach introduces two bosons to track both 𝑄 ⇌ 𝑄 ± 1 valence fluctuations and reduces to a single symmetric 𝑠 boson in the effective action, encoding all the high-energy atomic physics information in the boson's kinematics, while the low-energy part of the action remains unchanged across finite-𝑈, infinite-𝑈, and Kondo limits. We recover the infinite-𝑈- and Kondo-limit actions from our approach and show that the Kondo resonance already develops in the normal state when the slave boson has yet to condense. We show that the slave rotor and 𝑠 boson have the same algebraic structure, and we establish a unified functional integral framework connecting the 𝑠-boson and slave-rotor representations for the single-impurity Anderson model.

Anderson impurity model↗
DOE OSTI2024

Summary of Leach-Burn-Leach Round-Robin Test Results

Tristructural isotropic (TRISO) coated particle fuel is of interest to multiple member states associated with the Generation IV International Forum (GIF). Qualification and acceptance testing of TRISO fuel includes the determination of particle defect fractions and impurities in the as-fabricated TRISO fuel particle and/or composited fuel form, which are typically determined using burn-leach (BL) or leach-burn-leach (LBL) analysis methods. A round-robin test focusing on LBL analysis was undertaken by the GIF Very High Temperature Reactor System Fuel and Fuel Cycle Project Management Board. Participating in this LBL round-robin test were research staff from the Institute of Nuclear and New Energy Technology (INET) in the People’s Republic of China, the Korean Atomic Energy Research Institute (KAERI) in the Republic of Korea, and Oak Ridge National Laboratory (ORNL) in the United States of America. The motivation of this LBL benchmarking exercise was to explore consistency in results generated by LBL and to elucidate the potential influence of variations in the procedures used by the participating members on the accurate determination of the measured fuel properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Environmental Monitoring Report, Salmon, Mississippi, Site 2019

This report presents the monitoring results for groundwater and surface water samples collected in October 2019 by the U.S. Department of Energy (DOE) Office of Legacy Management for the Salmon, Mississippi, Site. The U.S. Atomic Energy Commission (AEC), a predecessor agency to DOE, conducted a series of underground detonations in the Tatum Salt Dome beneath the Salmon site between 1964 and 1970 to study seismic signatures. Groundwater sampling at the site is conducted to monitor shallow groundwater that became contaminated during AEC use of the site, contamination within the shot cavity, and contamination in the Aquifer 5 injection well. Surface water sampling is conducted to ensure no contaminants are leaving the site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Final state interactions for high energy scattering off atomic electrons

We consider the scattering of high energy leptons off bound atomic electrons focusing primarily on final state interactions between the outgoing energetic electron, and the heavy residual charged “debris” in the final state. These effects are inherently absent from calculations for a free electron at rest. Coulomb exchanges are enhanced by the large number of electrons in the atomic debris, and are unsuppressed by nonrelativistic velocities in the debris. We find that these exchanges can be resummed using operator methods, and cancel at the level of the cross section until at least 𝑂⁡(𝛼 3 ). Furthermore, we argue that both final and initial state Coulomb exchanges (enhanced by the number of electrons in the atom) do not affect the cross section until at least 𝑂⁡(𝛼 3 ). Transverse photon couplings to nonrelativistic electrons are proportional to their small velocities, and rotational invariance suppresses their contribution to 𝑂⁡(𝛼 3 ). Our results are relevant for precision experiments involving neutrinos, electrons, positrons, and muons scattering off of atomic electrons in a fixed target.

Atoms↗
DOE OSTI2020

Incorporating new evaluation on angular distributions and energy spectra of neutron-induced charged particle reactions into the next ENDF library release (Technical Progress Report)

This project aims to provide the improved Evaluated Nuclear Data File (ENDF) using the newly measured data as well as the latest nuclear reaction model for calculating angular distributions and energy spectra on neutron-induced charged particle reactions through the collaboration of Korea Atomic Energy Research Institute (KAERI) and Los Alamos National Laboratory (LANL). The LANL group will provide the experimental data for angular distributions and spectra of (n,p) and (n,α) on several structural materials such as Fe, Ni and Zn isotopes using the Low Energy Neutron-induced Charged-particle (Z) Chamber (LENZ) instrument at Los Alamos Neutron Science Center (LANSCE). The KAERI group will provide the improved evaluated nuclear library which is based on the LANL experimental data, and further will predict angular distributions and spectra of (n,p) and (n,α) reactions on unmeasured nuclides, such as Cr, Mn, Co, Cu and so on. For the first year of this project, we planned to analyze (n,p) and (n,α) reactions for 54,56 Fe and perform new measurements on those reactions for 58,60 Ni isotopes with the LENZ instrument at LANSCE. For improving our evaluation quality, we have studied reaction models to reproduce LANL’s experimental angular distributions and energy spectra using the full Hauser-Feshbach model code, CoH3 with no approximations used. As the first year’s deliverables, we provided the experimental (n,p) and (n,α) reaction cross sections for 54,56 Fe and incorporate new evaluation on angular distributions and energy spectra of neutron-induced charged particle reactions into the current ENDF/B-VIII.0.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

First atomic weapons following WWII were detonated 75 years ago Operation Crossroads kicks off era of testing Los Alamos-created weapons [Slides]

It was time to test. After the Los Alamos-created atomic bombs helped end World War II, the no-longer-secret Lab transitioned into an era of weapons testing, starting 75 years ago with Operation Crossroads with the goal of studying nuclear weapons’ effects on warships. In August 1945, U.S. Senator Brien McMahon, who later authored the Atomic Energy Act of 1946, said: “In order to test the destructive powers of the atomic bomb against naval vessels, I would like to see these (Japanese naval) ships taken to sea and an atomic bomb dropped on them. The resulting explosion should prove to us just how effective the atomic bomb is when used against the giant naval ships.”

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2016

Atomic Oxygen Erosion Yield Dependence Upon Texture Development in Polymers

The atomic oxygen erosion yield (volume of a polymer that is lost due to oxidation per incident atom) of polymers is typically assumed to be reasonably constant with increasing fluence. However polymers containing ash or inorganic pigments, tend to have erosion yields that decrease with fluence due to an increasing presence of protective particles on the polymer surface. This paper investigates two additional possible causes for erosion yields of polymers that are dependent upon atomic oxygen. These are the development of surface texture which can cause the erosion yield to change with fluence due to changes in the aspect ratio of the surface texture that develops and polymer specific atomic oxygen interaction parameters. The surface texture development under directed hyperthermal attack produces higher aspect ratio surface texture than isotropic thermal energy atomic oxygen attack. The fluence dependence of erosion yields is documented for low Kapton H (DuPont, Wilmington, DE) effective fluences for a variety of polymers under directed hyperthermal and isotropic thermal energy attack.

oxygen atoms↗
DOE OSTI2023

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
NASA NTRS1984

Atomic mean excitation energies for stopping powers from local plasma oscillator strengths

The stopping of a charged particle by isolated atoms is investigated theoretically using an 'atomic plasma' model in which atomic oscillator strengths are replaced by the plasma frequency spectrum. The plasma-frequency correction factor for individual electron motion proposed by Pines (1953) is incorporated, and atomic mean excitation energies are calculated for atoms through Sr. The results are compared in a graph with those obtained theoretically by Inokuti et al. (1978, 1981) and Dehmer et al. (1975) and with the experimental values compiled by Seltzer and Berger (1982): good agreement is shown.

Wilson, J. W.↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1997

Residual Stress Measurements After Proof and Flight: ETP-0403

The intent of this testing was to evaluate the residual stresses that occur in and around the attachment details of a case stiffener segment that has been subjected to flight/recovery followed by proof loading. Not measured in this test were stresses relieved at joint disassembly due to out-of-round and interference effects, and those released by cutting the specimens out of the case segment. The test article was lightweight case stiffener segment 1U50715, S/N L023 which was flown in the forward stiffener position on flight SRM 14A and in the aft position on flight SRM24A. Both of these flights were flown with the 3 stiffener ring configuration. Stiffener L023 had a stiffener ring installed only on the aft stub in its first flight, and it had both rings installed on its second flight. No significant post flight damage was found on either flight. Finally, the segment was used on the DM-8 static test motor in the forward position. No stiffener rings were installed. It had only one proof pressurization prior to assignment to its first use, and it was cleaned and proof tested after each flight. Thus, the segment had seen 3 proof tests, two flight pressurizations, and two low intensity water impacts prior to manufacturing for use on DM-8. On DM-8 it received one static firing pressurization in the horizontal configuration. Residual stresses at the surface and in depth were evaluated by both the x-ray diffraction and neutron beam diffraction methods. The x-ray diffraction evaluations were conducted by Technology for Energy Corporation (TEC) at their facilities in Knoxville, TN. The neutron beam evaluations were done by Atomic Energy of Canada Limited (AECL) at the Chalk River Nuclear Laboratories in Ontario. The results showed general agreement with relatively high compressive residual stresses on the surface and moderate to low subsurface tensile residual stresses.

Webster, Ronald L..↗
NASA NTRS1981

Spatial aspects of electron energy degradation in atomic oxygen

Spatial (radial and longitudinal) yield spectra for electron energy degradation in atomic oxygen have been obtained using a Monte Carlo method for 25 eV to 10 keV incident electrons. Four-dimensional yield spectra have been analytically represented in terms of a model containing three simple microplumes. We find that the scaled spatial yield spectra for O is approximately the same as for N2. This feature provides a basis for inferring yield spectra for any atmosphere gas or mixture of gases.

Singhal, R. P.↗