Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orbit correction”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Electronic Structure Predictions of the Energetic Properties of Tellurium Fluorides

The heats of formation, bond dissociation energies (BDEs), fluoride affinities (FA), fluorocation affinities (FCA), electron affinities (EA), and ionization energies (IP) of TeFn (n = 1 - 6) have been predicted using the Feller-Peterson-Dixon (FPD) approach. To benchmark the approach, the bond dissociation energies of Te 2 and TeO, the heats of formation of Te 2 , TeH 2 , TeO and TeO 2 , and the electron affinity for TeO and TeO 2 were calculated as there are experimental thermodynamic data available for these tellurium compounds, which allow confirmation of the heat of formation of Te gas as ΔH f,0K (Te) = 50.7 ± 0.6 kcal/mol. Spin-orbit corrections are required for good results and cannot be ignored. A comparison between fluoride affinities, fluorocation affinities, electron affinities, and ionization energies of TeF n and SeF n is reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Reactions of U + with H 2 , D 2 , and HD Studied by Guided Ion Beam Tandem Mass Spectrometry and Theory

The kinetic energy-dependent reactions of the atomic actinide uranium cation (U + ) with H 2 , D 2 , and HD were examined by guided ion beam tandem mass spectrometry. An average 0 K bond dissociation energy of D 0 (U + – H) = 2.48 ± 0.06 eV is obtained by analysis of the endothermic product ion cross sections. Quantum chemistry calculations were performed for comparison with experimental thermochemistry, including high-level CASSCF–CASPT2–RASSI calculations of the spin–orbit corrections. CCSD(T) and the CASSCF levels show excellent agreement with experiment, whereas B3LYP and PBE0 slightly overestimate and the M06 approach badly underestimates the bond energy for UH + . Theory was also used to investigate the electronic structures of the reaction intermediates and potential energy surfaces. The experimental product branching ratio for the reaction of U + with HD indicates that these reactions occur primarily via a direct reaction mechanism, despite the presence of a deep-well for UH 2 + formation according to theory. The reactivity and hydride bond energy for U + are compared with those for transition metal, lanthanide, and actinide cations, and periodic trends are discussed. Furthermore, these comparisons suggest that the 5f electrons on uranium are largely core and uninvolved in the reactive chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The Role of Spin–Orbit Coupling in the Linear Absorption Spectrum and Intersystem Crossing Rate Coefficients of Ruthenium Polypyridyl Dyes

The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g., triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using ab initio electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. Here, in this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of 10 13 s -1 for each dye molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2018

Coherent backscattering in quasi-ballistic ultra-high mobility GaAs/AlGaAs 2DES

Abstract A small and narrow negative-magnetoresistance (MR) effect that appears about null magnetic field over the interval −0.025 ≤ B ≤ 0.025 T in magnetotransport studies of the GaAs/AlGaAs 2D system with μ ≈ 10 7 cm 2 / Vs is experimentally examined as a function of the sample temperature, T . The temperature dependent magnetoresistance data were fit using the Hikami et al . theory, without including the spin-orbit correction, to extract the inelastic length, l i , which decreases rapidly with increasing temperature. It turns out that l i < l e , where l e is the elastic length, for all T . Thus, we measured the single particle lifetime, τ s , and the single particle mean free path l s = v F τ s . A comparison between l i and l s indicates that l i > l s . The results suggest that the observed small and narrow magnetoresistance effect about null magnetic field could be a manifestation of coherent backscattering due to small angle scattering from remote ionized donors in the high mobility GaAs/AlGaAs 2DES.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Commissioning of the soft x-ray undulator beamline at the Siam Photon Laboratory

The synchrotron radiation from the first undulator at the Siam Photon Laboratory was characterized with the photon beam position monitors (BPMs) and grating monochromator. The soft x-ray undulator beamline employs a varied line-spacing plane grating monochromator with three interchangeable gratings. Since 2010, the beamline has delivered photons with energy of 40-160 and 220-1040 eV at the resolving power of 10,000 for user services at the two end- stations that utilize the photoemission electron spectroscopy and microscopy techniques. The undulator power-density distributions measured by the 0.05-mm wire-scan BPM were in good agreement with those in simulation. The flux-density distributions were evaluated in the red-shift measurements, which identify the central cone of radiation and its distribution. Since 2014, the operation of the other insertion devices in the storage ring has started, and consequently bought about the increases in the emittance from 41 to 61 nm·rad and the coupling constant from 4 to 11%. The local electron-orbit correction greatly improved the alignment of the electron beam in the undulator section resulting in the improvements of the photon flux and harmonics peaks of the undulator radiation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ a many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt2 and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt2, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu 3d electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I 2 , so that experiment and theory are brought into imperfect but reasonable agreement within about 120 meV/atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Formation of ECR Plasma in a Dielectric Plasma Guide under Self-Excitation of a Standing Ion-Acoustic Wave

The formation of a spatially localized plasma with a high brightness has been experimentally observed in a dielectric plasma guide under the electron cyclotron resonance discharge at the excitation of a standing ion-acoustic wave. The results obtained show the possibility of designing compact high-intensity radiation sources with a spectrum determined by the working gas or gas mixture type, high-intensity chemically active particle flow sources, and plasma thrusters for correcting orbits of light spacecraft.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

STATUS OF THE JLEIC ION COLLIDER RING DESIGN

We present an update on the lattice design and beam dynamics study of the ion collider ring of JLEIC (Jefferson Lab Electron Ion Collider). The collider ring consists of two 261.7 degree arcs connected by two straight sections crossing each other. One of the straights houses an interaction region (IR) and is shaped to make a 50 mrad crossing angle with the electron beam at the interaction point (IP) to meet physics requirements. The forward acceptance requirements downstream of the IP in the ion direction lead to an asymmetric IR lattice design. The detector solenoid effects and the multipole fields of the IR magnets further complicate this picture. In this paper, compensation of the detector solenoid effects is considered together with orbit correction and multipole effects. We also study local compensation of the magnet multipoles using dedicated multipole correctors. And an optimization of the betatron tunes is also presented.

Wei, Guohui↗
DOE OSTI2021

Detector Solenoid Compensation for the Electron-Ion Collider

The central detector in the present EIC design includes a 4 m long solenoid with an integrated strength of up to 12 Tm. The electron beam passes on-axis through the solenoid, but the hadron beam has an angle of 25 mrad. Thus the solenoid couples horizontal and vertical betatron motion in both electron and hadron storage rings, and causes a vertical closed orbit excursion in the hadron ring. The solenoid also couples the transverse and longitudinal motions of both beams, when crab cavities are also considered. In this paper, we present schemes for closed orbit correction and coupling compensation at the IP, including crabbing.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Universal Injector at LERF – Baseline Design

A new electron injector complex has been proposed within the Low Energy Recirculating Facility (LERF) vault, which would serve both the 22 GeV CEBAF and the positron programs, while being compatible with the electron source required to produce positrons for Ce+BAF. The baseline design of a 3-pass recirculator features a main linac configured with three cryomodules, five isochronous arcs and three straight sections, allowing electron beam acceleration up to 650 MeV. The injector offers the flexibility to extract 1-pass and 2-pass energy electrons, as needed for positron production, with final reconfiguration to a 3-pass 650 MeV injector for 22 GeV CEBAF. This paper provides an overview of the baseline design of the Universal injector complex, including the 8 MeV injector and the recirculator. A comprehensive suite of beam dynamics studies to validate our design is being launched, starting with start-to-end tracking with Elegant followed by the orbit correction and error mitigation.

Bogacz, A. [Thomas Jefferson National Accelerator ↗
DOE OSTI2025

Virtual to Physical: Reinforcement Learning to Optimize SNS Particle Accelerator Controls

Complex accelerators must have control systems that can handle dynamic nonlinear environments. This makes traditional control methods unsuitable as they can struggle to adapt to these uncertainties. This provides an ideal environment for reinforcement learning algorithms as they are adaptable and generalizable. We present a reinforcement learning pipeline that can effectively handle the dynamics of a complex accelerator. We test and prove our pipelines capabilities on multiple environments including the Spallation Neutron Source (SNS) and the Beam Test Facility (BTF) at Oakridge National Lab (ORNL). Due to the limited time available to train an online algorithm like reinforcement learning on a real accelerator, we utilize a virtual twin accelerator (VIRAC) developed by ORNL to pretrain the policy and show its ability to converge in the virtual environment. We then test the adaptability of the pretrained RL model by applying it on the real accelerator and comparing the results. Utilizing our Scientific Optimization and Controls Toolkit (SOCT) and open-source standards such as Gymnasium we create and solve for a MEBT orbit correction problem in the SNS and an emittance maximization problem in the BTF. We show how Twin Delayed Deep Deterministic Policy Gradient (TD3) can solve this optimization environment in the virtual accelerator and transfer this policy onto the real accelerator for inference and model retraining. We show how reinforcement learning can be utilized as a control system for complex accelerators and provide a model pipeline for how an implementation performs and can be adapted to new accelerator control problems.

Kasparian, Armen [Thomas Jefferson National Accele↗
DOE OSTI2026

‘Universal Injector’ at LERF - Layout and Optics Architecture

Here, we present a next level design of a compact injector within the LERF vault, which would serve both the 22 GeV CEBAF and the positron program, while being compatible with the electron source required to produce positrons for Ce+BAF. The baseline design of a 3-pass recirculator features a main linac configured with three C-75 cryo-modules, five isochronous return arcs and three straight sections, facilitating electron beam acceleration up to 650 MeV. The Universal Injector offers flexibility of extraction 1-pass and 2-pass energy electrons, as needed for positron production, with a final 3-pass extraction required by the 650 MeV injector for 22 GeV CEBAF. This note provides an overview of the baseline optics design of the Universal Injector complex, including the 8 MeV injector merger and the recirculator racetrack, including individual extraction lines for all three passes. A comprehensive suite of beam dynamics studies to validate the design is under way; starting with the orbit correction scheme, followed by start-to-end tracking with lattice misalignments and magnet errors.

Bogacz, Alex [Thomas Jefferson National Accelerato↗